Решение линейного неоднородного дифференциального уравнения состоит из двух частей: решения соответствующего однородного уравнения и непосредственно решения неоднородного. Поэтому на первом шаге решается уравнение без учета правой части. Затем ищется решение в зависимости от вида правой части. И, наконец, эти решения складываются.
О решении однородного уравнения читайте в статье.
Правая часть вида Pn (x) — многочлен степени n
Решение ищется в виде
y~ =Qn (x)xs,
где Qn (x) - многочлен той же степени, что и в правой части. Показатель степени s равен количеству нулевых корней характеристического уравнения.
y′′−y′−2y=x2−3x+1
k2−k−2=0
k1 =−1;k2 =2
yo =C1 e−x+C2 e2x
Теперь рассмотрим правую часть. Она представляет собой многочлен второй степени, при этом нулевых корней у характеристического уравнения нет (s=0), следовательно, решение нелинейного уравнения ищем в виде
y~ =Ax2+Bx+C
Первая производная
y~ ′=2Ax+B
Вторая производная
y~ ′′=2A
Подставляем их в начальное уравнение
2A−(2Ax+B)−2(Ax2+Bx+C)=x2−3x+1
−2Ax2+(−2A−2B)x+(2A−B−2C)=x2−3x+1
Коэффициенты при равных степенях должны быть равны
⎩⎪⎨⎪⎧ −2A=1−2A−2B=−32A−B−2C=1
⎩⎪⎨⎪⎧ A=−21 B=23 −A=23 −(−21 )=2C=2−1+2A−B =−1
Получаем
y~ =−21 x2+2x−1
И окончательно
y=yo +y~ =C1 e−x+C2 e2x−21 x2+2x−1
Правая часть вида eαxPn
Решение ищется в виде
y~ =eαxQn (x)xs,
где Qn (x) - многочлен той же степени, что и в правой части. Показатель степени s равен количеству кратных корней характеристического уравнения, которые при это еще равны α.
y′′+3y′+2y=ex(x+1)
k1 =−2;k2 =−1
yo =C1 e−2x+C2 e−x
Теперь рассмотрим правую часть. Она содержит многочлен первой степени, при этом корней, равных единицу (в правой части α=1) у характеристического уравнения нет (s=0), следовательно, решение нелинейного уравнения ищем в виде
y~ =ex(Ax+B)
Первая производная
y~ ′=ex(Ax+B)+ex⋅A=ex(Ax+A+B)
Вторая производная
y~ ′′=ex(Ax+A+B)+ex⋅A=ex(Ax+2A+B)
Подставляем их в начальное уравнение
ex(Ax+2A+B)+3ex(Ax+A+B)+2ex(Ax+B)=ex(x+1)
(Ax+2A+B)+3(Ax+A+B)+2(Ax+B)=x+1
6Ax+5A+6B=x+1
Коэффициенты при равных степенях должны быть равны
{6A=15A+6B=1 {A=61 B=361
Получаем
y~ =ex(61 x+361 )
И окончательно
y=yo +y~ =C1 e−2x+C2 e−x+ex(61 x+361 )
Правая часть вида eαx(Pn cosβx+Qm sinβx)
Решение ищется в виде
y~ =eαx(Rl cosβx+Ql sinβx)xs
Многочлены, входящие в решение имеют максимальную степень из тех, что содержат многочлены правой части l=max(n;m). Далее рассматриваются корни характеристического уравнения, и если они комплексные и равны α±βi, где α и β из правой части, то s=1 (не забываем, что мы рассматриваем уравнения второй степени, поэтому больше 1 в данном случае s не может быть).
y′′+3y′+2y=exsin3x
k1 =−2;k2 =−1
yo =C1 e−2x+C2 e−x
Проанализировав правую часть, видим, что решение необходимо искать в виде
y~ =ex(Acos3x+Bsin3x)
Первая производная
y~ ′=ex(Acos3x+Bsin3x)+ex(−3Asin3x+3Bcos3x)=
=ex(−3Asin3x+3Bcos3x+Acos3x+Bsin3x)
Вторая производная
y~ ′′=ex(−3Asin3x+3Bcos3x+Acos3x+Bsin3x)+
+ex(−9Acos3x−9Bsin3x−3Asin3x+3Bcos3x)=
=ex(−6Asin3x+8Acos3x+8Bsin3x−6Bcos3x)
Подставляем их в начальное уравнение
−ex(−6Asin3x+8Acos3x+8Bsin3x−6Bcos3x)+
+3ex(−3Asin3x+3Bcos3x+Acos3x+Bsin3x)+
+2ex(Acos3x+Bsin3x)=exsin3x
Приводим подобные
−15Aexsin3x−3Aexcos3x−3Bexsin3x+15Bexcos3x=exsin3x
Коэффициенты при одинаковых функциях должны быть равны
{−15A−3B=1−3A+15B=0 {A=−785 B=−781
Получаем
y~ =ex(−785 cos3x−781 sin3x)
И окончательно
y=yo +y~ =C1 e−2x+C2 e−x+ex(−785 cos3x−781 sin3x)
На Студворк вы можете заказать статью по математике онлайн у профильных экспертов!
Комментарии