Laplace Transforms — Question 2

PDF ↗

Question 2

A sixth-order initial-value problem is (D2+1)3y=0,y(j)(0)=0(0≤j≤4),y(5)(0)=1,(D^2+1)^3y=0,\qquad y^{(j)}(0)=0\ (0\leq j\leq 4),\qquad y^{(5)}(0)=1, where D=d/dtD=d/dt and t≥0t\geq 0. To invert a repeated quadratic pole, introduce Ha(s)=1/(s2+a2)H_a(s)=1/(s^2+a^2) for a>0a>0.

Tasks

  1. Expand the differential operator and transform the problem. Explain why only one initial term survives.

  2. Starting with ℒ−1Ha=sin⁡(at)/a\mathcal L^{-1}H_a=\sin(at)/a, use differentiation with respect to aa twice to invert (s2+a2)−3(s^2+a^2)^{-3}. Then set a=1a=1.

  3. Verify the equation and initial state without relying on the transform calculation. Identify the first nonzero Taylor term.

  4. Describe the pole multiplicities and a convergence half-plane. Decide whether the solution is bounded and whether the final-value theorem can be used, giving a decisive argument for each claim.

Original worksheet page 1: question and worked solution for 7-5-002
Show solutionHide solution

Question 2 – Solution

Strategy. Parameter differentiation resolves repeated quadratic poles while a Taylor check tests the high-order initial data.

Step 1: Transform the sixth-order problem. The equation is y(6)+3y(4)+3y″+y=0y^{(6)}+3y^{(4)}+3y''+y=0. Only ℒ(y(6))=s6Y−1\mathcal L(y^{(6)})=s^6Y-1 has a nonzero initial subtraction, so Y(s)=1(s2+1)3.\boxed{Y(s)=\frac 1{(s^2+1)^3}.}

Step 2: Differentiate a transform family. Since ∂a(s2+a2)−m=−2ma(s2+a2)−m−1\partial_a(s^2+a^2)^{-m}=-2ma(s^2+a^2)^{-m-1}, h2(t)=sin⁡(at)−atcos⁡(at)2a3,h3(t)=−14a∂ah2(t)=3sin⁡(at)−3atcos⁡(at)−a2t2sin⁡(at)8a5.h_2(t)=\frac{\sin(at)-at\cos(at)}{2a^3},\qquad h_3(t)=-\frac 1{4a}\partial_a h_2(t) =\frac{3\sin(at)-3at\cos(at)-a^2t^2\sin(at)}{8a^5}. Differentiation under the integral is valid for aa in a compact positive interval and Re⁡s>0\operatorname{Re}s>0, since the derivatives have polynomial growth. Hence y(t)=3sin⁡t−3tcos⁡t−t2sin⁡t8.\boxed{y(t)=\frac{3\sin t-3t\cos t-t^2\sin t}{8}.}

Step 3: Verify independently. Each term is the real or imaginary part of eite^{it} times a polynomial of degree at most two; (D−i)3(D-i)^3 annihilates these complex modes, and the conjugate factor annihilates their conjugates. Thus (D2+1)3y=0(D^2+1)^3y=0. Expanding the displayed trigonometric expression gives y(t)=t5120−t71680+O(t9).y(t)=\frac{t^5}{120}-\frac{t^7}{1680}+O(t^9). This proves the five zero initial derivatives and y(5)(0)=1y^{(5)}(0)=1.

Step 4: Distinguish oscillation from boundedness. There are triple poles at s=±is=\pm i; Re⁡s>0\operatorname{Re}s>0 gives absolute convergence. Along tn=π/2+2πnt_n=\pi/2+2\pi n, y(tn)=(3−tn2)/8→−∞y(t_n)=(3-t_n^2)/8\to-\infty, so the response is unbounded. The poles of sYsY on the imaginary axis invalidate the final-value theorem. The formal value lim⁡s→0sY=0\lim_{s\to 0}sY=0 is not a time limit.

Original worksheet page 2: question and worked solution for 7-5-002

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