Experiment-5
Q1)
b) MATLAB Code:
clc
clear all
close all
% Original System
num=[1 7 24 24];
den=[1 10 35 50 24];
b=tf(num,den);
subplot(521)
step(b)
subplot(522)
impulse(b)
% 2nd order reduced system
n1=[0.2925 1];
d1=[0.4001 1.3758 1];
b1=tf(n1,d1);
subplot(525)
step(b1)
subplot(526)
impulse(b1)
% 3rd order reduced system
n2=[0.2133 0.59 1];
d2=[0.2558 0.9356 1.6733 1];
b2=tf(n2,d2);
subplot(5,2,9)
step(b2)
subplot(5,2,10)
impulse(b2)
Output:
c)
MATLAB Code:
clc
clear all
close all
% Original System
num=[1 7 24 24];
den=[1 10 35 50 24];
b=tf(num,den);
subplot(221)
bode(b)
% 2nd order reduced system
n1=[0.2925 1];
d1=[0.4001 1.3758 1];
b1=tf(n1,d1);
subplot(222)
bode(b1)
% 3rd order reduced system
n2=[0.2133 0.59 1];
d2=[0.2558 0.9356 1.6733 1];
b2=tf(n2,d2);
subplot(223)
bode(b2)
Output:
d)
MATLAB Code:
clc
clear all
close all
%Original System
num=[1 7 24 24];
den=[1 10 35 50 24];
b=tf(num,den)
pole(b)
% 2nd order reduced system
n1=[0.2925 1];
d1=[0.4001 1.3758 1];
b1=tf(n1,d1)
pole(b1)
% 3rd order reduced system
n2=[0.2133 0.59 1];
d2=[0.2558 0.9356 1.6733 1];
b2=tf(n2,d2)
pole(b2)
Output:
b =
s^3 + 7 s^2 + 24 s + 24
---------------------------------
s^4 + 10 s^3 + 35 s^2 + 50 s + 24
Continuous-time transfer function.
ans =
-4.0000
-3.0000
-2.0000
-1.0000
b1 =
0.2925 s + 1
------------------------
0.4001 s^2 + 1.376 s + 1
Continuous-time transfer function.
ans =
-2.3951
-1.0435
b2 =
0.2133 s^2 + 0.59 s + 1
-------------------------------------
0.2558 s^3 + 0.9356 s^2 + 1.673 s + 1
Continuous-time transfer function.
ans =
-1.3230 + 1.4542i
-1.3230 - 1.4542i
-1.0115 + 0.0000i