Solving IVP’s with Laplace Transforms — Question 7

PDF ↗

Question 7

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.

The data are prescribed at time one, not at zero: y″+y=t,y(1)=2,y′(1)=0,t≥1.y''+y=t,\qquad y(1)=2,\qquad y'(1)=0,\qquad t\ge 1. Use a new time origin before applying a one-sided transform.

Tasks

  1. Define τ=t−1\tau=t-1 and z(τ)=y(1+τ)z(\tau)=y(1+\tau). Derive the transformed IVP for zz.

  2. Find Z(s)=ℒτ{z}(s)Z(s)=\mathcal L_\tau\{z\}(s) and invert it.

  3. Convert back to tt and verify the equation and the data at time one.

  4. A student treats the supplied values as y(0),y′(0)y(0),y'(0) and transforms the forcing as 1/s21/s^2. Explain both changes needed to repair that calculation. State the exact domain of ZZ.

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

Question 7 – Solution

Strategy. Shift the independent variable and the forcing together. A Laplace transform beginning at the new origin must use data at that origin.

Step 1: Move the initial time. With τ=t−1\tau=t-1, the derivatives have the same form but the forcing becomes τ+1\tau+1. Thus z″+z=τ+1,z(0)=2,z′(0)=0.z''+z=\tau+1,\qquad z(0)=2,\qquad z'(0)=0. Writing Z=ℒτ{z}Z=\mathcal L_\tau\{z\} gives (s2+1)Z−2s=1s2+1s.(s^2+1)Z-2s=\frac 1{s^2}+\frac 1s.

Step 2: Resolve and invert. A useful decomposition is Z=2s+1/s2+1/ss2+1=1s2+1s+ss2+1−1s2+1.Z=\frac{2s+1/s^2+1/s}{s^2+1} =\frac 1{s^2}+\frac 1s+\frac{s}{s^2+1}-\frac 1{s^2+1}. Therefore z(τ)=τ+1+cos⁡τ−sin⁡τ\boxed{z(\tau)=\tau+1+\cos\tau-\sin\tau}. Transforming each term back checks the decomposition, including the constant part of the shifted forcing.

Step 3: Return to the original clock. The solution on the requested domain is y(t)=t+cos⁡(t−1)−sin⁡(t−1),t≥1.\boxed{y(t)=t+\cos(t-1)-\sin(t-1),\qquad t\ge 1.} Its first derivative is 1−sin⁡(t−1)−cos⁡(t−1)1-\sin(t-1)-\cos(t-1) and its second is −cos⁡(t−1)+sin⁡(t−1)-\cos(t-1)+\sin(t-1). Hence y″+y=ty''+y=t, while y(1)=2y(1)=2 and y′(1)=0y'(1)=0. These are the actual specified data, so uniqueness on [1,∞)[1,\infty) verifies the result.

Step 4: Diagnose the two clock errors. The values two and zero belong to z(0),z′(0)z(0),z'(0), not to values of yy at an unspecified earlier time. Once that clock change is made, the forcing is τ+1\tau+1 and its transform is 1/s2+1/s1/s^2+1/s, not just 1/s21/s^2. The transform being inverted is ZZ, defined using τ≥0\tau\ge 0. Since z(τ)=τ+O(1)z(\tau)=\tau+O(1) is eventually positive, it has exact domain s>0\boxed{s>0}. No extension of yy to times before one was needed.

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

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