Le 24/07/2020 à 13:11, Samuel Gougeon a écrit :
.../...
How to apply any absolute tolerance is quite clear (although unsafe).
IMHO the absolute tolerance must be set to 0 (so not keeping the
default one set to 1e-10), for the reason given hereabove.
I meant, in %r_simp(). In %r_clean(),
Le 20/03/2020 à 14:42, Stéphane Mottelet a écrit :
(num,den) calling style seems to be OK:
--> [a,b]=simp(prod(q.num),prod(q.den))
a =
2 3 4
3.432D-17 +1.230D-16s +3.079D-16s +5.709D-17s +3.432D-17s
b =
Stéphane,
In the meantime, a workaround is to toggle simplification off:
simp_mode(%f).
This inhibits simplification.
Regards,
Federico Miyara
On 20/03/2020 10:42, Stéphane Mottelet wrote:
(num,den) calling style seems to be OK:
--> [a,b]=simp(prod(q.num),prod(q.den))
a =
(num,den) calling style seems to be OK:
--> [a,b]=simp(prod(q.num),prod(q.den))
a =
2 3 4
3.432D-17 +1.230D-16s +3.079D-16s +5.709D-17s +3.432D-17s
b =
2 3 4
1.004
Hello Frederico,
The problem is in simp() :
---> rlist(prod(a.num),prod(a.den),a.dt)
ans =
2 3 4
3.432D-17 + 1.230D-16s + 3.079D-16s + 5.709D-17s + 3.432D-17s
--
* Federico
Miyara
*Envoyé :* mardi 17 mars 2020 10:31
*À :* Users mailing list for Scilab
*Objet :* [Scilab-users] Strange behaviour of prod on rationals
Dear all,
Look at this code (the coefficients are actually the result of pevious
calculations):
NUM = [5.858D-09 + 2.011D-08*%s + 4.884D
Here are examples of my process in Open loop (FTBO) or Close loop (FTBF)
Depending of managemat, D can have s14 …
De : Perrichon
Envoyé : mardi 17 mars 2020 10:49
À : 'Users mailing list for Scilab'
Objet : RE: [Scilab-users] Strange behaviour of prod on rationals
Hello
Objet : [Scilab-users] Strange behaviour of prod on rationals
Dear all,
Look at this code (the coefficients are actually the result of pevious
calculations):
NUM = [5.858D-09 + 2.011D-08*%s + 4.884D-08*%s^2 ...
5.858D-09 + 8.796D-10*%s + 7.028D-10*%s^2]
DEN = [0.1199597 + 7.2765093
Dear all,
Look at this code (the coefficients are actually the result of pevious
calculations):
NUM = [5.858D-09 + 2.011D-08*%s + 4.884D-08*%s^2 ...
5.858D-09 + 8.796D-10*%s + 7.028D-10*%s^2]
DEN = [0.1199597 + 7.2765093*%s + %s^2 ...
8.336136 + 7.0282601*%s + %s^2]
q = NUM./DEN