Solving IVP’s with Laplace Transforms — Question 4

PDF ↗

Question 4

Use one-sided Laplace transforms and retain all initial-value terms. Write Y(s)=ℒ{y}(s)Y(s)=\mathcal L\{y\}(s), with real ss sufficiently large during the transformation. Unless stated otherwise, solve on t≥0t\ge 0.

An undamped oscillator starts at rest: y″+4y=sin⁡2t,y(0)=0,y′(0)=0.y''+4y=\sin 2t,\qquad y(0)=0,\qquad y'(0)=0. Define its energy by E(t)=[y′(t)2+4y(t)2]/2E(t)=[y'(t)^2+4y(t)^2]/2.

Tasks

  1. Derive YY and invert it. You may use ℒ{tcos⁡2t}=(s2−4)/(s2+4)2\mathcal L\{t\cos 2t\}=(s^2-4)/(s^2+4)^2.

  2. Check the equation and both initial data directly.

  3. Prove that the response is unbounded by evaluating a suitable explicit time sequence. State its exact real transform domain.

  4. Derive the energy identity and calculate E(nπ)E(n\pi) for positive integers nn. Explain why bounded forcing can supply unbounded accumulated energy.

Original worksheet page 1: question and worked solution for 4-5-004
Show solutionHide solution

Question 4 – Solution

Strategy. Matching forcing and natural frequencies produce repeated transform poles. Explicit values and the work identity distinguish real growth from a mere upper bound.

Step 1: Transform and invert. Zero initial data give (s2+4)Y=2s2+4,Y=2(s2+4)2.(s^2+4)Y=\frac 2{s^2+4},\qquad Y=\frac 2{(s^2+4)^2}. Combining the sine pair with the supplied time-cosine pair yields y(t)=sin⁡2t−2tcos⁡2t8.\boxed{y(t)=\frac{\sin 2t-2t\cos 2t}{8}.} Indeed the forward numerator is (s2+4)/4−(s2−4)/4=2(s^2+4)/4-(s^2-4)/4=2.

Step 2: Verify the IVP. Direct differentiation gives y′=t2sin⁡2t,y″=12sin⁡2t+tcos⁡2t.y'=\tfrac t2\sin 2t,\qquad y''=\tfrac 12\sin 2t+t\cos 2t. Adding 4y=12sin⁡2t−tcos⁡2t4y=\tfrac 12\sin 2t-t\cos 2t gives sin⁡2t\sin 2t. Both initial values are zero, proving the claimed solution by uniqueness.

Step 3: Exhibit unbounded values and the domain. At t=nπt=n\pi, y(nπ)=−nπ/4→−∞.\boxed{y(n\pi)=-n\pi/4\longrightarrow-\infty.} The actual response is therefore unbounded; a growing envelope alone would not prove that. Its linear envelope gives absolute transform convergence for s>0s>0. At zero, its finite integral is [1−cos⁡2R−Rsin⁡2R]/8[1-\cos 2R-R\sin 2R]/8, which has no limit. Below zero, intervals around successive peaks of the tcos⁡2tt\cos 2t term violate the Cauchy criterion. Hence the exact domain is s>0\boxed{s>0}.

Step 4: Account for the energy input. Multiplying the equation by y′y' gives E′=y′(y″+4y)=y′sin⁡2t=t2sin⁡22t≥0.E'=y'(y''+4y)=y'\sin 2t=\tfrac t2\sin^2 2t\ge 0. At nπn\pi, the velocity is zero, so E(nπ)=2y(nπ)2=n2π2/8.\boxed{E(n\pi)=2y(n\pi)^2=n^2\pi^2/8.} Although the forcing amplitude is at most one, its work accumulates over an increasing time interval. Resonance makes that accumulated work unbounded. The graph shows the growing displacement oscillations, not the energy itself.

See the diagram in the original worksheet below.

Original worksheet page 2: question and worked solution for 4-5-004

Original worksheet layout. Use Enlarge or open the PDF for a closer view.