Skip to main content

Bode design: PM 55°, GM 10 dB, Kv = 5

Short answer

GD(z) = 40.78(z-0.8187)/(z-0.2608). This lead compensator achieves a static velocity error constant Kv = 5 s⁻¹, phase margin ≈55°, and gain margin ≈12.3 dB.

Design a digital controller for the system shown below. Use the Bode diagram approach in s plane. The design specifications are that the phase margin be 55°, the gain margin be at least 10 dB and the static velocity error constant be 5 sec⁻¹. The sampling period is specified as 0.1 sec. [Figure: unity-feedback loop: R(z) → summing junction → sampler → digital controller → ZOH → plant 1/s(s+2) → C(z)]

The design is done in the w-plane (bilinear transformation), where the Bode diagram can be used as in the s-plane. A lead compensator GD(w)=K1+τw1+ατwG_D(w) = K\dfrac{1+\tau w}{1+\alpha\tau w} is used.

Step 1: Pulse transfer function of the plant

With T=0.1T = 0.1 s and e−0.2=0.8187e^{-0.2} = 0.8187:

G(z)=(1−z−1) Z[1s2(s+2)]=0.004683z+0.004381(z−1)(z−0.8187)=0.004683(z+0.9355)(z−1)(z−0.8187)\begin{aligned} G(z) &= (1-z^{-1})\,\mathcal{Z}\left[\frac{1}{s^2(s+2)}\right]\\ &= \frac{0.004683z + 0.004381}{(z-1)(z-0.8187)} = \frac{0.004683(z+0.9355)}{(z-1)(z-0.8187)} \end{aligned}

Step 2: Transform to the w-plane

Substitute z=1+0.05w1−0.05wz = \dfrac{1+0.05w}{1-0.05w}:

G(w)=0.5(1−w20)(1+w600.4)w (1+0.5017w)G(w) = \frac{0.5\left(1-\frac{w}{20}\right)\left(1+\frac{w}{600.4}\right)}{w\,(1+0.5017w)}

Step 3: Gain for Kv=5K_v = 5

Kv=lim⁡w→0w GD(w)G(w)=K×0.5=5  ⇒  K=10K_v = \lim_{w\to0} w\,G_D(w)G(w) = K\times0.5 = 5 \;\Rightarrow\; K = 10

Step 4: Uncompensated Bode data (10 G(jν)10\,G(j\nu))

ν\nu (rad/s)12461020
Magnitude (dB)13.05.0−4.9−11.2−19.2−29.1
Phase (deg)−119.4−140.6−164.4−177.8−194.3−217.4

Gain crossover is at ν=2.88\nu = 2.88 rad/s and the phase margin is only 26.8°. About 28° or more of lead is needed. The w-plane zero at w=20w=20 adds extra lag at high frequency, so a larger lead is needed than in a continuous design.

Step 5: Lead compensator

Choose the compensator zero to cancel the plant pole at w=−1.993w = -1.993, so τ=0.5017\tau = 0.5017. Then reduce α\alpha until PM = 55°:

α\alphaNew crossover νc\nu_cPMGM
0.254.4748.6°12.3 dB
0.204.6552.3°12.3 dB
0.174.7655.0°12.3 dB
GD(w)=10 1+0.5017w1+0.0853wG_D(w) = 10\,\frac{1 + 0.5017w}{1 + 0.0853w}

Compensated loop: GD(w)G(w)=5(1−w/20)(1+w/600.4)w(1+0.0853w)G_D(w)G(w) = \dfrac{5(1-w/20)(1+w/600.4)}{w(1+0.0853w)}. Phase margin = 55°, phase crossover at ν=15.7\nu = 15.7 rad/s, gain margin = 12.3 dB (≥ 10 dB).

Step 6: Back to the z-plane

Substitute w=2Tz−1z+1=20z−1z+1w = \dfrac{2}{T}\dfrac{z-1}{z+1} = 20\dfrac{z-1}{z+1}:

GD(z)=10 1+10.033z−1z+11+1.706z−1z+1=40.78 z−0.8187z−0.2608\begin{aligned} G_D(z) &= 10\,\frac{1+10.033\frac{z-1}{z+1}}{1+1.706\frac{z-1}{z+1}} = 40.78\,\frac{z-0.8187}{z-0.2608} \end{aligned}

Check of KvK_v

Kv=1Tlim⁡z→1(z−1)GD(z)G(z)=10.1×40.78×0.18130.7392×0.0090630.1813=5.0 s−1K_v = \frac{1}{T}\lim_{z\to1}(z-1)G_D(z)G(z) = \frac{1}{0.1}\times40.78\times\frac{0.1813}{0.7392}\times\frac{0.009063}{0.1813} = 5.0\ \text{s}^{-1}

Answer: GD(z)=40.78 z−0.8187z−0.2608G_D(z) = 40.78\,\dfrac{z-0.8187}{z-0.2608}, giving PM ≈ 55°, GM ≈ 12.3 dB and Kv=5 s−1K_v = 5\ \text{s}^{-1}.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Digital Control System chapter 4 questions →