MATLAB实现部分分式展开式.doc

MATLAB实现部分分式展开式
1. 用部分分式法求F（s）的逆变换
>> format rat
>> num=[1 2];
>> den=[1 4 3 0];
>> [r,p]=residue(num,den)

r =

-1/6
-1/2
2/3

p =

-3
-1
0

>>
2. 用部分分式法求F（s）的逆变换

>> num=[1 -2];
>> a=conv([1 0],[1 1]);b=conv([1 1],[1 1]);
>> den=conv(a,b);
>> [r,p]=residue(num,den)

r =

2
2
3
-2

p =

-1
-1
-1
0
Or
>> num=[1 -2];
>> den=poly([0 -1 -1 -1]);
>> [r,p]=residue(num,den)

r =

2
2
3
-2

p =

-1
-1
-1
0

>>

3. 用部分分式法求F（s）的逆变换
>> num=[2 3 0 5];
>> den=conv([1 1],[1 1 2]);
>> [r,p,k]=residue(num,den)

r =

-2 + 2024/1785i
-2 - 2024/1785i
3

p =

-1/2 + 1012/765i
-1/2 - 1012/765i
-1

k =

2
Or
>> num=[2 3 0 5];
>> den=conv([1 1],[1 1 2]);
>> [r,p,k]=residue(num,den);
>> [angle,mag]=cart2pol(real(r),imag(r))

angle =

1972/751
-1972/751
0

mag =

7895/3434
7895/3434
3

>>