看完不会?让 AI 一步步解给你看
Calculus · real student question

用拉普拉斯变换求解初值问题 y'' − y = 4·delta(t − 2) + t²,其中 y(0) = 0,y'(0) = 2,delta 表示狄拉克 δ 函数。

Question

求解初值问题

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

其中 δ\delta 表示狄拉克 δ 函数。

Step-by-step solution

  1. 因为出现了冲激,所以选用拉普拉斯变换。 强迫项里含有 δ(t−2)\delta(t-2),它不是能用普通方法积分的函数;待定系数法和常数变易法都会失效。拉普拉斯变换把这个冲激变成普通的指数 e−2se^{-2s},于是整个问题化成了代数。

  2. 两边同时做变换,把初始条件一并折进去。 用 L{y′′}=s2Y(s)−sy(0)−y′(0)\mathcal{L}\{y''\} = s^2Y(s) - sy(0) - y'(0),代入 y(0)=0y(0)=0、y′(0)=2y'(0)=2,再用标准变换 L{δ(t−a)}=e−as\mathcal{L}\{\delta(t-a)\} = e^{-as} 和 L{t2}=2/s3\mathcal{L}\{t^2\} = 2/s^3:

    s2Y−2−Y=4e−2s+2s3s^2Y - 2 - Y = 4e^{-2s} + \frac{2}{s^3}

  3. 解出 Y(s)Y(s)。 把左边的 YY 提出来:

    Y(s)=2s2−1+4e−2ss2−1+2s3(s2−1)Y(s) = \frac{2}{s^2-1} + \frac{4e^{-2s}}{s^2-1} + \frac{2}{s^3(s^2-1)}

    这三块可以分别做反变换——这正是转到 ss 域里工作换来的好处。

  4. 用部分分式拆开难缠的第三项。 设 1s3(s2−1)=As+Bs2+Cs3+Ds+Es2−1\dfrac{1}{s^3(s^2-1)} = \dfrac{A}{s} + \dfrac{B}{s^2} + \dfrac{C}{s^3} + \dfrac{Ds+E}{s^2-1},比较系数得 A=−1A = -1、B=0B = 0、C=−1C = -1、D=1D = 1、E=0E = 0,于是

    2s3(s2−1)=−2s−2s3+2ss2−1\frac{2}{s^3(s^2-1)} = -\frac{2}{s} - \frac{2}{s^3} + \frac{2s}{s^2-1}

  5. 逐项做反变换。 用 L−1{1/(s2−1)}=sinh⁡t\mathcal{L}^{-1}\{1/(s^2-1)\} = \sinh t、L−1{s/(s2−1)}=cosh⁡t\mathcal{L}^{-1}\{s/(s^2-1)\} = \cosh t、L−1{2/s3}=t2\mathcal{L}^{-1}\{2/s^3\} = t^2,以及第二平移(延迟)定理 L−1{e−asF(s)}=u(t−a)f(t−a)\mathcal{L}^{-1}\{e^{-as}F(s)\} = u(t-a)f(t-a):

    y(t)=2sinh⁡t+2cosh⁡t−t2−2+4u(t−2)sinh⁡(t−2)y(t) = 2\sinh t + 2\cosh t - t^2 - 2 + 4u(t-2)\sinh(t-2)

  6. 化简并检验。 由于 2sinh⁡t+2cosh⁡t=2et2\sinh t + 2\cosh t = 2e^{t},

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

    检验一下:y(0)=2−0−2=0y(0) = 2 - 0 - 2 = 0,y′(0)=2−0=2y'(0) = 2 - 0 = 2,都对。当 t<2t<2 时阶跃项还没被激活,y′′−y=(2et−2)−(2et−t2−2)=t2y'' - y = (2e^t - 2) - (2e^t - t^2 - 2) = t^2,正合要求。到 t=2t=2 处,冲激让 y′y' 恰好跳变 44,而 yy 本身保持连续——这正是强度为 44 的单位冲激的物理特征。

Answer

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

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