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(ta)+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)}=eas/s\mathcal{L}\{u_a(t)\} = e^{-as}/s.

Now the model problem. Take a damped spring from the family of Chapter 2 — damping a=1a = 1, frequency b=3b = 3, 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)=10eas/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 etcos3te^{-t}\cos 3t and 13etsin3t\tfrac{1}{3} e^{-t}\sin 3t. So without its tag the response would be f(t)=1et(cos3t+13sin3t)f(t) = 1 - e^{-t}\left(\cos 3t + \tfrac{1}{3}\sin 3t\right), and the tag ease^{-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.

(2s)·(10+2·s+s2)1+s1\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, 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 etcos3te^{-t}\cos 3t, and what is left over is 13etsin3t\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 that lies flat at zero, then at the switch-on time springs upward, overshoots a dashed horizontal line at height 1, rings through a few shrinking oscillations, and settles onto the dashed line. The slider a slides the switch-on moment anywhere from t = 1 to t = 8, carrying the entire ringing pattern rigidly with it.Explore in 3D (opens in a new tab)
The rest-to-rest step response x(t)=ua(t)[1e(ta)(cos3(ta)+13sin3(ta))]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). The dashed 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 dashed 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 ease^{-as} made visible: a delay in time and nothing else.
  4. Observe. Read the first peak at any aa: height about 1.351.35, reached about 1.051.05 after switch-on. The second peak is barely above the dashed line, near 1.041.04. 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 t12t \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(ta)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+3·I,13·I-1 + 3 \cdot I,\; -1 - 3 \cdot 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 atba \le t \le b is ua(t)ub(t)u_a(t) - u_b(t), and by linearity the response is f(ta)ua(t)f(tb)ub(t)f(t-a)\,u_a(t) - f(t-b)\,u_b(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 τ=ta\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 ADA Title II. Every page is built against an automated accessibility gate, and the published editions will meet ADA Title II requirements when they release in late September 2026. If something is unusable, please tell us.