EXPERIMENT-1
OBJECT:
8 ( s + 3) ( s + 4 )
The transfer function of a system is given as, G ( s ) = s s + 2 s2 + 2s + 5 . Determine the poles and zeros
(
)(
)
and show the pole zero configurations in s-plane.
EXPLANATIONS:
8 ( s + 3) ( s + 4 )
8s 2 + 56s + 96
Given transfer function is G ( s ) = s s + 2 s 2 + 2s + 5 = s 4 + 4s3 + 9s 2 + 10s . Zeros are obtained by equating
(
)(
)
numerator equal to zero, and poles are obtained by equating denominator equal to zero.
Zeros: 8 ( s + 3) ( s + 4 ) = 0 s = 3,s = 4
2
Poles: s ( s + 2 ) ( s + 2s + 5) = 0 s = 0,s = 2,s = 1 + j1,s = 1 j1
MATLAB PROGRAMMING:
>> p=[8 56 96];
>> q=[1 4 9 10 0];
>> sys=tf(p,q)
Transfer function:
8 s^2 + 56 s + 96
-------------------------s^4 + 4 s^3 + 9 s^2 + 10 s
>> pzmap(sys)
Pole-Zero Map
2.5
poles
zeros
2
1.5
System: sys
Zero : -4
Damping: 1
Overshoot (%): 0
Frequency (rad/sec): 4
Imaginary Axis
1
0.5
System: sys
Pole : -1 + 2i
Damping: 0.447
Overshoot (%): 20.8
Frequency (rad/sec): 2.24
0
-0.5
System: sys
Pole : -1 - 2i
Damping: 0.447
Overshoot (%): 20.8
Frequency (rad/sec): 2.24
-1
-1.5
-2
-2.5
-4.5
-4
-3.5
-3
-2.5
-2
Real Axis
-1.5
-1
-0.5
EXPERIMENT-2
OBJECT:
A control system is shown in figure. Determine the transfer function.
R(s)
+ -
4/s(s+4)
+ -
s+1.2
s+0.8
EXPLANATIONS:
C(s)
Transfer function is given as, H ( s ) = R(s) = s 2 + 12s + 12 .
MATLAB PROGRAMMING:
>> n1=[4];
>> d1=[1 4 0];
>> printsys(n1,d1)
num/den =
4
--------s^2 + 4 s
>> n2=[1 1.2];
>> d2=[0 1];
>> printsys(n2,d2)
num/den =
s + 1.2
------1
>> n3=[1 0.8];
>> d3=[0 1];
>> printsys(n3,d3)
num/den =
s + 0.8
------1
>> [n4,d4]=cloop(n1,d1,-1)
n4 =
0
C(s)
d4 =
1
>> printsys(n4,d4)
num/den =
4
------------s^2 + 4 s + 4
>> [n5,d5]=parallel(n2,d2,n3,d3)
n5 =
0
d5 =
0
>> printsys(n5,d5)
num/den =
2s+2
------1
>> [n6,d6]=feedback(n4,d4,n5,d5,-1)
n6 =
0
12
12
d6 =
0
>> printsys(n6,d6)
num/den =
4
--------------s^2 + 12 s + 12