Differential Equations, Interactive EditionXYZ Homework Edition

⇩ Download ▾

3.1 Step Forcing

A furnace, a relay, a charging circuit — real inputs are switched. The cleanest model of a switch is the Heaviside step function ua(t)u_a(t), equal to 00 for t<at < a and to 11 for t>at > a: nothing, then suddenly something. (The figure below builds it arithmetically as 12(sign(t−a)+1)\tfrac{1}{2}(\mathrm{sign}(t - a) + 1), which is 00 before the switch and 11 after it.) An input with a jump defeats any method that wants one smooth formula for all time, so this section adopts a tool that ingests the whole history at once.

The Laplace transform trades a function of time for a function of a new variable ss,

F ( s ) = ℒ { f } ( s ) = ∫ 0 ∞ e − s t f ( t ) d t , F(s) = \mathcal{L}\{f\}(s) = \int_0^\infty e^{-st} f(t)\, dt,

and its worth rests on three dictionary entries. Differentiation becomes multiplication: ℒ{x′}=sX(s)−x(0)\mathcal{L}\{x'\} = sX(s) - x(0) and ℒ{x″}=s2X(s)−sx(0)−x′(0)\mathcal{L}\{x''\} = s^2 X(s) - s\, x(0) - x'(0). The constant 11 becomes 1/s1/s. And a delayed copy of a signal becomes an exponential tag — the shifting rule

ℒ { u a ( t ) f ( t − a ) } = e − a s F ( s ) , \mathcal{L}\{u_a(t)\, f(t - a)\} = e^{-as} F(s),

whose right side is the transform of ff multiplied by a factor that says "start this at t=at = a instead of t=0t = 0." In particular ℒ{ua(t)}=e−as/s\mathcal{L}\{u_a(t)\} = e^{-as}/s.

Now the model problem. Take a damped spring from the family of Chapter 2 — decay rate 11, ringing frequency 33, characteristic roots −1±3i-1 \pm 3i — sitting at rest until a constant force switches on:

x ″ + 2 x ′ + 10 x = 10 u a ( t ) , x ( 0 ) = x ′ ( 0 ) = 0 . x'' + 2x' + 10x = 10\, u_a(t), \qquad x(0) = x'(0) = 0.

Transforming every term (the initial conditions contribute nothing, being zero) gives (s2+2s+10)X(s)=10e−as/s(s^2 + 2s + 10)\, X(s) = 10\, e^{-as}/s, so

X ( s ) = e − a s · 10 s ( s 2 + 2 s + 10 ) = e − a s ( 1 s − s + 1 ( s + 1 ) 2 + 9 − 1 3 · 3 ( s + 1 ) 2 + 9 ) . X(s) = e^{-as} \cdot \frac{10}{s\,(s^2 + 2s + 10)} = e^{-as} \left( \frac{1}{s} - \frac{s + 1}{(s+1)^2 + 9} - \frac{1}{3} \cdot \frac{3}{(s+1)^2 + 9} \right).

The three pieces inside the parentheses are standard: 1/s1/s is the constant 11, and the two shifted fractions are the transforms of e−tcos3te^{-t}\cos 3t and 13e−tsin3t\tfrac{1}{3} e^{-t}\sin 3t. So without its tag the response would be f(t)=1−e−t(cos3t+13sin3t)f(t) = 1 - e^{-t}\left(\cos 3t + \tfrac{1}{3}\sin 3t\right), and the tag e−ase^{-as} delays it wholesale:

x ( t ) = u a ( t ) [ 1 − e − ( t − a ) ( cos 3 ( t − a ) + 1 3 sin 3 ( t − a ) ) ] . x(t) = u_a(t)\left[\, 1 - e^{-(t-a)} \left( \cos 3(t-a) + \tfrac{1}{3} \sin 3(t-a) \right) \right].

Nothing happens until t=at = a; then the mass climbs toward the new equilibrium x=1x = 1, overshoots it, and rings its way down onto it at the free-vibration frequency 33.

(−2−s)·(10+2·s+s2)−1+s−1\left(-2 - s\right) \cdot \left(10 + 2 \cdot s + s^{2}\right)^{-1} + s^{-1}
The decomposition, done by machine ✓ Computed · mojocas 0.1.0 ✓ Agrees with the text The decomposition, done by machine, computed exactly by mojocas 0.1.0, and confirmed to agree with the result stated in the text.

The one algebraic step between X(s)X(s) and the answer, and the step the section names in its own objective. The engine decomposes 10/(s(s2+2s+10))10/\bigl(s(s^2+2s+10)\bigr) independently and confirms the three pieces above are that fraction.

Read the two forms side by side, because they do not look alike. The engine returns the canonical decomposition, displayed with negative powers as (−2−s)(10+2s+s2)−1+s−1(-2 - s)(10 + 2s + s^2)^{-1} + s^{-1} — that is, 1/s−(s+2)/(s2+2s+10)1/s - (s+2)/(s^2+2s+10): two terms, no completed square. The text splits that second term in two and shifts both halves, giving three terms once the untouched 1/s1/s is counted back in, because that form is the one you can read a transform off: (s+1)(s+1) over a shifted denominator is e−tcos3te^{-t}\cos 3t, and what is left over is 13e−tsin3t\tfrac{1}{3} e^{-t}\sin 3t. Algebra prefers the first form; a table lookup needs the second. The certificate says they are the same fraction, which is what makes the rewrite legitimate rather than merely traditional.

The decomposition is also self-verifying at the engine boundary: it reassembles the pieces and requires the difference from the input to normalize to zero before returning anything at all.

A red curve on a time axis, over a grid of unit squares, that lies flat at zero, then at the switch-on time springs upward, overshoots an amber horizontal line at height 1 by about a third of a unit, rings through a few shrinking oscillations, and settles onto the amber line. The slider a slides the switch-on moment anywhere from time 1 to time 8, carrying the entire ringing pattern rigidly with it.Explore in 3D (opens in a new tab)

A red curve on a time axis, over a grid of unit squares, that lies flat at zero, then at the switch-on time springs upward, overshoots an amber horizontal line at height 1 by about a third of a unit, rings through a few shrinking oscillations, and settles onto the amber line. The slider a slides the switch-on moment anywhere from time 1 to time 8, carrying the entire ringing pattern rigidly with it.

The rest-to-rest step response x(t)=ua(t)[1−e−(t−a)(cos3(t−a)+13sin3(t−a))]x(t) = u_a(t)\left[1 - e^{-(t-a)}(\cos 3(t-a) + \tfrac{1}{3}\sin 3(t-a))\right] of x″+2x′+10x=10ua(t)x'' + 2x' + 10x = 10\,u_a(t). Time tt runs along the horizontal axis from 00 to 1212, and xx is measured upward (grid squares one unit). The amber line is the new equilibrium x=1x = 1; the slider aa sets the switch-on time and carries the entire ringing pattern rigidly with it.

Explore

  1. Observe. With a=3a = 3, the curve before the switch is not merely small — it is exactly zero. Explain from the equation why rest plus zero forcing means zero, identically. Then count how many times the curve visibly crosses the amber line before settling onto it.
  2. Predict. You are about to drag aa from 33 to 66. Commit first: will the response change its shape, its height, or only its starting point?
  3. Verify. Drag. The whole pattern translates rigidly — the climb, the first peak, and each ring sit at the same distances after the switch, wherever the switch is. That rigidity is the factor e−ase^{-as} made visible: a delay in time and nothing else.
  4. Observe. Estimate the first peak at any aa against the amber line, which stands one unit above the axis: the peak clears it by about a third of that — height about 1.351.35, reached about 1.051.05 after switch-on. The second peak is barely above the amber line, near 1.041.04 — too close to separate by eye, which is itself the point. The overshoot has collapsed to roughly an eighth of its first value in a single ring.
  5. Predict, then verify. Push aa to 88. Which features of the response still fit in the window t≤12t \le 12, and which have fallen off the right edge? What would a switch at t=12t = 12 show?

Connect

The response has two personalities, and the partial fractions found them mechanically. The term 1/s1/s inverts to the steady state x=1x = 1, where the new constant force 1010 exactly balances the spring force 10x10x. The remaining terms invert to the transient — a free damped vibration straight out of Chapter 2, born at the switching instant and dying like e−(t−a)e^{-(t-a)}. No new behavior was invented at the switch: the denominator (s+1)2+9(s+1)^2 + 9 is the characteristic polynomial s2+2s+10s^2 + 2s + 10 wearing its roots −1±3i-1 \pm 3i on its sleeve, so the ringing frequency and decay rate after the jump are the ones the free system always had.

−1+3i,−1−3i-1 + 3 i,\; -1 - 3 i
The roots behind the ringing ✓ Computed · mojocas 0.1.0 ✓ Agrees with the text The roots behind the ringing, computed exactly by mojocas 0.1.0, and confirmed to agree with the result stated in the text.

The characteristic polynomial the denominator was wearing, solved exactly: one conjugate pair, decay rate −1-1, ringing frequency 33 — certified against the roots the prose reads off its sleeve.

The shifting rule also handles switching off. A force applied only for a≤t≤ca \le t \le c is ua(t)−uc(t)u_a(t) - u_c(t), and by linearity the response is f(t−a)ua(t)−f(t−c)uc(t)f(t-a)\,u_a(t) - f(t-c)\,u_c(t): turn on the response, then turn on its negative. A whole chain of switches — on, off, on, off — costs one tagged term per switch, bookkeeping that the next section compresses into a geometric series.

−(−3·sin(3·t)+cos(3·t))·exp(−t)+(13·sin(3·t)+cos(3·t))·exp(−t)-\left(-3 \cdot \sin\left(3 \cdot t\right) + \cos\left(3 \cdot t\right)\right) \cdot \exp\left(-t\right) + \left(\frac{1}{3} \cdot \sin\left(3 \cdot t\right) + \cos\left(3 \cdot t\right)\right) \cdot \exp\left(-t\right)
The slope of the step response ✓ Computed · mojocas 0.1.0 ✓ Agrees with the text The slope of the step response, computed exactly by mojocas 0.1.0, and confirmed to agree with the result stated in the text.

The overshoot argument's one computation, certified. The engine shows its derivative in raw, uncollapsed form — and the check proves that expression equal to the 103e−τsin3τ\tfrac{10}{3} e^{-\tau} \sin 3\tau the example states, which is exactly the claim that the cosine terms cancel. (Written in tt here; the example's τ=t−a\tau = t - a shift changes nothing, as the example itself notes.)

An original work of XYZ Homework, built around interactive XYZ 3D figures, following the standard first-course sequence (first-order equations through systems and the phase plane). Distinct from the LibreTexts-sourced Differential Equations for Engineers (Lebl) edition, which has its own attribution. License: CC-BY-NC-SA-4.0.

These eBooks are a prerelease and are not yet certified conformant with WCAG 2.1 AA or any other accessibility standard. Every page is built against an automated accessibility gate, and our target for the published editions is WCAG 2.2 Level AA. If something is unusable, please tell us.