Calculus · real student question

ラプラス変換を使って、積分微分方程式 y′(t) − 2∫₀ᵗ y(τ) sin(t − τ) dτ = 1 を初期条件 y(0) = −1 のもとで解け。

Question

次の積分微分方程式を解け。

y(t)20ty(τ)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τ=(yg)(t)\int_0^t y(\tau)g(t-\tau)\,d\tau=(y*g)(t) の形をしており、g(t)=sintg(t)=\sin t である。だからこそラプラス変換が適切な道具になる。畳み込み定理により L{yg}=Y(s)G(s)\mathcal{L}\{y*g\}=Y(s)G(s) なので、面倒な積分が普通の積に変わる。L{sint}=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)+12Y(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=1s1=1ss.Y(s)\cdot\frac{s(s^2+1)-2}{s^2+1}=\frac{1}{s}-1=\frac{1-s}{s}.

    3次式は s3+s2=(s1)(s2+s+2)s^3+s-2=(s-1)(s^2+s+2) と因数分解でき、1s=(s1)1-s=-(s-1) なので (s1)(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=1A+C=0A+C=0A+B=1A+B=1、すなわち A=12A=\tfrac12B=12B=\tfrac12C=12C=-\tfrac12 である。よって

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

  5. 平方完成して移動定理の形に合わせる。 2次式は虚数解をもつので因数分解せず、s2+s+2=(s+12)2+74s^2+s+2=\left(s+\tfrac12\right)^2+\tfrac74 と書き、分子も同じずらし方に合わせて s1=(s+12)32s-1=\left(s+\tfrac12\right)-\tfrac32 と表す:

    Y(s)=12s12s+12(s+12)2+(72)2+341(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. 項ごとに逆変換する。 3つの項はそれぞれ定数、減衰する余弦、減衰する正弦の変換である(最後の項には因子 1ω=27\tfrac{1}{\omega}=\tfrac{2}{\sqrt7} が必要なので 3427=327\tfrac34\cdot\tfrac{2}{\sqrt7}=\tfrac{3}{2\sqrt7} となる):

    y(t)=1212et/2cos ⁣(72t)+327et/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)=1212=1y(0)=-\tfrac12-\tfrac12=-1 となり、初期条件と一致する。閉じた形を代入して y(t)20ty(τ)sin(tτ)dτy'(t)-2\int_0^t y(\tau)\sin(t-\tau)\,d\tau をいくつかの tt で数値的に評価すると、毎回 1.0000001.000000 が返る。したがってこの答えは初期条件だけでなく、もとの方程式そのものを満たしている。

Answer

y(t)=1212et/2cos ⁣(72t)+327et/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 →