Skip to main content

Root-locus design: ζ = 0.5, ts = 2 s, plant 1/(s(s+2))

Short answer

Desired closed-loop poles: |z| = e^{−ζωnT} = 0.6703, angle ωdT = 39.69°. Plant G(z) = 0.01758(z+0.8753)/[(z−1)(z−0.6703)] has angle deficiency +51.26°. The lead controller zero cancels the plant pole at 0.6703; pole placed at 0.2543. Gain K = 12.67 gives G_D(z) = 12.67(z−0.6703)/(z−0.2543) with static velocity error constant K_v = 2.80 s⁻¹.

Design a digital controller such that the compensated system should have damping ratio ζ=0.5 and settling time tₛ=2 s using root locus approach. [Figure: unity-feedback loop: R(s) → summing junction → E(s) → sampler (T = 0.2) → digital controller → ZOH → 1/s(s+2) → C(s)]

Step 1: Desired closed-loop poles

From the specifications ζ=0.5\zeta = 0.5 and ts=2t_s = 2 s (2% criterion, ts=4/ζωnt_s = 4/\zeta\omega_n):

ζωn=4ts=2,ωn=4 rad/sωd=ωn1−ζ2=3.464 rad/s\begin{aligned} \zeta\omega_n &= \frac{4}{t_s} = 2, \qquad \omega_n = 4\ \text{rad/s} \\ \omega_d &= \omega_n\sqrt{1-\zeta^2} = 3.464\ \text{rad/s} \end{aligned}

With T=0.2T = 0.2 s, the samples per cycle are ωs/ωd=31.42/3.464≈9\omega_s/\omega_d = 31.42/3.464 \approx 9.

∣z∣=e−ζωnT=e−0.4=0.6703∠z=ωdT=0.6928 rad=39.69∘z1,2=0.6703∠±39.69∘=0.5158±j0.4281\begin{aligned} |z| &= e^{-\zeta\omega_nT} = e^{-0.4} = 0.6703 \\ \angle z &= \omega_dT = 0.6928\ \text{rad} = 39.69^\circ \\ z_{1,2} &= 0.6703\angle \pm 39.69^\circ = 0.5158 \pm j0.4281 \end{aligned}

Step 2: Plant pulse transfer function

G(z)=(1−z−1) Z[1s2(s+2)]=(2T−1+e−2T)z+(1−e−2T−2Te−2T)4(z−1)(z−e−2T)=0.01758(z+0.8753)(z−1)(z−0.6703)\begin{aligned} G(z) &= (1-z^{-1})\,\mathcal{Z}\left[\frac{1}{s^2(s+2)}\right] = \frac{(2T-1+e^{-2T})z + (1-e^{-2T}-2Te^{-2T})}{4(z-1)(z-e^{-2T})} \\ &= \frac{0.01758(z + 0.8753)}{(z-1)(z-0.6703)} \end{aligned}

Step 3: Angle deficiency at z1=0.5158+j0.4281z_1 = 0.5158 + j0.4281

FactorAngle at z1z_1
zero at −0.8753-0.8753+17.11∘+17.11^\circ
pole at 11−138.52∘-138.52^\circ
pole at 0.67030.6703−109.85∘-109.85^\circ
Total ∠G(z1)\angle G(z_1)−231.26∘-231.26^\circ

The angle condition requires −180∘-180^\circ, so the controller must add +51.26∘+51.26^\circ of phase lead.

Step 4: Lead controller

Take GD(z)=Kz−αz−βG_D(z) = K\dfrac{z-\alpha}{z-\beta}. Put the zero on the plant pole, α=0.6703\alpha = 0.6703, to cancel it. This zero adds +109.85∘+109.85^\circ. The pole must then contribute 109.85∘−51.26∘=58.59∘109.85^\circ - 51.26^\circ = 58.59^\circ:

β=0.5158−0.4281tan⁡58.59∘=0.2543\beta = 0.5158 - \frac{0.4281}{\tan 58.59^\circ} = 0.2543

Step 5: Gain from the magnitude condition

After cancellation, the open-loop function is 0.01758K(z+0.8753)(z−1)(z−0.2543)\dfrac{0.01758K(z+0.8753)}{(z-1)(z-0.2543)}, and ∣GDG∣=1|G_DG| = 1 at z1z_1:

K=∣z1−1∣ ∣z1−0.2543∣0.01758 ∣z1+0.8753∣=0.6464×0.50170.01758×1.4555=12.67K = \frac{|z_1-1|\,|z_1-0.2543|}{0.01758\,|z_1+0.8753|} = \frac{0.6464 \times 0.5017}{0.01758 \times 1.4555} = 12.67 GD(z)=12.67 z−0.6703z−0.2543G_D(z) = 12.67\,\frac{z - 0.6703}{z - 0.2543}

Step 6: Check

GD(z)G(z)=0.2228(z+0.8753)(z−1)(z−0.2543)Characteristic eq.:z2−1.0316z+0.4493=0\begin{aligned} G_D(z)G(z) &= \frac{0.2228(z + 0.8753)}{(z-1)(z-0.2543)} \\ \text{Characteristic eq.:}\quad z^2 &- 1.0316z + 0.4493 = 0 \end{aligned}

The roots are z=0.5158±j0.4281z = 0.5158 \pm j0.4281, the desired poles. ✔ This is a second-order closed loop, so the poles are exactly dominant.

Static velocity error constant:

Kv=lim⁡z→1(1−z−1)GDGT=10.2⋅0.2228×1.87530.7457=2.80 s−1K_v = \lim_{z\to1}\frac{(1-z^{-1})G_DG}{T} = \frac{1}{0.2}\cdot\frac{0.2228 \times 1.8753}{0.7457} = 2.80\ \text{s}^{-1}

Answer: GD(z)=12.67z−0.6703z−0.2543G_D(z) = 12.67\dfrac{z-0.6703}{z-0.2543}, a lead compensator, with Kv=2.80 s−1K_v = 2.80\ \text{s}^{-1}.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Digital Control System chapter 4 questions →