Keyboard shortcuts

Press ← or → to navigate between chapters

Press S or / to search in the book

Press ? to show this help

Press Esc to hide this help

Forcing and resonance

Chapter 8 - Separate natural motion from response

A nonhomogeneous equation adds an input:

ay′′+by′+cy=g(t). ay''+by'+cy=g(t).

Every solution splits into

y=yh+yp, \boxed{y=y_h+y_p},

where yhy_h is homogeneous natural motion and ypy_p is one particular response to the forcing. Natural motion is a decaying transient only when all homogeneous modes decay; undamped oscillations persist, and unstable modes can grow.

Match the input's shape

For constant coefficients, undetermined coefficients proposes a particular solution with the same family of shapes as g(t)g(t).

Polynomialg(t)=t²+1

Try a general polynomial of the same degree.

Exponentialg(t)=e^{at}

Try AeatAe^{at}.

Sinusoidg(t)=cos(ωt)

Try Acos⁡(ωt)+Bsin⁡(ωt)A\cos(ωt)+B\sin(ωt).

Productsg(t)=te^{at}

Multiply the corresponding trial families.

If the trial already belongs to yhy_h, multiply it by enough powers of tt to make it independent. This is the algebraic signature of resonance.

Worked example: ordinary forcing

Solve y′′+y=3cos⁡2ty''+y=3\cos2tshape matching
  1. Natural motion. r2+1=0r^2+1=0 gives yh=C1cos⁡t+C2sin⁡ty_h=C_1\cos t+C_2\sin t.
  2. Choose a trial. Since frequency 2 does not appear in yhy_h, try yp=Acos⁡2t+Bsin⁡2ty_p=A\cos2t+B\sin2t.
  3. Substitute. yp′′+yp=(−4A+A)cos⁡2t+(−4B+B)sin⁡2ty_p''+y_p=(-4A+A)\cos2t+(-4B+B)\sin2t.
  4. Match coefficients. −3A=3-3A=3 and −3B=0-3B=0, so A=−1A=-1, B=0B=0.
  5. Combine. y=C1cos⁡t+C2sin⁡t−cos⁡2ty=C_1\cos t+C_2\sin t-\cos2t.
  6. Interpret. The output contains natural frequency 1 and forced frequency 2. Initial conditions would set the natural-motion amplitudes.

Lab: sweep through resonance

Forced oscillatorx′′+cx′+x=cos⁡(ωt)x''+cx'+x=\cos(\omega t)

time responsesteady-state amplitude

Reduce damping, then move ω\omega near 1. The response curve rises sharply because the input frequency approaches the natural frequency. Damping keeps the steady-state amplitude finite and shifts the exact peak slightly below 1.

Resonant forcing

For

y′′+y=cos⁡t, y''+y=\cos t,

the ordinary trial Acos⁡t+Bsin⁡tA\cos t+B\sin t duplicates yhy_h. Multiply by tt:

yp=t(Acos⁡t+Bsin⁡t). y_p=t(A\cos t+B\sin t).

Substitution gives one valid choice yp=12tsin⁡ty_p=\tfrac12 t\sin t. Its envelope grows linearly because an undamped input adds energy in phase with the natural motion.

Trial collision. For y′′−2y′+y=ety''-2y'+y=e^t, which particular trial is appropriate?

The root r=1r=1 repeats twice in the homogeneous equation, so multiply the exponential trial by t2t^2.

Fade the response

Find a particular solution of y′′+3y′+2y=4e−ty''+3y'+2y=4e^{-t}.

  1. Characteristic roots: ‾\underline{\hspace{5em}}.
  2. Collision? The input e−te^{-t} [is / is not] a natural mode.
  3. Correct trial: yp=‾y_p=\underline{\hspace{6em}}.
  4. Substitute and find the coefficient.
Check the response

The roots are −1,−2-1,-2, so e−te^{-t} collides once. Try yp=Ate−ty_p=Ate^{-t}. Substitution yields Ae−t=4e−tA e^{-t}=4e^{-t}, so A=4A=4. One particular solution is 4te−t4te^{-t}.

Transient versus steady state

In a damped stable system, yhy_h decays. What remains is the steady-state response ypy_p. This split explains why two experiments with different initial conditions can eventually show the same driven behavior.

Error clinic: applying initial conditions to ypy_p alone

The particular solution only supplies one response to the input. Initial conditions select constants after yh+ypy_h+y_p has been assembled.

Practice

  1. Practice - polynomial input

    8.1 Find the general solution of y′′+y=ty''+y=t.

    Detailed solution

  2. Practice - sinusoidal input

    8.2 Solve y′′+4y=cos⁡ty''+4y=\cos t, y(0)=0y(0)=0, y′(0)=0y'(0)=0.

    Detailed solution

  3. Integrated - collision

    8.3 For each input to y′′+4y=g(t)y''+4y= g(t), state the trial without solving: (a) e−te^{-t}, (b) cos⁡2t\cos2t, (c) t2t^2.

    Detailed solution

  4. Challenge - beating

    8.4 Use cos⁡at−cos⁡bt=−2sin⁡((a+b)t/2)sin⁡((a−b)t/2)\cos at-\cos bt=-2\sin((a+b)t/2)\sin((a-b)t/2) to explain the slow amplitude modulation when an undamped oscillator is forced near, but not at, its natural frequency.

    Detailed solution

  5. Practice - extension

    8.5 Solve y′′−2y′+y=ety''-2y'+y=e^t, y(0)=0y(0)=0, y′(0)=0y'(0)=0. Explain the required power of tt in the particular trial, and verify the complete solution using the substitution y=etuy=e^t u.

    Detailed solution

  6. Challenge - transfer

    8.6 Solve y′′+2y′+2y=cos⁡ty''+2y'+2y=\cos t, y(0)=0y(0)=0, y′(0)=0y'(0)=0. Find the steady-state amplitude, identify the transient, and explain why two solutions with different initial conditions approach the same periodic response.

    Detailed solution

Retrieve the design logic

Explain without formulas:

  1. why a particular solution should resemble the input;
  2. why a collision with natural motion requires a factor of tt; and
  3. why damping changes a resonance peak.

Next: Read coupled states as vector fields, eigen-directions, and phase portraits.