Линейные неоднородные дифференциальные уравнения второго порядка с постоянными коэффициентами.

Линейное неоднородное дифференциальное уравнение второго порядка с постоянными коэффициентами имеет вид

y'' + py' +qy = f(x)

Для нахождения общего решения этого уравнения нужно найти общее решение однородного уравнения, как это было показано в предыдущем параграфе, и какое-либо частное решение у* неоднородного уравнения. Их сумма будет общим решением данного неоднородного уравнения:

у = + у*.

Рассмотрим один из методов нахождения частного решения – метод неопределенных коэффициентов.

Суть метода заключается в следующем: если правая часть уравнения неоднородного уравнения имеет вид

f(x)=e^αx[P_n(x)cosβx + Q_m(x)sinβх],

где α и β – действительные числа, а P_n(x) и Q_m(x) – многочлены соответственно n - йи m – й степени с действительными коэффициентами, то частное решение у* ищется в виде

y* = x^l e^αx[M_s(x)cosβx+N_s(x)sinβx],

где M_s(x) и N_s(x) – многочлены степени s = max(n,m) с неопределенными буквенными коэффициентами, а l – кратность, с которой α + βi входит в число корней характеристического уравнения.

Следует отметить, что в общем виде многочлены соответствующей степени, имеют вид:

- многочлен 0-ой степени - А

- многочлен 1-ой степени - Ах+В

- многочлен 2-ой степени - Ах²+Вх+С

- многочлен 3-ей степени - Ах³+Вх²+Сх+D и т.д.

А, В, С, D, … - неопределенные коэффициенты.

Для того, чтобы найти неопределенные коэффициенты, частное решение у*, его производные у^*' и у^*'' подставляют в левую часть неоднородного уравнения и производят соответствующие упрощения; затем в полученном тождестве приравнивают коэффициенты при подобных членах в левой и правой частях, что дает систему уравнений относительно искомых неопределенных коэффициентов, решив которую находят эти коэффициенты.

Замечание. Если правая часть неоднородного уравнения есть сумма функций вида

f(x)=f₁(x)+f₂(x)+ … +f_n(x),

 

нужно предварительно найти частные решения y₁*, y₂*, …,y_n*, соответствующие функциям f₁(x, f₂(x), …,f_n(x). Тогда частное решение данного уравнения запишутся в виде

y*= y₁*+ y₂*+, …,+ y_n*.

 

Более общим методом решения уравнений является метод вариации произвольных постоянных.

Пусть у₁ и у₂ – линейно независимые частные решения однородного уравнения. Тогда общее решение неоднородного уравнения следует искать в виде

у=С₁(х)у₁+С₂(х)у₂

где функции С₁(х) и С₂(х) определяются из системы уравнений

 

 

Решая эту систему, получим

 

,

где - определитель Вронского.

Интегрируя С^'₁(х) и С^'₂(х) получаем

 

,

откуда, подставляя найденные функции в функцию решения, найдем общее решение линейного неоднородного уравнения.

Пример 10.12. Найти общее решение дифференциального уравнения

y"+4y'=-2xe^-4х.

Решение. Общее решение неоднородного уравнения равно сумме общего решения однородного уравнения и какого-либо частного решения у* неоднородного уравнения:

у = + у*.

Найдем общее решение однородного уравнения

y"+4y'=0.

Составим характеристическое уравнение и решим его

k² + 4k = 0, k₁=0, k₂=-4.

Общее решение однородного уравнения имеет вид

= С₁ +C₂e^-4х.

Перейдем к отысканию частного решения у* данного уравнения. Здесь правая часть неоднородного уравнения f(x)=-2xe^-4x содержит произведение многочлена первой степени и показательной функции, а также коэффициент в показателе степени совпадает с одни из корней характеристического уравнения, следовательно, необходимо в частное решение добавить множитель х и частное решение будт иметь вид

y*= x(Ax+B)e^-4x=(Ax²+Bx)e^-4x.

Вычислим производные функции у*

y*'= 2Axe^-4x-4e^-4x(Ax²+Bx) = (-4Ax²+2Ax-4Bx+B)e^-4x

y*" = (-8Ax+2A-4B) e^-4x-4e^-4x(-4Ax²+2Ax -4Bx+B)=

= (16Ax²-16Ax +16Bx+2A – 8B)e^-4x.

Подставим найденные производные в неоднородное уравнение

(16Ax²-16Ax +16Bx+2A – 8B)e^-4x+4(-4Ax²+2Ax-4Bx+B)e^-4x=-2xe^-4x

(16Ax²-16Ax +16Bx+2A – 8B+16Ax²+8Ax-16Bx+4B)e^-4x=-2xe^-4x.

Приведем подобные члены и сократим равенство на е^-4х

-8Ax+2A -4B =-2x.

Приравняем коэффициенты левой и правой части равенства, стоящие перед х в одинаковой степени, получим систему уравнений, решив которую найдем значения неопределенных коэффициентов А и В:

Получили частное решение неоднородного уравнения

y*=

Теперь можно записать общее решение данного неоднородного уравнения

y = С₁ +C₂e^-4х+ .

Пример 10.13. Найти общее решение дифференциального уравнения

y"+4y=3xcosx.

Решение. Общее решение неоднородного уравнения равно сумме общего решения однородного уравнения и какого-либо частного решения у* неоднородного уравнения:

у = + у*.

Найдем общее решение однородного уравнения

y"+4y=0.

Составим характеристическое уравнение и решим его

k² + 4 = 0, k_1,2=±2i (α=0, β=2)

Общее решение однородного уравнения имеет вид

= С₁cos2x +C₂ sin2x.

Перейдем к отысканию частного решения у* данного уравнения. Здесь правая часть неоднородного уравнения f(x)=3xcosx содержит произведение многочлена первой степени и тригонометрической функции (α =0, β =1), коэффициенты в правой части неоднородного уравнения не совпадают с корнями характеристического уравнения и частное решение будет иметь вид

y*= (Ax+B)cosx+(Cx+D)sinx,

где A,B,C,D – неопределенные коэффициенты.

Вычислим производные функции у*:

y*'=Acosx - (Ax+B)sinx + Csinx + (Cx+D)cosx=(Cx + D + A)cosx +

+(-Ax +C - B)sinx

y*" =Ccosx - (Cx + D +A)sinx – Asinx+(-Ax+C-B)cosx=(-Ax + 2C – B)cosx + +(-Cx – 2A – D)sinx

Подставим найденные производные в неоднородное уравнение, приведем подобные члены и сгруппируем члены при cosx и sinx

(3Ax + 3B + 2C)cosx + (3Cx + 3D - 2A)sinx = 3xcosx

3Axcosx + 3Cxsinx + (3B+2C)cosx + (3D-2A)sinx = 3xcosx

Приравняем коэффициенты левой и правой части равенства, стоящие перед xcosx, xsinx, cosx и sinx получим систему уравнений, решив которую найдем значения неопределенных коэффициентов A,B,C,D:

 

 

Таким образом, частное решение имеет вид

y*=xcosx+ sinx.

Окончательно получаем общее решение неоднородного дифференциального уравнения

у = С₁cos2x +C₂ sin2x+ xcosx+ sinx.

Задачи контрольной работы

В задачах 10.1- 10.20 решить дифференциальные уравнения первого порядка.

10.1. ; 10. 2. ; 10.3. ;

10.4. ; 10. 5. ; 10.6. ;

10.7. ; 10.8. ; 10.9. ;

10.10. ; 10.11. ; 10.12. ;

10.13. ; 10.14. ; 10.15. ;

16. ; 10.17. ; 10.18. ;

10.19. ; 10.20. .

 

В задачах 10.21 – 10.40 найти частное решение однородного дифференциального уравнения второго порядка с постоянными коэффициентами, удовлетворяющее заданным начальным условиям.

 

10.21. y'' – 7y' + 10y = 0; y(0) = 2; y'(0) = -1.

10.22. y'' + 2y' + 10y = 0; y() = 0; y'() = 1.

10.23. y'' – 6y' + 9y = 0; y(0) = 1; y' (0) = 0.

10.24. y'' + 8y' + 7y = 0; y(0) = 2; y'(0) = 1.

10.25. y'' + 9y = 0; y(π) = 0; y'(π) = -1.

10.26. y'' – 7y' + 12y = 0; y(0) = 2; у'(0)=-2.

10.27. y'' + 9y' = 0; y(0) = 1; y'(0) = -3.

10.28. y'' – 3y' + 2y = 0; y(0) = 0; y'(0) = 1.

10.29. y'' – 5y' + 6y = 0; y(0) = 5; y'(0) = 0.

10.30. y'' – 2y' + 5y = 0; y(0) = -1; y'(0) = 0.

10.31. y'' + 16y = 0; y(π) = -1; y'(π) = 0.

10.32. y'' + 10y' + 25y = 0; y(0) = 1; y'(0) = 1.

10.33. y'' – 6y' = 0; y(0) = 2; y'(0) = -2.

10.34. y'' – 4y' + 4y = 0; y(0) = 1; y'(0) = 3.

10.35. y'' – 8y' + 15y = 0; y(0) = 1; y'(0) = -2.

10.36. y'' – 4y' + 17y = 0; y() = 0; y'() = 1.

10.37. y'' – 2y' + y = 0; y(1) = 0; у'(1)=2

10.38. y'' + y = 0; y(π) = -1; y'(π) = -4.

10.39. y'' – 7y' + 6y = 0; y(0) = 2; y'(0) = 0.

10.40. y'' + 8y' + 16y = 0; y(0) = 1; y'(0) = 0.

 

В задачах 10.41 – 10.60 найти общее решение однородного дифференциального уравнения второго порядка с постоянными коэффициентами.

 

10.41. у'' - 2у' = 3х² + 1 10.51. у'' + 2у' = х² -3х +1

10.42. у'' - 5у' + 6у = 2хе^-х 10.52. у'' - 5у' - 24у = (2х + 3)е^х

10.43. у'' + 8у' = (х-1)е^2х 10.53. у'' - 2у' - 3у = 8е^3х

10.44. у'' - 6у' + 8у = 3е^4х 10.54. у'' + 2у' - 3у =- 2е^3х

10.45. у'' - 2у' - 3у = хе^-х 10.55. у'' + 8у' = (х²+1) е^-х

10.46. у'' + у' - 2у = (х = 2)е^-2х 10.56. у'' + 4у' + 3у = -хе^-х

10.47. у'' + 2у' - 8у = (3х+1)е^2х 10.57. у'' - 2у' - 3у =(х + 2)е^-х

10.48. у'' + 7у' = 2х² + х 10.58. у'' + у' + 6у = 2(х – 1)е^2х

10.49. у'' - у' = 8х² е^х 10.59. у'' - 4у' = 2х² – 3х +1

10.50. у'' + 3у' -10у = 2х²е^х 10.60. у'' - 5у' + 6у = 2хе^3х