Дифференциальные уравнения n-го порядка

ОГЛАВЛЕНИЕ

1.ВВЕДЕНИЕ 2

2.ОБЫКНОВЕННЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ N-ГО ПОРЯДКА. 3

3.ЛИНЕЙНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ N-ГО ПОРЯДКА. 6

3.1.ОДНОРОДНЫЕ ЛИНЕЙНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ N-ГО ПОРЯДКА. 6

3.2.ЛИНЕЙНЫЕ ОДНОРОДНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ. 7

3.3.ЛИНЕЙНЫЕ НЕОДНОРОДНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ. 11

3.4.МЕТОД ВАРИАЦИИ ПРОИЗВОЛЬНЫХ ПОСТОЯННЫХ РЕШЕНИЯ ЗАДАЧИ КОШИ ДЛЯ ЛИНЕЙНОГО НЕОДНОРОДНОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ 15

3.5.ЛИНЕЙНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ N-ГО ПОРЯДКА С ПЕРЕМЕННЫМИ КОЭФФИЦИЕНТАМИ. 17

3.6.ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ, ДОПУСКАЮЩИЕ ПОНИЖЕНИЕ ПОРЯДКА 18

Уравнение, не содержащее в явном виде неизвестную функцию и её младшие производные. 19

Уравнение, не содержащее в явном виде независимую переменную x. 20

4.НЕЛИНЕЙНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ. 22

   4.1.ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ В ЧАСТНЫХ ПРОИЗВОДНЫХ. 23

5.ЗАКЛЮЧЕНИЕ 25

6.   ИСПОЛЬЗОВАННЫЕ ИСТОЧНИКИ 26

Введение

Эффективное решение задач  во многих областях  невозможно без  описания их дифференциальными уравнениями  различного порядка. Для получения  результатов необходимы доступные  и проверенные методы решения  дифференциальных уравнений. Особое место  занимают дифференциальные уравнения n-го порядка, поэтому тема курсовой работы представляет интерес и является актуальной.

Целью данной работы является описание методов решения дифференциальных уравнений n-го порядка с предварительной классификацией по видам уравнений.

Для достижения поставленной цели в рамках работы необходимо решить следующие задачи:

 

1. Провести анализ существующих источников и привести определения и теоретические сведения по дифференциальным уравнениям n-го порядка.
2. На основе проведенного анализа описать методы решений уравнений с наглядными примерами решений.

 

 

 

 

 

 

 

 

 

 

 

 

 

ОБЫКНОВЕННЫЕ  ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ N-ГО ПОРЯДКА.

Обыкновенным дифференциальным уравнением n-го порядка называется уравнение вида 
F(x, y, y', y'', ..., y⁽ⁿ⁾) = 0,  
где F - известная функция (n+2) переменных, определенная в области D М Rⁿ⁺², x - независимая переменная из интервала (a, b), y = y(x) - неизвестная функция, n - порядок уравнения.

Пример 1. Уравнение движения материальной точки.

Движение материальной точки  массы m под действием внешних сил F описывается вторым законом Ньютона ma = F .  
Пусть точка движется по оси 0x , тогда функция x = x(t) — абсцисса точки в момент времени t , удовлетворяет обыкновенному дифференциальному уравнению 2-го порядка mx'' = F(t, x, x'). 
Например, уравнение 
x'' + w²x=0 
— уравнение гармонического осциллятора, описывает периодические колебания материальной точки с периодом T=2p/w.

Пример 2. Уравнение изменения объема производства в замкнутой экономической системе.

Изменение объема производства в некоторой замкнутой экономической  системе описывает дифференциальное уравнение второго порядка 
y''+2|k|y'+w²y=0.  
В замкнутой экономической системе нет экспорта, импорта и притока капитала извне. Уравнение описывает поведение разности 
y(x)=Y(x)-G/s 
между объемом производства Y(x) и фиксированной величиной G/s отношения правительственных расходов к предельной склонности населения к сбережению. 
Ниже приведен график решения уравнения при k=0.25, w²=0.25. 
Колебания решения уравнения около нуля соответствуют периодам спада и подъема в экономике.  

В дальнейшем будем рассматривать  обыкновенные дифференциальные уравнения, разрешенные относительно старшей  производной - уравнения, записанные в нормальной форме: 
y⁽ⁿ⁾) = f(x, y, y', y'', ..., y^(n-1)).

Функция y(x) называется решением дифференциального уравнения n-го порядка, если она n раз непрерывно дифференцируема на промежутке (a, b) и удовлетворяет  уравнению для всех  
xÎ (a, b).

Пример 3. Интегральная кривая для уравнения затухающих колебаний.

Уравнение второго порядка 
x'' + 2ax' + bx = 0 
при a²<b описывает затухающие (ангармонические) колебания.  
Ниже приведен график решения уравнения для a=0.1, b=1.

Дифференциальное  уравнение n-го порядка имеет, вообще говоря, бесконечное множество решений. Чтобы выделить единственное решение уравнения достаточно определить начальные условия: 
y(x₀) = y₀ , y'(x₀) = y_0,1 , y''(x₀) = y_0,2 , ..., y^(n-1)(x₀) = y_0,n-1.

При определенных ограничениях на правую часть уравнения эта задача, она  называется задачей Коши, имеет единственное решение.

Пример4. Решения задачи Коши для уравнения изменения объема производства в замкнутой экономической системе при различных начальных условиях.

Изменение объема производства в некоторой замкнутой экономической  системе описывает дифференциальное уравнение второго порядка 
y'' +2|k|y' + w²y = 0. 
В замкнутой экономической системе нет экспорта, импорта и притока капитала извне. Уравнение описывает поведение разности y(x)=Y(x)-G/s 
между объемом производства Y(x) и фиксированной величиной G/s отношения правительственных расходов к предельной склонности населения к сбережению. 
Ниже приведены графики решений уравнения при k=0.25, w²=0.25 при различных начальных условиях. Видно, что колебания решений около нуля — периоды спада и подъема в экономике — зависят от начального состояния системы. 

 

Справедлива следующая теорема существования и единственности решения задачи Коши.  
Если правая часть уравнения  
y⁽ⁿ⁾) = f(x, y, y', y'', ..., y^(n-1)) 
и ее частные производные по переменным y, y', y'', ..., y^(n-1) непрерывны в области G Rⁿ⁺¹, то для любой точки (x₀, y₀, y_0,1, y_0,2, ..., y_0,n-1) из G на некотором интервале (x₀-h, x₀+h) существует единственное решение y(x) уравнения, удовлетворяющее начальным условиям 
y(x₀) = y₀ , y'(x₀) = y_0,1 , y''(x₀) = y_0,2 , ..., y^(n-1)(x₀) = y_0,n-1.

Численное решение задачи Коши 
y⁽ⁿ⁾) = f(x, y, y', y'', ..., y^(n-1)), 
y(x₀) = y₀ , y'(x₀) = y_0,1 , y''(x₀) = y_0,2 , ..., y^(n-1)(x₀) = y_0,n-1  
состоит в построении таблицы приближенных значений y_i решения y=y(x) в точках x₁, x₂, ..., x_i, ... .

Задача о численном решении  дифференциального уравнения порядка  выше первого чаще всего сводится к численному решению решению задачи Коши для нормальной системы дифференциальных уравнений.  
Обозначив  
y(x)=y₁(x), y'(x)=y₂(x), y''(x)=y₃(x), ..., y^(n-1)(x)=y_n (x),  
получим задачу Коши для системы n дифференциальных уравнений 1-го порядка 
y₁'=y₂ , y₂'=y₃ , ..., y_n' =f(x, y₁, y₂ , ..., y_n ),  
y₁(x₀ )=y₀, y₂(x₀)=y_1,0 , ..., y_n-1(x₀)=y_n-1,0,  
которая в векторной форме имеет вид

`Y '= `F(x,`Y), `Y(x₀) =`Y₀,

`Y (x)=(y<sub>1</sub>(x), y<sub>2</sub>(x), ..., y<sub>n</sub>(x)), `Y '(x)=(y₁'(x), y₂'(x), ..., y_n'(x)),

                                  `F(x,`Y)= (y₂, y₃, ..., y_n, f(x, y₁, y₂ , ..., y_n )).

 
Численное решение задачи Коши для этой системы состоит в построении таблицы приближенных значений y_i,1 , y_i,2 , ..., y_i,N компонент y_i(x_j) вектора решения в точках x₁ , x₂ , ..., x_N. Чтобы получить расчетные формулы метода Рунге—Кутты для систем дифференциальных уравнений, достаточно в расчетных формулах для уравнений первого порядка заменить 
y, f(x, y), k₁, k₂, k₃, k₄ на `Y, `F(x,`Y), `k₁, `k<sub>2</sub>, `k₃, `k₄ .

 

 

 

 

 

 

 

 

ЛИНЕЙНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ N-ГО ПОРЯДКА.

Линейным дифференциальным уравнением n -го порядка называется уравнение вида  
y⁽ⁿ⁾ + a₁(x) y^(n-1) + ... + a_n-1 (x) y' + a_n(x) y = f(x),  
где y = y(x) — неизвестная функция,  
a₁(x), a₂(x), ..., a_n-1(x), a_n(x),  f(x) — известные функции, которые будем полагать непрерывными на промежутке (a, b).

Выражение в левой части уравнения называется линейным дифференциальным оператором n -го порядка:  
L(y) = y⁽ⁿ⁾ + a₁(x) y^(n-1) + ... + a_n-1 (x) y' + a_n(x) y .  
Уравнения  
y⁽ⁿ⁾ + a₁(x) y^(n-1) + ... + a_n-1 (x) y' + a_n(x) y = 0 и  
y⁽ⁿ⁾ + a₁(x) y^(n-1) + ... + a_n-1 (x) y' + a_n(x) y = f(x),  f(x) № 0,  
называются соответственно однородным и неоднородным линейным дифференциальным уравнением n -го порядка.

Будем записывать однородное и неоднородное линейные дифференциальные уравнения  в виде:  
L(y) = 0 и L(y) = f(x).

Принцип суперпозиции основан на следующих свойствах решений линейных уравнений:  
а) Если  y₁(x) и y₂(x) — два решения однородного линейного уравнения L(y)=0, то их линейная комбинация y(x) = c₁ y₁(x) + c₂ y₂(x) 
при любых постоянных c₁, c₂   является решением однородного уравнения.  
б) Если y₁(x) и  y₂(x) — два решения неоднородного линейного уравнения  
L(y) = f(x), то их разность y(x) = y₁(x) - y₂(x) 
является решением однородного уравнения L(y) = 0.  
в) Любое решение неоднородного линейного уравнения L(y) = f(x) есть сумма частного (фиксированного) решения неоднородного уравнения и некоторого решения однородного уравнения.

Принцип суперпозиции: 
Если y₁(x) и  y₂(x) — решения неоднородных линейных уравнений  
L(y) =  f₁(x) и L(y) =  f₂(x), то их сумма y(x) = y₁(x) + y₂(x) является решением уравнения 
L(y) =  f₁(x) +  f₂(x).  

Пример  5. Проверка принципа суперпозиции для линейного дифференциального уравнения 2-го порядка.

Рассмотрим  два линейных дифференциальных уравнения 2-го порядка: однородное и неоднородное уравнения.  
Функции y₁(x) = lnx и y₂(x) = x — два решения однородного уравнения 
 
а функция 
 
— решение неоднородного уравнения 
.  
Подстановкой в уравнения легко проверить, что функция 
y(x) = c₁ y₁(x) + c₂ y₂(x)  
является решением однородного уравнения при любых значения констант c₁, c₂, а функция  
y(x) = c₁ y₁(x) + c₂ y₂(x) + y₃(x) — решение приведенного выше неоднородного уравнения.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЛИНЕЙНЫЕ ОДНОРОДНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ С ПОСТОЯННЫМИ  КОЭФФИЦИЕНТАМИ

Для линейного однородного  дифференциального уравнения n-го порядка

y⁽ⁿ⁾ + a₁(x) y^(n-1) + ... + a_n-1 (x) y' + a_n(x) y = 0, 

где y = y(x) — неизвестная функция, a₁(x), a₂(x), ..., a_n-1(x), a_n(x) — известные, непрерывные, справедливо:  
1) существуют n линейно независимых решений уравнения  
y₁(x), y₂(x), ..., y_n(x); 
2) при любых значениях констант c₁, c₂, ..., c_n функция 
y(x)= c₁ y₁(x) +  c₂ y₂(x) + ... + c_n y_n(x) 
является решением уравнения; 
3) для любых начальных значений x₀,  y₀,   y_0,1, ..., y_0,n-1 существуют такие значения c*₁, c*_n, ..., c*_n, что решение  
y*(x)=c*₁ y₁(x) + c*₂ y₂(x) + ... + c*_n y_n (x) 
удовлетворяет при x = x₀ начальным условиям 
y*(x₀)=y₀, (y*)'(x₀)=y_0,1 , ...,(y*)^(n-1)(x₀)=y_0,n-1.

Выражение  y(x)= c₁ y₁(x) +  c₂ y₂(x) + ... + c_n y_n(x) называется общим решением линейного однородного дифференциального уравнения n-го порядка.

Совокупность n линейно независимых решений линейного однородного дифференциального уравнения n-го порядка y₁(x), y₂(x), ..., y_n(x) называется фундаментальной системой решений уравнения.

Для линейного однородного  дифференциального уравнения с постоянными коэффициентами существует простой алгоритм построения фундаментальной системы решений. Будем искать решение уравнения в виде y(x) = exp(lx): 
exp(lx)⁽ⁿ⁾ + a₁exp(lx)^(n-1) + ... + a_n-1exp(lx)' + a_nexp(lx)= 
=  (lⁿ + a₁l^n-1 + ... + a_n-1l + a_n)exp(lx) = 0,  
т.е. число l является корнем характеристического уравнения 
lⁿ + a₁l^n-1 + ... + a_n-1l + a_n = 0. 
Левая часть характеристического уравнения называется характеристическим многочленом  линейного дифференциального уравнения: 
P(l) = lⁿ + a₁l^n-1 + ... + a_n-1l + a_n. 
Таким образом, задача о решении линейного однородного уравнения n -го порядка с постоянными коэффициентами сводится к решению алгебраического уравнения.

Если характеристическое уравнение имеет n различных действительных корней  
l₁, l₂ ,... , l_n,  
то фундаментальная система решений состоит из функций    
y₁(x) = exp(l₁x), y₂(x) = exp(l₂x), ..., y_n(x) = exp(l_nx),  
и общее решение однородного уравнения имеет вид: 
y(x)= c₁ exp(l₁x) +  c₂ exp(l₂x) + ... + c_n exp(l_nx).

Пример 6. Фундаментальная система решений и общее решение для случая простых действительных корней

Фундаментальная система  решений и общее решение для  случая простых действительных корней

Рассмотрим уравнение  y'' - 3y' + 2y = 0.  
Его характеристическое уравнение l² - 3l + 2 = 0 имеет два различных действительных корня  l₁ =1 и l₂ =2.  
Фундаментальная система решений уравнения: 
y₁ = exp(l₁x)=exp(x) и y₂ = exp(l₂x)=exp(2x)  
Общее решение уравнения: y(x) = c₁exp(x) + c₂exp(2x)

 

Если какой-либо из действительных корней характеристического уравнения  повторяется r раз (r-кратный корень), то в фундаментальной системе решений ему отвечают r функций; если  
l_k=l_k+1 = ... = l_k+r-1,  
то в фундаментальную систему решений уравнения входят r функций: 
y_k(x) = exp(l_kx), 
y_k+1(x) = xexp(l_kx),  
y_k+2(x) = x²exp(l_kx), ...,  
y_k+r-1(x) =x^r-1 exp(l_nx).

Пример 7. Фундаментальная система решений и общее решение для случая кратных действительных корней.

Фундаментальная система  решений и общее решение для  случая кратных действительных корней

Рассмотрим уравнение y''- 2y' + y = 0.  
Его характеристическое уравнение l² - 2l + 1 = 0

имеет один кратный действительный корень l ₁ = l ₂ = 1.  
Фундаментальная система решений уравнения:  y₁ = exp(x) и y₂= xexp(x) 
Общее решение уравнения: y(x) = c₁exp(x) + c₂xexp(x).

Если характеристическое уравнение имеет комплексные корни, то каждой паре простых (имеющих кратность 1 ) комплексных корней 
l_k,k+1=a_k ± ib_k  
в фундаментальной системе решений отвечает пара функций 
y_k(x) = exp(a_kx)cos(b_kx), y_k+1(x) = exp(a_kx)sin(b_kx).

Пример 8. Фундаментальная система решений и общее решение для случая п простых комплексных корней.

Фундаментальная система  решений и общее решение для  случая простых комплексных корней

Рассмотрим уравнение y'' - 2y' + 5y = 0.  
Его характеристическое уравнение  l² - 2l + 5 = 0 
имеет пару комплексно сопряженных корней  l₁ = 1-2i, l₂ = 1+ 2i. 
Фундаментальная система решений уравнения: exp(x)cos2x, exp(x)sin2x. 
Общее решение уравнения: y(x) = c₁exp(x)cos2x + c₂exp(x)sin2x.

ПРИМЕР 4. Фундаментальная система решений и общее решение для случая простых комплексных корней. Мнимые корни.

Фундаментальная система  решений и общее решение для  случая простых комплексных корней. Мнимые корни

Рассмотрим уравнение y'' + y = 0. 
Его характеристическое уравнение l² + 1 = 0 
имеет пару комплексно сопряженных корней l₁ =i, l₂ = -i. 
Фундаментальная система решений уравнения: cosx, sinx 
Общее решение уравнения: y(x) = c₁cosx + c₂sinx.

Если же комплексная пара корней имеет кратность r, то  такой паре  
l_k=l_k+1 = ... = l_2k+2r-1=a_k ± ib_k,  
в фундаментальной системе решений отвечают функции  
    exp(a_kx)cos(b_kx),          exp(a_kx)sin(b_kx), 
   xexp(a_kx)cos(b_kx),        xexp(a_kx)sin(b_kx), 
  x²exp(a_kx)cos(b_kx),      x²exp(a_kx)sin(b_kx), 
................ 
x^r-1exp(a_kx)cos(b_kx),    x^r-1exp(a_kx)sin(b_kx).

ПРИМЕР 5. Фундаментальная система решений и общее решение для случая кратных комплексных корней.

Фундаментальная система  решений и общее решение для  случая кратных  комплексных корней

Рассмотрим уравнение y''''- 4y''' + 14y'' - 20y' + 25y = 0.  
Его характеристическое уравнение l⁴- 4l³ + 14l² - 20l + 25 = 0 
имеет пару кратных комплексно сопряженных корней 
l_1,2 =1- 2i, l_3,4 = 1 + 2 i.  
Фундаментальная система решений уравнения: 
y₁ = exp(x)cos2x, y₂= exp(x)sin2x, y₃ = xexp(x)cos2x и y₄ = xexp(x)sin2x.  
Общее решение уравнения:  
y(x) = c₁exp(x)cos2x + c₂exp(x)sin2x + c₃xexp(x)cos2x + c₄xexp(x)sin2x.

Таким образом, для отыскания общего решения линейного однородного дифференциального уравнения с постоянными коэффициентами следует: 
записать характеристическое уравнение; 
найти все корни характеристического уравнения l₁, l₂, ... , l_n; 
записать фундаментальную систему решений y₁(x), y₂(x), ..., y_n(x); 
записать выражение для общего решения y(x)= c₁ y₁(x) +  c₂ y₂(x) + ... + c_n y_n(x). 
Для решения задачи Коши нужно подставить выражение для общего решения в начальные условия и определить значения постоянных c₁,..., c_n, которые являются решениями системы линейных алгебраических уравнений   
c₁ y₁(x₀) +  c₂ y₂(x₀) + ... + c_n y_n(x₀) = y₀, 
c₁ y'₁(x₀) +  c₂ y'₂(x₀) + ... + c_n y'_n(x₀)   =y_0,1, 
......... , 
c₁ y₁^(n-1)(x₀) +  c₂ y₂^(n-1)(x₀) + ... + c_n y_n^(n-1)(x₀) = y_0,n-1

ПРИМЕР 6. Решение задачи Коши.

Решение задачи Коши

Рассмотрим задачу Коши для  однородного дифференциального  уравнения 
y'' + 2y' + 3y = 0, y(0) = 1, y'(0) = 1.  
Его характеристическое уравнение l² + 2l + 3 = 0 
имеет пару комплексно сопряженных корней l₁ = -1- i, l₂ = -1 + i. 
Фундаментальная система решений содержит два решения  
exp(-x)cos x, y=exp(-x)sin x, 
его общее решение имеет вид  
y(x) = c₁exp(-x)cos x + c₂exp(-x)sin x. 
Решение задачи Коши y(0)=1, y'(0)=1 находим из условий 
y(0) = c₁exp(0)cos(0) + c₁exp(0)sin(0) = c₁ =1,  
y'(0) = -c₁exp(0)cos(0) -c₁ exp(0)sin(0) - c₂exp(0)sin(0) + c₂ exp(0)cos(0) =  
= - c₁ + c₂ =1,  откуда c₁ = 1 и c₂ = . Подставив константы в выражение для общего решения получим решение задачи Коши  
y(x) = exp(-x)cos x + exp(-x)sin x.

ЛИНЕЙНЫЕ НЕОДНОРОДНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ.

Для линейного неоднородного дифференциального уравнения n-го порядка

y⁽ⁿ⁾ + a₁(x) y^(n-1) + ... + a_n-1 (x) y' + a_n(x) y = f(x), 

где y = y(x) — неизвестная функция, a₁(x), a₂(x), ..., a_n-1(x), a_n(x), f(x) — известные, непрерывные, справедливо:  
1) если y₁(x) и y₂(x) — два решения неоднородного уравнения, то функция 
y(x) = y₁(x) - y₂(x) — решение соответствующего  однородного уравнения; 
2) если y₁(x) решение неоднородного уравнения, а  y₂(x) — решение соответствующего однородного уравнения, то функция 
y(x) = y₁(x) + y₂(x) — решение неоднородного уравнения; 
3) если  y₁(x), y₂(x), ...,  y_n(x) — n линейно независимых решений однородного уравнения, а y_ч(x) — произвольное решение неоднородного уравнения,  
то для любых начальных значений  
x₀,  y₀,   y_0,1, ..., y_0,n-1  
существуют такие значения  
c*₁, c*_n, ..., c*_n, что решение  
y*(x)=c*₁ y₁(x) + c*₂ y₂(x) + ... + c*_n y_n (x) + y_ч(x) 
удовлетворяет при x = x₀ начальным условиям 
y*(x₀)=y₀, (y*)'(x₀)=y_0,1 , ...,(y*)^(n-1)(x₀)=y_0,n-1.

Выражение  
y(x)= c₁ y₁(x) +  c₂ y₂(x) + ... + c_n y_n(x)   + y_ч(x)  
называется общим решением линейного неоднородного дифференциального уравнения n-го порядка.

Для отыскания частных  решений неоднородных дифференциальных уравнений с постоянными коэффициентами с правыми частями вида:  
P_k(x)exp(ax)cos(bx) + Q_m(x)exp(ax)sin(bx),  
где P_k(x), Q_m(x) — многочлены степени k и m соответственно, существует простой алгоритм построения частного решения, называемый методом подбора.

Метод подбора, или метод  неопределенных коэффициентов, состоит  в следующем.  
Искомое решение уравнения записывается в виде: 
(P_r(x)exp(ax)cos(bx) + Q_r(x)exp(ax)sin(bx))x^s,  
где P_r(x), Q_r(x) — многочлены степени r = max(k, m) с неизвестными коэффициентами 
p_r , p_r-1, ..., p₁, p₀, q_r, q_r-1, ..., q₁, q₀.  
Сомножитель x^s называют резонансным сомножителем. Резонанс имеет место в случаях, когда среди корней характеристического уравнения есть корень 
l =a ± ib  кратности s. 
Т.е. если среди корней характеристического уравнения соответствующего однородного уравнения есть такой, что его действительная часть совпадает с коэффициентом в показателе степени экспоненты, а мнимая — с коэффициентом в аргументе тригонометрической функции в правой части уравнения, и кратность этого корня s, то в искомом частном решении присутствует резонансный сомножитель x^s. Если же такого совпадения нет (s=0),  то резонансный сомножитель отсутствует.

Подставив выражение для  частного решения  в левую часть уравнения, получим обобщенный многочлен того же вида, что и многочлен в правой части уравнения, коэффициенты которого неизвестны.  
Два обобщенных многочлена равны тогда и только тогда, когда равны коэффициенты при сомножителях вида x^texp(ax)sin(bx), x^texp(ax)cos(bx)  с одинаковыми степенями t.  
Приравняв коэффициенты при таких сомножителях, получим систему 2(r+1) линейных алгебраических уравнений относительно 2(r+1) неизвестных. Можно показать, что такая система совместна и имеет единственное решение.

Пример 8.

Решить  уравнение

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

общее решение  соответствующего однородного уравнения,   частное решение.

Решим соответствующее однородное уравнение 

Характеристическое  уравнение:

 

 

Частное решение ищем в виде:

 

Подставим эти производные в данное уравнение, получим:

 

 

 

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

• найти общее решение соответствующего однородного уравнения (записать характеристическое уравнение, найти все корни характеристического уравнения l₁, l₂, ... , l_n, записать фундаментальную систему решений y₁(x), y₂(x), ..., y_n(x));
• найти любое частное решение неоднородного уравнения y_ч(x);
• записать выражение для общего решения

y(x)= c₁ y₁(x) + c₂ y₂(x) + ... + c_n y_n(x) + y_ч(x);

Для решения задачи Коши нужно подставить выражение для общего решения в начальные условия и определить значения постоянных c₁,..., c_n, которые являются решениями системы линейных алгебраических уравнений

c₁ y₁(x₀) + c₂ y₂(x₀) + ... + c_n y_n(x₀) + y_ч(x₀)= y₀,

c₁ y'₁(x₀) + c₂ y'₂(x₀) + ... + c_n y'_n(x₀) + y_ч(x₀)=y_0,1,

......... ,

c₁ y₁^(n-1)(x₀) + c₂ y₂^(n-1)(x₀) + ... + c_n y_n^(n-1)(x₀) + y_ч(x₀)= y_0,n-1

Пример 9.

Решить  задачу Коши 

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

общее решение  соответствующего однородного уравнения,   частное решение.

Решим соответствующее однородное уравнение 

Характеристическое  уравнение 

 

 

Подставим в исходное уравнение, получим:

 

 

Найдем  и ,  используя начальные условия:

 

МЕТОД ВАРИАЦИИ ПРОИЗВОЛЬНЫХ ПОСТОЯННЫХ РЕШЕНИЯ ЗАДАЧИ КОШИ ДЛЯ  ЛИНЕЙНОГО НЕОДНОРОДНОГО ДИФФЕРЕНЦИАЛЬНОГО  УРАВНЕНИЯ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ

Доказано, что для линейного неоднородного дифференциального уравнения  
y⁽ⁿ⁾ + a₁ y^(n-1) + ... + a_n-1 y' + a_n y = f(x) 
при непрерывной правой части f(x), для любых начальных значений  
x₀,  y₀,   y_0,1, ..., y_0,n-1  
существует и единственно решение задачи Коши  
y(x₀)=y₀, (y)'(x₀)=y_0,1 , ...,(y)^(n-1)(x₀)=y_0,n-1.

Решение задачи Коши для  неоднородного дифференциального  уравнения с постоянными коэффициентами можно найти методом вариации произвольных постоянных (методом Лагранжа), который состоит в следующем:  
Записываем искомое решение задачи Коши для неоднородного уравнения  в виде  
y(x)= c₁(x) y₁(x) +  c₂(x) y₂(x) + ... + c_n(x) y_n(x), 
где  y₁(x),  y₂(x), ..., y_n(x)  — линейно независимые решения соответствующего однородного уравнения, и находим неизвестные функции 
c₁(x) ,  c₂(x), ...,  c_n(x), 
такие, чтобы функция y = y(x) удовлетворяла неоднородному уравнению и заданным начальным условиям.

Опишем алгоритм решения задачи Коши для уравнения  второго порядка 
y'' + a₁ y' + a₂ y = f(x), y(x₀)=y₀, (y)'(x₀)=y_0,1.  
Будем искать решение задачи в виде 
y(x)= c₁(x) y₁(x) +  c₂(x) y₂(x),  
где  y₁(x),  y₂(x) — линейно независимые решения однородного уравнения 
y'' + a₁ y' + a₂ y = 0. 
Вычислим  y'(x), y''(x)  и подставим полученные выражения в уравнение.  
Вычислим первую производную 
y'(x)= (c₁'(x) y₁(x) +  c₂(x)' y₂(x)) + (c₁(x) y₁'(x) +  c₂(x) y₂'(x)),  
положим  
c₁'(x) y₁(x) +  c₂(x)' y₂(x)  = 0 
и тогда  
y'(x)= c₁(x) y₁'(x) +  c₂(x) y₂'(x), 
y''(x)= (y'(x))'= (c₁(x) y₁'(x) +  c₂(x) y₂'(x))'= 
=c₁'(x) y₁'(x) +  c₂(x)' y₂'(x) + c₁(x) y₁''(x) +  c₂(x) y₂''(x). 

 

Подставив y(x) и ее производные в уравнение, получим: 
y'' + a₁ y' + a₂ y = 
= c₁'(x) y₁'(x) +  c₂(x)' y₂'(x) + c₁(x) y₁''(x) +  c₂(x) y₂''(x) +  
+ a₁(c₁(x) y₁'(x)+c₂(x) y₂'(x)) + a₂(c₁(x) y₁(x)+c₂(x) y₂(x)) = 
= c₁(x)( y₁''(x)+a₁ y₁'(x)+a₂ y₁(x)) + c₂(x)( y₂''(x)+a₁ y₂'(x)+a₂ y₂(x)) + 
+ c₁'(x) y₁'(x) +  c₂(x)' y₂'(x) = 0 + 0 + c₁'(x) y₁'(x) +  c₂(x)' y₂'(x) = f(x), 
при условии c₁'(x) y₁(x) +  c₂(x)' y₂(x)  = 0.