Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Дифференциальные уравнения.docx
Скачиваний:
2
Добавлен:
01.07.2025
Размер:
122 Кб
Скачать
☆

1. Однородные линейные уравнения

Нахождение ФСР однородного уравнения начинается с решения характеристического уравнения. Характеристическое уравнение имеет вид:

http://www.cyberforum.ru/cgi-bin/latex.cgi?<br /> k^{n}+p_1k^{n-1}+p_2k^{n-2}+...+p_{n-1}=0 т.е. попросту заменяем производные i-ого порядка на http://www.cyberforum.ru/cgi-bin/latex.cgi?<br /> k^i .

Решив это уравнение получим корни http://www.cyberforum.ru/cgi-bin/latex.cgi?k_1,\;k_2,...,\;k_n . Тут возможны следующие варианты:

  • Все корни - действительные числа, кратность каждого корня - единица.

  • Все корни - действительные числа, но есть корни с кратностью большей, чем 1.

  • Среди корней есть комплексные.

***

Первый случай довольно тривиален. Базисные функции будут иметь вид http://www.cyberforum.ru/cgi-bin/latex.cgi?\varphi_i(t)=e^{k_it} .

Общее решение уравнения http://www.cyberforum.ru/cgi-bin/latex.cgi?y=C_1e^{k_1t}+C_2e^{k_2t}+...+C_ne^{k_nt}

Пример:

http://www.cyberforum.ru/cgi-bin/latex.cgi?<br /> y''-4y'+3y=0,\;\;k^2-4k+3=0,\;\;k_1=1, \;k_2=3,\;\Rightarrow y=C_1e^t+C_2e^{3t}

***

Второй уже посложнее. Для всех однократных корней соответствующие им базисные функции будут иметь вид, указанный выше. Корень кратности m даст нам m базисных функций вида http://www.cyberforum.ru/cgi-bin/latex.cgi?\varphi_i(t)=t^{p}e^{k_it}, \;p=\bar{0, \;m-1} .

Примеры:

http://www.cyberforum.ru/cgi-bin/latex.cgi?<br /> 1. \;y''-4y'+4y=0,\;\;k^2-4k+4=0,\;\;k_1=k_2=2,\;\Rightarrow y=C_1e^{2t}+C_2te^{2t}\\<br /> 2.\;y'''-5y''+8y'-4y=0,\;\;k^3-5k^2+8k-4=0,\;\;k_1=k_2=2,\;k_3=1\;\Rightarrow y=C_1e^{2t}+C_2te^{2t}+C_3e^t

***

Если же есть комплексный корень http://www.cyberforum.ru/cgi-bin/latex.cgi?<br /> k_i , то, очевидно, среди корней есть сопряженный к нему http://www.cyberforum.ru/cgi-bin/latex.cgi?<br /> k_{i+1}=\bar{k_i} . Пара таких корней даст две базисные функции http://www.cyberforum.ru/cgi-bin/latex.cgi?\varphi_i(t)=e^{\operatorname{Re}(k_i)t}\cos (\operatorname{Im}(k_i)t), \;\varphi_{i+1}(t)=e^{\operatorname{Re}(k_i)t}\sin (\operatorname{Im}(k_i)t)

Если пара имеет кратность m>1, то, аналогично предыдущему случаю, получим m пар базисных функций http://www.cyberforum.ru/cgi-bin/latex.cgi?\varphi_i(t)=t^pe^{\operatorname{Re}(k_i)t}\cos (\operatorname{Im}(k_i)t), \;\varphi_{i+1}(t)=t^pe^{\operatorname{Re}(k_i)t}\sin (\operatorname{Im}(k_i)t), \;p=\bar{0, \;m-1}

Пример:

http://www.cyberforum.ru/cgi-bin/latex.cgi?y''-4y'+13y=0,\;k^2-4k+13=0,\;k_1=2+3i,\;k_2=2-3i,\;y=C_1e^{2t}\cos 3t+C_2e^{2t}\sin 3t

2. Неоднородные уравнения со специальной правой частью.

Уравнения со специальной правой частью - это неоднородные линейные уравнения, в которых функция f(t) имеет вид:

http://www.cyberforum.ru/cgi-bin/latex.cgi?e^{\alpha t}(P_n(t)\cos \beta t + Q_m(t)\sin \beta t)

где http://www.cyberforum.ru/cgi-bin/latex.cgi?P_n(t),\;Q_m(t) - полиномы степеней n и m соответственно.

Возможны следующие варианты:

  • http://www.cyberforum.ru/cgi-bin/latex.cgi?\lambda_{1,2}=\alpha \pm i\beta - не являются корнями характеристического уравнения.

  • - пара характеристических корней кратности p;

В первом случае частное решение будет иметь вид: http://www.cyberforum.ru/cgi-bin/latex.cgi?e^{\alpha t}(M_n(t)\cos \beta t+N_m(t)\sin \beta t) , где http://www.cyberforum.ru/cgi-bin/latex.cgi?M_n(t),\;N_m(t) - некоторые полиномы степеней n и m соответственно. Коэффициенты полиномов при степенях t еще требуется определить. Сделать это можно так: надо подставить наше частное решение в уравнение и привести подобные. Затем получившиеся коэффициенты приравнять к нулю.

Примеры:

http://www.cyberforum.ru/cgi-bin/latex.cgi?<br /> 1)\;y''-2y'+1=(t^2+1)e^{3t},\; k^2-2k+1=0,\;k_1=k_2=1,\;\lambda = 3\\\;\;\;\bar{y}=(At^2+Bt+C)e^{3t},\;\bar{y}'=e^{3t}(3At^2+(3B+2A)t+3C+B)\\\;\;\;\bar{y}''=e^{3t}(9At^2+(9B+12A)t+9C+6B+2A)\\\;\;\;\bar{y}''-2\bar{y}'+\bar{y}=e^{3t}(9At^2+(9B+12A)t+9C+6B+2A-2(3At^2+(3B+2A)t+3C+B)+At^2+Bt+C)=\\\;\;\;=e^{3t}(4At^2+(4B+8A)t+4C+4B+2A)=(t^2+1)e^{3t}\\\;\;\;4At^2-t^2+(4B+8A)t+4C+4B+2A-1=0\\\;\;\;x^2:\;4A-1=0\;\;\;x:\;4B+8A=0\;\;\;x^0:\;4C+4B+2A-1=0A=0.25,\;B=-0.5,\;C=0.625\Rightarrow \bar{y}=(0.25t^2-0.5t+0.625)e^{3t}\\\;\;\;y=C_1e^{t}+C_2te^{t}+(0.25t^2-0.5t+0.625)e^{3t}

http://www.cyberforum.ru/cgi-bin/latex.cgi?<br /> 2)\;y''+y=e^t\cos t,\;k^2+1=0,\;k_1=i,\;k_2=-i,\;\lambda=1\pm i\\\;\;\;\bar{y}=e^t(A\cos t +B\sin t),\;\;\bar{y}''=e^t(2B\cos t - 2A\sin t)\\\;\;\;\bar{y}''+\bar{y}=e^t(A\cos t +B\sin t+2B\cos t - 2A\sin t)=e^t\cos t,\;A\cos t +B\sin t+2B\cos t - 2A\sin t-\cos t = 0\\\;\;\;\cos t: \;A+2B-1=0\\\;\;\;\sin t:\;B-2A=0\\\;\;\;A=0.2,\;B=0.4\;\bar{y}=e^t(0.2\cos t+0.4\sin t)\;\;\;y=C_1\cos t + C_2 \sin t + e^t(0.2\cos t+0.4\sin t)

***

Второй случай приводит к несколько другому частному решению:

http://www.cyberforum.ru/cgi-bin/latex.cgi?t^pe^{\alpha t}(M_n(t)\cos \beta t+N_m(t)\sin \beta t) , т.е. надо еще умножить на t^p, p - кратность корня. Почему так получается, можете посмотреть в книгах Краснова, Камке и др.

Поиск коэффициентов полиномов производится также, как и в первом случае.