Застряли на этом? Решите с помощью ИИ
Calculus · real student question

Решите с помощью преобразования Лапласа интегро-дифференциальное уравнение y′(t) − 2∫₀ᵗ y(τ) sin(t − τ) dτ = 1 при начальном условии y(0) = −1.

Question

Solve the integro-differential equation

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, что совпадает с начальными условиями. Подстановка найденного выражения и вычисление y′(t)−2∫0ty(τ)sin⁡(t−τ) dτy'(t)-2\int_0^t y(\tau)\sin(t-\tau)\,d\tau при разных tt даёт 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 →