Control Systems
Second-Order System Response
Natural frequency, damping ratio, and the shape of every step response
A second-order system response is the transient behavior of any plant with two energy-storage elements, governed by the standard transfer function G(s)=ωn²/(s²+2ζωn s+ωn²). Just two numbers decide everything: the undamped natural frequency ωn, in radians per second, which sets how fast the system reacts, and the dimensionless damping ratio ζ, which sets the shape. When ζ is below 1 the system is underdamped and overshoots — a ζ of 0.5 overshoots by about 16% — when ζ equals 1 it is critically damped and reaches the target as fast as possible with no overshoot, and when ζ exceeds 1 it is overdamped and slow. From ζ and ωn alone you can predict percent overshoot, peak time, rise time, and the 2% settling time Ts≈4/(ζωn) in closed form. The same equations describe an RLC circuit, a spring-mass-damper, a phase-locked loop, and virtually every well-designed feedback controller.
- Standard formωn²/(s²+2ζωn s+ωn²)
- Poless = −ζωn ± jωn√(1−ζ²)
- Overshoot%OS = 100·e^(−ζπ/√(1−ζ²))
- 2% settlingTs ≈ 4 / (ζωn)
- Peak timeTp = π / ωn√(1−ζ²)
- Design sweet spotζ ≈ 0.7 → ~4.3% OS
Interactive visualization
Press play, or step through manually. The visualization is yours to drive — try it before reading on.
Watch the 60-second explainer
A condensed visual walkthrough — narrated, captioned, under a minute.
Why second-order response matters
Almost every dynamic system an engineer designs boils down, near its dominant behavior, to a second-order model. The reason is physical: a second-order system is one that stores energy in two independent places and trades that energy back and forth. A mass stores kinetic energy while a spring stores potential energy; an inductor stores energy in its magnetic field while a capacitor stores it in its electric field. Whenever two such reservoirs are coupled, the system can oscillate, and the mathematics is identical whether the medium is mechanical, electrical, thermal, or hydraulic.
- Control loops. A well-tuned closed loop is deliberately shaped to look second-order with ζ≈0.5–0.8, so overshoot and settling time are predictable.
- RLC circuits. A series RLC has ωn = 1/√(LC) and ζ = (R/2)·√(C/L); it is the electrical twin of the spring-mass-damper.
- Vehicle suspension. A quarter-car spring and shock absorber target ζ≈0.25–0.35 for ride comfort, trading a little overshoot for a soft response.
- Servo and motion systems. Position loops are tuned so the closed-loop poles form a dominant second-order pair setting bandwidth and overshoot.
- Phase-locked loops. A PLL's loop filter makes the phase error a second-order system, with ζ and ωn setting lock time and jitter.
- Instrumentation. Accelerometers, galvanometers, and pressure sensors are second-order; their damping sets bandwidth and ringing.
Because one model covers so much ground, the closed-form formulas for overshoot, rise time, and settling time are among the most-used equations in all of engineering. Learn them once and you can specify a suspension, a filter, and a robot joint with the same arithmetic.
How it works, step by step
The standard second-order transfer function relating output Y(s) to input X(s) is
G(s) = Y(s)/X(s) = ωn² / (s² + 2ζωn·s + ωn²)
where the symbols are:
- s — the Laplace variable, in rad/s (a complex frequency).
- ωn — undamped natural frequency, in rad/s. It is the frequency at which the system would ring if there were no damping.
- ζ (zeta) — damping ratio, dimensionless. The ratio of actual damping to the critical damping that just prevents oscillation.
The DC gain here is 1 (set numerator to ωn² and s→0), so the output settles to the input value; a physical gain K simply multiplies the whole thing. The character of the response comes from the roots of the denominator — the poles:
s = −ζωn ± jωn·√(1 − ζ²)
Read those poles geometrically. Their distance from the origin is always ωn. Their horizontal position (real part) is σ = −ζωn, which controls decay. Their vertical position (imaginary part) is the damped natural frequency ωd = ωn·√(1−ζ²), the frequency the system actually rings at. The angle β = arccos(ζ) that the pole makes with the negative real axis is set entirely by ζ. So ζ steers the pole angle and ωn its radius.
- Apply a unit step. The input jumps from 0 to 1 and stays there. The system must move from its old equilibrium to the new one.
- The poles decide the path. Complex poles (ζ<1) mean the output overshoots and rings at ωd inside a decaying envelope e^(−ζωn·t). Real poles (ζ≥1) mean a smooth, non-oscillating climb.
- The envelope sets settling. The oscillation amplitude is squeezed by e^(−ζωn·t); once it drops below the tolerance band, the system is "settled." Since σ = ζωn, farther-left poles settle faster.
- ζ sets overshoot, ωn sets speed. Change ωn and the whole waveform stretches or compresses in time but keeps its shape and overshoot. Change ζ and the shape itself changes.
For the underdamped case, the closed-form unit-step response is y(t) = 1 − e^(−ζωn·t)·[cos(ωd·t) + (ζ/√(1−ζ²))·sin(ωd·t)], which peaks above 1 and rings down to 1.
The three damping regimes
The single parameter ζ divides all second-order behavior into distinct regimes. The table below lists the pole structure and step-response character for each, with representative time-domain metrics for a fixed ωn.
| Regime | ζ range | Poles | Step response | Overshoot |
|---|---|---|---|---|
| Undamped | ζ = 0 | ±jωn on imaginary axis | Rings forever, never settles | 100% |
| Underdamped | 0 < ζ < 1 | Complex conjugate pair | Overshoots, then rings down | 0–100% |
| Critically damped | ζ = 1 | Repeated real pole at −ωn | Fastest with no overshoot | 0% |
| Overdamped | ζ > 1 | Two distinct real poles | Slow, sluggish, no overshoot | 0% |
The most important design insight hides in this table: critically damped is the fastest response that does not overshoot. Beginners often assume "more damping is safer, so crank ζ up." But once ζ passes 1, a slow real pole appears and dominates, so the overdamped system is slower than the critically damped one while gaining nothing. Meanwhile a slightly underdamped design at ζ≈0.7 rises faster than either and pays only about 4% overshoot for it.
Transient metrics and where the numbers come from
For an underdamped second-order system responding to a step, four metrics fully characterize the transient. Each follows directly from ζ and ωn:
- Percent overshoot: %OS = 100·e^(−ζπ/√(1−ζ²)). Depends on ζ only — not on ωn.
- Peak time: Tp = π/ωd = π/(ωn√(1−ζ²)). Time to the first peak.
- Rise time (0–100%): Tr = (π − β)/ωd, with β = arccos(ζ). A handy 10–90% approximation is Tr ≈ 1.8/ωn near ζ=0.5.
- Settling time (2%): Ts ≈ 4/(ζωn); the 5% band gives Ts ≈ 3/(ζωn).
Notice the clean division of labor. Overshoot is a function of ζ alone. Settling time depends on the product ζωn, which is exactly σ, the real part of the poles. Peak and rise time scale as 1/ωn. This is why the design recipe is simple: pick ζ from the overshoot you can tolerate, then pick ωn from the speed you need. The table below shows the overshoot that different damping ratios produce.
| Damping ratio ζ | Percent overshoot | Regime / note |
|---|---|---|
| 0.1 | ≈ 73% | Very lightly damped, long ringing |
| 0.3 | ≈ 37% | Typical vehicle suspension range |
| 0.5 | ≈ 16% | Fast but bouncy |
| 0.6 | ≈ 9.5% | Common servo target |
| 0.707 | ≈ 4.3% | Butterworth pole angle, classic sweet spot |
| 0.8 | ≈ 1.5% | Well damped, slower rise |
| 1.0 | 0% | Critically damped |
Worked example: designing a positioning stage
Suppose a linear positioning stage must move to a commanded position, overshoot no more than 5%, and settle to within 2% of the target in 0.5 s. Design the target closed-loop second-order dynamics.
- Pick ζ from the overshoot spec. Invert the overshoot formula: ζ = −ln(0.05)/√(π² + ln²(0.05)). Here ln(0.05) = −2.996, so ζ = 2.996/√(9.87 + 8.98) = 2.996/4.34 ≈ 0.69. Round to ζ = 0.7, giving about 4.6% overshoot — safely under 5%.
- Pick ωn from the settling spec. Require Ts = 4/(ζωn) ≤ 0.5 s. Solve ωn ≥ 4/(ζ·Ts) = 4/(0.7·0.5) = 11.4 rad/s. Choose ωn = 12 rad/s for margin.
- Check the other metrics. ωd = ωn√(1−ζ²) = 12·√(1−0.49) = 12·0.714 = 8.57 rad/s. Peak time Tp = π/8.57 = 0.37 s. Rise time Tr = (π − arccos(0.7))/8.57 = (3.1416 − 0.795)/8.57 = 0.27 s.
- Place the poles. s = −ζωn ± jωd = −8.4 ± j8.57. Any controller — a PD, a lead compensator, or a state-feedback gain — that drives the dominant closed-loop poles to this location meets the spec.
This four-line procedure is the everyday workflow of control design: translate time-domain specs into a target ζ and ωn, convert those into a desired pole location, then choose gains that put the poles there. Tools such as the root locus and Bode plot exist precisely to help you steer the dominant poles to that spot.
Common misconceptions and failure modes
- "More damping is always safer." Past ζ=1 the response only gets slower, not steadier. Overdamped ≠ better.
- "Overshoot depends on how fast you drive it." No — overshoot is a function of ζ alone; ωn changes only timing, not the overshoot percentage.
- "ωn is the frequency I'll observe." The ringing frequency is ωd = ωn√(1−ζ²), always below ωn; only for ζ=0 do they coincide.
- "A nearby zero doesn't matter." A left-half-plane zero close to the poles increases overshoot and speeds the rise; the pole-only formulas underestimate it.
- "Any high-order system reduces to second order." Only if the non-dominant poles/zeros are roughly 5–10× farther left; otherwise the approximation breaks down.
- "Settling time scales with ωn." It scales with the product ζωn (the real pole part σ), so a high-ωn but very lightly damped system can settle slowly.
Frequently asked questions
What is a second-order system response?
It is the output that a system with two energy-storage elements produces when disturbed, governed by the standard transfer function G(s)=ωn²/(s²+2ζωn s+ωn²). The response is fully determined by two parameters: the natural frequency ωn in rad/s, which sets how fast it reacts, and the damping ratio ζ, which sets its shape. A spring-mass-damper, an RLC circuit, and a tuned servo loop all share this same second-order form, so the same overshoot, rise-time, and settling-time formulas apply to every one of them.
What do the damping ratio and natural frequency mean?
The natural frequency ωn is the angular frequency at which the system would oscillate with no damping, in rad/s; it scales the whole time axis, so doubling ωn halves rise time and settling time without changing the shape. The damping ratio ζ is dimensionless and sets the shape: ζ=0 is a pure undamped oscillator, 0<ζ<1 rings and decays, ζ=1 is the fastest non-overshooting response, and ζ>1 is sluggish with no ring. For an underdamped system the poles sit at s=−ζωn ± jωn√(1−ζ²), so ζ fixes their angle off the imaginary axis and ωn fixes their distance from the origin.
How do you calculate percent overshoot from the damping ratio?
For an underdamped step response, percent overshoot depends only on ζ, not on ωn: %OS = 100·exp(−ζπ/√(1−ζ²)). A damping ratio of ζ=0.5 gives about 16% overshoot, ζ=0.707 gives about 4.3%, and ζ=0.6 gives about 9.5%. Overshoot reaches zero only as ζ approaches 1. Inverting the formula lets you pick the ζ needed for a target overshoot: ζ = −ln(%OS/100) / √(π² + ln²(%OS/100)).
What is the difference between underdamped, critically damped, and overdamped?
Underdamped (0<ζ<1) has complex-conjugate poles, so it overshoots the target and rings before settling. Critically damped (ζ=1) has a repeated real pole at −ωn and gives the fastest possible approach with no overshoot at all. Overdamped (ζ>1) has two distinct real poles; it never overshoots but is slower than the critical case because the slower pole dominates. Undamped (ζ=0) has poles on the imaginary axis and oscillates forever. In control design ζ≈0.7 is a common sweet spot, trading a small 4-5% overshoot for a fast rise time.
How do you find the settling time of a second-order system?
For an underdamped system the response is bounded by an envelope that decays as exp(−ζωn t), so the settling time is governed by the real part of the poles, σ = ζωn. The 2% settling time is Ts ≈ 4/(ζωn) and the 5% settling time is Ts ≈ 3/(ζωn). Because σ = ζωn is the horizontal distance of the poles from the imaginary axis, pushing the poles farther left speeds up settling. Note that settling time depends on the product ζωn, whereas overshoot depends on ζ alone.
What are the rise time and peak time formulas?
The peak time, when the underdamped response first reaches its maximum, is Tp = π/ωd where ωd = ωn√(1−ζ²) is the damped natural frequency. The 0-to-100% rise time is Tr = (π − β)/ωd where β = arccos(ζ). A common engineering approximation for the 10-90% rise time near ζ=0.5 is Tr ≈ 1.8/ωn. All three quantities scale inversely with ωn, so raising the natural frequency makes the system faster without altering its overshoot.
Why do second-order approximations work for higher-order systems?
Most real plants are higher order, but their step response is often dominated by a pair of complex poles closest to the imaginary axis. If the remaining poles and any zeros are at least about 5-10 times farther left, they decay quickly and contribute little, so the system behaves like a second-order one with the dominant-pole ζ and ωn. This dominant-pole approximation lets engineers use the simple overshoot and settling-time formulas for design. A nearby zero, however, can increase overshoot and speed rise time, and a right-half-plane zero causes non-minimum-phase undershoot that the standard formulas miss.