April 5, 2014 FAHEM_HW 06
Problem -1- Solution:
a)
Transfer frequency function is :
(
b)
Bode diagram :
clc
close
s=tf('s')
T=8;
G=1/(T*s+1);
bode(G)
Bode Diagram
0
-5
System: G
Frequency (rad/s): 0.125
Magnitude (dB): -3.01
Magnitude (dB)
-10
-15
-20
-25
-30
-35
-40
0
Phase (deg)
-45
System: G
Frequency (rad/s): 0.125
Phase (deg): -45
-90
-3
10
-2
10
-1
10
10
10
Frequency (rad/s)
Page 1 of 10
April 5, 2014 FAHEM_HW 06
)
From the bode diagram the cutoff frequency is (
Recall to the lecture:
Discussion: I can get the
from the figure by moving the cursor on the
intersection point between the 45 deg. And the phase curve. Also this is the
same result when I calculated this value theoretically. Also, at this point the low
and high frequencies are equal.
c)
Bode Diagram
0
-10
Magnitude (dB)
-20
System: G
Frequency (rad/s): 10
Magnitude (dB): -38.1
-30
-40
System: G
Frequency (rad/s): 100
Magnitude (dB): -58.1
-50
-60
-70
Phase (deg)
-80
-30
-60
-90
-1
10
10
10
10
10
Frequency (rad/s)
Form the bode diagram we choose two point in high frequency:
Point -1- has:
Point -2- has:
Slope=
Page 2 of 10
April 5, 2014 FAHEM_HW 06
Theoretically:
The magnitude of first order system is :
| |
| |
| |
(
| |
| |
Then :
| |
| |
d)
Then:
( )
( )
Page 3 of 10
April 5, 2014 FAHEM_HW 06
Bode Diagram
0
G
System: G
Frequency (rad/s): 50.1
Magnitude (dB): -3
-5
Magnitude (dB)
-10
-15
-20
-25
Phase (deg)
-30
0
-45
System: G
Frequency (rad/s): 50.3
Phase (deg): -45
-90
0
10
10
10
10
Frequency (rad/s)
Page 4 of 10
April 5, 2014 FAHEM_HW 06
Problem -2- Solution:
a) Transfer frequency function is :
(
Bode diagram :
clc
close
s=tf('s')
T=5;
G=2/(T*s+1);
bode(G)
Bode Diagram
10
G
5
Magnitude (dB)
System: G
Frequency (rad/s): 0.2
Magnitude (dB): 3.01
-5
-10
-15
-20
-25
Phase (deg)
-30
0
-45
System: G
Frequency (rad/s): 0.2
Phase (deg): -45
-90
-2
10
10
-1
10
10
Frequency (rad/s)
From the figure the system bandwidth is :
( )
Page 5 of 10
April 5, 2014 FAHEM_HW 06
b)
Step Response
2
1.8
1.6
1.4
Amplitude
1.2
1
0.8
0.6
0.4
0.2
0
10
15
20
25
30
35
40
45
50
Time (seconds)
s =
RiseTime: 11.5109
SettlingTime:
SettlingMin:
SettlingMax:
Overshoot:
Undershoot:
Peak:
PeakTime:
19.5490
1.7999
1.9999
0
0
1.9999
50
clc
close
s=tf('s')
T1=5;
t=0:0.001:50;
figure(1)
G1=2/(T1*s+1);
step(G1,t);
x=step(G1,t);
S=stepinfo(x,t,'RiseTimeLimits',[0 0.9])
figure(2)
bode(G1)
hold on
T2=.5;
G2=2/(T*s+1);
bode(G2)
Page 6 of 10
April 5, 2014 FAHEM_HW 06
c)
Then:
( )
( )
Bode Diagram
20
10
Magnitude (dB)
0
-10
System: G1
Frequency (rad/s): 0.2
Magnitude (dB): 3.01
T=0.5
System: G2
Frequency (rad/s): 50.3
Magnitude (dB): 3.01
-20
T=5
-30
-40
-50
-60
0
Phase (deg)
G1
G2
-45
System: G1
Frequency (rad/s): 0.2
Phase (deg): -45
-90
-2
10
10
-1
System: G2
Frequency (rad/s): 50.3
Phase (deg): -45
10
10
10
Frequency (rad/s)
Page 7 of 10
10
April 5, 2014 FAHEM_HW 06
Problem -3- Solution :
a)
( )
b) From the characteristic equation
, we can find :
Bode Diagram
20
Magnitude (dB)
-20
-40
-60
-80
0
G
Phase (deg)
-45
-90
-135
-180
-1
10
10
10
10
Frequency (rad/s)
Page 8 of 10
April 5, 2014 FAHEM_HW 06
c)
Bode Diagram
System: G
Frequency (rad/s): 2.74
Magnitude (dB): 9.17
Magnitude (dB)
10
System: G
Frequency (rad/s): 2.83
Magnitude (dB): 8.93
System: G
Frequency (rad/s): 4.3
Magnitude (dB): -3
0
-5
-10
-15
-45
Phase (deg)
System: G
Frequency (rad/s): 2.74
Phase (deg): -79.4
-90
System: G
Frequency (rad/s): 2.83
Phase (deg): -90
-135
System: G
Frequency (rad/s): 4.3
Phase (deg): -157
-180
10
0.4
10
0.5
10
0.6
10
0.7
Frequency (rad/s)
From the bode diagram
d)
the cutoff frequency from bode diagram is :
Page 9 of 10
April 5, 2014 FAHEM_HW 06
e) From the characteristic equation
calculate :
, and use FRF, we can
And:
Discussion:
The result I got from these two methods (c and e), are very close.
So, the frequency value and amplitude show me the real resonance ( )
accrued before theoretical resonance states, also the system at real
resonance state has bigger value of amplitude then theoretical value. As
a result, the designer should give the
good tension, because the
systems maybe go to failure before it arrives to theoretical value of
resonance.
clc
close
s=tf('s')
z=0.1768;
wn=8^0.5
G=wn^2/(s^2+2*z*wn*s+wn^2);
bode(G)
Page 10 of 10
