Frequency Response
- Steady-state response to a sinusoidal input for a stable system
- The frequency of input is varied over a range of interests
- Can be determined by experiments without modeling
- Popular in industry
Consider a stable system G(s).
Let the input be
r(t) = R sin ωt
with
The resulting steady-state response is
yss(t) = Y sin(ωt + φ)
where
A: yss(t) is also a sinusoidal where
G(s) = 1s + 2
Suppose r(t) = √2 sin 2t
R = ? ω = ? rad/s
The above results imply that we can compute the steady-state response to any sinusoidal input, as long as we have
|G(iω)|and∠G(iω)
Bode plot is the plot of log10 |G(iω)| and ∠G(iω) for varying log10 ω
Q: why log ?
A: useful to illustrate both of high-frequency and low-frequency responses in a single axis
Derive the magnitude and phase for ω > 0.
G(iω) = g(iω)n(g > 0, ω > 0)
|G(iω)| = g|iω|n = gωn
Find g and n. Then sketch magnitude and phase.
G(s) = τs + 1 (τ > 0)
Complete the table together.
| ωτ | ω | G(iω) = iωτ + 1 | |G(iω)| | ∠G(iω) |
|---|---|---|---|---|
| ≪ 1 | ||||
| ≈ 1 | ||||
| ≫ 1 |
Substitute s = iω: G(iω) = 1 + iωτ
Magnitude: √(real² + imaginary²)
|G(iω)| = √(1² + (ωτ)²)
Phase: tan⁻¹(imaginary / real)
∠G(iω) = tan−1(ωτ/1)
| ωτ | ω | G(iω) | |G(iω)| | ∠G(iω) |
|---|---|---|---|---|
| ≪ 1 | ||||
| = 1 | ||||
| ≫ 1 |
Find the corner. Sketch magnitude and phase.
Standard rules derived from G(iω), summarized for sketching.
Step 0 Rearrange the transfer function G(s) as a product of the following basic terms:
Step 1 Find the corner frequencies 1τ and ωn, and mark them. The horizontal axis for ω is divided into several sections by these corner frequencies.
Step 2 (Low freq. section) Using the factor g sn, draw asymptotes at low frequencies where ω is less than the smallest corner frequency.
| Magnitude plot | a line with slope n, passing through (ω = 1, |G(iω)| = g) |
|---|---|
| Phase plot | 90°n |
Step 3 (Repeat Step 3 and 4 for each section) Draw asymptotes until the next corner frequency using the following table, and repeat this for each corner frequency.
| Magnitude plot | increase the slope by 1 at the corner frequency 1τ |
|---|---|
| Phase plot | increase the phase by 90° at the corner frequency 1τ |
| Magnitude plot | decrease the slope by 1 at the corner frequency 1τ |
|---|---|
| Phase plot | decrease the phase by 90° at the corner frequency 1τ |
| Magnitude plot | increase the slope by 2 at the corner frequency ωn |
|---|---|
| Phase plot | increase the phase by 180° at the corner frequency ωn |
| Magnitude plot | decrease the slope by 2 at the corner frequency ωn |
|---|---|
| Phase plot | decrease the phase by 180° at the corner frequency ωn |
Step 4 Sketch approximate magnitude / phase curves by the following guide.
| Magnitude plot | increase the asymptote value by a factor of 1.4 at the corner frequency 1τ |
|---|---|
| Phase plot | increase (decrease) the asymptote value by 5.7° at ω = 110τ (at ω = 10τ) |
| Magnitude plot | decrease the asymptote value by a factor of 0.7 at the corner frequency 1τ |
|---|---|
| Phase plot | decrease (increase) the asymptote value by 5.7° at ω = 110τ (at ω = 10τ) |
| Magnitude plot | multiply the asymptote value by a factor of 2ζ at the corner frequency ωn |
|---|---|
| Phase plot | the curve is close to the asymptote when ζ is small |
| Magnitude plot | multiply the asymptote value by a factor of 12ζ at the corner frequency ωn |
|---|---|
| Phase plot | the curve is close to the asymptote when ζ is small |
Rewrite the system using the standard factors in the table.
G(s) = 2000(s + 0.5)s(s + 10)(s + 50)
Match each factor. Complete one row at a time.
Start with slope n
and phase 90°n.
At 1/τ: slope +1,
phase +90°.
At 1/τ: slope −1,
phase −90°.
| Our factor | Match to table | Corner (rad/s) | Slope | Phase |
|---|---|---|---|---|
| 2s−1 | ||||
| 1 + 2s | ||||
| 11 + 0.1s | ||||
| 11 + 0.02s |
Start with 2/s. At each corner, change the current slope.
| Frequency | Current slope | Magnitude sketch |
|---|---|---|
| 0 < ω < 0.5 | ||
| 0.5 < ω < 10 | ||
| 10 < ω < 50 | ||
| ω > 50 |
Find the starting slope before drawing.
Start at −90°. Add the phase change from each factor.
| Frequency | Table rule | Phase asymptote |
|---|---|---|
| 0 < ω < 0.5 | ||
| 0.5 < ω < 10 | ||
| 10 < ω < 50 | ||
| ω > 50 |
Use the Phase plot row in the reference table.
r(t) = R sin ωt ⇒ yss(t) = R|G(iω)| sin(ωt + ∠G(iω))
Maximum |G(iω)|.
The height of the peak.
For the prototype second-order system,
G(s) = ωn2s2 + 2ζωns + ωn2
The frequency-domain specifications are given by
| 0 < ζ ≤ 1/√2 ≈ 0.7071 | ζ > 1/√2 | |
|---|---|---|
| Mr | 12ζ√(1 − ζ2) | 1 |
| ωr | ωn√(1 − 2ζ2) | 0 |
BW = ωn [(1 − 2ζ2) + √(4ζ4 − 4ζ2 + 2)]1/2
Resonant peak Mrζ ↑ → Mr ↓ to 1.Independent of ωn.
Bandwidth BWAt fixed ωn: ζ ↑ → BW ↓.At fixed ζ: ωn ↑ → BW ↑.
G(s) = ωn2s2 + 2ζωns + ωn2
ωn > 0, ζ > 0
| ωωn | ω | G(iω) | |G(iω)| | ∠G(iω) |
|---|---|---|---|---|
| ≪ 1 | ||||
| = 1 | ||||
| ≫ 1 |