Skip to main content

Controller design for peak time ≤ 1.2 s (ζ = 0.5)

Short answer

Using lead compensation with zero cancelling the plant pole at 0.8187, the controller is G_D(z)=7.340(z-0.8187)/(z-0.4177). The resulting closed-loop dominant poles are 0.6401±j0.3696, giving damping ratio 0.5 and peak time 1.2 s.

Design a digital controller for the system shown in figure below so that dominant closed loop poles of the system will have damping ratio of 0.5. It is required that the peak time of the response tₚ≤ 1.2 sec. Assume sampling time T = 0.2 sec. [Figure: unity-feedback loop: R(z) → summing junction → sampler → digital controller → m*(kT) → ZOH → 1/s(s+1) → C(z)]

A lead-type digital controller GD(z)=Kz−αz−βG_D(z) = K\dfrac{z-\alpha}{z-\beta} is designed by the z-plane root-locus method. Its zero cancels the slow plant pole.

Plant pulse transfer function (T=0.2T = 0.2 s)

With e−0.2=0.8187e^{-0.2} = 0.8187:

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

Desired dominant poles

tp=πωd=1.2⇒ωd=2.618 rad/sωn=ωd1−0.52=3.023 rad/s,ζωn=1.511∣z∣=e−ζωnT=e−0.3023=0.7391,∠z=ωdT=0.5236 rad=30∘z1=0.6401+j0.3696\begin{aligned} t_p &= \frac{\pi}{\omega_d} = 1.2 \Rightarrow \omega_d = 2.618\ \text{rad/s}\\ \omega_n &= \frac{\omega_d}{\sqrt{1-0.5^2}} = 3.023\ \text{rad/s},\quad \zeta\omega_n = 1.511\\ |z| &= e^{-\zeta\omega_nT} = e^{-0.3023} = 0.7391, \quad \angle z = \omega_dT = 0.5236\ \text{rad} = 30^\circ\\ z_1 &= 0.6401 + j0.3696 \end{aligned}

That is 12 samples per cycle of damped oscillation, which is adequate.

Angle deficiency

FactorAngle at z1z_1
zero at −0.9355-0.9355+13.20∘+13.20^\circ
pole at 11−134.24∘-134.24^\circ
pole at 0.81870.8187−115.80∘-115.80^\circ
Total−236.84∘-236.84^\circ

The controller must add +56.84∘+56.84^\circ.

Controller zero and pole

  • Zero at α=0.8187\alpha = 0.8187 cancels the plant pole and removes its −115.80∘-115.80^\circ.
  • Then the controller pole must give 13.20∘−134.24∘−θβ=−180∘13.20^\circ - 134.24^\circ - \theta_\beta = -180^\circ, so θβ=58.96∘\theta_\beta = 58.96^\circ:
β=0.6401−0.3696tan⁡58.96∘=0.4177\beta = 0.6401 - \frac{0.3696}{\tan 58.96^\circ} = 0.4177

Gain

K=∣(z1−1)(z1−0.4177)0.01873(z1+0.9355)∣=7.340K = \left|\frac{(z_1-1)(z_1-0.4177)}{0.01873(z_1+0.9355)}\right| = 7.340 GD(z)=7.340 z−0.8187z−0.4177G_D(z) = 7.340\,\frac{z-0.8187}{z-0.4177}

Check

  • Characteristic equation: z2−1.2802z+0.5463=0z^2 - 1.2802z + 0.5463 = 0, roots 0.6401±j0.36960.6401 \pm j0.3696, which is the design point.
  • Step response: c(kT)=0, 0.137, 0.442, 0.757, 0.994, 1.125, 1.163, 1.141, …c(kT) = 0,\ 0.137,\ 0.442,\ 0.757,\ 0.994,\ 1.125,\ 1.163,\ 1.141,\ \ldots. The peak is 1.163 at t=6T=1.2t = 6T = 1.2 s, so tp≤1.2t_p \le 1.2 s is met.
  • Velocity error constant:
Kv=1Tlim⁡z→1(z−1)GD(z)G(z)=7.340×0.036250.2×(1−0.4177)=2.28 s−1K_v = \frac{1}{T}\lim_{z\to1}(z-1)G_D(z)G(z) = \frac{7.340\times0.03625}{0.2\times(1-0.4177)} = 2.28\ \text{s}^{-1}

Answer: GD(z)=7.340 z−0.8187z−0.4177G_D(z) = 7.340\,\dfrac{z-0.8187}{z-0.4177}, giving closed-loop poles 0.6401±j0.36960.6401\pm j0.3696 (ζ=0.5\zeta=0.5, tp=1.2t_p=1.2 s).

✓ Answer checked by Claude Opus 5.5 (high effort)

All Digital Control System chapter 4 questions →