Calculus · real student question

Use the Laplace transform to solve the initial value problem y'' - y = 4*delta(t - 2) + t^2 with y(0) = 0 and y'(0) = 2, where delta is the Dirac delta function.

Question

Solve the initial value problem

yy=4δ(t2)+t2,y(0)=0,y(0)=2y'' - y = 4\delta(t-2) + t^2, \qquad y(0)=0,\quad y'(0)=2

where δ\delta denotes the Dirac delta function.

Step-by-step solution

  1. Choose the Laplace transform because of the impulse. The forcing term contains δ(t2)\delta(t-2), which is not a function that can be integrated by ordinary methods; undetermined coefficients and variation of parameters both fail. The Laplace transform turns the impulse into the ordinary exponential e2se^{-2s}, so the whole problem becomes algebra.

  2. Transform both sides, folding in the initial conditions. Using L{y}=s2Y(s)sy(0)y(0)\mathcal{L}\{y''\} = s^2Y(s) - sy(0) - y'(0) with y(0)=0y(0)=0, y(0)=2y'(0)=2, and the standard transforms L{δ(ta)}=eas\mathcal{L}\{\delta(t-a)\} = e^{-as} and L{t2}=2/s3\mathcal{L}\{t^2\} = 2/s^3:

    s2Y2Y=4e2s+2s3s^2Y - 2 - Y = 4e^{-2s} + \frac{2}{s^3}

  3. Solve for Y(s)Y(s). Factoring YY out of the left side:

    Y(s)=2s21+4e2ss21+2s3(s21)Y(s) = \frac{2}{s^2-1} + \frac{4e^{-2s}}{s^2-1} + \frac{2}{s^3(s^2-1)}

    The three pieces are inverted separately — that is the payoff of working in the ss-domain.

  4. Break the awkward third term with partial fractions. Writing 1s3(s21)=As+Bs2+Cs3+Ds+Es21\dfrac{1}{s^3(s^2-1)} = \dfrac{A}{s} + \dfrac{B}{s^2} + \dfrac{C}{s^3} + \dfrac{Ds+E}{s^2-1} and matching coefficients gives A=1A = -1, B=0B = 0, C=1C = -1, D=1D = 1, E=0E = 0, so

    2s3(s21)=2s2s3+2ss21\frac{2}{s^3(s^2-1)} = -\frac{2}{s} - \frac{2}{s^3} + \frac{2s}{s^2-1}

  5. Invert each term. Using L1{1/(s21)}=sinht\mathcal{L}^{-1}\{1/(s^2-1)\} = \sinh t, L1{s/(s21)}=cosht\mathcal{L}^{-1}\{s/(s^2-1)\} = \cosh t, L1{2/s3}=t2\mathcal{L}^{-1}\{2/s^3\} = t^2, and the second shifting theorem L1{easF(s)}=u(ta)f(ta)\mathcal{L}^{-1}\{e^{-as}F(s)\} = u(t-a)f(t-a):

    y(t)=2sinht+2coshtt22+4u(t2)sinh(t2)y(t) = 2\sinh t + 2\cosh t - t^2 - 2 + 4u(t-2)\sinh(t-2)

  6. Simplify and check. Since 2sinht+2cosht=2et2\sinh t + 2\cosh t = 2e^{t},

    y(t)=2ett22+4u(t2)sinh(t2)y(t) = 2e^{t} - t^2 - 2 + 4u(t-2)\sinh(t-2)

    Check it: y(0)=202=0y(0) = 2 - 0 - 2 = 0 and y(0)=20=2y'(0) = 2 - 0 = 2, both correct. For t<2t<2 the Heaviside term is silent and yy=(2et2)(2ett22)=t2y'' - y = (2e^t - 2) - (2e^t - t^2 - 2) = t^2, as required. At t=2t=2 the impulse makes yy' jump by exactly 44 while yy itself stays continuous — the physical signature of a unit impulse of strength 44.

Answer

y(t)=2ett22+4u(t2)sinh(t2)y(t) = 2e^{t} - t^2 - 2 + 4u(t-2)\sinh(t-2)

Need to solve a different problem like this? Open the solver →