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

Материал из Википедии — свободной энциклопедии

Перейти к: навигация, поиск

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

\sum^{n}_{k=0} {a_k y^{(k)}(x)}=a_n y^{(n)} + a_{n-1} y^{(n-1)}+\dots+ a_1 y' + a_0 y=f(t)

где

  • y = y(t) — искомая функция,
  • y(k) = y(k)(t) — её k-тая производная,
  • a_0, a_1, a_2,\dots a_n — фиксированые числа,
  • f(t) — заданая функция (когда f(t)\equiv 0, имеем линейное однородное уравнение, иначе — линейное неоднородное уравнение).

Содержание

[править] Однородное уравнение

[править] Уравнение порядка n

Однородное уравнение:

a_n y^{(n)} + a_{n-1} y^{(n-1)} + \dots + a_1 y' + a_0 y=0

интегрируется следующим образом:

Пусть \lambda_1, \dots,\lambda_k — все различные корни характеристического многочлена, являющегося левой частью характеристического уравнения

a_n \lambda^n+a_{n-1}\lambda^{n-1}+\dots+a_1 \lambda + a_0 = 0

кратностей m_1,m_2,\dots, m_k, соответственно, m_1+m_2+\dots+m_k=n.

Тогда функции

t^\nu e^{\lambda_jt},\  \ 1\le j\le k,\ \  0\le \nu \le m_j-1

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

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

Воспользовавшись формулой Эйлера для пар комплексно сопряженных корней \lambda_j = \alpha_j \pm i\beta_j,\  \ 1\le j\le k можно заменить соответствующие пары комплексных функций в фундаментальной системе решений парами вещественных функций вида

t^\nu e^{\alpha_j t}\cos (\beta_j t),\  \ t^\nu e^{\alpha_j t}\sin (\beta_j t),\  \ j\in \overline{1 \dots k},\ \  0\le \nu \le m_j-1

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

[править] Уравнение второго порядка

Однородное уравнение второго порядка:

a2y'' + a1y' + a0y = 0

интегрируется следующим образом:

Пусть λ12 — корни характеристического уравнения.

a2λ2 + a1λ + a0 = 0,

являющегося квадратным уравнением.

Вид общего решения однородного уравнения зависит от значения дискриминанта \Delta=a_1^2 - 4a_2a_0:

  • при Δ > 0 уравнение имеет два различных вещественных корня
\lambda_{1,2} =\alpha_{1,2} = \frac{-a_1 \pm \sqrt{\Delta}}{2a_2}.

Общее решение имеет вид:

y(t) = c_1e^{\alpha_1 t} + c_2e^{\alpha_2 t}
  • при Δ = 0 — двa совпадающих вещественных корня
\lambda_1 = \lambda_2 = \alpha = \frac{-a_1}{2a_2}.

Общее решение имеет вид:

y(t) = c1eαt + c2teαt
\lambda_{1,2} =\alpha \pm i\beta = \frac{-a_1}{2a_2} \pm i\frac{ \sqrt{|\Delta|}}{2a_2}.

Общее решение имеет вид:

y(t) = c1eαtcos(βt) + c2eαtsin(βt)

[править] Неоднородное уравнение

Неоднородное уравнение интегрируется методом вариации произвольных постоянных (Метод Лагранжа).

[править] Вид общего решения неоднородного уравнения

Если дано частное решение неоднородного уравнения y0(t), и y_1(t),\ldots, y_n(t) — фундаментальная система решений соответствующего однородного уравнения, то общее решение уравнения задается формулой

y(t) = c_1y_1(t) + \ldots + c_ny_n(t) + y_0(t),

где c_1,\dots,c_n — произвольные постоянные.

[править] Принцип суперпозиции

Как в общем случае линейных уравнений, имеет место принцип суперпозиции, используемый в разных формулировках принципа суперпозиции в физике.

В случае, когда функция в правой части состоит из суммы двух функций

f(t) = f1(t) + f2(t),

частное решение неоднородного уравнения тоже состоит из суммы двух функций

y0(t) = y01(t) + y02(t),

где y_{0j}(t),\ \  j\in \overline{1,2} являются решениями неоднородного уравнения с правыми частями f_j(t),\ \  j\in \overline{1,2}, соответственно.

[править] Частный случай: квазимногочлен

В случае, когда f(t) — квазимногочлен, то есть

f(t) = p(t)eαtcos(βt) + q(t)eαtsin(βt)

где p(t),\ q(t) — многочлены, частное решение уравнения ищется в виде

y0(t) = (P(t)eαtcos(βt) + Q(t)eαtsin(βt))ts

где

  • P(t),\ Q(t) многочлены, deg(P)=deg(Q)=Max(deg(p),\ deg(q)), коэффициенты которых находятся подстановкой y0(t) в уравнение и вычисление методом неопределенных коэффициентов.
  • s является кратностью комплексного числа w = α + iβ, как корня характеристического уравнения однородного уравнения.

В частности, когда

f(t) = p(t)eαt

где p(t) — многочлен, частное решение уравнения ищется в виде

y0(t) = P(t)eαtts

Здесь P(t) — многочлен, deg(P) = deg(p), с неопределенными коэффициентами, которые находятся подстановкой y0(t) в уравнение. s является кратностью α, как корня характеристического уравнения однородного уравнения.

Когда же

f(t) = p(t)

где p(t) — многочлен, частное решение уравнения ищется в виде

y0(t) = P(t)ts

Здесь P(t) — многочлен, deg(P) = deg(p), а s является кратностью нуля, как корня характеристического уравнения однородного уравнения.

[править] Уравнение Коши — Эйлера

Уравнение Коши — Эйлера является частным случаем линейного дифференциального уравнения вида:


\sum^{n}_{k=1} {a_k(\alpha x + \beta )^k y^{(k)}(x)}=
a_n(\alpha x + \beta )^n y^{(n)}(x) + ... +
a_2(\alpha x + \beta )^2 y''(x) + a_1(\alpha x + \beta ) y'(x) + a_0y(x) = f(x)
,

приводимым к линейному дифференциальному уравнению с постоянными коэффициентами подстановкой вида x + β) = et.