Skip to main content

z-plane design with Kv, plant 1/(s(s+2))

Short answer

Desired closed-loop poles are z=0.5158±j0.4281 with ζ=0.5 and ts=2 s. Plant with ZOH: G(z)=0.01758(z+0.8753)/[(z-1)(z-0.6703)]. Controller: G_D(z)=12.67(z-0.6703)/(z-0.2543). Static velocity error constant Kv=2.80 s^-1.

Consider the digital control system shown in figure below. In the z plane, design a digital controller such that the dominant closed-loop poles have a damping ratio of 0.5 and a settling time of 2 seconds. The sampling period is assumed to be 0.2 sec. Also, obtain the static velocity error constant of the system. [Figure: unity-feedback loop: r(t) → summing junction → e(t) → sampler → e(kT) → digital controller → ZOH → 1/s(s+2) → C(t)]

Step 1: Desired closed-loop poles

With ζ=0.5\zeta = 0.5 and ts=4/(ζωn)=2t_s = 4/(\zeta\omega_n) = 2 s:

ζωn=2,ωn=4,ωd=3.464 rad/s\zeta\omega_n = 2, \quad \omega_n = 4, \quad \omega_d = 3.464\ \text{rad/s}

With T=0.2T = 0.2 s:

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

Step 2: ZOH + plant 1/[s(s+2)]1/[s(s+2)]

G(z)=(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)G(z) = \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)}

Step 3: Angle deficiency

At z1z_1: the zero gives +17.11∘+17.11^\circ, the pole at 1 gives −138.52∘-138.52^\circ, and the pole at 0.6703 gives −109.85∘-109.85^\circ. The total is −231.26∘-231.26^\circ. The controller must supply +51.26∘+51.26^\circ (lead).

Step 4: Lead controller GD(z)=Kz−αz−βG_D(z) = K\dfrac{z-\alpha}{z-\beta}

Choose α=0.6703\alpha = 0.6703 to cancel the plant pole. This contributes +109.85∘+109.85^\circ. The controller 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

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

GDG=0.2228(z+0.8753)(z−1)(z−0.2543)G_DG = \frac{0.2228(z + 0.8753)}{(z-1)(z-0.2543)}

The characteristic equation z2−1.0316z+0.4493=0z^2 - 1.0316z + 0.4493 = 0 gives z=0.5158±j0.4281z = 0.5158 \pm j0.4281. ✔

Step 7: Static velocity error constant

Kv=lim⁡z→1(1−z−1)GD(z)G(z)T=10.2⋅0.2228(1+0.8753)1−0.2543=5×0.41780.7457=2.80 s−1\begin{aligned} K_v &= \lim_{z\to1}\frac{(1-z^{-1})G_D(z)G(z)}{T} = \frac{1}{0.2}\cdot\frac{0.2228(1 + 0.8753)}{1 - 0.2543} \\ &= 5 \times \frac{0.4178}{0.7457} = 2.80\ \text{s}^{-1} \end{aligned}

Answer: GD(z)=12.67z−0.6703z−0.2543G_D(z) = 12.67\dfrac{z-0.6703}{z-0.2543} and Kv=2.80 s−1K_v = 2.80\ \text{s}^{-1}. The steady-state error for a unit ramp is 1/Kv=0.3571/K_v = 0.357.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Digital Control System chapter 4 questions →