Основы решения линейных неоднородных дифференциальных уравнений второго порядка (ЛНДУ-2) с постоянными коэффициентами (ПК)
ЛНДУ 2-го порядка с постоянными коэффициентами $p$ и $q$ имеет вид $y''+p\cdot y'+q\cdot y=f\left(x\right)$, где $f\left(x\right)$ - непрерывная функция.
В отношении ЛНДУ 2-го с ПК справедливы два следующих утверждения.
Предположим, что некоторая функция $U$ является произвольным частным решением неоднородного дифференциального уравнения. Предположим также, что некоторая функция $Y$ является общим решением (ОР) соответствующего линейного однородного дифференциального уравнения (ЛОДУ) $y''+p\cdot y'+q\cdot y=0$. Тогда ОР ЛНДУ-2 равно сумме указанных частного и общего решений, то есть $y=U+Y$.
Если правая часть ЛНДУ 2-го порядка представляет собой сумму функций, то есть $f\left(x\right)=f_{1} \left(x\right)+f_{2} \left(x\right)+...+f_{r} \left(x\right)$, то сначала можно найти ЧР $U_{1} ,U_{2} ,...,U_{r} $, которые соответствуют каждой из функций $f_{1} \left(x\right),f_{2} \left(x\right),...,f_{r} \left(x\right)$, а уже после этого записать ЧР ЛНДУ-2 в виде $U=U_{1} +U_{2} +...+U_{r} $.
Решение ЛНДУ 2-го порядка с ПК
Очевидно, что вид того или иного ЧР $U$ данного ЛНДУ-2 зависит от конкретного вида его правой части $f\left(x\right)$. Простейшие случаи поиска ЧР ЛНДУ-2 сформулированы в виде четырех следующих правил.
Правило № 1.
Правая часть ЛНДУ-2 имеет вид $f\left(x\right)=P_{n} \left(x\right)$, где $P_{n} \left(x\right)=a_{0} \cdot x^{n} +a_{1} \cdot x^{n-1} +...+a_{n-1} \cdot x+a_{n} $, то есть называется многочленом степени $n$. Тогда его ЧР $U$ ищут в виде $U=Q_{n} \left(x\right)\cdot x^{r} $, где $Q_{n} \left(x\right)$ - другой многочлен той же степени, что и $P_{n} \left(x\right)$, а $r$ - количество корней характеристического уравнения соответствующего ЛОДУ-2, равных нулю. Коэффициенты многочлена $Q_{n} \left(x\right)$ находят методом неопределенных коэффициентов (НК).
Правило № 2.
Правая часть ЛНДУ-2 имеет вид $f\left(x\right)=e^{\alpha \cdot x} \cdot P_{n} \left(x\right)$, где $P_{n} \left(x\right)$ представляет собой многочлен степени $n$. Тогда его ЧР $U$ ищут в виде $U=Q_{n} \left(x\right)\cdot x^{r} \cdot e^{\alpha \cdot x} $, где $Q_{n} \left(x\right)$ - другой многочлен той же степени, что и $P_{n} \left(x\right)$, а $r$ - количество корней характеристического уравнения соответствующего ЛОДУ-2, равных $\alpha $. Коэффициенты многочлена $Q_{n} \left(x\right)$ находят методом НК.
Правило № 3.
Правая часть ЛНДУ-2 имеет вид $f\left(x\right)=a\cdot \cos \left(\beta \cdot x\right)+b\cdot \sin \left(\beta \cdot x\right)$, где $a$, $b$ и $\beta $ - известные числа. Тогда его ЧР $U$ ищут в виде $U=\left(A\cdot \cos \left(\beta \cdot x\right)+B\cdot \sin \left(\beta \cdot x\right)\right)\cdot x^{r} $, где $A$ и $B$ - неизвестные коэффициенты, а $r$ - количество корней характеристического уравнения соответствующего ЛОДУ-2, равных $i\cdot \beta $. Коэффициенты $A$ и $B$ находят методом НК.
Правило № 4.
Правая часть ЛНДУ-2 имеет вид $f\left(x\right)=e^{\alpha \cdot x} \cdot \left[P_{n} \left(x\right)\cdot \cos \left(\beta \cdot x\right)+P_{m} \left(x\right)\cdot \sin \left(\beta \cdot x\right)\right]$, где $P_{n} \left(x\right)$ - многочлен степени $n$, а $P_{m} \left(x\right)$ - многочлен степени $m$. Тогда его ЧР $U$ ищут в виде $U=e^{\alpha \cdot x} \cdot \left[Q_{s} \left(x\right)\cdot \cos \left(\beta \cdot x\right)+R_{s} \left(x\right)\cdot \sin \left(\beta \cdot x\right)\right]\cdot x^{r} $, где $Q_{s} \left(x\right)$ и $R_{s} \left(x\right)$ - многочлены степени $s$, число $s$ - максимальное из двух чисел $n$ и $m$, а $r$ - количество корней характеристического уравнения соответствующего ЛОДУ-2, равных $\alpha +i\cdot \beta $. Коэффициенты многочленов $Q_{s} \left(x\right)$ и $R_{s} \left(x\right)$ находят методом НК.
Метод НК состоит в применении следующего правила. Для того чтобы найти неизвестные коэффициенты многочлена, которые входят в состав частного решения неоднородного дифференциального уравнения ЛНДУ-2, необходимо:
- подставить ЧР $U$, записанное в общем виде, в левую часть ЛНДУ-2;
- в левой части ЛНДУ-2 выполнить упрощения и сгруппировать члены с одинаковыми степенями $x$;
- в полученном тождестве приравнять коэффициенты при членах с одинаковыми степенями $x$ левой и правой частей;
- решить полученную систему линейных уравнений относительно неизвестных коэффициентов.
Задача: найти ОР ЛНДУ-2 $y''-3\cdot y'-18\cdot y=\left(36\cdot x+12\right)\cdot e^{3\cdot x} $. Найти также ЧР, удовлетворяющее начальным условиям $y=6$ при $x=0$ и $y'=1$ при $x=0$.
Записываем соответствующее ЛОДУ-2: $y''-3\cdot y'-18\cdot y=0$.
Характеристическое уравнение: $k^{2} -3\cdot k-18=0$. Корни характеристического уравнения: $k_{1} =-3$, $k_{2} =6$. Эти корни действительны и различны. Таким образом, ОР соответствующего ЛОДУ-2 имеет вид: $Y=C_{1} \cdot e^{-3\cdot x} +C_{2} \cdot e^{6\cdot x} $.
Правая часть данного ЛНДУ-2 имеет вид $\left(36\cdot x+12\right)\cdot e^{3\cdot x} $. В ней необходимо рассматривать коэффициент показателя степени экспоненты $\alpha =3$. Этот коэффициент не совпадает ни с одним из корней характеристического уравнения. Поэтому ЧР данного ЛНДУ-2 имеет вид $U=\left(A\cdot x+B\right)\cdot e^{3\cdot x} $.
Будем искать коэффициенты $A$, $B$ методом НК.
Находим первую производную ЧР:
$U'=\left(A\cdot x+B\right)^{{'} } \cdot e^{3\cdot x} +\left(A\cdot x+B\right)\cdot \left(e^{3\cdot x} \right)^{{'} } =$
$=A\cdot e^{3\cdot x} +\left(A\cdot x+B\right)\cdot 3\cdot e^{3\cdot x} =\left(A+3\cdot A\cdot x+3\cdot B\right)\cdot e^{3\cdot x} .$
Находим вторую производную ЧР:
$U''=\left(A+3\cdot A\cdot x+3\cdot B\right)^{{'} } \cdot e^{3\cdot x} +\left(A+3\cdot A\cdot x+3\cdot B\right)\cdot \left(e^{3\cdot x} \right)^{{'} } =$
$=3\cdot A\cdot e^{3\cdot x} +\left(A+3\cdot A\cdot x+3\cdot B\right)\cdot 3\cdot e^{3\cdot x} =\left(6\cdot A+9\cdot A\cdot x+9\cdot B\right)\cdot e^{3\cdot x} .$
Подставляем функции $U''$, $U'$ и $U$ вместо $y''$, $y'$ и $y$ в данное ЛНДУ-2 $y''-3\cdot y'-18\cdot y=\left(36\cdot x+12\right)\cdot e^{3\cdot x}. $ При этом, поскольку экспонента $e^{3\cdot x} $ входит как множитель во все составляющие, то её можно опустить. Получаем:
$6\cdot A+9\cdot A\cdot x+9\cdot B-3\cdot \left(A+3\cdot A\cdot x+3\cdot B\right)-18\cdot \left(A\cdot x+B\right)=36\cdot x+12.$
Выполняем действия в левой части полученного равенства:
$-18\cdot A\cdot x+3\cdot A-18\cdot B=36\cdot x+12.$
Применяем метод НК. Получаем систему линейных уравнений с двумя неизвестными:
$-18\cdot A=36;$
$3\cdot A-18\cdot B=12.$
Решение этой системы таково: $A=-2$, $B=-1$.
ЧР $U=\left(A\cdot x+B\right)\cdot e^{3\cdot x} $ для нашей задачи выглядит следующим образом: $U=\left(-2\cdot x-1\right)\cdot e^{3\cdot x} $.
ОР $y=Y+U$ для нашей задачи выглядит следующим образом: $y=C_{1} \cdot e^{-3\cdot x} +C_{2} \cdot e^{6\cdot x} +\left(-2\cdot x-1\right)\cdot e^{3\cdot x} $.
С целью поиска ЧР, удовлетворяющего заданным начальным условиям, находим производную $y'$ ОР:
$y'=-3\cdot C_{1} \cdot e^{-3\cdot x} +6\cdot C_{2} \cdot e^{6\cdot x} -2\cdot e^{3\cdot x} +\left(-2\cdot x-1\right)\cdot 3\cdot e^{3\cdot x} .$
Подставляем в $y$ и $y'$ начальные условия $y=6$ при $x=0$ и $y'=1$ при $x=0$:
$6=C_{1} +C_{2} -1; $
$1=-3\cdot C_{1} +6\cdot C_{2} -2-3=-3\cdot C_{1} +6\cdot C_{2} -5.$
Получили систему уравнений:
$C_{1} +C_{2} =7;$
$-3\cdot C_{1} +6\cdot C_{2} =6.$
Решаем её. Находим $C_{1} $ по формуле Крамера, а $C_{2} $ определяем из первого уравнения:
$C_{1} =\frac{\left|\begin{array}{cc} {7} & {1} \\ {6} & {6} \end{array}\right|}{\left|\begin{array}{cc} {1} & {1} \\ {-3} & {6} \end{array}\right|} =\frac{7\cdot 6-6\cdot 1}{1\cdot 6-\left(-3\right)\cdot 1} =\frac{36}{9} =4; C_{2} =7-C_{1} =7-4=3.$
Таким образом, ЧР данного дифференциального уравнения имеет вид: $y=4\cdot e^{-3\cdot x} +3\cdot e^{6\cdot x} +\left(-2\cdot x-1\right)\cdot e^{3\cdot x} $.