Part 2-4: Frequency-Response Analysis: T U T U
Part 2-4: Frequency-Response Analysis: T U T U
• The nature of the linear system [i.e., the TF G(s) ] will determine the
amplitude ratio (Y0/U0) and phase shift (f ), and these quantities will
be a function of the input frequency
Imaginary
sin wt
1
wt Real
0
cos wt
Result:
(( jw ) 3
+ 4( jw ) 2 + 12 jw + 26)Y ( jw )e jwt = (2 jw + 3)U ( jw )e jwt
Frequency Transfer Function (4)
• Repeat last equation:
( ( jω )3
+ 4( jω ) 2
+12 jω + 26 ) Y ( jω )e jω t
= ( 2 jω + 3) U( jω )e jω t
Y ( jw ) 2 jw + 3
= = G ( jw )
U ( jw ) ( jw ) + 4( jw ) + 12 jw + 26
3 2
Y ( jw ) 2 jw + 3
= = G ( jw )
U ( jw ) ( jw ) + 4( jw ) + 12 jw + 26
3 2
Polar form G ( jw ) = G ( jw ) e jf
G(jw)
Im
f
0 Re Real axis
G ( jw ) = Re 2 + Im 2 Magnitude of G(jw)
⎛ Im ⎞
−1
ϕ ( jω ) = arg[G( jω )] = tan ⎜ ⎟ Angle (or argument) of G(jw)
⎝ Re ⎠
Derivation of Frequency Response
• Next, we wish to derive an equation for the frequency response,
which is the steady-state output of a system driven by a sinusoidal
input
• However, we can sub for the sinusoidal TF using its polar form
G ( j w ) = G ( j w ) e jf
Derivation of Frequency Response (2)
yss (t ) = G ( jw ) U 0 cos(wt + f )
ϕ ( jω ) = arg[G( jω )] = phase shift between the output sinusoid and input sinusoid
Frequency Response:
Example
= 1.2 W
RL system:
I (s) 1 1 current
G (s) = = = =
Ein ( s ) Ls + R 0.02 s + 1.2 source voltage
w = 50 rad/s
Frequency Response:
Example (2)
Frequency transfer function: G ( jw ) = 1
0.02 jw + 1.2
12 + 0 2
Magnitude: G ( jw ) =
1.2 2 + (0.02w )
2
⎛0⎞ ⎛ 0.02ω ⎞
ϕ ( jω ) = arg [G( jω )] = tan −1 ⎜ ⎟ − tan −1 ⎜ ⎟ à f = -0.6947 rad
⎝1⎠ ⎝ 1.2 ⎠
Frequency Response:
Example (3)
Put it all together to obtain the frequency response:
Set amplitude = 2
and frequency = 50 rad/s
Frequency Response:
Example (2)
20 log(0.124) = -18.13 dB
6
G ( jw ) =
jw + 4
62 + 02
Magnitude G ( jw ) =
42 + w 2
−1 ⎛0⎞ −1 ⎛ ω ⎞
Phase angle ϕ ( jω ) = arg [G( jω )] = tan ⎜ ⎟ − tan ⎜ ⎟
⎝6⎠ ⎝4⎠
Bode Diagram: 1st-Order System (3)
DC gain
Bode diagram
1st-Order System Bode Diagram (5)
Note the 2 asymptotes for the Magnitude Bode plot:
“break” or “corner” frequency = 4 rad/s
10 Low-frequency
G (j w ), dB
0
asymptote (flat
line at DC gain)
-10
Magnitude of
-20
High-frequency
-30
10
-1
10
0
10
1
10
2 asymptote:
Frequency, w , rad/s -20 dB/decade
0
Phase of G (j w ), deg
-15
-30
-45
-60
-75
-90
-1 0 1 2
10 10 10 10
Frequency, w , rad/s
Bode Diagram via MATLAB
• True Bode diagram of the first-order system (from MATLAB):
>> numG = 6;
>> denG = [ 1 4 ];
Bode Diagram
>> sysG = tf(numG,denG) 5
>> bode(sysG) 0
Magnitude (dB)
-5
-10
-15
6
-25
G (s) =
0
s+4
Phase (deg)
-45
-90
-1 0 1 2
10 10 10 10
Frequency (rad/sec)
How to Use the Bode Diagram
6
• Suppose the input is u(t) = 7 sin 2t (recall that G ( s) = )
s+4
Magnitude = 2.5527 dB
Bode Diagram Amplitude ratio is
5
0 G ( j 2) = 2.5527 dB
Magnitude (dB)
-5
-20 w = 2 rad/s
-25 Frequency response
0
yss = 9.392 sin(2t – 0.464)
Phase = -26.57 deg
Phase (deg)
-45
7 x 1.3416 Phase
(in rad)
-90
-1 0 1 2
10 10 10 10
Frequency (rad/sec)
Bode Diagram:
Example
• Consider the simple RC circuit driven by a sinusoidal voltage
Problem: Use the Bode diagram to find input frequency w where the
magnitude of the output voltage eC drops to 1.2 V (one-half the input)
Example (2)
Parameter values:
R = 4 W and C = 0.003 F
EC ( s ) 1 Time constant
à Transfer function is G ( s ) = = t = RC = 0.012 s
Ein ( s ) RCs + 1
à The Bode diagram will show a corner frequency at wc = 1/t = 83.333 rad/s
(note DC gain = 1)
Example (3)
EC ( s ) 1 Time constant
Transfer function G (s) = =
Ein ( s ) RCs + 1 t = RC = 0.012 s
à We simply read the -6 dB drop in magnitude from the Bode diagram (next)
Low-Pass Filter: Bode Diagram
Example (4)
10
Magnitude, dB
-10 ½ magnitude:
-20
6 dB
G ( jw ) = 0.5
w = 144 rad/s
-30
0 1 2 3 or
10 10 10 10
Frequency, w , rad/s G ( jw ) dB = 20 log10 (0.5) = -6.02 dB
0
-15
Therefore, output voltage
Phase, deg
-30
-45 eC is reduced to one-half
-60 mag. at input frequency
-75 w = 144 rad/s
-90
0 1 2 3
(or, 23 Hz)
10 10 10 10
Frequency, w , rad/s
RC circuit is a
Corner frequency wc = 83.33 rad/s “low-pass filter”
Low-Pass Filter: Noisy Input
Example
• Consider again the RC circuit driven by a sinusoidal voltage
Note: voltage signal eS is the “clean” desired signal (w = 10 rad/s or 1.6 Hz)
The filter output voltage eC nearly matches the desired signal (noise is filtered out)
Kw n2
G (s) = 2
s + 2zw n s + w n2
100
G (s) = 2
s + 4 s + 100
Low-frequency High-frequency
asymptote asymptote
(DC gain = 1) -40 dB/decade
Asymptotes intersect at
corner frequency wn = 10 rad/s
100
G (s) = 2
s + 4 s + 100 High-frequency
phase à -180 deg
Bode Diagram of 2nd-Order System (3)
100
G (s) = 2
s + 20zs + 100
1. When z = 0, wr = wn
w r = w n 1 - 2z 2 Resonant frequency
2. Resonant frequency is
always < wn
3. If z > 0.7071, there is no
resonant peak
Bandwidth
• Cutoff frequency (wB ) is defined as the maximum input
frequency at which the output of the system will track the input
sinusoid in a satisfactory manner
15 N input 6.4 mm
force poppet
response
Note that the steady-state amplitude of the poppet valve is about 6.4 mm, as predicted
by the Bode diagram; also, the output is essentially “in phase” with the input. The Bode
diagram (previous page) shows a very small phase lag of 3 deg when w1 = 20 rad/s.
Frequency Response of a Poppet Valve:
Resonant Input w2 = 148.6 rad/s (23.7 Hz)
Consider again the Bode diagram of the
poppet valve, where the input force is
f = 15sin148.6t N. Therefore, the input
frequency is the resonant frequency
-59 dB
-78.2 deg
15.9 mm
15 N input
poppet
force
response
Note that the steady-state amplitude of the poppet valve is about 15.9 mm, as predicted
by the Bode diagram; also, the output lags the input. The phase lag translates into a
time lag = (1.36 rad)/(148.6 rad/s) = 0.0092 sec, as shown on the time response
Frequency Response of a Poppet Valve:
High-Frequency Input w3 = 600 rad/s (95.5 Hz)
Finally, assume the input force is f = 15sin600t N. Therefore, the
input (forcing) frequency is the high frequency, w3 = 600 rad/s (95.5 Hz)
w3 = 600 rad/s
xss(t) = 0.44sin(600t – 3.031) mm
Frequency response
-173.7 deg
Frequency Response of a Poppet Valve:
High-Frequency Input w3 = 600 rad/s
Simulink response of input force at high frequency, f = 15sin600t N
Note that the steady-state amplitude of the poppet valve is about 0.44 mm, as
predicted by the Bode diagram; also, the output is ~180 deg out-of-phase with the
input. The phase lag translates into a time lag = (3.03 rad)/(600 rad/s) = 0.005 sec,
as shown on the time response
Example: Vibration Isolation
• Vibrations in mechanical systems are typically an
undesirable phenomena (transmit forces, generate noise)
(a) (b)
X (s) bs + k (b / m) s + (k / m)
G(s) = = 2 = 2
X b ( s ) ms + bs + k s + (b / m) s + (k / m)
Transmissibility (3)
• The transmitted displacement will depend on the magnitude of
transfer function G(s)
(b / m) jw + (k / m)
G ( jw ) =
- w 2 + (b / m) jw + (k / m)
• Substitute k / m = w 2
n , b / m = 2zwn and after some algebra, we get
1 + j (2zw / w n )
G ( jw ) =
1 - (w 2 / w n2 ) + j (2zw / w n )
Transmissibility (4)
• Finally, we can compute the amplitude ratio
12 + (2zw / wn ) 2
G ( jw ) =
[1 - (w 2
]2
/ wn2 ) + (2zw / wn ) 2
x 12 + (2zb ) 2
TR = = G ( jw ) =
xb [1 - b ] + (2zb )
2 2 2
Transmissibility
Transmissibility vs. b
Notes:
1. TR depends on damping
ratio z and frequency ratio b
2. When z » 0, TR ® ¥ when
b = 1 (or, input is resonance)
b = w / wn
Transmissibility vs. b (2)
Notes:
4. When b is very small (stiff
spring), TR ~ 1 for all z
b = w / wn
TR < 1 occurs only when b > 2
Transmissibility: Summary
1. For minimum transmitted
force or motion, choose
vibration isolation mounts
with undamped natural
frequency (i.e., stiffness)
much lower than the
expected input vibration
frequency (therefore, b = w / w n
is high).
b = w / wn
Transmissibility: Summary (2)
3. Many times the disturbing
input frequency may not be
constant but may be in a
known range. The isolator
design must consider a
range of b
b = w / wn