[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #34266] Four issues with "residue" and residue
From: |
Bernard Grung |
Subject: |
[Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez" |
Date: |
Sun, 11 Sep 2011 02:19:46 +0000 |
User-agent: |
Mozilla/5.0 (X11; Linux i686; rv:6.0) Gecko/20100101 Firefox/6.0 |
Follow-up Comment #5, bug #34266 (project octave):
Sorry, I forget to remove the "bg_" from the code in the last comment.
Here is a better way -- just run the following in your simulator:
z1=+7.0372976777e+06;
p1=-3.1415926536e+09;
p2=-4.9964813512e+08;
r1=-(1+z1/p1)/(1-p1/p2)/p2/p1 ;
r2=-(1+z1/p2)/(1-p2/p1)/p2/p1 ;
r3=(1+(p2+p1)/p2/p1*z1)/p2/p1 ;
r4=z1/p2/p1 ;
fi= 1.05e+08 ;
T =1/fi ;
p1=exp(p1*T);
p2=exp(p2*T);
p3=1;
p4=1;
p =[p1;p2;p3;p4];
r4=r4*T;
r3=r3-r4;
r =[r1;r2;r3;r4];
m=[ 1 , 1 , 1 , 1 ;
-(2+p2) , -(2+p1) , -(1+p1+p2) , -(p1+p2) ;
+(1+2*p2) , +(1+2*p1) , (p1*p2+p1+p2) , +p1*p2 ;
-p2 , -p1 , -p1*p2 , 0 ];
a1=1.00;
a2=-(2+p1+p2);
a3=1+2*(p1+p2)+p1*p2;
a4=-(p1+p2+2*p1*p2);
a5=p2*p1;
a =[a1,a2,a3,a4,a5];
b=[m*r]';
k=[];
e=[1;1;1;2];
[rx,px,kx,ex] = bg_residuez(b,a);
assert ((abs (rx - r) < 1e-8
&& abs (px - p) < 1e-8
&& isempty (kx)
&& ex == e));
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?34266>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Bernard Grung, 2011/09/10
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Jordi Gutiérrez Hermoso, 2011/09/10
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Bernard Grung, 2011/09/10
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Jordi Gutiérrez Hermoso, 2011/09/10
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Bernard Grung, 2011/09/10
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Jordi Gutiérrez Hermoso, 2011/09/10
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez",
Bernard Grung <=
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Bernard Grung, 2011/09/10
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Bernard Grung, 2011/09/11
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Jordi Gutiérrez Hermoso, 2011/09/11
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Jordi Gutiérrez Hermoso, 2011/09/11
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Jordi Gutiérrez Hermoso, 2011/09/11
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Bernard Grung, 2011/09/14
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Bernard Grung, 2011/09/16
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Jordi Gutiérrez Hermoso, 2011/09/16
- [Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez", Jordi Gutiérrez Hermoso, 2011/09/29