Skip to main content

Lead compensator for k/[s(1+0.1s)(1+0.001s)], Kv = 1000

Short answer

Set K = 1000. Using Bode design, the required phase lead is 55°, α = 0.1, and ω_m is placed at 176 rad/s where uncompensated gain is −10 dB. The lead compensator is (1+0.018s)/(1+0.0018s) with amplifier gain 10, yielding PM ≈ 48° and Kv = 1000 s⁻¹.

Design a suitable lead compensator for a system having open loop TF G(s) = K/[s(1+0.1s)(1+0.001s)] such that the compensated system should have phase margin of at least 45° and static velocity error constant of at least 1000.

A lead compensator Gc(s)=1+Ts1+αTsG_c(s)=\dfrac{1+Ts}{1+\alpha Ts} (α<1\alpha<1) is designed on the Bode plot: it supplies its maximum lead ϕm\phi_m (sin⁡ϕm=1−α1+α\sin\phi_m=\frac{1-\alpha}{1+\alpha}) at the new gain crossover ωm=1Tα\omega_m=\frac{1}{T\sqrt\alpha}, where it adds 10log⁡(1/α)10\log(1/\alpha) dB of gain.

Step 1: K from the velocity error constant

Kv=lim⁡s→0sG(s)=K=1000 s−1K_v=\lim_{s\to0}sG(s)=K=1000\ \text{s}^{-1} G(s)=1000s(1+0.1s)(1+0.001s)G(s)=\frac{1000}{s(1+0.1s)(1+0.001s)}

Step 2: Uncompensated Bode plot

Corners at 10 and 1000 rad/s. Asymptotic magnitude: 60 dB at ω=1\omega=1, −20-20 dB/dec to 10 (40 dB), −40-40 dB/dec to 1000 (−40-40 dB), then −60-60 dB/dec.

ϕ(ω)=−90∘−tan⁡−10.1ω−tan⁡−10.001ω\phi(\omega)=-90^\circ-\tan^{-1}0.1\omega-\tan^{-1}0.001\omega
ω\omega11055.61001765561000
$G$ (dB)60.037.010.1−0.1-0.1−10.0-10.0
ϕ\phi−95.8∘-95.8^\circ−135.6∘-135.6^\circ−173∘-173^\circ−180∘-180^\circ−186.7∘-186.7^\circ−208∘-208^\circ−224.4∘-224.4^\circ

From ∣G∣=1|G|=1: ωgc=99.5\omega_{gc}=99.5 rad/s, ϕ=−179.94∘\phi=-179.94^\circ, so

PMunc≈0∘PM_{unc}\approx0^\circ

(and ωpc=100\omega_{pc}=100 rad/s, GM ≈ 0 dB). The system is on the verge of instability.

Step 3: Required phase lead

ϕm=PMd−PMunc+ϵ=45∘−0∘+10∘=55∘\phi_m=PM_d-PM_{unc}+\epsilon=45^\circ-0^\circ+10^\circ=55^\circ

A margin of ϵ=10∘\epsilon=10^\circ is taken because the crossover will move to a higher frequency where the plant phase is lower (a 5° margin gives only 44° and fails the check).

Step 4: α

α=1−sin⁡55∘1+sin⁡55∘=1−0.8191+0.819=0.0994≈0.1\alpha=\frac{1-\sin55^\circ}{1+\sin55^\circ}=\frac{1-0.819}{1+0.819}=0.0994\approx0.1

With α=0.1\alpha=0.1: ϕm=sin⁡−1(0.9/1.1)=54.9∘\phi_m=\sin^{-1}(0.9/1.1)=54.9^\circ.

Step 5: New gain crossover frequency

The lead adds 10log⁡(1/0.1)=1010\log(1/0.1)=10 dB at ωm\omega_m, so place ωm\omega_m where ∣G(jω)∣=−10|G(j\omega)|=-10 dB:

1000ω1+0.01ω21+10−6ω2=0.316⇒ωm≈176 rad/s\frac{1000}{\omega\sqrt{1+0.01\omega^2}\sqrt{1+10^{-6}\omega^2}}=0.316\Rightarrow\omega_m\approx176\ \text{rad/s}

Step 6: Time constants

T=1ωmα=1176×0.316=0.018 sαT=0.0018 s\begin{aligned} T&=\frac{1}{\omega_m\sqrt\alpha}=\frac{1}{176\times0.316}=0.018\ \text{s}\\ \alpha T&=0.0018\ \text{s} \end{aligned}

Zero at 1/T=55.61/T=55.6 rad/s, pole at 1/αT=5561/\alpha T=556 rad/s.

Gc(s)=1+0.018s1+0.0018sG_c(s)=\frac{1+0.018s}{1+0.0018s}

An amplifier of gain 1/α=101/\alpha=10 makes up the lead network's attenuation, so KvK_v stays 1000.

Step 7: Compensated system and check

Gc(s)G(s)=1000(1+0.018s)s(1+0.1s)(1+0.001s)(1+0.0018s)G_c(s)G(s)=\frac{1000(1+0.018s)}{s(1+0.1s)(1+0.001s)(1+0.0018s)}
ω\omega1055.6100176.8556696
$G_cG$ (dB)37.113.06.10
ϕ\phi−126∘-126^\circ−134∘-134^\circ−129∘-129^\circ−131.9∘-131.9^\circ−169∘-169^\circ−180∘-180^\circ

At ωgc′=176.8\omega_{gc}'=176.8 rad/s:

ϕ=−90∘−tan⁡−117.68−tan⁡−10.177+tan⁡−13.18−tan⁡−10.318=−90∘−86.8∘−10.0∘+72.5∘−17.6∘=−131.9∘PM=180∘−131.9∘=48.1∘ (≥45∘)\begin{aligned} \phi&=-90^\circ-\tan^{-1}17.68-\tan^{-1}0.177+\tan^{-1}3.18-\tan^{-1}0.318\\ &=-90^\circ-86.8^\circ-10.0^\circ+72.5^\circ-17.6^\circ=-131.9^\circ\\ PM&=180^\circ-131.9^\circ=48.1^\circ\ (\ge45^\circ) \end{aligned}

GM ≈ 17.5 dB at ωpc≈696\omega_{pc}\approx696 rad/s.

 dB 60 |\  uncompensated  -20
    37 |  \__ 10    compensated
     0 |------\----x----------- w
       |    100 \  176.8
       |         \__  -40/-60
 lead adds phase between 55.6 and 556 rad/s
ItemBeforeAfter
KvK_v10001000
ωgc\omega_{gc}99.5 rad/s176.8 rad/s
PM≈ 0°48.1°
GM≈ 0 dB17.5 dB

Answer: Gc(s)=1+0.018s1+0.0018sG_c(s)=\dfrac{1+0.018s}{1+0.0018s} (α = 0.1, T = 0.018 s) with K = 1000 gives Kv = 1000 s⁻¹ and PM ≈ 48° (≥ 45°).

✓ Answer checked by Claude Opus 5.5 (high effort)

All Control Systems chapter 7 questions →