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

用拉普拉斯变换求解积分微分方程 y′(t) − 2∫₀ᵗ y(τ) sin(t − τ) dτ = 1,初始条件 y(0) = −1。

Question

求解积分微分方程

y′(t)−2∫0ty(τ)sin⁡(t−τ) dτ=1,y(0)=−1.y'(t)-2\int_0^t y(\tau)\sin(t-\tau)\,d\tau = 1,\qquad y(0)=-1.

Step-by-step solution

  1. 先认出那个积分其实是卷积。 它恰好长成 ∫0ty(τ)g(t−τ) dτ=(y∗g)(t)\int_0^t y(\tau)g(t-\tau)\,d\tau=(y*g)(t) 的形状,其中 g(t)=sin⁡tg(t)=\sin t。这正是拉普拉斯变换成为趁手工具的原因:卷积定理说 L{y∗g}=Y(s)G(s)\mathcal{L}\{y*g\}=Y(s)G(s),一个麻烦的积分就此变成普通乘积。由 L{sin⁡t}=1s2+1\mathcal{L}\{\sin t\}=\tfrac{1}{s^2+1},积分那一项变换后就是 2Y(s)s2+1\tfrac{2Y(s)}{s^2+1}。

  2. 对整个方程做变换。 用 L{y′}=sY(s)−y(0)\mathcal{L}\{y'\}=sY(s)-y(0) 以及 y(0)=−1y(0)=-1,再加上 L{1}=1s\mathcal{L}\{1\}=\tfrac1s:

    sY(s)+1−2Y(s)s2+1=1s.sY(s)+1-\frac{2Y(s)}{s^2+1}=\frac{1}{s}.

  3. 用代数解出 Y(s)Y(s)。 把含 YY 的项通分合并,常数移到右边:

    Y(s)⋅s(s2+1)−2s2+1=1s−1=1−ss.Y(s)\cdot\frac{s(s^2+1)-2}{s^2+1}=\frac{1}{s}-1=\frac{1-s}{s}.

    那个三次式分解为 s3+s−2=(s−1)(s2+s+2)s^3+s-2=(s-1)(s^2+s+2),而 1−s=−(s−1)1-s=-(s-1),于是 (s−1)(s-1) 正好约掉 —— 这是运算走在正轨上的有力信号:

    Y(s)=−s2+1s (s2+s+2).Y(s)=-\frac{s^2+1}{s\,(s^2+s+2)}.

  4. 做部分分式分解。 写成 s2+1s(s2+s+2)=As+Bs+Cs2+s+2\dfrac{s^2+1}{s(s^2+s+2)}=\dfrac{A}{s}+\dfrac{Bs+C}{s^2+s+2}。去分母得到 (A+B)s2+(A+C)s+2A=s2+1(A+B)s^2+(A+C)s+2A=s^2+1,于是 2A=12A=1、A+C=0A+C=0、A+B=1A+B=1,即 A=12A=\tfrac12、B=12B=\tfrac12、C=−12C=-\tfrac12。因此

    Y(s)=−12s−s−12(s2+s+2).Y(s)=-\frac{1}{2s}-\frac{s-1}{2(s^2+s+2)}.

  5. 配方,好套上平移定理。 这个二次式有复根,所以不去分解它,而是写成 s2+s+2=(s+12)2+74s^2+s+2=\left(s+\tfrac12\right)^2+\tfrac74,并把分子也按同样的平移改写,s−1=(s+12)−32s-1=\left(s+\tfrac12\right)-\tfrac32:

    Y(s)=−12s−12⋅s+12(s+12)2+(72)2+34⋅1(s+12)2+(72)2.Y(s)=-\frac{1}{2s}-\frac12\cdot\frac{s+\tfrac12}{\left(s+\tfrac12\right)^2+\left(\tfrac{\sqrt7}{2}\right)^2}+\frac34\cdot\frac{1}{\left(s+\tfrac12\right)^2+\left(\tfrac{\sqrt7}{2}\right)^2}.

  6. 逐项做逆变换。 这三块分别对应常数、衰减余弦和衰减正弦(最后一块还要配上因子 1ω=27\tfrac{1}{\omega}=\tfrac{2}{\sqrt7},于是 34⋅27=327\tfrac34\cdot\tfrac{2}{\sqrt7}=\tfrac{3}{2\sqrt7}):

    y(t)=−12−12e−t/2cos⁡ ⁣(72t)+327e−t/2sin⁡ ⁣(72t).y(t)=-\frac12-\frac12e^{-t/2}\cos\!\left(\frac{\sqrt7}{2}t\right)+\frac{3}{2\sqrt7}e^{-t/2}\sin\!\left(\frac{\sqrt7}{2}t\right).

  7. 检验这个解。 在 t=0t=0 处正弦项为零,y(0)=−12−12=−1y(0)=-\tfrac12-\tfrac12=-1,与初始条件吻合。把闭式解代回去,在若干个 tt 处数值计算 y′(t)−2∫0ty(τ)sin⁡(t−τ) dτy'(t)-2\int_0^t y(\tau)\sin(t-\tau)\,d\tau,每次都返回 1.0000001.000000,说明这个答案满足的是原方程,而不只是初始条件。

Answer

y(t)=−12−12e−t/2cos⁡ ⁣(72t)+327e−t/2sin⁡ ⁣(72t)y(t)=-\frac{1}{2}-\frac{1}{2}e^{-t/2}\cos\!\left(\frac{\sqrt{7}}{2}t\right)+\frac{3}{2\sqrt{7}}e^{-t/2}\sin\!\left(\frac{\sqrt{7}}{2}t\right)

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