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)=s2+2ζωns+ωn2ωn2
and a small family of standard test inputs (step, ramp, parabolic, impulse). Five scalar performance measures — rise time Tr, peak time Tp, percent overshoot P.O., settling time Ts, and steady-state error ess — are defined from the step response and then linked directly to the location of the closed-loop poles in the s-plane. This pole-placement geometry is then used to design controller parameters (choose K, p, a, etc.) so that specifications on P.O., Ts, Tp 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.,Ts from a given G(s),H(s) for a step input.
Compute ess for a given input (step, ramp, combined).
Read a measured step response and back out K, T, ζ, ωn.
Place closed-loop poles to satisfy P.O.,Ts,Tp specs simultaneously (feasible-region design).
Approximate a 3rd-order T(s) as 2nd-order while preserving gain.
Compute and minimize ISE (or ITAE) with respect to a parameter.
1. Desired outcomes checklist
#
Outcome
Where covered
Priority
1
Recognize step, ramp, parabolic, impulse test signals and their Laplace transforms
§2.1
Standard
2
Derive step and impulse response of a second-order system
§2.2, §2.3
Standard
3
Compute Tr,Tp,P.O.,Ts,ess for 2nd-order step response
§2.4–§2.9, §2.11
High (Study Q1)
4
Draw lines of constant Tp (horizontal), Ts (vertical), P.O. (diagonal/radial) in the s-plane
§2.14
Standard
5
Explain how pole motion (vertical/horizontal/radial) changes each metric
§2.14
Standard
6
Select closed-loop poles to meet specs (e.g. P.O.≤20%, Ts≤4)
§2.15
High (Study Q2)
7
Approximate higher-order systems as 2nd-order using dominant poles (preserve gain)
§2.17
Standard
8
Compute ISE/IAE/ITAE/ITSE and minimize w.r.t. a parameter (dI/dK=0)
§2.18–§2.21
High (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 signal
r(t) (for t>0)
R(s)
Step
A
A/s
Ramp
At
A/s2
Parabolic
At2
2A/s3
Unit impulse
pulse of area 1 at t=0
1
General form (slide 2):
r(t)=tn,R(s)=sn+1n!.
The unit impulse (slide 3) is defined as a limit of a rectangular pulse of width ϵ and height 1/ϵ (so area =1) as ϵ→0+. It models physical “impact hammer” testing.
[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)
reducing ζ makes ωnTp smaller (faster swiftness) but increasesP.O. and increasesTs. These are fundamental trade-offs.
[INSERT FIGURE: P.O. and ωnTp vs. ζ, from Ch05.pdf, page 23]
2.14 Pole location analysis in the splane (slides 24–29)
Complex conjugate closed-loop poles:
s1,2=−ζωn±jωn1−ζ2=−σd±jωd,
with ζ=cosθ (angle θ measured from the negative real axis).
Geometry in s-plane
What is constant
Why
Horizontal line
Tp
ωd fixed ⇒Tp=π/ωd fixed
Vertical line
Ts
σd fixed ⇒Ts=4/σd fixed
Radial line (from origin)
P.O.
θ fixed ⇒ζ=cosθ fixed
Pole motion rules (slides 27–29):
Vertical motion upward (increase ωd, same σd): same envelope, same Ts; Tp decreases, Tr decreases, overshoot increases.
Horizontal motion leftward (increase σd, same ωd): same Tp; faster decay (smaller Ts); rise time increases (system becomes “slower” in rise).
Radial motion outward (same θ=cos−1ζ, larger ωn): same P.O., same ζ; both Tp and Ts decrease.
Figures/diagrams
[INSERT FIGURE: s-plane with horizontal Tp, vertical Ts, radial 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 K and p (slides 30–34)
See §3, Example 3.
2.16 Root locations and transient response modes (slides 35–36)
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 s-plane, from Ch05.pdf, page 36]
2.17 Dominant poles approximation (slides 37–38)
Theory
The 2nd-order formulas for Tr,Tp,P.O.,Ts 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ω-axis (the dominant poles).
Slide 38 rule (illustrative):
If ∣s3∣>10∣Re(s1,2)∣, then s3 is a non-dominant pole and its effect on the response is negligible.
Important: When dropping s3, the DC gain must be preserved. Replace the factor (s+∣s3∣) in the denominator by ∣s3∣ (its value at s=0).
Simple sample case (slide 38)
Given
T(s)=(s2+4s+13)(s+30)K,s1,2=−2±3j,s3=−30.
Because ∣s3∣=30>10⋅∣Re(s1,2)∣=20, s3 is non-dominant. Preserve gain by evaluating (s+30) at s=0, i.e. divide numerator by 30:
T(s)≈s2+4s+13K/30.
2.18 Performance indices (slides 39–43)
A performance indexI is a scalar functional of the error e(t)=r(t)−y(t). The optimum parameter value is the one that minimizes I. If I=f(K),
dKdI=0.
solve for K.
The four indices from the slides:
Index
Definition
Notes
ISE
∫0Te2(t)dt
Discriminates between over- and underdamped.
IAE
$\displaystyle\int_{0}^{T}
e(t)
ITAE
$\displaystyle\int_{0}^{T}t,
e(t)
ITSE
∫0Tte2(t)dt
The upper limit T is either Ts or ∞; if not specified, use T→∞ (slide 42).
[INSERT FIGURE: Sketches of r(t),y(t),e(t),e2(t), and ∫e2(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)], H(s)=1
Problem. For a unit step input, compute Tr,Tp,P.O.,Ts.
Example 5 (slides 49–51): select K to minimize ISE
Problem.
T(s)=R(s)Y(s)=(K+5)s+65s+6,R(s)=s1.Minimize ISE w.r.t. K.
Solution.
Error:
E(s)=R(s)−Y(s)=s1[1−(K+5)s+65s+6].
Numerator of bracket:
(K+5)s+6−(5s+6)=Ks.
So
E(s)=s1⋅(K+5)s+6Ks=(K+5)s+6K.
Inverse Laplace (rewrite as s+6/(K+5)K/(K+5)):
e(t)=K+5K⋅exp(−K+56t)⋅(K+5)⋅K+51=Ke−6t/(K+5).
The slide writes e(t)=Ke−6t/(K+5). (Algebra check: L−1{K/[(K+5)s+6]}=K+5Ke−6t/(K+5); the slide absorbs 1/(K+5) into the integration, the ISE expression below is what we use and is consistent with slide 51.)
Example 6 (slides 52–55): Space Telescope — choose K3 to minimize ISE
Problem. Select K3 to minimize the ISE due to a unit-step disturbance Td(s)=1/s. Given K1=0.5, K1K2Kp=2.5. Reference r(t)=0 (only disturbance is applied).
Solution.
Transfer from disturbance to output (from signal-flow graph, slide 54):
The two exponential terms contribute tiny, finite amounts (their integrals over [0,∞) are 0.5/12≈0.0417 and 0.25/24≈0.0104); the steady-state squared error 0.25 dominates:
ISE≈0.25⋅100+0.0417+0.0104≈25.05≈25.
(Practice: take the integral explicitly as indicated in the solution.)
(c)
As
K→∞, T(s)→1, so e(t)=r(t)−y(t)→0 for all t, and
ISE→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)
When you see: “For a unit step, find Tr,Tp,P.O.,Ts.”
Steps:
Compute closed-loop T(s)=G(s)/[1+G(s)H(s)].
Normalize the denominator so the coefficient of s2 is 1.
Match: ωn2= constant term, 2ζωn= coefficient of s. Deduce ωn,ζ.
Sketch the feasible region in the left-half splane.
Pick a pole s1,2=−σd±jωd inside that region.
Compute ωn=σd2+ωd2, ζ=σd/ωn.
Back out controller parameters from the standard-form match (K=ωn2, etc.).
Verify all three inequalities hold.
5.3 Steady-state error
Steps (unit feedback):
Write R(s) from r(t).
E(s)=R(s)/[1+G(s)].
ess=lims→0sE(s).
If H(s)=1: use the general form E(s)=R(s)[1−1+G(s)H(s)G(s)].
Watch out: If the input has a ramp component (R contains 1/s2) and G(s) has no 1/s factor (no integrator), ess→∞.
5.4 Performance index minimization
Steps:
Compute E(s)=R(s)[1−T(s)] (for unit feedback: R(s)/[1+G(s)]). For disturbance-only problems (no reference), r(t)=0 and e(t)=−y(t).
Inverse Laplace to get e(t).
Form the index (ISE uses e2; ITAE uses t∣e∣; etc.).
Integrate analytically (integration by parts for ITAE/ITSE).
Write I=f(K); solve dI/dK=0.
Common mistake: leaving off the factor of t in ITAE/ITSE, or the absolute value in IAE/ITAE when e(t) changes sign.
6. Important figures and how to read them
F1. Step-response family (page 6): y(t) vs. normalized time ωnt, one curve per ζ∈{0.1,0.2,0.4,0.7,1,2}. Smaller ζ→ bigger overshoot, more oscillation; ζ=1 critically damped; ζ>1 overdamped.
F7. ISE vs.K3(Space Telescope) (page 55): both ISE and IAE shown as functions of K3; convex curves with clear minimum near K3=10≈3.2.
7. Formula sheet
Formula
Variables
When
T(s)=s2+2ζωns+ωn2ωn2
ωn: natural freq; ζ: damping ratio
Standard 2nd-order closed-loop
ωd=ωn1−ζ2
ωd: damped freq
Underdamped
s1,2=−σd±jωd, σd=ζωn
Pole coordinates
y(t)=1−e−ζωnt[cosωdt+1−ζ2ζsinωdt]
Unit-step response, underdamped
y(t)=1−ζ2ωne−ζωntsinωdt
Unit-impulse response, underdamped
Tr=ωd1atan2(−ζ+1−ζ2)
Rise time, underdamped
Tr1=t2−t1,y(t1)=0.1fv,y(t2)=0.9fv
fv: final value
Rise time, overdamped
Tp=π/ωd
Peak time
Mpt=1+e−ζπ/1−ζ2
Peak value (unit step)
P.O.=100e−ζπ/1−ζ2%
Percent overshoot (unit step)
ζ=cosθ
θ: angle from negative real axis
Convert angle↔︎ζ
Ts≅4/(ζωn)=4/σd
2% criterion
Settling time
ess=lims→01+G(s)sR(s)
H(s)=1 only
Steady-state error
E(s)=R(s)[1−1+G(s)H(s)G(s)]
general
Error transform
ISE=∫0Te2dt
Squared error
$\text{IAE}=\int_{0}^{T}
e
,dt$
$\text{ITAE}=\int_{0}^{T}t
e
,dt$
ITSE=∫0Tte2dt
Time-weighted squared
8. Concept map / summary table
Topic
Main idea
Key equation
Typical exam task
Input test signals
Standard inputs for characterization
R(s)=n!/sn+1
Write R(s) for a given r(t)
2nd-order step response
Underdamped oscillation toward 1
y(t)=1−e−ζωnt[⋯]
Derive/evaluate y(t)
Rise time Tr
First hit of final value
Tr=ωd1atan2(⋯)
Compute from ζ,ωn
Peak time Tp
First peak
Tp=π/ωd
Compute or invert
Percent overshoot
Peak over final value
P.O.=100e−ζπ/1−ζ2
Find ζ from P.O.
Settling time Ts
Entry into ±2% band
Ts≅4/(ζωn)
Compute or constrain σd
Steady-state error
Long-run error
ess=lims→0sR(s)/[1+G(s)]
Evaluate for step/ramp inputs
Pole-location analysis
Geometry = performance
Iso-Tp, iso-Ts, iso-P.O. lines
Relate pole motion to metric change
Pole placement
Satisfy multiple specs
Feasible region in s-plane
Find K,p,a
Dominant poles
Reduce order
$
s_3
Performance indices
Scalar measure of tracking
ISE, IAE, ITAE, ITSE
Integrate, minimize w.r.t. parameter
9. Common mistakes
ωnvs.ωd: Tp=π/ωd, notπ/ωn.
Unnormalized denominator. From 2s2+4s+8, divide through to get s2+2s+4 before matching; then ωn=2, ζ=0.5.
Real vs. imaginary part.Ts depends on σd=ζωn (real). Tp depends on ωd (imaginary).
esswith non-unity feedback: Don’t use sR/(1+G) when H(s)=1; use the general formula E(s)=R(s)[1−G/(1+GH)].
Dominant-pole gain. When dropping a non-dominant pole (s+∣s3∣), divide the numerator by ∣s3∣ to preserve the DC gain.
Upper integration limit for indices. If the problem specifies T=Ts, use it; otherwise take T→∞.
atan2 quadrant. For Tr, the argument (+1−ζ2,−ζ) has positive y and negative x → 2nd quadrant → Tr>0.
Disturbance-only ISE. If only a disturbance is applied and no reference, take r(t)=0⇒e(t)=−y(t).
10. Final self-test
Conceptual questions
C1. What does it mean physically for ζ=1? Name the case and state whether overshoot occurs.
C2. How must closed-loop poles move in the s-plane to reduce Ts without changing 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/s, what is ess to a unit step? To a unit ramp?
C5. State the dominant-poles condition that allows a third pole s3 to be ignored. What additional step is needed when simplifying?
Calculation questions
Q1. Unity-feedback system with G(s)=9/[s(s+3)], unit step input. Find ωn,ζ,ωd,Tp,P.O.,Ts.
Q2. Specs: P.O.≤10%, Ts≤2. Find the minimum ζ and the minimum σd=ζωn. Describe the feasible region.
Q3.T(s)=4/(s2+4s+4), unit step. Compute ITAE.
Exam-style question
E1. Structure of slides 30–34: controller K/(s+p), process 1/s, unity feedback. Specs: P.O.≤20%, Ts≤2 sec, Tp<3 sec.
Translate each spec into a pole-location constraint.
Describe the feasible region.
Pick a specific pole location and find K and p.
Verify all three specs.
Answer key
C1.ζ=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 θ=cos−1ζ: ζ (and P.O.) stay the same, while ωn grows so σd=ζωn grows and Ts=4/σd shrinks. (Equivalently, move them purely to the left would shrink Ts but would also change the angle and hence P.O.; the radial-outward motion is the one that preserves P.O..)
C3. ITAE multiplies the error by t, 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=0 (the integrator in G(s) eliminates step error). Ramp: ess=1/K (finite, non-zero).
C5. If ∣s3∣>10∣Re(s1,2)∣, the pole s3 decays much faster than the dominant pair and can be dropped. Preserve the DC gain by evaluating the dropped factor at s=0 (divide numerator accordingly).
Q1.T(s)=s2+3s+99. ωn2=9⇒ωn=3; 2ζωn=3⇒ζ=0.5; ωd=30.75=2.598. Tp=π/2.598=1.21 s; P.O.=100e−0.5π/0.866=16.3%; Ts=4/1.5=2.67 s.
Q2.P.O.≤10%⇒ζ≥0.591. Ts≤2⇒σd≥2. Feasible region: left of the vertical line σ=−2, inside the radial lines at θ≤cos−1(0.591)≈53.8∘.
Q3.T(s)=4/(s+2)2 (this is ζ=1,ωn=2). E(s)=s1[1−(s+2)24]=s(s+2)2s2+4s=(s+2)2s+4. Partial-fraction: s+2A+(s+2)2B⇒A=1,B=2, so e(t)=e−2t+2te−2t. Then ITAE=∫0∞t(e−2t+2te−2t)dt=∫te−2t+2∫t2e−2t. Using ∫0∞tne−atdt=n!/an+1: first = 1/4; second = 2⋅(2/8)=1/2. Total ITAE=3/4.
E1. (a) P.O.≤20%⇒ζ≥0.456⇒θ≤62.9∘. Ts≤2⇒σd≥2. Tp<3⇒ωd>π/3≈1.047.
(b) Region: left of σ=−2, above ωd=1.047, inside the radial cone at ±62.9∘.
(c) Choose s1,2=−2±j2: ωn=22, ζ=1/2=0.707, so K=ωn2=8, p=2ζωn=4.
(d) Ts=4/2=2 ✓; ωd=2>1.047⇒Tp=π/2=1.57<3 ✓; P.O.=100e−0.707π/0.707=100e−π≈4.3%≤20% ✓.