求取留数（续） +5 [r,p,k]=residue(num,den A(s)s-P1 +++k) s-P2 S-P 对于重根部分 + + s-Pi(s-P:)2 +…+6-p,严 B(S)_s2+2s+3 r= 1 0 2 A阿s+3s+3s+打s++s++s++0 1.0000 0.0000 2.0000 p= -1.0000 num=[123小： -1.0000 den=[1331, -1.0000 [r.p,k]=residue(num,den) k= School of Mechanical Engineering ME369-Lecture 4.5 Shanghai Jiao Tong University Fall 2015 求取留数（续） B=5+5++5-+k6) [r,p,k]=residue(num,den) A(s)s-P1 s-P2 s-Pn 对于重根部分 +…+ s-Pi (s-Pi) (s-p:)m [num,den]=residue(r,p,k) B(9=1+ 0 2 52+25+3 A)s+1(s++(s+厅 s2+3s2+3s+1 =102 num p=[-1-1-1]: 123 k=0: den= [num,den]=residue(r,p,k) 1331 School of Mechanical Engineering ME369-Lecture 4.5 Shanghai Jiao Tong University Fall 2015 44 ME369-Lecture 4.5 Fall 2015 School of Mechanical Engineering Shanghai Jiao Tong University 0 (s 1) 2 (s 1) 0 s 1 1 s 3s 3s 1 s 2s 3 A(s) B(s) 3 2 2 3 2              m j j m 2 j j 1 j j (s p ) r (s p ) r s p r          k(s) s p r s p r s p r A(s) B(s) n n 2 2 1 1          [r,p,k]=residue(num,den) num=[ 1 2 3]; den=[1 3 3 1]; [r,p,k]=residue(num,den) r = 1.0000 0.0000 2.0000 p = -1.0000 -1.0000 -1.0000 k = [] 求取留数（续） 对于重根部分 ME369-Lecture 4.5 Fall 2015 School of Mechanical Engineering Shanghai Jiao Tong University [num,den]=residue(r,p,k) m j j m 2 j j 1 j j (s p ) r (s p ) r s p r          对于重根部分 k(s) s p r s p r s p r A(s) B(s) n n 2 2 1 1          [r,p,k]=residue(num,den) 2 2 3 3 2 ( ) 1 0 2 2 3 ( ) 1 ( 1) ( 1) 3 3 1 B s s s A s s s s s s s             r=[ 1 0 2]; p=[-1 -1 -1]; k=[]; [num,den]=residue(r,p,k) num = 1 2 3 den = 1 3 3 1 求取留数（续）