Skip to main content

Controller with integral action, plant 10/((s+1)(s+5))

Short answer

Digital controller: GD(z)=2.204(z-0.8187)(z-0.3679)/[(z-1)(z-0.2009)] with integral action. Static velocity error constant Kv=3.16 s^-1.

Design a suitable digital controller that includes an integral control action. The design specifications are that the damping ratio of the closed loop poles be 0.5 and there will be at least eight samples per cycle of the damped sinusoidal oscillations. The sampling time is assumed to be 0.2 sec. [Figure: unity-feedback loop: R(z) → summing junction → sampler (T) → digital controller → (1-e⁻ᵀˢ)/s → 10/((s+1)(s+5)) → C(z)]

Step 1: Plant with ZOH (T=0.2T = 0.2 s)

10s(s+1)(s+5)=2s−2.5s+1+0.5s+5G(z)=(1−z−1)[2zz−1−2.5zz−e−0.2+0.5zz−e−1]=0.1371z+0.0921(z−0.8187)(z−0.3679)=0.1371(z+0.6714)(z−0.8187)(z−0.3679)\begin{aligned} \frac{10}{s(s+1)(s+5)} &= \frac{2}{s} - \frac{2.5}{s+1} + \frac{0.5}{s+5} \\ G(z) &= (1-z^{-1})\left[\frac{2z}{z-1} - \frac{2.5z}{z-e^{-0.2}} + \frac{0.5z}{z-e^{-1}}\right] \\ &= \frac{0.1371z + 0.0921}{(z-0.8187)(z-0.3679)} = \frac{0.1371(z + 0.6714)}{(z-0.8187)(z-0.3679)} \end{aligned}

Step 2: Desired poles

With ζ=0.5\zeta = 0.5 and 8 samples per cycle:

∣z∣=exp⁡(−2π(0.5)0.866⋅18)=0.6354,∠z=45∘,z1=0.4493+j0.4493|z| = \exp\left(-\frac{2\pi(0.5)}{0.866}\cdot\frac18\right) = 0.6354, \quad \angle z = 45^\circ, \quad z_1 = 0.4493 + j0.4493

Step 3: Try PI alone

Integral action requires a controller pole at z=1z = 1. With a PI controller K(z−0.8187)/(z−1)K(z-0.8187)/(z-1) cancelling the slow plant pole, the angle at z1z_1 is:

+21.85∘⏟zero −0.6714−140.79∘⏟pole 1−79.73∘⏟pole 0.3679=−198.67∘\underbrace{+21.85^\circ}_{\text{zero } -0.6714} \underbrace{-140.79^\circ}_{\text{pole } 1} \underbrace{-79.73^\circ}_{\text{pole } 0.3679} = -198.67^\circ

This is not −180∘-180^\circ. There is 18.67∘18.67^\circ too much lag, so PI alone cannot place the poles. A lead part is added, which makes the controller PID-like.

Step 4: Controller with integral + lead action

GD(z)=K(z−0.8187)(z−0.3679)(z−1)(z−β)G_D(z) = K\frac{(z - 0.8187)(z - 0.3679)}{(z-1)(z-\beta)}

Both plant poles are cancelled. The angle condition at z1z_1 is:

21.85∘−140.79∘−∠(z1−β)=−180∘  ⇒  ∠(z1−β)=61.06∘21.85^\circ - 140.79^\circ - \angle(z_1 - \beta) = -180^\circ \;\Rightarrow\; \angle(z_1-\beta) = 61.06^\circ β=0.4493−0.4493tan⁡61.06∘=0.2009\beta = 0.4493 - \frac{0.4493}{\tan 61.06^\circ} = 0.2009

Step 5: Gain

K=∣z1−1∣ ∣z1−0.2009∣0.1371 ∣z1+0.6714∣=0.7107×0.51340.1371×1.2074=2.204K = \frac{|z_1-1|\,|z_1-0.2009|}{0.1371\,|z_1+0.6714|} = \frac{0.7107 \times 0.5134}{0.1371 \times 1.2074} = 2.204 GD(z)=2.204 (z−0.8187)(z−0.3679)(z−1)(z−0.2009)G_D(z) = 2.204\,\frac{(z-0.8187)(z-0.3679)}{(z-1)(z-0.2009)}

Step 6: Check

GDG=0.3022(z+0.6714)(z−1)(z−0.2009)G_DG = \frac{0.3022(z + 0.6714)}{(z-1)(z-0.2009)}

The characteristic equation is z2−0.8986z+0.4038=0z^2 - 0.8986z + 0.4038 = 0, with roots z=0.4493±j0.4493z = 0.4493 \pm j0.4493. ✔ That is ζ=0.5\zeta = 0.5 with exactly 8 samples per cycle.

Kv=1T⋅0.3022×1.67141−0.2009=3.16 s−1K_v = \frac{1}{T}\cdot\frac{0.3022 \times 1.6714}{1 - 0.2009} = 3.16\ \text{s}^{-1}

The pole at z=1z = 1 gives zero steady-state error for a step input.

Answer: GD(z)=2.204(z−0.8187)(z−0.3679)(z−1)(z−0.2009)G_D(z) = 2.204\dfrac{(z-0.8187)(z-0.3679)}{(z-1)(z-0.2009)}, with Kv=3.16 s−1K_v = 3.16\ \text{s}^{-1}.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Digital Control System chapter 4 questions →