Login
📚 Differential Equations, Interactive Edition
Chapters ▾
⇩ 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), equal to 0 for t<a and to 1 for t>a: nothing, then suddenly something. (The figure below builds it arithmetically as 12(sign(ta)+1), which is 0 before the switch and 1 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 s,

F ( s ) = { f } ( s ) = 0 e s t f ( t ) d t ,

and its worth rests on three dictionary entries. Differentiation becomes multiplication: {x}=sX(s)x(0) and {x}=s2X(s)sx(0)x(0). The constant 1 becomes 1/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 ) ,

whose right side is the transform of f multiplied by a factor that says "start this at t=a instead of t=0." In particular {ua(t)}=eas/s.

Now the model problem. Take a damped spring from the family of Chapter 2 — damping a=1, frequency b=3, characteristic roots 1±3i — sitting at rest until a constant force switches on:

x + 2 x + 10 x = 10 u a ( t ) , x ( 0 ) = x ( 0 ) = 0 .

Transforming every term (the initial conditions contribute nothing, being zero) gives (s2+2s+10)X(s)=10eas/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 ) .

The three pieces inside the parentheses are standard: 1/s is the constant 1, and the two shifted fractions are the transforms of etcos3t and 13etsin3t. So without its tag the response would be f(t)=1et(cos3t+13sin3t), and the tag eas delays it wholesale:

x ( t ) = u a ( t ) [ 1 e ( t a ) ( cos 3 ( t a ) + 1 3 sin 3 ( t a ) ) ] .

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

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))] of x+2x+10x=10ua(t). The dashed line is the new equilibrium x=1; the slider a sets the switch-on time and carries the entire ringing pattern rigidly with it.

Explore

  1. Observe. With a=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 a from 3 to 6. 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 eas made visible: a delay in time and nothing else.
  4. Observe. Read the first peak at any a: height about 1.35, reached about 1.05 after switch-on. The second peak is barely above the dashed line, near 1.04. The overshoot has collapsed to roughly an eighth of its first value in a single ring.
  5. Predict, then verify. Push a to 8. Which features of the response still fit in the window t12, and which have fallen off the right edge? What would a switch at t=12 show?

Connect

The response has two personalities, and the partial fractions found them mechanically. The term 1/s inverts to the steady state x=1, where the new constant force 10 exactly balances the spring force 10x. 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). No new behavior was invented at the switch: the denominator (s+1)2+9 is the characteristic polynomial s2+2s+10 wearing its roots 1±3i on its sleeve, so the ringing frequency and decay rate after the jump are the ones the free system always had.

The shifting rule also handles switching off. A force applied only for atb is ua(t)ub(t), and by linearity the response is f(ta)ua(t)f(tb)ub(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.

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.