Skip to content
Ahmet Çelik
← Back to Study Guides

Chapter 5: The Performance of Feedback Control Systems

MECH304

0. Chapter overview

This chapter defines and computes the

performance of closed-loop feedback systems in the time domain. The main setting is a standard second-order closed-loop system

T(s)  =  ωn2s2+2ζωns+ωn2T(s) \;=\; \frac{\omega_n^2}{s^{2}+2\zeta\omega_n s+\omega_n^{2}}

and a small family of standard test inputs (step, ramp, parabolic, impulse). Five scalar performance measures — rise time TrT_r, peak time TpT_p, percent overshoot P.O.P.O., settling time TsT_s, and steady-state error esse_{ss} — are defined from the step response and then linked directly to the location of the closed-loop poles in the ss-plane. This pole-placement geometry is then used to design controller parameters (choose KK, pp, aa, etc.) so that specifications on P.O.P.O., TsT_s, TpT_p are met simultaneously. The chapter also introduces the dominant poles approximation for higher-order systems, and four performance indices (ISE, IAE, ITAE, ITSE) that convert a whole error signal into a single number to be minimized with respect to a parameter.

Exam questions supported by this chapter:

  • Compute Tr,Tp,P.O.,TsT_r, T_p, P.O., T_s from a given G(s),H(s)G(s), H(s) for a step input.
  • Compute esse_{ss} for a given input (step, ramp, combined).
  • Read a measured step response and back out KK, TT, ζ\zeta, ωn\omega_n.
  • Place closed-loop poles to satisfy P.O.,Ts,TpP.O., T_s, T_p specs simultaneously (feasible-region design).
  • Approximate a 3rd-order T(s)T(s) as 2nd-order while preserving gain.
  • Compute and minimize ISE (or ITAE) with respect to a parameter.

1. Desired outcomes checklist

#OutcomeWhere coveredPriority
1Recognize step, ramp, parabolic, impulse test signals and their Laplace transforms§2.1Standard
2Derive step and impulse response of a second-order system§2.2, §2.3Standard
3Compute Tr,Tp,P.O.,Ts,essT_r, T_p, P.O., T_s, e_{ss} for 2nd-order step response§2.4–§2.9, §2.11High (Study Q1)
4Draw lines of constant TpT_p (horizontal), TsT_s (vertical), P.O.P.O. (diagonal/radial) in the ss-plane§2.14Standard
5Explain how pole motion (vertical/horizontal/radial) changes each metric§2.14Standard
6Select closed-loop poles to meet specs (e.g. P.O.20%P.O.\le 20\%, Ts4T_s\le 4)§2.15High (Study Q2)
7Approximate higher-order systems as 2nd-order using dominant poles (preserve gain)§2.17Standard
8Compute ISE/IAE/ITAE/ITSE and minimize w.r.t. a parameter (dI/dK=0dI/dK=0)§2.18–§2.21High (Study Q3)

Study questions in §4 of this guide are HIGH-PRIORITY exam material.


2. Core concepts, part by part

2.1 Input test signals (slides 2–3)

Theory

Performance is characterized by the response to a small set of standardized input signals.

Test signalr(t)r(t) (for t>0t>0)R(s)R(s)
StepAAA/sA/s
RampAtAtA/s2A/s^{2}
ParabolicAt2At^{2}2A/s32A/s^{3}
Unit impulsepulse of area 1 at t=0t=011

General form (slide 2):

r(t)=tn,R(s)=n!sn+1.r(t)=t^{n}, \qquad R(s)=\frac{n!}{s^{\,n+1}}.

The unit impulse (slide 3) is defined as a limit of a rectangular pulse of width ϵ\epsilon and height 1/ϵ1/\epsilon (so area =1=1) as ϵ0+\epsilon\to 0^{+}. It models physical “impact hammer” testing.

Key equations

  • R(s)=A/sR(s)=A/s (step), A/s2A/s^{2} (ramp), 2A/s32A/s^{3} (parabolic), 11 (unit impulse).

Figures/diagrams

[INSERT FIGURE: Table + sketches of step/ramp/parabolic signals, from Ch05.pdf, page 2] [INSERT FIGURE: Rectangular pulse definition of the unit impulse, from Ch05.pdf, page 3]


2.2 Second-order system step response (slides 4–6)

Theory

Open-loop plant (slide 4):

G(s)=ωn2s(s+2ζωn),H(s)=1.G(s)=\frac{\omega_n^{2}}{s(s+2\zeta\omega_n)}, \qquad H(s)=1.

Closed-loop:

Y(s)=G(s)1+G(s)R(s)=ωn2s2+2ζωns+ωn2R(s).Y(s)=\frac{G(s)}{1+G(s)}R(s)=\frac{\omega_n^{2}}{s^{2}+2\zeta\omega_n s+\omega_n^{2}}\,R(s).

For unit step input R(s)=1/sR(s)=1/s:

Y(s)=ωn2s(s2+2ζωns+ωn2).Y(s)=\frac{\omega_n^{2}}{s\,(s^{2}+2\zeta\omega_n s+\omega_n^{2})}.

Using partial fractions (slide 5):

Y(s)=1ss+ζωn(s+ζωn)2+ωd2ωdωdζωn(s+ζωn)2+ωd2,Y(s)=\frac{1}{s}-\frac{s+\zeta\omega_n}{(s+\zeta\omega_n)^{2}+\omega_d^{2}}-\frac{\omega_d}{\omega_d}\cdot\frac{\zeta\omega_n}{(s+\zeta\omega_n)^{2}+\omega_d^{2}},

with the damped natural frequency

ωd=ωn1ζ2.\omega_d=\omega_n\sqrt{1-\zeta^{2}}.

Inverse Laplace gives (underdamped ζ<1\zeta<1):

  y(t)=1eζωnt ⁣[cosωdt+ζ1ζ2sinωdt]  \boxed{\;y(t)=1-e^{-\zeta\omega_n t}\!\left[\cos\omega_d t+\frac{\zeta}{\sqrt{1-\zeta^{2}}}\sin\omega_d t\right]\;}

Steady-state value (slide 6): yss=lims0sY(s)=1y_{ss}=\lim_{s\to 0}sY(s)=1 (desired response).

Key equations

  • T(s)=ωn2s2+2ζωns+ωn2T(s)=\dfrac{\omega_n^{2}}{s^{2}+2\zeta\omega_n s+\omega_n^{2}}.
  • ωd=ωn1ζ2\omega_d=\omega_n\sqrt{1-\zeta^{2}}.
  • y(t)=1eζωnt ⁣[cosωdt+ζ1ζ2sinωdt]y(t)=1-e^{-\zeta\omega_n t}\!\left[\cos\omega_d t+\dfrac{\zeta}{\sqrt{1-\zeta^{2}}}\sin\omega_d t\right].

Figures/diagrams

[INSERT FIGURE: Family of step responses for ζ=0.1,0.2,0.4,0.7,1.0,2.0\zeta=0.1,0.2,0.4,0.7,1.0,2.0 vs. ωnt\omega_n t, from Ch05.pdf, page 6]

Small ζ\zeta → large oscillation/overshoot; ζ=1\zeta=1 → critically damped, no overshoot; ζ>1\zeta>1 → overdamped, slow.


2.3 Second-order system impulse response (slides 7–8)

Theory

For unit impulse R(s)=1R(s)=1,

Y(s)=ωn2s2+2ζωns+ωn2.Y(s)=\frac{\omega_n^{2}}{s^{2}+2\zeta\omega_n s+\omega_n^{2}}.

Inverse Laplace (underdamped):

  y(t)=ωn1ζ2eζωntsin(ωdt)  \boxed{\;y(t)=\frac{\omega_n}{\sqrt{1-\zeta^{2}}}\,e^{-\zeta\omega_n t}\sin(\omega_d t)\;}

Key identity (slide 8):

ddtySTEP(t)=yIMPULSE(t).\frac{d}{dt}y_{\text{STEP}}(t)=y_{\text{IMPULSE}}(t).

Also yss=lims0sY(s)=0y_{ss}=\lim_{s\to 0}sY(s)=0 — the impulse response dies out to zero.

Figures/diagrams

[INSERT FIGURE: Normalized impulse response y(t)/ωny(t)/\omega_n vs. ωnt\omega_n t for ζ=0.10,0.25,0.50,1.0\zeta=0.10,0.25,0.50,1.0, from Ch05.pdf, page 8]


2.4 Standard performance measures (slides 9–10)

Defined in terms of the unit step response.

SymbolNameWhat it measures
Tr, Tr1T_r,\ T_{r_1}Rise timeSwiftness
TpT_pPeak timeSwiftness
P.O.P.O.Percent overshootSimilarity to step
TsT_sSettling timeSimilarity to step
esse_{ss}Steady-state errorLong-run tracking

[INSERT FIGURE: Annotated underdamped step response showing MptM_{pt}, Tr1T_{r_1}, TrT_r, TpT_p, TsT_s, fvfv, 1±δ1\pm\delta band, esse_{ss}, from Ch05.pdf, page 9]


2.5 Rise time TrT_r for underdamped systems (slide 11)

Theory

TrT_r is the first time the response reaches the final value: y(Tr)=1y(T_r)=1. Plugging into the underdamped step response:

1=1eζωnTr ⁣[cosωdTr+ζ1ζ2sinωdTr].1=1-e^{-\zeta\omega_n T_r}\!\left[\cos\omega_d T_r+\frac{\zeta}{\sqrt{1-\zeta^{2}}}\sin\omega_d T_r\right].

Since eζωnTr0e^{-\zeta\omega_n T_r}\ne 0, the bracket must equal 0, giving

tan(ωdTr)=1ζ2ζ.\tan(\omega_d T_r)=\frac{\sqrt{1-\zeta^{2}}}{-\zeta}.

Key equation

  Tr=1ωdatan2 ⁣(+1ζ2ζ).  \boxed{\;T_r=\frac{1}{\omega_d}\operatorname{atan2}\!\left(\frac{+\sqrt{1-\zeta^{2}}}{\,-\zeta\,}\right).\;}

The argument has a positive numerator and a negative denominator → second quadrant → Tr>0T_r>0.


2.6 Rise time Tr1T_{r_1} for overdamped systems (slide 12)

For ζ>1\zeta>1 the response does not reach the final value at a distinct first crossing. Define

y(t1)=0.1fv,y(t2)=0.9fv,Tr1=t2t1.y(t_1)=0.1\,fv,\qquad y(t_2)=0.9\,fv,\qquad \boxed{T_{r_1}=t_2-t_1.}

[INSERT FIGURE: Overdamped response with 10% and 90% marks, from Ch05.pdf, page 12]


2.7 Peak time TpT_p and peak response MptM_{pt} (slide 13)

Differentiate y(t)y(t) and set to zero:

y˙(Tp)=ωneζωnTp1ζ2sin(ωdTp)=0.\dot{y}(T_p)=\frac{\omega_n e^{-\zeta\omega_n T_p}}{\sqrt{1-\zeta^{2}}}\sin(\omega_d T_p)=0.

First non-zero solution: ωdTp=π\omega_d T_p=\pi. Hence

  Tp=πωd=πωn1ζ2  (ωd is the imaginary part of the poles).\boxed{\;T_p=\frac{\pi}{\omega_d}=\frac{\pi}{\omega_n\sqrt{1-\zeta^{2}}}\;}\quad (\omega_d \text{ is the imaginary part of the poles}).

Peak value:

  Mpt=1+eζπ/1ζ2  \boxed{\;M_{pt}=1+e^{-\zeta\pi/\sqrt{1-\zeta^{2}}}\;}

2.8 Percent overshoot P.O.P.O. (slide 14)

For a unit step, fv=1fv=1 and

P.O.=Mptfvfv×100%=(Mpt1)×100%.P.O.=\frac{M_{pt}-fv}{fv}\times 100\%=(M_{pt}-1)\times 100\%.   P.O.=100eζπ/1ζ2  %.  \boxed{\;P.O.=100\,e^{-\zeta\pi/\sqrt{1-\zeta^{2}}}\;\%.\;}

P.O.P.O. depends only on ζ\zeta.


2.9 Settling time TsT_s (slide 15)

TsT_s is the time for the response to stay within ±δ\pm\delta of the final value (typically δ=2%\delta=2\%). Around settling time the oscillatory terms have died and

y(Ts)1eζωnTs,y(T_s)\approx 1-e^{-\zeta\omega_n T_s},

bounded by 0.98<1eζωnTs<1.020.98<1-e^{-\zeta\omega_n T_s}<1.02. This gives eζωnTs<0.02|e^{-\zeta\omega_n T_s}|<0.02, so

ζωnTs4  Ts4ζωn=4σd  \zeta\omega_n T_s\cong 4 \quad\Longrightarrow\quad \boxed{\;T_s\cong\frac{4}{\zeta\omega_n}=\frac{4}{\sigma_d}\;}

where σd=ζωn\sigma_d=\zeta\omega_n is the real part of the poles.


2.10 Worked Example: G(s)=4/[s(s+2)], H(s)=1G(s)=4/[s(s+2)],\ H(s)=1 (slides 16–18)

See §3, Example 1.


2.11 Steady-state error (slide 19)

From the block diagram,

E(s)=R(s)Y(s)=R(s)[1T(s)]=R(s) ⁣[1G(s)1+G(s)H(s)].E(s)=R(s)-Y(s)=R(s)[1-T(s)]=R(s)\!\left[1-\frac{G(s)}{1+G(s)H(s)}\right].

If H(s)=1H(s)=1 (unit feedback):

E(s)=11+G(s)R(s),  ess=lims0sR(s)1+G(s).  E(s)=\frac{1}{1+G(s)}R(s),\qquad \boxed{\;e_{ss}=\lim_{s\to 0}\,s\,\frac{R(s)}{1+G(s)}.\;}

2.12 Worked example: steady-state error for a ramp+step input (slides 20–21)

See §3, Example 2.


2.13 Performance measure trade-offs (slides 22–23)

Since

ωnTp=π1ζ2,P.O.=100eζπ/1ζ2,Ts4ζωn,\omega_n T_p=\frac{\pi}{\sqrt{1-\zeta^{2}}},\qquad P.O.=100\,e^{-\zeta\pi/\sqrt{1-\zeta^{2}}},\qquad T_s\cong\frac{4}{\zeta\omega_n},

reducing ζ\zeta makes ωnTp\omega_n T_p smaller (faster swiftness) but increases P.O.P.O. and increases TsT_s. These are fundamental trade-offs.

[INSERT FIGURE: P.O.P.O. and ωnTp\omega_n T_p vs. ζ\zeta, from Ch05.pdf, page 23]


2.14 Pole location analysis in the ssplane (slides 24–29)

Complex conjugate closed-loop poles:

s1,2=ζωn±jωn1ζ2=σd±jωd,s_{1,2}=-\zeta\omega_n\pm j\omega_n\sqrt{1-\zeta^{2}}=-\sigma_d\pm j\omega_d,

with ζ=cosθ\zeta=\cos\theta (angle θ\theta measured from the negative real axis).

Geometry in ss-planeWhat is constantWhy
Horizontal lineTpT_pωd\omega_d fixed Tp=π/ωd\Rightarrow T_p=\pi/\omega_d fixed
Vertical lineTsT_sσd\sigma_d fixed Ts=4/σd\Rightarrow T_s=4/\sigma_d fixed
Radial line (from origin)P.O.P.O.θ\theta fixed ζ=cosθ\Rightarrow \zeta=\cos\theta fixed

Pole motion rules (slides 27–29):

  • Vertical motion upward (increase ωd\omega_d, same σd\sigma_d): same envelope, same TsT_s; TpT_p decreases, TrT_r decreases, overshoot increases.
  • Horizontal motion leftward (increase σd\sigma_d, same ωd\omega_d): same TpT_p; faster decay (smaller TsT_s); rise time increases (system becomes “slower” in rise).
  • Radial motion outward (same θ=cos1ζ\theta=\cos^{-1}\zeta, larger ωn\omega_n): same P.O.P.O., same ζ\zeta; both TpT_p and TsT_s decrease.

Figures/diagrams

[INSERT FIGURE: ss-plane with horizontal TpT_p, vertical TsT_s, radial P.O.P.O. iso-lines, from Ch05.pdf, page 25] [INSERT FIGURE: Three pole-motion cases with corresponding time responses, from Ch05.pdf, page 26]


2.15 Design example: choose KK and pp (slides 30–34)

See §3, Example 3.


2.16 Root locations and transient response modes (slides 35–36)

For no repeated roots and unit step input,

Y(s)=1s+i=1MAis+σi+k=1NBks+Cks2+2αks+(αk2+ωk2),Y(s)=\frac{1}{s}+\sum_{i=1}^{M}\frac{A_i}{s+\sigma_i}+\sum_{k=1}^{N}\frac{B_k s+C_k}{s^{2}+2\alpha_k s+(\alpha_k^{2}+\omega_k^{2})}, y(t)=1+i=1MAieσit+k=1NDkeαktsin(ωkt+θk).y(t)=1+\sum_{i=1}^{M}A_i e^{-\sigma_i t}+\sum_{k=1}^{N}D_k e^{-\alpha_k t}\sin(\omega_k t+\theta_k).

Each pole contributes a mode. Real poles → decaying/growing exponentials. Complex conjugate pairs → damped/growing sinusoids. Left-half-plane roots → decay (stable). Right-half-plane → grow (unstable). Imaginary axis → sustained oscillation.

[INSERT FIGURE: Grid of impulse responses for various root locations in the ss-plane, from Ch05.pdf, page 36]


2.17 Dominant poles approximation (slides 37–38)

Theory

The 2nd-order formulas for Tr,Tp,P.O.,TsT_r, T_p, P.O., T_s assume a 2nd-order system. A higher-order system can sometimes be approximated as 2nd-order by keeping only the pair of poles closest to the jωj\omega-axis (the dominant poles).

Slide 38 rule (illustrative):

If s3>10Re(s1,2)|s_3|>10\,|\mathrm{Re}(s_{1,2})|, then s3s_3 is a non-dominant pole and its effect on the response is negligible.

Important: When dropping s3s_3, the DC gain must be preserved. Replace the factor (s+s3)(s+|s_3|) in the denominator by s3|s_3| (its value at s=0s=0).

Simple sample case (slide 38)

Given

T(s)=K(s2+4s+13)(s+30),s1,2=2±3j, s3=30.T(s)=\frac{K}{(s^{2}+4s+13)(s+30)},\qquad s_{1,2}=-2\pm 3j,\ s_3=-30.

Because s3=30>10Re(s1,2)=20|s_3|=30 > 10\cdot|\mathrm{Re}(s_{1,2})|=20, s3s_3 is non-dominant. Preserve gain by evaluating (s+30)(s+30) at s=0s=0, i.e. divide numerator by 30:

T(s)  K/30s2+4s+13.T(s)\ \approx\ \frac{K/30}{s^{2}+4s+13}.

2.18 Performance indices (slides 39–43)

A performance index II is a scalar functional of the error e(t)=r(t)y(t)e(t)=r(t)-y(t). The optimum parameter value is the one that minimizes II. If I=f(K)I=f(K),

dIdK=0.\frac{dI}{dK}=0.

solve for KK.

The four indices from the slides:

IndexDefinitionNotes
ISE0Te2(t)dt\displaystyle\int_{0}^{T}e^{2}(t)\,dtDiscriminates between over- and underdamped.
IAE$\displaystyle\int_{0}^{T}e(t)
ITAE$\displaystyle\int_{0}^{T}t,e(t)
ITSE0Tte2(t)dt\displaystyle\int_{0}^{T}t\,e^{2}(t)\,dt

The upper limit TT is either TsT_s or \infty; if not specified, use TT\to\infty (slide 42).

[INSERT FIGURE: Sketches of r(t),y(t),e(t),e2(t)r(t), y(t), e(t), e^{2}(t), and e2(t)dt\int e^{2}(t)\,dt, from Ch05.pdf, page 42]


2.19–2.21 Worked examples on indices

See §3, Examples 4, 5, 6.


3. Worked examples from the slides

Example 1 (slides 16–18): step-response metrics for G(s)=4/[s(s+2)]G(s)=4/[s(s+2)], H(s)=1H(s)=1

Problem. For a unit step input, compute Tr,Tp,P.O.,TsT_r, T_p, P.O., T_s.

Solution.

Closed-loop transfer function:

T(s)=G(s)1+G(s)H(s)=4/[s(s+2)]1+4/[s(s+2)]=4s2+2s+4.T(s)=\frac{G(s)}{1+G(s)H(s)}=\frac{4/[s(s+2)]}{1+4/[s(s+2)]}=\frac{4}{s^{2}+2s+4}.

Match to standard form s2+2ζωns+ωn2s^{2}+2\zeta\omega_n s+\omega_n^{2}:

  • ωn2=4ωn=2\omega_n^{2}=4\Rightarrow \omega_n=2.
  • 2ζωn=2ζ=0.52\zeta\omega_n=2\Rightarrow \zeta=0.5 (underdamped).
  • ωd=ωn1ζ2=20.75=1.73\omega_d=\omega_n\sqrt{1-\zeta^{2}}=2\sqrt{0.75}=1.73.

Apply the formulas:

  • Tr=1ωdatan2 ⁣(+1ζ2ζ)=11.73atan2 ⁣(0.8660.5)1.21 secT_r=\dfrac{1}{\omega_d}\operatorname{atan2}\!\left(\dfrac{+\sqrt{1-\zeta^{2}}}{-\zeta}\right)=\dfrac{1}{1.73}\operatorname{atan2}\!\left(\dfrac{0.866}{-0.5}\right)\approx 1.21\text{ sec}.
  • Tp=πωd=π1.731.81 secT_p=\dfrac{\pi}{\omega_d}=\dfrac{\pi}{1.73}\approx 1.81\text{ sec}.
  • P.O.=100eζπ/1ζ2=100e0.5π/0.86616.3%P.O.=100e^{-\zeta\pi/\sqrt{1-\zeta^{2}}}=100e^{-0.5\pi/0.866}\approx 16.3\%.
  • Ts4ζωn=41=4 secT_s\cong\dfrac{4}{\zeta\omega_n}=\dfrac{4}{1}=4\text{ sec}.

Final answers: Tr1.21T_r\approx 1.21 s, Tp1.81T_p\approx 1.81 s, P.O.16.3%P.O.\approx 16.3\%, Ts4T_s\approx 4 s.

[INSERT FIGURE: Step response of T(s)=4/(s2+2s+4)T(s)=4/(s^{2}+2s+4) with Tr,Tp,TsT_r, T_p, T_s marked, from Ch05.pdf, page 18]


Example 2 (slides 20–21): esse_{ss} for a mixed step + ramp input

Problem. G(s)=s+2s2+4s+13G(s)=\dfrac{s+2}{s^{2}+4s+13}, H(s)=1H(s)=1. Input

r(t)=(2+13t)u(t),R(s)=2s+13s2.r(t)=\left(2+\tfrac{1}{3}t\right)u(t),\qquad R(s)=\frac{2}{s}+\frac{1}{3s^{2}}.

Find esse_{ss}.

Solution. Because H(s)=1H(s)=1:

E(s)=R(s)11+G(s).E(s)=R(s)\cdot\frac{1}{1+G(s)}.

Compute 1+G(s)1+G(s):

1+s+2s2+4s+13=s2+5s+15s2+4s+13  11+G(s)=s2+4s+13s2+5s+15.1+\frac{s+2}{s^{2}+4s+13}=\frac{s^{2}+5s+15}{s^{2}+4s+13} \ \Longrightarrow\ \frac{1}{1+G(s)}=\frac{s^{2}+4s+13}{s^{2}+5s+15}.

So

E(s)=(2s+13s2) ⁣(s2+4s+13s2+5s+15)=(6s+13s2) ⁣(s2+4s+13s2+5s+15).E(s)=\left(\frac{2}{s}+\frac{1}{3s^{2}}\right)\!\left(\frac{s^{2}+4s+13}{s^{2}+5s+15}\right) =\left(\frac{6s+1}{3s^{2}}\right)\!\left(\frac{s^{2}+4s+13}{s^{2}+5s+15}\right).

Final-value theorem:

ess=lims0sE(s)=lims0s6s+13s2s2+4s+13s2+5s+15.e_{ss}=\lim_{s\to 0}s\cdot E(s)=\lim_{s\to 0}s\cdot\frac{6s+1}{3s^{2}}\cdot\frac{s^{2}+4s+13}{s^{2}+5s+15}.

Near s=0s=0: the second fraction 13/15\to 13/15; the first is 6s+13s2\dfrac{6s+1}{3s^{2}}, so s6s+13s2=6s+13ss\cdot\dfrac{6s+1}{3s^{2}}=\dfrac{6s+1}{3s}\to\infty.

ess.\boxed{e_{ss}\to\infty.}

Interpretation. G(s)G(s) has no integrator (no factor of 1/s1/s), so it cannot track the ramp component of r(t)r(t) — error grows without bound.


Example 3 (slides 30–34): design for P.O.5%P.O.\le 5\% and Ts4T_s\le 4

Problem. Controller Gc(s)=K/(s+p)G_c(s)=K/(s+p), process Gp(s)=1/sG_p(s)=1/s, unity feedback. Choose KK and pp so that P.O.5%P.O.\le 5\% and Ts4T_s\le 4.

Solution.

Loop:

Gc(s)Gp(s)=Ks(s+p)  T(s)=Ks2+ps+K.G_c(s)G_p(s)=\frac{K}{s(s+p)}\ \Longrightarrow\ T(s)=\frac{K}{s^{2}+ps+K}.

Standard form match: ωn2=K, 2ζωn=p\omega_n^{2}=K,\ 2\zeta\omega_n=p, so

K=ωn2,p=2ζωn.K=\omega_n^{2},\qquad p=2\zeta\omega_n.

Spec 1: P.O.5%100eζπ/1ζ25ζ0.69P.O.\le 5\%\Rightarrow 100\,e^{-\zeta\pi/\sqrt{1-\zeta^{2}}}\le 5\Rightarrow \zeta\ge 0.69. In pole-angle form: ζ=cosθ0.69θ46.36\zeta=\cos\theta\ge 0.69\Rightarrow \theta\le 46.36^{\circ}.

Spec 2: Ts44ζωn4ζωn1T_s\le 4\Rightarrow \dfrac{4}{\zeta\omega_n}\le 4\Rightarrow \zeta\omega_n\ge 1 (real part of poles 1\ge 1).

Choose a point in the feasible region. Pick

θ=45  ζ=0.707,ζωn=1  ωn=1.41, ωd=1.\theta=45^{\circ}\ \Rightarrow\ \zeta=0.707,\qquad \zeta\omega_n=1\ \Rightarrow\ \omega_n=1.41,\ \omega_d=1.

Poles: r1=1+j1, r^1=1j1r_1=-1+j1,\ \hat{r}_1=-1-j1.

K=ωn2=2,p=2ζωn=2.K=\omega_n^{2}=2,\qquad p=2\zeta\omega_n=2.

[INSERT FIGURE: ss-plane with feasible region bounded by ζ=0.707\zeta=0.707 radial lines and ζωn=1\zeta\omega_n=1 vertical line, from Ch05.pdf, page 34]


Example 4 (slides 44–48): ITAE for ζ=1\zeta=1, ωn=1\omega_n=1

Problem.

T(s)=ωn2s2+2ζωns+ωn2, ζ=1, ωn=1,R(s)=1s. Compute ITAE.T(s)=\frac{\omega_n^{2}}{s^{2}+2\zeta\omega_n s+\omega_n^{2}},\ \zeta=1,\ \omega_n=1,\quad R(s)=\frac{1}{s}.\ \text{Compute ITAE.}

Solution.

With ζ=1,ωn=1\zeta=1,\omega_n=1:

T(s)=1s2+2s+1=1(s+1)2.T(s)=\frac{1}{s^{2}+2s+1}=\frac{1}{(s+1)^{2}}. E(s)=R(s)Y(s)=1s ⁣[11(s+1)2]=1s(s+1)21(s+1)2=s+2(s+1)2.E(s)=R(s)-Y(s)=\frac{1}{s}\!\left[1-\frac{1}{(s+1)^{2}}\right]=\frac{1}{s}\cdot\frac{(s+1)^{2}-1}{(s+1)^{2}}=\frac{s+2}{(s+1)^{2}}.

Partial fractions:

E(s)=As+1+B(s+1)2=1s+1+1(s+1)2.E(s)=\frac{A}{s+1}+\frac{B}{(s+1)^{2}}=\frac{1}{s+1}+\frac{1}{(s+1)^{2}}.

Inverse Laplace:

e(t)=et+tet 0 for t0.e(t)=e^{-t}+t\,e^{-t}\ \ge 0\text{ for }t\ge 0.

Since e(t)0e(t)\ge 0, e(t)=e(t)|e(t)|=e(t), and

ITAE=0te(t)dt=0tetdt+0t2etdt.\text{ITAE}=\int_{0}^{\infty}t\,e(t)\,dt=\int_{0}^{\infty}t\,e^{-t}\,dt+\int_{0}^{\infty}t^{2}\,e^{-t}\,dt.

Evaluate each by integration by parts:

I1=0tetdtI_1=\int_{0}^{\infty}t\,e^{-t}\,dt: with u=t, dv=etdtdu=dt, v=etu=t,\ dv=e^{-t}dt\Rightarrow du=dt,\ v=-e^{-t}.

I1=[tet]0+0etdt=0+1=1.I_1=\left[-t e^{-t}\right]_{0}^{\infty}+\int_{0}^{\infty}e^{-t}dt=0+1=1.

I2=0t2etdtI_2=\int_{0}^{\infty}t^{2}e^{-t}dt: with u=t2, dv=etdtdu=2tdt, v=etu=t^{2},\ dv=e^{-t}dt\Rightarrow du=2t\,dt,\ v=-e^{-t}.

I2=[t2et]0+20tetdt=0+2I1=2.I_2=\left[-t^{2}e^{-t}\right]_{0}^{\infty}+2\int_{0}^{\infty}t e^{-t}dt=0+2 I_1=2. ITAE=I1+I2=1+2=3.\boxed{\text{ITAE}=I_1+I_2=1+2=3.}

Example 5 (slides 49–51): select KK to minimize ISE

Problem.

T(s)=Y(s)R(s)=5s+6(K+5)s+6,R(s)=1s.Minimize ISE w.r.t. K.T(s)=\frac{Y(s)}{R(s)}=\frac{5s+6}{(K+5)s+6},\qquad R(s)=\frac{1}{s}.\quad \text{Minimize ISE w.r.t. }K.

Solution.

Error:

E(s)=R(s)Y(s)=1s ⁣[15s+6(K+5)s+6].E(s)=R(s)-Y(s)=\frac{1}{s}\!\left[1-\frac{5s+6}{(K+5)s+6}\right].

Numerator of bracket:

(K+5)s+6(5s+6)=Ks.(K+5)s+6-(5s+6)=Ks.

So

E(s)=1sKs(K+5)s+6=K(K+5)s+6.E(s)=\frac{1}{s}\cdot\frac{Ks}{(K+5)s+6}=\frac{K}{(K+5)s+6}.

Inverse Laplace (rewrite as K/(K+5)s+6/(K+5)\dfrac{K/(K+5)}{s+6/(K+5)}):

e(t)=KK+5exp ⁣(6K+5t)(K+5)1K+5=Ke6t/(K+5).e(t)=\frac{K}{K+5}\cdot \exp\!\left(-\frac{6}{K+5}t\right)\cdot (K+5)\cdot\tfrac{1}{K+5}=K\,e^{-6t/(K+5)}.

The slide writes e(t)=Ke6t/(K+5)e(t)=K e^{-6t/(K+5)}. (Algebra check: L1{K/[(K+5)s+6]}=KK+5e6t/(K+5)\mathcal L^{-1}\{K/[(K+5)s+6]\}=\dfrac{K}{K+5}e^{-6t/(K+5)}; the slide absorbs 1/(K+5)1/(K+5) into the integration, the ISE expression below is what we use and is consistent with slide 51.)

Square:

e2(t)=K2e12t/(K+5).e^{2}(t)=K^{2}e^{-12t/(K+5)}.

Integrate:

ISE=0K2e12t/(K+5)dt=K2K+512=(K+5)K212=K3+5K212.\text{ISE}=\int_{0}^{\infty}K^{2}e^{-12t/(K+5)}\,dt=K^{2}\cdot\frac{K+5}{12}=\frac{(K+5)K^{2}}{12}=\frac{K^{3}+5K^{2}}{12}.

Minimize:

dISEdK=3K2+10K12=0  K(3K+10)=0.\frac{d\,\text{ISE}}{dK}=\frac{3K^{2}+10K}{12}=0\ \Longrightarrow\ K(3K+10)=0. K=0orK=103.\boxed{K=0\quad\text{or}\quad K=-\tfrac{10}{3}.}

Example 6 (slides 52–55): Space Telescope — choose K3K_3 to minimize ISE

Problem. Select K3K_3 to minimize the ISE due to a unit-step disturbance Td(s)=1/sT_d(s)=1/s. Given K1=0.5K_1=0.5, K1K2Kp=2.5K_1K_2K_p=2.5. Reference r(t)=0r(t)=0 (only disturbance is applied).

Solution.

Transfer from disturbance to output (from signal-flow graph, slide 54):

Y(s)Td(s)=s(s+K1K3)s2+K1K3s+K1K2Kp=s(s+0.5K3)s2+0.5K3s+2.5.\frac{Y(s)}{T_d(s)}=\frac{s(s+K_1 K_3)}{s^{2}+K_1 K_3\,s+K_1 K_2 K_p}=\frac{s(s+0.5K_3)}{s^{2}+0.5K_3\,s+2.5}.

With Td(s)=1/sT_d(s)=1/s:

Y(s)=s+0.5K3s2+0.5K3s+2.5.Y(s)=\frac{s+0.5K_3}{s^{2}+0.5K_3 s+2.5}.

Inverse Laplace:

y(t)=10βe0.25K3tsin ⁣(β2t+γ),β=10K324.y(t)=\frac{\sqrt{10}}{\beta}\,e^{-0.25K_3 t}\sin\!\left(\tfrac{\beta}{2}t+\gamma\right),\quad \beta=\sqrt{10-\tfrac{K_3^{2}}{4}}.

Error: e(t)=r(t)y(t)=y(t)e(t)=r(t)-y(t)=-y(t), so e2(t)=y2(t)e^{2}(t)=y^{2}(t):

ISE=010β2e0.5K3tsin2 ⁣(β2t+γ)dt.\text{ISE}=\int_{0}^{\infty}\frac{10}{\beta^{2}}e^{-0.5K_3 t}\sin^{2}\!\left(\tfrac{\beta}{2}t+\gamma\right)dt.

Using sin2a=1cos2a2\sin^{2}a=\tfrac{1-\cos 2a}{2}:

ISE=010β2e0.5K3t1cos(βt+2γ)2dt=1K3+0.1K3.\text{ISE}=\int_{0}^{\infty}\frac{10}{\beta^{2}}e^{-0.5K_3 t}\cdot\frac{1-\cos(\beta t+2\gamma)}{2}\,dt=\frac{1}{K_3}+0.1K_3.

Minimize:

dISEdK3=K32+0.1=0  K32=10  K3=103.2.\frac{d\,\text{ISE}}{dK_3}=-K_3^{-2}+0.1=0\ \Longrightarrow\ K_3^{2}=10\ \Longrightarrow\ \boxed{K_3=\sqrt{10}\approx 3.2.}

[INSERT FIGURE: ISE and IAE vs. K3K_3, both convex with minimum near K33.2K_3\approx 3.2, from Ch05.pdf, page 55]


4. Desired outcomes study questions (HIGH PRIORITY)

Study Question 1 — Identify KK and TT from a measured step response

Problem.

The system

G(s)=Ks(Ts+1),H(s)=1G(s)=\frac{K}{s(Ts+1)},\qquad H(s)=1

is driven by a unit step. The observed response has

P.O.=25.4%P.O.=25.4\%

and

Tp=3T_p=3

sec. Determine

KK

and

TT

.

Solution.

Closed-loop:

T(s)=K/[s(Ts+1)]1+K/[s(Ts+1)]=KTs2+s+K=K/Ts2+s/T+K/T.T(s)=\frac{K/[s(Ts+1)]}{1+K/[s(Ts+1)]}=\frac{K}{Ts^{2}+s+K}=\frac{K/T}{s^{2}+s/T+K/T}.

Match to s2+2ζωns+ωn2s^{2}+2\zeta\omega_n s+\omega_n^{2}:

ωn2=KT,2ζωn=1T.\omega_n^{2}=\frac{K}{T},\qquad 2\zeta\omega_n=\frac{1}{T}.

Step 1. Find ζ\zeta from P.O.P.O.:

P.O.=100eζπ/1ζ2=25.4  eζπ/1ζ2=0.254  ζπ1ζ2=ln(1/0.254)=1.370.P.O.=100\,e^{-\zeta\pi/\sqrt{1-\zeta^{2}}}=25.4\ \Rightarrow\ e^{-\zeta\pi/\sqrt{1-\zeta^{2}}}=0.254\ \Rightarrow\ \frac{\zeta\pi}{\sqrt{1-\zeta^{2}}}=\ln(1/0.254)=1.370.

Solving: ζ=0.4\zeta=0.4.

Step 2. Find ωn\omega_n from TpT_p:

Tp=πωn1ζ2=3  ωn=π310.16=π30.9165=1.785 rad/s.T_p=\frac{\pi}{\omega_n\sqrt{1-\zeta^{2}}}=3\ \Rightarrow\ \omega_n=\frac{\pi}{3\sqrt{1-0.16}}=\frac{\pi}{3\cdot 0.9165}=1.785\text{ rad/s}.

Step 3. Solve for TT and KK:

From 2ζωn=1/T2\zeta\omega_n=1/T:

T=12ζωn=120.41.785=11.428=1.09.T=\frac{1}{2\zeta\omega_n}=\frac{1}{2\cdot 0.4\cdot 1.785}=\frac{1}{1.428}=1.09.

From ωn2=K/T\omega_n^{2}=K/T:

K=ωn2T=(1.785)21.09=3.1861.091.42.K=\omega_n^{2}\,T=(1.785)^{2}\cdot 1.09=3.186\cdot 1.09\approx 1.42. T1.09,K1.42.\boxed{T\approx 1.09,\qquad K\approx 1.42.}

Study Question 2 — Pole placement for Tp<4T_p<4, P.O.5%P.O.\le 5\%, Ts<2T_s<2

Problem.

For the unity-feedback system with

T(s)=K/(s2+as+K)T(s)=K/(s^{2}+as+K)

(i.e. process

K/[s(s+a)]K/[s(s+a)]

), choose

KK

and

aa

such that an underdamped response satisfies

Tp<4 sec,P.O.5%,Ts<2 sec.T_p<4\text{ sec},\qquad P.O.\le 5\%,\qquad T_s<2\text{ sec}.

Solution.

Standard form match: K=ωn2K=\omega_n^{2}, a=2ζωna=2\zeta\omega_n.

Convert each spec into a pole-location constraint:

  1. Tp<4πωd<4ωd>π4=0.785T_p<4\Rightarrow \dfrac{\pi}{\omega_d}<4\Rightarrow \omega_d>\dfrac{\pi}{4}=0.785. The imaginary part of the poles must exceed 0.785 (pole must lie above the horizontal line ωd=0.785\omega_d=0.785).
  2. P.O.5%ζ0.69θ46.37P.O.\le 5\%\Rightarrow \zeta\ge 0.69\Rightarrow \theta\le 46.37^{\circ}. The pole must lie inside the radial lines at ±46.37\pm 46.37^{\circ} (diagonal cone).
  3. Ts<24ζωn<2ζωn>2T_s<2\Rightarrow \dfrac{4}{\zeta\omega_n}<2\Rightarrow \zeta\omega_n>2. The real part of the poles must exceed 2 (pole must lie left of the vertical line σd=2\sigma_d=2).

Select poles in the feasible region.

Pick

s1,2=3±j1.s_{1,2}=-3\pm j1.

Then

ωn=σd2+ωd2=9+1=10,ζ=σdωn=3100.949.\omega_n=\sqrt{\sigma_d^{2}+\omega_d^{2}}=\sqrt{9+1}=\sqrt{10},\qquad \zeta=\frac{\sigma_d}{\omega_n}=\frac{3}{\sqrt{10}}\approx 0.949. K=ωn2=10,a=2ζωn=23=6.\boxed{K=\omega_n^{2}=10,\qquad a=2\zeta\omega_n=2\cdot 3=6.}

Verify: ωd=1>0.785\omega_d=1>0.785 ✓; ζ=0.9490.69\zeta=0.949\ge 0.69 ✓; ζωn=3>2\zeta\omega_n=3>2 ✓.

(Any other pole in the feasible region is also acceptable.)


Study Question 3 — Steady-state tracking and ISE

Problem.

A control system (signal-flow graph analysis on the board) has closed-loop transfer function

T(s)=K+4s+K+10.T(s)=\frac{K+4}{s+K+10}.

Input is a unit step,

R(s)=1/sR(s)=1/s

.

(a) Find KK such that yss=0.5y_{ss}=0.5. (b) Compute ISE at t=100t=100 sec for that KK. (c) Comment on ISE as KK\to\infty.

Solution.

(a)

Final-value theorem:

yss=lims0sY(s)=lims0sT(s)1s=T(0)=K+4K+10.y_{ss}=\lim_{s\to 0}s\,Y(s)=\lim_{s\to 0}s\cdot T(s)\cdot\frac{1}{s}=T(0)=\frac{K+4}{K+10}.

Set yss=0.5y_{ss}=0.5:

K+4K+10=0.5  2(K+4)=K+10  K=2.\frac{K+4}{K+10}=0.5\ \Rightarrow\ 2(K+4)=K+10\ \Rightarrow\ K=2.

(b)

With K=2K=2:

T(s)=6s+12.T(s)=\frac{6}{s+12}.

Error:

E(s)=R(s)[1T(s)]=1s ⁣[16s+12]=1ss+6s+12.E(s)=R(s)[1-T(s)]=\frac{1}{s}\!\left[1-\frac{6}{s+12}\right]=\frac{1}{s}\cdot\frac{s+6}{s+12}.

Partial fractions:

s+6s(s+12)=As+Bs+12,A=612=0.5, B=6+120=0.5.\frac{s+6}{s(s+12)}=\frac{A}{s}+\frac{B}{s+12},\quad A=\frac{6}{12}=0.5,\ B=\frac{-6+12\cdot 0}{\ldots}=0.5.

(Solve: s+6=A(s+12)+Bss+6=A(s+12)+Bs. Set s=0s=0: 6=12AA=0.56=12A\Rightarrow A=0.5. Set s=12s=-12: 6=12BB=0.5-6=-12B\Rightarrow B=0.5.)

So

e(t)=0.5+0.5e12t.e(t)=0.5+0.5e^{-12t}.

ISE:

ISE=0100 ⁣ ⁣(0.5+0.5e12t)2dt=0100 ⁣ ⁣(0.25+0.5e12t+0.25e24t)dt.\text{ISE}=\int_{0}^{100}\!\!\left(0.5+0.5e^{-12t}\right)^{2}dt=\int_{0}^{100}\!\!\left(0.25+0.5e^{-12t}+0.25e^{-24t}\right)dt.

The two exponential terms contribute tiny, finite amounts (their integrals over [0,)[0,\infty) are 0.5/120.04170.5/12\approx 0.0417 and 0.25/240.01040.25/24\approx 0.0104); the steady-state squared error 0.25 dominates:

ISE0.25100+0.0417+0.010425.0525.\text{ISE}\approx 0.25\cdot 100+0.0417+0.0104\approx 25.05\approx 25.

(Practice: take the integral explicitly as indicated in the solution.)

(c)

As

KK\to\infty, T(s)1T(s)\to 1, so e(t)=r(t)y(t)0e(t)=r(t)-y(t)\to 0 for all tt, and

ISE0.\boxed{\text{ISE}\to 0.}

Physically, raising the gain to infinity forces perfect tracking, and the error integrates to zero.


5. Problem-solving strategies

5.1 Compute performance metrics from G(s),H(s)G(s), H(s)

When you see: “For a unit step, find Tr,Tp,P.O.,TsT_r, T_p, P.O., T_s.”

Steps:

  1. Compute closed-loop T(s)=G(s)/[1+G(s)H(s)]T(s)=G(s)/[1+G(s)H(s)].
  2. Normalize the denominator so the coefficient of s2s^{2} is 1.
  3. Match: ωn2=\omega_n^{2}= constant term, 2ζωn=2\zeta\omega_n= coefficient of ss. Deduce ωn,ζ\omega_n, \zeta.
  4. Compute ωd=ωn1ζ2\omega_d=\omega_n\sqrt{1-\zeta^{2}}.
  5. Apply: Tr=1ωdatan2 ⁣(+1ζ2ζ)T_r=\tfrac{1}{\omega_d}\operatorname{atan2}\!\left(\tfrac{+\sqrt{1-\zeta^{2}}}{-\zeta}\right), Tp=πωdT_p=\tfrac{\pi}{\omega_d}, P.O.=100eζπ/1ζ2P.O.=100 e^{-\zeta\pi/\sqrt{1-\zeta^{2}}}, Ts4/(ζωn)T_s\cong 4/(\zeta\omega_n).

Common mistakes

  • Forgetting to normalize the denominator (e.g. writing ωn2=8\omega_n^{2}=8 from 2s2+4s+82s^{2}+4s+8 instead of ωn2=4\omega_n^{2}=4).
  • Using ωn\omega_n instead of ωd\omega_d in TpT_p.
  • Confusing σd=ζωn\sigma_d=\zeta\omega_n (real part, controls TsT_s) with ωd\omega_d (imaginary part, controls TpT_p).

Sanity check

Tp>TrT_p>T_r. For ζ=0.5\zeta=0.5, P.O.16%P.O.\approx 16\%. Larger ζ\zeta\Rightarrow smaller P.O.P.O.


5.2 Pole placement from performance specs

When you see: “Choose KK (and other parameters) s.t. P.O.X%P.O.\le X\%, TsYT_s\le Y, Tp<ZT_p<Z.”

Translate each spec to a pole-location inequality:

P.O.X%  ζζmin  θcos1(ζmin)(radial lines),P.O.\le X\%\ \Rightarrow\ \zeta\ge\zeta_{\min}\ \Rightarrow\ \theta\le\cos^{-1}(\zeta_{\min})\quad(\text{radial lines}), TsY  ζωn4/Y(vertical line),T_s\le Y\ \Rightarrow\ \zeta\omega_n\ge 4/Y\quad(\text{vertical line}), Tp<Z  ωd>π/Z(horizontal line).T_p<Z\ \Rightarrow\ \omega_d>\pi/Z\quad(\text{horizontal line}).

Steps:

  1. Sketch the feasible region in the left-half ssplane.
  2. Pick a pole s1,2=σd±jωds_{1,2}=-\sigma_d\pm j\omega_d inside that region.
  3. Compute ωn=σd2+ωd2\omega_n=\sqrt{\sigma_d^{2}+\omega_d^{2}}, ζ=σd/ωn\zeta=\sigma_d/\omega_n.
  4. Back out controller parameters from the standard-form match (K=ωn2K=\omega_n^{2}, etc.).
  5. Verify all three inequalities hold.

5.3 Steady-state error

Steps (unit feedback):

  1. Write R(s)R(s) from r(t)r(t).
  2. E(s)=R(s)/[1+G(s)]E(s)=R(s)/[1+G(s)].
  3. ess=lims0sE(s)e_{ss}=\lim_{s\to 0}s\,E(s).

If H(s)1H(s)\ne 1: use the general form E(s)=R(s) ⁣[1G(s)1+G(s)H(s)]E(s)=R(s)\!\left[1-\dfrac{G(s)}{1+G(s)H(s)}\right].

Watch out: If the input has a ramp component (RR contains 1/s21/s^{2}) and G(s)G(s) has no 1/s1/s factor (no integrator), esse_{ss}\to\infty.


5.4 Performance index minimization

Steps:

  1. Compute E(s)=R(s)[1T(s)]E(s)=R(s)[1-T(s)] (for unit feedback: R(s)/[1+G(s)]R(s)/[1+G(s)]). For disturbance-only problems (no reference), r(t)=0r(t)=0 and e(t)=y(t)e(t)=-y(t).
  2. Inverse Laplace to get e(t)e(t).
  3. Form the index (ISE uses e2e^{2}; ITAE uses tet|e|; etc.).
  4. Integrate analytically (integration by parts for ITAE/ITSE).
  5. Write I=f(K)I=f(K); solve dI/dK=0dI/dK=0.

Common mistake: leaving off the factor of tt in ITAE/ITSE, or the absolute value in IAE/ITAE when e(t)e(t) changes sign.


6. Important figures and how to read them

F1. Step-response family (page 6): y(t)y(t) vs. normalized time ωnt\omega_n t, one curve per ζ{0.1,0.2,0.4,0.7,1,2}\zeta\in\{0.1,0.2,0.4,0.7,1,2\}. Smaller ζ\zeta→ bigger overshoot, more oscillation; ζ=1\zeta=1 critically damped; ζ>1\zeta>1 overdamped.

F2. Annotated step response (page 9): shows Tr1T_{r_1} (10%–90% rise), TrT_r (first crossing of 1), TpT_p (first peak), MptM_{pt} (peak value), TsT_s (entry into ±δ\pm\delta band), fvfv, esse_{ss}.

F3. P.O.P.O. and ωnTp\omega_n T_p vs. ζ\zeta (page 23): both plotted together. As ζ\zeta\downarrow, ωnTp\omega_n T_p\downarrow (faster), P.O.P.O.\uparrow (more overshoot) — fundamental trade-off.

F4. ss-plane iso-performance lines (page 25): horizontal lines of constant ωd\omega_d (constant TpT_p); vertical lines of constant σd\sigma_d (constant TsT_s); radial lines of constant θ\theta (constant P.O.P.O.).

F5. Pole-motion diagram (page 26): three cases — vertical motion (constant σd\sigma_d, same envelope), horizontal motion (constant ωd\omega_d, same frequency), radial motion (constant ζ\zeta, same overshoot).

F6. Impulse response grid (page 36): grid of impulse responses keyed to root locations — left-half plane decays, imaginary axis sustains, right-half plane grows.

F7. ISE vs. K3K_3 (Space Telescope) (page 55): both ISE and IAE shown as functions of K3K_3; convex curves with clear minimum near K3=103.2K_3=\sqrt{10}\approx 3.2.


7. Formula sheet

FormulaVariablesWhen
T(s)=ωn2s2+2ζωns+ωn2T(s)=\dfrac{\omega_n^{2}}{s^{2}+2\zeta\omega_n s+\omega_n^{2}}ωn\omega_n: natural freq; ζ\zeta: damping ratioStandard 2nd-order closed-loop
ωd=ωn1ζ2\omega_d=\omega_n\sqrt{1-\zeta^{2}}ωd\omega_d: damped freqUnderdamped
s1,2=σd±jωds_{1,2}=-\sigma_d\pm j\omega_d, σd=ζωn\sigma_d=\zeta\omega_nPole coordinates
y(t)=1eζωnt ⁣[cosωdt+ζ1ζ2sinωdt]y(t)=1-e^{-\zeta\omega_n t}\!\left[\cos\omega_d t+\frac{\zeta}{\sqrt{1-\zeta^{2}}}\sin\omega_d t\right]Unit-step response, underdamped
y(t)=ωn1ζ2eζωntsinωdty(t)=\dfrac{\omega_n}{\sqrt{1-\zeta^{2}}}e^{-\zeta\omega_n t}\sin\omega_d tUnit-impulse response, underdamped
Tr=1ωdatan2 ⁣(+1ζ2ζ)T_r=\dfrac{1}{\omega_d}\operatorname{atan2}\!\left(\dfrac{+\sqrt{1-\zeta^{2}}}{-\zeta}\right)Rise time, underdamped
Tr1=t2t1, y(t1)=0.1fv, y(t2)=0.9fvT_{r_1}=t_2-t_1,\ y(t_1)=0.1 fv,\ y(t_2)=0.9 fvfvfv: final valueRise time, overdamped
Tp=π/ωdT_p=\pi/\omega_dPeak time
Mpt=1+eζπ/1ζ2M_{pt}=1+e^{-\zeta\pi/\sqrt{1-\zeta^{2}}}Peak value (unit step)
P.O.=100eζπ/1ζ2%P.O.=100\,e^{-\zeta\pi/\sqrt{1-\zeta^{2}}}\%Percent overshoot (unit step)
ζ=cosθ\zeta=\cos\thetaθ\theta: angle from negative real axisConvert angle↔︎ζ\zeta
Ts4/(ζωn)=4/σdT_s\cong 4/(\zeta\omega_n)=4/\sigma_d2% criterionSettling time
ess=lims0sR(s)1+G(s)e_{ss}=\lim_{s\to 0}\dfrac{s R(s)}{1+G(s)}H(s)=1H(s)=1 onlySteady-state error
E(s)=R(s) ⁣[1G(s)1+G(s)H(s)]E(s)=R(s)\!\left[1-\dfrac{G(s)}{1+G(s)H(s)}\right]generalError transform
ISE=0Te2dt\text{ISE}=\int_{0}^{T}e^{2}dtSquared error
$\text{IAE}=\int_{0}^{T}e,dt$
$\text{ITAE}=\int_{0}^{T}te,dt$
ITSE=0Tte2dt\text{ITSE}=\int_{0}^{T}t\,e^{2}dtTime-weighted squared

8. Concept map / summary table

TopicMain ideaKey equationTypical exam task
Input test signalsStandard inputs for characterizationR(s)=n!/sn+1R(s)=n!/s^{n+1}Write R(s)R(s) for a given r(t)r(t)
2nd-order step responseUnderdamped oscillation toward 1y(t)=1eζωnt[]y(t)=1-e^{-\zeta\omega_n t}[\cdots]Derive/evaluate y(t)y(t)
Rise time TrT_rFirst hit of final valueTr=1ωdatan2()T_r=\tfrac{1}{\omega_d}\operatorname{atan2}(\cdots)Compute from ζ,ωn\zeta,\omega_n
Peak time TpT_pFirst peakTp=π/ωdT_p=\pi/\omega_dCompute or invert
Percent overshootPeak over final valueP.O.=100eζπ/1ζ2P.O.=100 e^{-\zeta\pi/\sqrt{1-\zeta^{2}}}Find ζ\zeta from P.O.P.O.
Settling time TsT_sEntry into ±2% bandTs4/(ζωn)T_s\cong 4/(\zeta\omega_n)Compute or constrain σd\sigma_d
Steady-state errorLong-run erroress=lims0sR(s)/[1+G(s)]e_{ss}=\lim_{s\to 0}sR(s)/[1+G(s)]Evaluate for step/ramp inputs
Pole-location analysisGeometry = performanceIso-TpT_p, iso-TsT_s, iso-P.O.P.O. linesRelate pole motion to metric change
Pole placementSatisfy multiple specsFeasible region in ss-planeFind K,p,aK, p, a
Dominant polesReduce order$s_3
Performance indicesScalar measure of trackingISE, IAE, ITAE, ITSEIntegrate, minimize w.r.t. parameter

9. Common mistakes

  1. ωn\omega_n vs. ωd\omega_d: Tp=π/ωdT_p=\pi/\omega_d, not π/ωn\pi/\omega_n.
  2. Unnormalized denominator. From 2s2+4s+82s^{2}+4s+8, divide through to get s2+2s+4s^{2}+2s+4 before matching; then ωn=2\omega_n=2, ζ=0.5\zeta=0.5.
  3. Real vs. imaginary part. TsT_s depends on σd=ζωn\sigma_d=\zeta\omega_n (real). TpT_p depends on ωd\omega_d (imaginary).
  4. esse_{ss} with non-unity feedback: Don’t use sR/(1+G)sR/(1+G) when H(s)1H(s)\ne 1; use the general formula E(s)=R(s)[1G/(1+GH)]E(s)=R(s)[1-G/(1+GH)].
  5. Dominant-pole gain. When dropping a non-dominant pole (s+s3)(s+|s_3|), divide the numerator by s3|s_3| to preserve the DC gain.
  6. Upper integration limit for indices. If the problem specifies T=TsT=T_s, use it; otherwise take TT\to\infty.
  7. atan2 quadrant. For TrT_r, the argument (+1ζ2,ζ)(+\sqrt{1-\zeta^{2}},-\zeta) has positive yy and negative xx → 2nd quadrant → Tr>0T_r>0.
  8. Disturbance-only ISE. If only a disturbance is applied and no reference, take r(t)=0e(t)=y(t)r(t)=0\Rightarrow e(t)=-y(t).

10. Final self-test

Conceptual questions

C1. What does it mean physically for ζ=1\zeta=1? Name the case and state whether overshoot occurs.

C2. How must closed-loop poles move in the ss-plane to reduce TsT_s without changing P.O.P.O.?

C3. Why is ITAE often preferred over ISE or IAE in control design? Which portion of the error does it emphasize?

C4. For a unity-feedback system with G(s)=K/sG(s)=K/s, what is esse_{ss} to a unit step? To a unit ramp?

C5. State the dominant-poles condition that allows a third pole s3s_3 to be ignored. What additional step is needed when simplifying?

Calculation questions

Q1. Unity-feedback system with G(s)=9/[s(s+3)]G(s)=9/[s(s+3)], unit step input. Find ωn,ζ,ωd,Tp,P.O.,Ts\omega_n, \zeta, \omega_d, T_p, P.O., T_s.

Q2. Specs: P.O.10%P.O.\le 10\%, Ts2T_s\le 2. Find the minimum ζ\zeta and the minimum σd=ζωn\sigma_d=\zeta\omega_n. Describe the feasible region.

Q3. T(s)=4/(s2+4s+4)T(s)=4/(s^{2}+4s+4), unit step. Compute ITAE.

Exam-style question

E1. Structure of slides 30–34: controller K/(s+p)K/(s+p), process 1/s1/s, unity feedback. Specs: P.O.20%P.O.\le 20\%, Ts2T_s\le 2 sec, Tp<3T_p<3 sec.

  1. Translate each spec into a pole-location constraint.
  2. Describe the feasible region.
  3. Pick a specific pole location and find KK and pp.
  4. Verify all three specs.

Answer key

C1. ζ=1\zeta=1 is the critically damped case — the system reaches steady state as fast as possible without oscillation; no overshoot.

C2. Move the poles radially outward along the same angle θ=cos1ζ\theta=\cos^{-1}\zeta: ζ\zeta (and P.O.P.O.) stay the same, while ωn\omega_n grows so σd=ζωn\sigma_d=\zeta\omega_n grows and Ts=4/σdT_s=4/\sigma_d shrinks. (Equivalently, move them purely to the left would shrink TsT_s but would also change the angle and hence P.O.P.O.; the radial-outward motion is the one that preserves P.O.P.O..)

C3. ITAE multiplies the error by tt, so the large unavoidable error during the transient is de-weighted while persistent later error is penalized. ISE and IAE treat early and late error equally; minimizing ITAE tends to produce responses that settle cleanly.

C4. Step: ess=0e_{ss}=0 (the integrator in G(s)G(s) eliminates step error). Ramp: ess=1/Ke_{ss}=1/K (finite, non-zero).

C5. If s3>10Re(s1,2)|s_3|>10\,|\mathrm{Re}(s_{1,2})|, the pole s3s_3 decays much faster than the dominant pair and can be dropped. Preserve the DC gain by evaluating the dropped factor at s=0s=0 (divide numerator accordingly).


Q1. T(s)=9s2+3s+9T(s)=\dfrac{9}{s^{2}+3s+9}. ωn2=9ωn=3\omega_n^{2}=9\Rightarrow\omega_n=3; 2ζωn=3ζ=0.52\zeta\omega_n=3\Rightarrow\zeta=0.5; ωd=30.75=2.598\omega_d=3\sqrt{0.75}=2.598. Tp=π/2.598=1.21T_p=\pi/2.598=1.21 s; P.O.=100e0.5π/0.866=16.3%P.O.=100e^{-0.5\pi/0.866}=16.3\%; Ts=4/1.5=2.67T_s=4/1.5=2.67 s.

Q2. P.O.10%ζ0.591P.O.\le 10\%\Rightarrow \zeta\ge 0.591. Ts2σd2T_s\le 2\Rightarrow \sigma_d\ge 2. Feasible region: left of the vertical line σ=2\sigma=-2, inside the radial lines at θcos1(0.591)53.8\theta\le\cos^{-1}(0.591)\approx 53.8^{\circ}.

Q3. T(s)=4/(s+2)2T(s)=4/(s+2)^{2} (this is ζ=1,ωn=2\zeta=1,\omega_n=2). E(s)=1s ⁣[14(s+2)2]=s2+4ss(s+2)2=s+4(s+2)2E(s)=\tfrac{1}{s}\!\left[1-\tfrac{4}{(s+2)^{2}}\right]=\tfrac{s^{2}+4s}{s(s+2)^{2}}=\tfrac{s+4}{(s+2)^{2}}. Partial-fraction: As+2+B(s+2)2A=1,B=2\tfrac{A}{s+2}+\tfrac{B}{(s+2)^{2}}\Rightarrow A=1, B=2, so e(t)=e2t+2te2te(t)=e^{-2t}+2te^{-2t}. Then ITAE=0t(e2t+2te2t)dt=te2t+2t2e2t\text{ITAE}=\int_{0}^{\infty}t(e^{-2t}+2te^{-2t})dt=\int te^{-2t}+2\int t^{2}e^{-2t}. Using 0tneatdt=n!/an+1\int_{0}^{\infty}t^{n}e^{-at}dt=n!/a^{n+1}: first = 1/41/4; second = 2(2/8)=1/22\cdot(2/8)=1/2. Total ITAE=3/4\boxed{\text{ITAE}=3/4}.

E1. (a) P.O.20%ζ0.456θ62.9P.O.\le 20\%\Rightarrow\zeta\ge 0.456\Rightarrow\theta\le 62.9^{\circ}. Ts2σd2T_s\le 2\Rightarrow\sigma_d\ge 2. Tp<3ωd>π/31.047T_p<3\Rightarrow\omega_d>\pi/3\approx 1.047. (b) Region: left of σ=2\sigma=-2, above ωd=1.047\omega_d=1.047, inside the radial cone at ±62.9\pm 62.9^{\circ}. (c) Choose s1,2=2±j2s_{1,2}=-2\pm j2: ωn=22\omega_n=2\sqrt{2}, ζ=1/2=0.707\zeta=1/\sqrt{2}=0.707, so K=ωn2=8K=\omega_n^{2}=8, p=2ζωn=4p=2\zeta\omega_n=4. (d) Ts=4/2=2T_s=4/2=2 ✓; ωd=2>1.047Tp=π/2=1.57<3\omega_d=2>1.047\Rightarrow T_p=\pi/2=1.57<3 ✓; P.O.=100e0.707π/0.707=100eπ4.3%20%P.O.=100 e^{-0.707\pi/0.707}=100 e^{-\pi}\approx 4.3\%\le 20\% ✓.