blob: c1ab83342244a4c22258cf30b5bf4f2e1f8a3523 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
//Example 3.14
//Partial fraction expansion for distinct real roots
clear;
clc;
//------------------------------------------------------------------
// Transfer function
s=%s;
num=2;
p1=(s+1);
p2=(s+2);
p3=(s+4);
sys=syslin('c',num/(p1*p2*p3))
//Partial fraction expansion is: sys= r1/p1 + r2/p2 + r3/p3
//residue calculation
r1=residu(num,p1,(p2*p3))
r2=residu(num,p2,(p1*p3))
r3=residu(num,p3,(p1*p2))
disp([r1 r2 r3]',"Residues of the poles p1, p2 and p3 are")
disp([roots(p1), roots(p2), roots(p3)]',"Poles p1, p2 and p3 are at")
disp('k=[]')
//------------------------------------------------------------------
|
