Transient and Steady-State Response

Set up the classic lab demonstration: a mass on a spring, a signal generator, a shaker, an oscilloscope. Switch on a sinusoidal drive and watch the trace. It does not immediately oscillate at a fixed size. For the first seconds the motion writhes — two rhythms fighting, the amplitude swelling and shrinking as if the system cannot make up its mind. Then, gradually, the struggle fades and the trace locks into a clean, constant sinusoid at the drive frequency, and stays there forever. Run it again with a completely different launch — pull the mass aside first, or fling it — and the early writhing changes completely, but the final oscillation is identical. That two-act structure is the subject of this article. The general solution of the driven, damped oscillator is the sum of two pieces with utterly different jobs: a steady state that oscillates at the drive frequency with amplitude and phase the system imposes, and a transient — a free, damped oscillation at the system's own ringing frequency — that decays as e−βte^{-\beta t} and carries every trace of the initial conditions to the grave. The steady state is where the motion is going; the transient is how it gets there, and it is the only part of the solution you are allowed to adjust when matching how the motion started. Mastering the split buys you real laboratory power: you will compute exactly how a resonator builds up from rest, why the startup shows beat-like interference, and why a high-QQ oscillator — a quartz fork, a superconducting readout cavity — takes roughly Q/πQ/\pi cycles to settle, which is the hidden speed limit on every precision resonant measurement.

Key takeaways
  • The general solution of x¨+2βx˙+ω02x=(F0/m)cos⁡ωt\ddot x + 2\beta\dot x + \omega_0^2 x = (F_0/m)\cos\omega t is x=xp+xhx = x_{\text{p}} + x_{\text{h}}: a particular (steady-state) piece plus a homogeneous (transient) piece.
  • The steady state xp=Acos⁡(ωt−δ)x_{\text{p}} = A\cos(\omega t - \delta) is fixed entirely by the system and the drive — its A(ω)A(\omega) and δ(ω)\delta(\omega) are derived in Driven Oscillations and contain no free constants.
  • The transient xh=e−βt(C1cos⁡ω1t+C2sin⁡ω1t)x_{\text{h}} = e^{-\beta t}(C_1\cos\omega_1 t + C_2\sin\omega_1 t) decays as e−βte^{-\beta t} and its two constants C1,C2C_1, C_2 absorb the initial position and velocity — never touch AA or δ\delta to fit initial conditions.
  • Started from rest, the amplitude builds toward the steady state as 1−e−βt1 - e^{-\beta t} on resonance; off resonance the transient (at ω1\omega_1) interferes with the drive (at ω\omega), producing decaying beat-like ringing.
  • Settling takes a time ∼1/β\sim 1/\beta, which is about Q/πQ/\pi oscillation cycles to reach 63% — high-QQ systems respond sharply but slowly, and every resonant instrument must wait out its own memory.
  • Every launch converges to the same steady state: dissipation erases initial-condition memory exponentially.
  • Key idea: Linearity splits the driven motion into transient + steady state: x(t)=xh(t)+xp(t)x(t)=x_h(t)+x_p(t), where xhx_h decays and absorbs the initial conditions, and xpx_p is the persistent drive-locked response whose amplitude and phase are independent of how the motion began.

Where this lives today

The settling rule N≈Q/πN \approx Q/\pi is a live engineering constraint at the frontier of quantum computing. Reading out a superconducting qubit means driving a microwave resonator coupled to it and watching the reflected steady state; the resonator cannot report faster than it can settle, so designers deliberately spoil its quality factor — internal QQ values above a million are throttled to a few thousand of external coupling — trading sensitivity for a ring-up of tens of nanoseconds. Every 2026-era processor's measurement chain is a negotiated compromise on exactly equation (8): too high a readout QQ and the computer idles waiting for photons to build up; too low and the signal drowns. The same waiting game runs at every scale. A lock-in amplifier hunting a sharp mechanical resonance must dwell ∼Q\sim Q cycles per frequency point, which is why characterizing a Q=107Q = 10^7 microresonator can take hours of automated sweeping. Concert-hall acousticians measure the room's transient directly — the reverberation time RT60 is the ring-down of the hall's modes — and tune absorbers (damping) so music's steady states bloom and decay on musical timescales. And car suspension engineers still buy Kelvin's bargain in reverse: damping near critical, Q∼1Q \sim 1, so the transient from each pothole is forgotten within a single bounce.

The physical picture

Think of the steady state as a moving target that the motion is falling toward. The drive and the damping between them define one — and only one — self-consistent perpetual motion: the oscillation at the drive frequency in which the energy fed in per cycle exactly equals the energy friction destroys per cycle. Start the system anywhere else and it is not on that motion; the difference between where it is and where it must end up is precisely the transient. And the transient is nothing new — it is a free, damped oscillation, exactly the ring-down studied in the Damped Harmonic Oscillator article, obeying the homogeneous equation because the drive is already fully "spent" sustaining the steady state. Free damped oscillations die as e−βte^{-\beta t}; therefore memory dies as e−βte^{-\beta t}.

Figure
Startup from rest on resonance — the amplitude climbs the envelope toward the steady state
steady-state amplitude A t = 1/β — envelope at 63% of A driven at ω = ω₀ from rest (Q = 10) envelope A(1 − e^(−βt)) — the transient cancels the steady state at t = 0, then dies time x(0) = 0, v(0) = 0 — yet the final amplitude forgets this completely

The figure above is the cleanest startup there is: drive switched on at t=0t = 0 with the oscillator at rest, drive frequency right on ω0\omega_0. At t=0t = 0 the steady state alone would say the mass is already moving — so the transient must start with exactly the opposite motion to cancel it. As the transient decays, the cancellation fails at the rate e−βte^{-\beta t}, and what is left grows along the envelope A(1−e−βt)A(1 - e^{-\beta t}): the oscillator does not so much "build up" as "stop cancelling." One damping time 1/β1/\beta brings the amplitude to 63% of its final value; three damping times bring it to 95%. Off resonance the startup is more dramatic, because the two pieces now run at different frequencies — the transient rings at the damped natural frequency ω1\omega_1, the steady state at the drive's ω\omega — and their sum shows interference:

Figure
The anatomy of a startup — steady state plus decaying transient equals the observed motion
steady state x_p — constant amplitude, drive frequency ω = 1.6 ω₀ transient x_h — rings at ω₁ ≈ ω₀, decays as e^(−βt), sized by the initial conditions observed motion x = x_p + x_h — messy start, clean finish rest start, ω = 1.6 ω₀, Q = 4 — by t ≈ 3/β only the pure steady state survives

Read the third trace against the first two. Early on, the observed motion is genuinely two-frequency: it is not "the oscillator hesitating," it is the arithmetic sum of a clean drive-frequency sinusoid and a clean decaying ring at ω1\omega_1. Where the two agree in sign the motion is big; where they oppose it is small. All the apparent complexity of startup is this superposition — and since one partner is doomed, the complexity is temporary. The initial conditions choose only C1C_1 and C2C_2, the size and phase of the doomed partner.

The equations, assembled

The system is the same driven, damped oscillator whose steady state is derived in Driven Oscillations; this article is about the complete solution, so let us fix notation once:

x¨+2βx˙+ω02x=F0mcos⁡ωt(1)\ddot x + 2\beta\dot x + \omega_0^2 x = \frac{F_0}{m}\cos\omega t \qquad (1)

where: ω0=k/m\omega_0 = \sqrt{k/m} — natural frequency; β=b/2m\beta = b/2m — damping parameter; F0,ωF_0, \omega — drive amplitude and frequency. (Several references — and this article's provenance note — use γ≡2β\gamma \equiv 2\beta, so their e−γt/2e^{-\gamma t/2} decay is our e−βte^{-\beta t}: same physics, different packing of the 2.) We assume the underdamped case β<ω0\beta < \omega_0, the one with visible ringing; the overdamped and critical cases swap the transient's sinusoids for decaying exponentials but change nothing structural. The theory of linear equations dictates the shape of every solution before any physics is done. One particular solution plus the general homogeneous solution is the general solution:

x(t)=xp(t)+xh(t)(2)x(t) = x_{\text{p}}(t) + x_{\text{h}}(t) \qquad (2)

where: xpx_{\text{p}} — any one motion that satisfies the driven equation (we use the steady state); xhx_{\text{h}} — the general solution of the free equation, carrying the two arbitrary constants a second-order equation owes us. This bland-looking line is the whole conceptual content of the article: the driven problem = the free problem + one forced motion. The homogeneous piece is the free damped oscillator, solved in its own article; for β<ω0\beta < \omega_0 it rings at the damped frequency and dies exponentially:

xh(t)=e−βt(C1cos⁡ω1t+C2sin⁡ω1t),ω1=ω02−β2(3)x_{\text{h}}(t) = e^{-\beta t}\big(C_1\cos\omega_1 t + C_2\sin\omega_1 t\big), \qquad \omega_1 = \sqrt{\omega_0^2 - \beta^2} \qquad (3)

where: ω1\omega_1 — the damped natural frequency, a bit below ω0\omega_0; C1,C2C_1, C_2 — the two free constants, soon to be fixed by initial conditions. Note the three distinct frequencies now in play: the drive's ω\omega (chosen by you), the natural ω0\omega_0 (property of the system), and the ringing ω1\omega_1 (property of the system with its damping). The transient rings at ω1\omega_1 — neither of the other two. The particular piece we take straight from the Driven Oscillations article, where the complex-exponential method produces it in five lines:

xp(t)=Acos⁡(ωt−δ),A=F0/m(ω02−ω2)2+(2βω)2,tan⁡δ=2βωω02−ω2(4)x_{\text{p}}(t) = A\cos(\omega t - \delta), \qquad A = \frac{F_0/m}{\sqrt{(\omega_0^2-\omega^2)^2 + (2\beta\omega)^2}}, \qquad \tan\delta = \frac{2\beta\omega}{\omega_0^2 - \omega^2} \qquad (4)

where: A(ω),δ(ω)A(\omega), \delta(\omega) — the steady-state amplitude and phase lag, functions fixed by the system and drive with no adjustable freedom whatsoever (0≤δ≤π0 \le \delta \le \pi). This is the piece that survives, and the peak it develops near ω0\omega_0 is the Resonance article's territory. Here its only role is: the destination.

Derivation AWhy every solution is transient plus steady state

We are establishing the fundamental decomposition, equation (2), with its three claims: that the sum solves the driven equation, that every solution has this form, and that all solutions converge to the same steady state — the system forgets its initial conditions. The strategy is pure linearity: work with the operator L[x]=x¨+2βx˙+ω02xL[x] = \ddot x + 2\beta\dot x + \omega_0^2 x and let superposition do everything; no formula-grinding is needed. The subtle claim is the second one — not just that transient-plus-steady-state is a solution, but that there are no others — and the trick for it is to look at the difference of two solutions. Watch that step; it is where uniqueness of the long-time motion comes from.

Step 1. Check the sum solves the driven equation. With L[xp]=(F0/m)cos⁡ωtL[x_{\text{p}}] = (F_0/m)\cos\omega t and L[xh]=0L[x_{\text{h}}] = 0, linearity gives

L[xp+xh]=L[xp]+L[xh]=F0mcos⁡ωt+0L[x_{\text{p}} + x_{\text{h}}] = L[x_{\text{p}}] + L[x_{\text{h}}] = \frac{F_0}{m}\cos\omega t + 0

Why linearity holds: every term of LL is proportional to xx or a derivative of xx — no squares, no products — so LL distributes over sums. This is the same structural fact that powers the superposition-of-drives result in the Driven Oscillations article.

Step 2. Show there are no other solutions. Let xx be any solution of the driven equation and consider its difference from our particular one, u=x−xpu = x - x_{\text{p}}. Then

L[u]=L[x]−L[xp]=F0mcos⁡ωt−F0mcos⁡ωt=0L[u] = L[x] - L[x_{\text{p}}] = \frac{F_0}{m}\cos\omega t - \frac{F_0}{m}\cos\omega t = 0

so uu solves the free equation, i.e. uu is one of the motions (3). Every solution is therefore xpx_{\text{p}} plus some homogeneous motion — the form (2) is exhaustive. Why this matters: it forbids exotic alternatives. There is no special launch that locks the oscillator into some other perpetual motion; whatever you do, you are on the steady state plus a doomed free ring.

Step 3. Count the freedom. The constants C1,C2C_1, C_2 in (3) are exactly two — the number a second-order equation needs to match an arbitrary initial position and velocity. The steady state contributes none: AA and δ\delta in (4) are already determined functions of ω\omega. Why the bookkeeping is satisfying: the drive fixes where you end up; the two-dimensional space of initial conditions maps one-to-one onto the two-dimensional space of transients. Nothing is over- or under-determined.

Step 4. Take the long-time limit. Since ∣xh∣≤e−βtC12+C22→0|x_{\text{h}}| \le e^{-\beta t}\sqrt{C_1^2 + C_2^2} \to 0 for any C1,C2C_1, C_2,

lim⁡t→∞[x(t)−xp(t)]=0\lim_{t\to\infty}\big[x(t) - x_{\text{p}}(t)\big] = 0

— every solution, from every launch, converges to the same steady state. Why this is memory loss: the only place the initial data lived was in C1,C2C_1, C_2, and both ride the dying exponential. Dissipation is the eraser; without it (β=0\beta = 0) the free part never dies and the system remembers its launch forever. *Verified: setting F0=0F_0 = 0 collapses (2) to the free damped oscillator, and setting t≫1/βt \gg 1/\beta collapses it to the steady state of the Driven Oscillations article; structure as in Morin Ch. 4 and Taylor Ch. 5.*

Matching initial conditions — the transient's one job

Here is where students lose the most points, so let us do it completely and then numerically. Suppose at t=0t = 0 you know the position x(0)=x0x(0) = x_0 and velocity x˙(0)=v0\dot x(0) = v_0. The temptation is to bend AA or δ\delta to fit — resist it; they are not yours. Evaluate the full solution (2) and its derivative at t=0t = 0:

x0=Acos⁡δ+C1,v0=Aωsin⁡δ−βC1+ω1C2(5)x_0 = A\cos\delta + C_1, \qquad v_0 = A\omega\sin\delta - \beta C_1 + \omega_1 C_2 \qquad (5)

where: the AA-terms are the steady state's position and velocity at t=0t = 0 — already fixed numbers once you know ω\omega — and C1,C2C_1, C_2 are the unknowns. Two linear equations, two unknowns, always solvable: C1=x0−Acos⁡δC_1 = x_0 - A\cos\delta, then C2=(v0+βC1−Aωsin⁡δ)/ω1C_2 = (v_0 + \beta C_1 - A\omega\sin\delta)/\omega_1. The transient is literally the difference between how you started and how the steady state wanted you to start, packaged as a free damped oscillation.

Go deeper · Derivation B · 2 minThe start from rest and the shape of the ringing
Derivation BThe start from rest and the shape of the ringing

We are deriving the explicit motion for the most common laboratory launch — oscillator at rest, drive switched on at t=0t = 0 — and reading off why the startup shows beat-like ringing when ω≠ω0\omega \ne \omega_0. The strategy is to feed x0=0,v0=0x_0 = 0, v_0 = 0 into the matching equations (5), solve for C1,C2C_1, C_2, and assemble the full x(t)x(t); the interference interpretation then falls out by comparing the two pieces' frequencies and initial sizes. Watch the signs: both constants come out negative multiples of the steady state's coefficients, which is the algebra saying the transient begins life opposing the steady state exactly.

Step 1. Impose x(0)=0x(0) = 0 in (5): C1=−Acos⁡δC_1 = -A\cos\delta. Why the minus sign: at t=0t = 0 the steady state alone sits at Acos⁡δA\cos\delta, but the mass is actually at zero — the transient must start exactly opposite to cancel it.

Step 2. Impose x˙(0)=0\dot x(0) = 0: 0=Aωsin⁡δ−βC1+ω1C20 = A\omega\sin\delta - \beta C_1 + \omega_1 C_2, so with C1C_1 from Step 1,

C2=−A (βcos⁡δ+ωsin⁡δ)ω1C_2 = -\frac{A\,(\beta\cos\delta + \omega\sin\delta)}{\omega_1}

Why β\beta appears here: the transient's velocity has two parts — the ring's own swing and the envelope's shrinkage — and the matching must cancel the steady state's initial velocity with both included.

Step 3. Assemble the complete rest-start motion:

x(t)=Acos⁡(ωt−δ)−Ae−βt[cos⁡δ cos⁡ω1t+βcos⁡δ+ωsin⁡δω1 sin⁡ω1t](6)x(t) = A\cos(\omega t - \delta) - A e^{-\beta t}\left[\cos\delta\,\cos\omega_1 t + \frac{\beta\cos\delta + \omega\sin\delta}{\omega_1}\,\sin\omega_1 t\right] \qquad (6)

where: the first term is the destination, the bracket is the memory. Both terms have amplitude of order AA at t=0t = 0 — the transient starts as big as the steady state itself — and the bracket rides the dying exponential.

Step 4. Read off the interference. For light damping the two pieces are comparable in size for the first few 1/β1/\beta, one oscillating at ω\omega, the other at ω1≈ω0\omega_1 \approx \omega_0. Their sum drifts periodically in and out of phase at the difference frequency: the envelope swells and shrinks with period 2π/∣ω−ω1∣2\pi/|\omega - \omega_1|, fading as e−βte^{-\beta t}. Why this is "beats with a doomed partner": the arithmetic is identical to two-tone beats — the topic of the Beats article, where both tones persist — but here one tone decays, so the interference pattern is transient ringing rather than a steady throb. When ω\omega is close to ω0\omega_0 the beat period is long and the startup looks like slow, dying surges; far from ω0\omega_0 the wiggles are fast and vanish within a few drive cycles. *Verified: direct substitution confirms x(0)=0x(0) = 0 exactly, and x˙(0)=0\dot x(0) = 0 using ddt\frac{\mathrm d}{\mathrm dt} of each term at t=0t=0; equation (6) reduces to the textbook rest-start response (Taylor §5.6) with γ=2β\gamma = 2\beta.*

Figure
Transient ringing — the dying transient beats against the drive before the steady state wins
envelope surges at the difference frequency rest start, ω = 1.25 ω₀, light damping — interference of x_p (at ω) with x_h (at ω₁) the surges fade as e^(−βt); only the drive tone survives

The settling time — why high Q means slow

How long must you wait for the steady state? Equation (6) answers with one number: the transient's envelope is e−βte^{-\beta t}, so the only timescale of forgetting is

τ=1β(7)\tau = \frac{1}{\beta} \qquad (7)

where: τ\tau — the settling (memory) time: amplitude within 63% of final after τ\tau, 95% after 3τ3\tau, 99% after ln⁡(100) τ≈4.6τ\ln(100)\,\tau \approx 4.6\tau. Note what τ\tau does not depend on: the drive. Loud or soft, on-resonance or off, the approach rate to the steady state is the system's own β\beta — the drive sets the destination, damping sets the travel time. Counting in cycles instead of seconds turns this into the practical rule of resonant instrumentation. Using Q=ω0/2βQ = \omega_0/2\beta (taken as a given here; the definitions live in the Quality Factor article), the number of oscillation periods in one settling time is

Nsettle=τT=ω02πβ=Qπ(8)N_{\text{settle}} = \frac{\tau}{T} = \frac{\omega_0}{2\pi\beta} = \frac{Q}{\pi} \qquad (8)

where: NsettleN_{\text{settle}} — cycles to reach 63% of the final amplitude; roughly 1.5 Q1.5\,Q cycles gets you to 99%. This is the flip side of everything that makes high QQ desirable: the same weak damping that gives a tall, narrow response makes the system slow to change its mind. A resonator that rings for QQ cycles must also be excited for QQ cycles.

Go deeper · Derivation C · 2 minThe resonant build-up envelope and the Q-cycle settling rule
Derivation CThe resonant build-up envelope and the Q-cycle settling rule

We are deriving the clean on-resonance limit of the rest-start solution — the envelope x≈A(1−e−βt)sin⁡ω0tx \approx A(1 - e^{-\beta t})\sin\omega_0 t behind the hero figure — and from it the settling rule of equations (7) and (8). The strategy is to specialize equation (6) to ω=ω0\omega = \omega_0, where the phase lag is exactly π/2\pi/2, and then use light damping (β≪ω0\beta \ll \omega_0) to identify ω1\omega_1 with ω0\omega_0. The step to watch is the ratio ω0/ω1\omega_0/\omega_1 multiplying the transient: it is what keeps the cancellation at t=0t = 0 exact before any approximation is made, and only after the envelope is extracted do we set it to one.

Step 1. Specialize to resonance. At ω=ω0\omega = \omega_0 the steady state has δ=π/2\delta = \pi/2 (from equation (4): the mismatch ω02−ω2\omega_0^2 - \omega^2 vanishes), so cos⁡δ=0\cos\delta = 0, sin⁡δ=1\sin\delta = 1, and the rest-start constants collapse to C1=0C_1 = 0, C2=−A ω0/ω1C_2 = -A\,\omega_0/\omega_1. Why C1C_1 dies: with a 90°90° lag the steady state passes through zero displacement at t=0t = 0 — the position needs no cancelling, only the velocity does.

Step 2. Write the motion with these constants:

x(t)=Asin⁡ω0t−A ω0ω1 e−βtsin⁡ω1tx(t) = A\sin\omega_0 t - A\,\frac{\omega_0}{\omega_1}\,e^{-\beta t}\sin\omega_1 t

Why both terms are sines: on resonance the response lags the cosine drive by a quarter cycle, and the transient must be a sine too, to cancel the response's initial velocity with its own.

Step 3. Apply light damping, β≪ω0\beta \ll \omega_0: then ω1=ω01−β2/ω02≈ω0\omega_1 = \omega_0\sqrt{1 - \beta^2/\omega_0^2} \approx \omega_0 (the difference is second order in β/ω0\beta/\omega_0), and the two sines merge:

x(t)≈A(1−e−βt)sin⁡ω0t(9)x(t) \approx A\big(1 - e^{-\beta t}\big)\sin\omega_0 t \qquad (9)

where: the oscillation is a pure tone at ω0\omega_0 from the first cycle; only its envelope evolves, rising as 1−e−βt1 - e^{-\beta t}. This is the trace on the oscilloscope in every driven-resonance lab: no beats on resonance (the two frequencies coincide), just a smooth climb.

Step 4. Extract the settling rule. The envelope reaches 1−e−1≈63%1 - e^{-1} \approx 63\% at t=1/βt = 1/\beta, which spans t/T=ω0/2πβ=Q/πt/T = \omega_0/2\pi\beta = Q/\pi oscillation periods — equation (8) — and 99%99\% at 4.6/β≈1.5 Q4.6/\beta \approx 1.5\,Q periods. Why the rule is worth memorizing: it converts a quoted QQ directly into waiting time. It is also the energy story in disguise: near the start the drive pumps at its full rate while friction, which bills as amplitude squared, barely collects — the gap closes only as the amplitude nears the balance point. *Verified: equation (9) satisfies x(0)=0x(0) = 0, x˙(0)=0\dot x(0) = 0 to first order in β/ω0\beta/\omega_0; differentiating and substituting into (1) at ω=ω0\omega = \omega_0 leaves only terms of order β2/ω02\beta^2/\omega_0^2; the t→∞t \to \infty limit reproduces A(ω0)A(\omega_0) and the β→0\beta \to 0 limit recovers linear secular growth x→(F0/2mω0) tsin⁡ω0tx \to (F_0/2m\omega_0)\,t\sin\omega_0 t, matching the undamped-catastrophe result quoted in the Resonance article.*

Figure
Settling is measured in cycles — and the price of a cycle is set by Q
final amplitude A 63% of A Q = 25: 63% in 8 cycles Q = 200: 63% in 64 cycles build-up envelope vs cycle number n: 1 − e^(−πn/Q) cycle number n amplitude / A 0

Limiting cases

CaseResultMeaning
long times, t≫1/βt \gg 1/\betax→Acos⁡(ωt−δ)x \to A\cos(\omega t - \delta)the transient is gone; motion is the pure steady state, independent of launch
high QQ (weak damping)τ=1/β\tau = 1/\beta long; settle ≈Q/π\approx Q/\pi cyclesslow forgetting: long beat-like ringing off resonance, slow build-up on resonance
low QQ (strong damping)transient dies within ∼\sim1 cyclesteady state established almost immediately; phasor analysis valid at once
undamped, β=0\beta = 0, ω≠ω0\omega \ne \omega_0transient never decaysfree and forced parts coexist forever as a permanent two-tone superposition — no true "steady state"
undamped on resonance, β=0\beta = 0, ω=ω0\omega = \omega_0x∝tsin⁡ω0tx \propto t\sin\omega_0 tsecular growth, no settling at all — the limit the Resonance article dissects

Worked examples

1 · Matching a launch, number by number (confidence builder)

A cart of mass m=0.50 kgm = 0.50\ \mathrm{kg} on a spring k=200 N/mk = 200\ \mathrm{N/m} with damping β=2.0 s−1\beta = 2.0\ \mathrm{s^{-1}} is driven by F(t)=(2.0 N)cos⁡(10 t)F(t) = (2.0\ \mathrm{N})\cos(10\,t). At t=0t = 0 the cart is held at x0=20 mmx_0 = 20\ \mathrm{mm} and released from rest. Find the complete motion.

Step 1. The fixed part first. ω0=200/0.50=20.0 rad/s\omega_0 = \sqrt{200/0.50} = 20.0\ \mathrm{rad/s}, ω1=ω02−β2=400−4=19.90 rad/s\omega_1 = \sqrt{\omega_0^2 - \beta^2} = \sqrt{400 - 4} = 19.90\ \mathrm{rad/s}. From equation (4) at ω=10 rad/s\omega = 10\ \mathrm{rad/s}: mismatch =400−100=300= 400 - 100 = 300, friction term =2βω=40= 2\beta\omega = 40, so A=4.0/3002+402=13.2 mmA = 4.0/\sqrt{300^2 + 40^2} = 13.2\ \mathrm{mm} and δ=arctan⁡(40/300)=0.133 rad\delta = \arctan(40/300) = 0.133\ \mathrm{rad} (7.6°7.6°). These numbers are now frozen — the launch cannot touch them.

Step 2. First matching equation (5): C1=x0−Acos⁡δ=20.0−13.2cos⁡(0.133)=20.0−13.1=6.9 mmC_1 = x_0 - A\cos\delta = 20.0 - 13.2\cos(0.133) = 20.0 - 13.1 = \boxed{6.9\ \mathrm{mm}}. The steady state already wants to start at 13.1 mm; the transient supplies the missing 6.9.

Step 3. Second matching equation. The steady state's initial velocity is Aωsin⁡δ=13.2×10×0.1325=17.5 mm/sA\omega\sin\delta = 13.2 \times 10 \times 0.1325 = 17.5\ \mathrm{mm/s}, so with v0=0v_0 = 0: C2=(v0+βC1−Aωsin⁡δ)/ω1=(0+2.0×6.9−17.5)/19.90=(13.8−17.5)/19.90=−0.18 mmC_2 = (v_0 + \beta C_1 - A\omega\sin\delta)/\omega_1 = (0 + 2.0 \times 6.9 - 17.5)/19.90 = (13.8 - 17.5)/19.90 = \boxed{-0.18\ \mathrm{mm}}.

Step 4. Assemble and check:

x(t)=13.2cos⁡(10t−0.133)+e−2.0t[6.9cos⁡(19.90 t)−0.18sin⁡(19.90 t)] mmx(t) = 13.2\cos(10t - 0.133) + e^{-2.0t}\big[6.9\cos(19.90\,t) - 0.18\sin(19.90\,t)\big] \ \mathrm{mm}

At t=0t = 0: 13.1+6.9=20.0 mm13.1 + 6.9 = 20.0\ \mathrm{mm} ✓, and the initial velocity terms cancel by construction ✓. The transient starts at 6.92+0.182=6.9 mm\sqrt{6.9^2 + 0.18^2} = 6.9\ \mathrm{mm} and is down to 5% of that by t=3/β=1.5 st = 3/\beta = 1.5\ \mathrm{s} — after which the cart oscillates at 10 rad/s10\ \mathrm{rad/s} with 13.2 mm amplitude, all memory of the 20 mm launch erased.

2 · How long before you can trust a resonant instrument? (exam level)

A quartz tuning-fork resonator runs at f0=32,768 Hzf_0 = 32{,}768\ \mathrm{Hz} with Q=1.0×105Q = 1.0\times10^5. It is driven on resonance from rest. How long, in seconds and in cycles, until its amplitude reaches 63% of the steady state? Until 99%? Compare with a loudspeaker cone (f0=50 Hzf_0 = 50\ \mathrm{Hz}, Q=0.7Q = 0.7).

Step 1. Convert QQ to a damping rate. ω0=2πf0=2.06×105 rad/s\omega_0 = 2\pi f_0 = 2.06\times10^5\ \mathrm{rad/s} and β=ω0/2Q=2.06×105/(2×105)=1.03 s−1\beta = \omega_0/2Q = 2.06\times10^5/(2\times10^5) = 1.03\ \mathrm{s^{-1}}. Remarkable in itself: a resonator vibrating two hundred thousand radians per second forgets at roughly one per second.

Step 2. Settling time from equation (7): τ=1/β=0.97 s\tau = 1/\beta = \boxed{0.97\ \mathrm{s}} to 63%. In cycles, equation (8): N=Q/π=1.0×105/π≈32,000 cyclesN = Q/\pi = 1.0\times10^5/\pi \approx \boxed{32{,}000\ \mathrm{cycles}}.

Step 3. The 99% mark needs ln⁡(100) τ=4.6×0.97=4.5 s\ln(100)\,\tau = 4.6 \times 0.97 = \boxed{4.5\ \mathrm{s}} — about 147,000147{,}000 cycles, or 1.5 Q1.5\,Q. Any measurement scheme that reads this resonator's amplitude must wait seconds per data point; this, not electronics, is the bandwidth bottleneck of high-QQ sensors, and it is why a lock-in measurement of a sharp resonance is inherently slow.

Step 4. The loudspeaker: β=ω0/2Q=2π(50)/1.4=224 s−1\beta = \omega_0/2Q = 2\pi(50)/1.4 = 224\ \mathrm{s^{-1}}, so τ=4.5 ms\tau = 4.5\ \mathrm{ms} and 99% arrives in 20 ms20\ \mathrm{ms} — about one cycle at 50 Hz. The design moral is the contrast: Q=0.7Q = 0.7 is chosen for loudspeakers precisely so the cone settles within a cycle and can track music's transients; Q=105Q = 10^5 is chosen for the fork precisely so it forgets nothing and keeps time. Same equations, opposite corners.

Origins · 1 minLord Kelvin and how this was discovered
LK
Lord Kelvin
1824–1907 · Glasgow

William Thomson — Lord Kelvin from 1892 — met the transient response of a driven system as a matter of intercontinental commerce. Consulting for the Atlantic telegraph in the 1850s, he showed that a pulse fed into a long undersea cable does not arrive as a pulse: the cable's distributed capacitance and resistance smear each symbol into a slow rise toward its steady value, so the operator must wait out a settling time before the next symbol is readable — his "law of squares" predicted that this waiting time grows as the cable length squared. His mirror galvanometer, sensitive enough to read the transient's earliest rise rather than its final level, is what made the 1866 cable commercially workable. Kelvin lived in the world this article describes: steady states you must wait for, transients that carry the message, and instruments designed around the settling time. His harmonic analyser — built with his brother James — even mechanized Fourier superposition, summing sinusoids by machine to predict tides decades before electronics.

History
1855
— William Thomson derives the telegraph equation for the Atlantic cable and the "law of squares": the settling time of a cable's response grows as the square of its length.
1858
— The first transatlantic cable carries faint, transient-smeared pulses; Thomson's mirror galvanometer reads the rising edge of each symbol before the steady level arrives.
1866
— The Great Eastern lays the first permanently successful Atlantic cable with Thomson aboard; he is knighted for making the transient readable.
1877
— Rayleigh's Theory of Sound formalizes the general solution of forced vibration as decaying free motion plus persistent forced motion.
1893
— Charles Steinmetz's phasor method lets engineers compute only the steady state — legitimate precisely because transients die exponentially.

Going further

References. Morin, Introduction to Classical Mechanics, Ch. 4 (Oscillations) · Goldstein, Poole & Safko, Classical Mechanics, 3rd ed., §3.x and small-oscillation chapter · Landau & Lifshitz, Mechanics (Vol. 1), §22, §26 (Forced oscillations) · Taylor, Classical Mechanics, Ch. 5 (Oscillations)

Provenance. Verify the amplitude A(ω)=(F0/m)/(ω02−ω2)2+γ2ω2A(\omega)=(F_0/m)/\sqrt{(\omega_0^2-\omega^2)^2+\gamma^2\omega^2}, phase tan⁡δ=γω/(ω02−ω2)\tan\delta=\gamma\omega/(\omega_0^2-\omega^2), damped frequency ω1=ω02−γ2/4\omega_1=\sqrt{\omega_0^2-\gamma^2/4}, amplitude-resonance frequency ωres=ω02−γ2/2\omega_{\text{res}}=\sqrt{\omega_0^2-\gamma^2/2}, and the power-balance result against Morin Ch. 4, Taylor Ch. 5, Goldstein small-oscillations chapter, and Landau–Lifshitz §26. The rest-start solution can be checked by direct substitution into x(0)=x˙(0)=0x(0)=\dot x(0)=0. Quality factor Q=ω0/γQ=\omega_0/\gamma is the standard convention (some texts use Q=ω0/ΔωQ=\omega_0/\Delta\omega which agrees at high QQ). This page writes β=γ/2\beta = \gamma/2 throughout, matching its sibling articles; all formulas map under γ=2β\gamma = 2\beta. Worked-example numbers (A=13.2A = 13.2 mm, δ=0.133\delta = 0.133 rad, C1=6.9C_1 = 6.9 mm, C2=−0.18C_2 = -0.18 mm; fork β=1.03 s−1\beta = 1.03\ \mathrm{s^{-1}}, τ=0.97\tau = 0.97 s, Q/π≈32,000Q/\pi \approx 32{,}000 cycles, 99% at 4.5 s) are illustrative arithmetic recomputed from equations (4)–(8).

Builds on: Driven Oscillations · Damped Harmonic Oscillator · Underdamped, Critically Damped, and Overdamped Regimes · Complex Exponential Representation of Oscillations · Simple Harmonic Motion Leads to: Resonance · Amplitude, Phase, and Angular Frequency · Energy in Simple Harmonic Motion

The transient amplitude decays as e^(−βt), but I have also seen e^(−2βt). Which is right?

Both — they track different quantities. The amplitude (and any displacement or velocity coordinate) decays as e−βte^{-\beta t}: that is the factor multiplying the transient in equation (3). The energy is quadratic in amplitude, so it decays as (e−βt)2=e−2βt(e^{-\beta t})^2 = e^{-2\beta t} — twice the rate. Confusing the two silently costs factors of 2 in settling estimates: in a ring-down the amplitude takes 1/β1/\beta to fall to 37%, while the stored energy gets there in half that time — and a power meter and a displacement sensor will report different-looking curves for the same event. When quoting a settling time, always say which quantity — and if you use the γ=2β\gamma = 2\beta convention, keep track of which symbol your textbook decays with.

Is the transient a real motion, or just a mathematical correction term?

Completely real — for the first few damping times it can dominate what you see. In the rest-start solution (6) the transient begins exactly as large as the steady state; the wild early oscilloscope trace, the surging envelope, the first violent swing of a pendulum you start pushing — all of that is mostly transient. Engineers fight real transients daily: inrush currents when circuits switch on are transients, the lurch of a car hitting a pothole is a transient, the "thump" a loudspeaker makes at power-on is a transient. What makes the transient special is not unreality but mortality: it is the only part of the motion with an expiry date, and the only part that knows how the motion started.

Why does the transient ring at ω₁ rather than at the natural frequency ω₀ or the drive frequency ω?

Because the transient is a solution of the free damped equation — the drive does not appear in it at all, so ω\omega is irrelevant to its rhythm; and damping slows a free oscillator's swing slightly below ω0\omega_0, to ω1=ω02−β2\omega_1 = \sqrt{\omega_0^2 - \beta^2}. Physically, friction steals a little energy each swing, so the oscillator takes marginally longer per cycle than it would undamped. For light damping the correction is tiny — β=0.05 ω0\beta = 0.05\,\omega_0 shifts the frequency by barely 0.1% — which is why one loosely says the transient "rings at the natural frequency." In precision work the distinction matters: the beat frequency in a startup measurement is ∣ω−ω1∣|\omega - \omega_1|, not ∣ω−ω0∣|\omega - \omega_0|, and fitting ring-downs with ω0\omega_0 instead of ω1\omega_1 biases damping estimates.

The startup wiggles look exactly like beats. Are they the same phenomenon as two-tuning-fork beats?

The interference arithmetic is identical; the cast is different. Classic beats — the throbbing of two tuning forks slightly out of tune, covered in the Beats article — involve two persistent tones, so the envelope oscillates at the difference frequency forever. In a driven startup the two tones are the steady state at ω\omega and the transient at ω1\omega_1, and the second is dying as e−βte^{-\beta t}: you get the same difference-frequency envelope, but each surge is weaker than the last, and after a few 1/β1/\beta the pattern is gone. A useful diagnostic reading: if you observe decaying beats after switching on a drive, you are watching the system forget; the beat rate tells you ∣ω−ω1∣|\omega - \omega_1| and the decay of the surges tells you β\beta — two system parameters from one scope trace.

Why exactly do high-Q systems take so long to settle? Is it just a formula, or is there a physical reason?

There is a clean energy reason. To sit at the steady-state amplitude, a resonator must hold a certain stored energy, and the drive can deliver energy only so fast — at a rate limited by the force and the (growing) velocity. A high-QQ system on resonance settles at a large amplitude (QQ times the static response, as the Resonance article shows), so there is a great deal of energy to pour in, while the drive's per-cycle deposit stays modest. Equivalently, in the frequency domain: building a response that distinguishes ω0\omega_0 from frequencies Δω∼ω0/Q\Delta\omega \sim \omega_0/Q away requires observing the drive for at least ∼Q\sim Q cycles — you cannot resolve a narrow line quickly. Sharpness in frequency and sluggishness in time are the same property, seen from two sides.

Does it matter at what phase of the force cycle I switch the drive on?

For the destination, no; for the journey, yes. Switching on at the force's peak versus its zero is equivalent to changing the drive's phase constant, which shifts the steady state's phase accordingly — but its amplitude AA and lag δ\delta relative to the force are untouched. What changes is the mismatch between the actual initial conditions and what the steady state "wants" at the switch-on instant, i.e. the constants C1,C2C_1, C_2: a different switch-on phase means a different-sized, differently-phased transient and a visibly different startup wiggle. The settling rate is always β\beta regardless. Pulse-sequence designers exploit this: choosing the switch-on phase can minimize (or deliberately maximize) the transient — a trick used to ring up resonators faster than 1−e−βt1 - e^{-\beta t} by briefly overdriving.

What does the transient look like if the system is critically damped or overdamped?

The decomposition survives untouched; only the transient's costume changes. For overdamping (β>ω0\beta > \omega_0) the homogeneous solution is a sum of two decaying exponentials, C1e−(β−β2−ω02)t+C2e−(β+β2−ω02)tC_1 e^{-(\beta - \sqrt{\beta^2 - \omega_0^2})t} + C_2 e^{-(\beta + \sqrt{\beta^2-\omega_0^2})t} — no ringing, just a slouch toward the steady state; at critical damping it is (C1+C2t)e−βt(C_1 + C_2 t)e^{-\beta t}. Two constants remain in every case, so initial conditions are matched the same way. One subtlety: overdamping does not settle faster — the slower of the two exponentials has a decay rate below β\beta, so heavy damping actually prolongs the memory. The fastest possible forgetting is at critical damping, which is exactly why measuring instruments and vehicle suspensions are tuned near it.

During the transient, is the energy balance of the steady state violated?

Not violated — just not yet in balance, and that is the whole point. In steady state the drive's average power input equals the dissipated power exactly, which is why the amplitude holds constant. During startup from rest the amplitude is below the balance point: friction (billing as amplitude squared) collects less than the drive deposits, and the surplus accumulates as growing stored energy. As the amplitude nears the steady value the books close and growth stops. Run it the other way — switch the drive off — and the input drops to zero while dissipation continues: the stored energy drains and you get the mirror-image transient, a free ring-down at ω1\omega_1. The detailed power bookkeeping of the settled state lives in the Power Absorption and Phase Lag article.

Before this

Where this leads