Материал: Diff_lektsii__1

удовлетворяющая системе

(эта система соответствует системе (3.32)), будет обладать следующим свойством: вектор-функция

Z(x) = Φ(x)α(x)

есть решение системы (3.15). Поэтому если для вектора α(x) с компонентами α1(x), ..., αn(x) выполнено равенство (3.38), то функция (3.37) является решением уравнения (3.14). Выпишем равенство (3.38) покомпонентно:

В случае отыскания решения неоднородного линейного дифференциального уравнения (3.14) обычно используют именно систему (3.39) для

функций α01(x), ..., α0n(x).

Пример. Найти общее решение уравнения

(0, +∞). Будем рассматривать уравнение (3.40) на (0, +∞). Соответствующее однородное уравнение

y00 − 2y0 + y = 0

– это уравнение с постоянными коэффициентами. Его характеристический полином λ2 −2λ+ 1 имеет корень λ = 1 кратности 2, ф.с.р. состоит из двух функций y1(x) = ex, y2(x) = xex. Система (3.39) запишется в виде

α01(x)ex + α02(x)xex = 0,α01(x)ex + α02(x)(x + 1)ex = exx .

50

Отсюда

α01(x) = −1, α02(x) = x1 .

Возьмем первообразные

α1(x) = −x, α2(x) = ln(x).

Тогда (3.37) имеет вид

ye(x) = −xex + (ln(x))xex.

Общее решение уравнения (3.40):

y = c1ex + c2xex − xex + (ln(x))xex.

Переобозначая c1 − 1 через c1 (понятно, что c1 пробегает множество всех вещественных чисел тогда и только тогда, когда c1 − 1 пробегает то же множество), можем переписать общее решение так:

y= ex(c1 + c2x + x ln(x)).

Вслучае, когда функции p1(x), ..., pn(x) в уравнении (3.14) постоянные, а q(x) имеет специальный вид, существует еще один метод нахождения решения уравнения (3.14).

Сформулируем (без доказательства) теорему, описывающую этот ме-

тод.

где Pk(x), Ql(x) – полиномы от x степеней k, l. Введем числа m = max(k, l),

α = a + bi,

где Pem(x), Qem(x) – полиномы от x степени m.

51

Теорему 3.7 применяют так: в выражении (3.43) для ye(x) неизвестными являются коэффициенты полиномов Pem(x), Qem(x). Запишем выражение вида (3.43) с неопределенными коэффициентами этих полиномов, подставим в (3.41) и приравняем коэффициенты при функциях вида xreax cos(bx), xreax sin(bx), r = 0, 1, ..., s + m. Теорема 3.7 гарантирует, что полученная алгебраическая система относительно коэффициентов полиномов Pem(x), Qem(x) имеет решение, и притом единственное.

Отметим, что в выражении (3.42) представление функции q(x) в указанном виде может быть не единственным. Например, функцию

q(x) = xex

можно представить как

q(x) = ex(x cos(0 · x) + 1 sin(0 · x)),

т. е. взять P1(x) = x, Q0(x) = 1. В этом случае m = max(1, 0) = 1. Эту же функцию можно представить как

q(x) = ex(x cos(0 · x) + (x2 − 1) sin(0 · x)),

т. е. взять P1(x) = 1, Q2(x) = x2 − 1. В этом случае m = max(1, 2) = 2. На практике следует выбирать представление (3.42) так, чтобы число

m было возможно наименьшим.

Это уравнение вида (3.41). Функцию q(x) = x в правой части можно представить в виде (3.42):

x = e0·x(x cos(0 · x) + 1 sin(0 · x))

с a = 0, b = 0, Pk(x) = x, k = 1, Ql(x) = 1, l = 0. Число α = a + bi = 0

– не корень полинома λ2 − 2λ + 1. Поэтому s = 0, m = max(1, 0) = 1. По теореме 3.7 у уравнения (3.44) есть единственное частное решение вида

где Pe1(x), Qe1(x) – полиномы степени 1, т. е.

ye(x) = Pe1(x) = Ax + B

(записываем неизвестный полином Pe1(x) степени 1 с неопределенными коэффициентами A, B). Дифференцируем ye(x):

ye0(x) = A, ye00(x) = 0

52

и подставляем в уравнение (3.44):

0 − 2A + Ax + B = x.

Приравниваем коэффициенты при x1, x0:

A = 1, −2A + B = 0.

Отсюда B = 2 и частное решение ye(x) имеет вид

ye(x) = x + 2.

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

y(x) = c1ex + c2xex.

Поэтому общее решение уравнения (3.44) запишется так:

y(x) = c1ex + c2xex + x + 2.

Сформулируем одно простое утверждение, связанное с нахождением решений неоднородных линейных систем и неоднородных линейных дифференциальных уравнений (так называемый принцип суперпозиции).

Пусть Z1(x) – решение системы

Тогда Z(x) = Z1(x) + Z2(x) – решение системы

Y0 = P (x)Y + Q1(x) + Q2(x).

Для доказательства этого утверждения подставим Z1(x) в (3.45), Z2(x)

– в (3.46) и сложим полученные выражения:

Z0 = Z01 + Z02 = P (x)Z1 + Q1(x) + P (x)Z2 + Q2(x) =

= P (x)(Z1 + Z2) + Q1(x) + Q2(x) = P (x)Z + Q1(x) + Q2(x).

Аналогичное утверждение для линейных дифференциальных уравнений порядка n формулируется так: если ye1(x) – решение уравнения

y(n) + p1(x)y(n−1) + ... + pn(x)y = q1(x),

ye2(x) – решение уравнения

y(n) + p1(x)y(n−1) + ... + pn(x)y = q2(x),

53

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

ye(x) = −xex + (ln(x))xex + x + 2

является решением уравнения (3.47), а общее решение этого уравнения име-

ет вид

y = ex(c1 + c2x + x ln(x)) + x + 2.

Как уже отмечалось, при решении линейных дифференциальных уравнений можно использовать операционный метод (иногда называемый методом преобразования Лапласа). Преобразование Лапласа применяется к так называемым функциям-оригиналам. Дадим вначале необходимые определения.

Рассмотрим функцию f(x) с комплексными значениями, определенную на промежутке a < x < b. Можно записать

f(x) = ϕ(x) + iψ(x),

где ϕ и ψ – вещественнозначные функции. Если существуют интегралы

b b

ZZ

ϕ(x) dx и ψ(x) dx, полагаем

Для функций с комплексными значениями верны аналоги всех основных теорем интегрального исчисления: формула Ньютона–Лейбница, формулы замены переменных и интегрирования по частям, а также неравенство

ZZ

f(x) dx ≤ |f(x)| dx.

(Здесь предполагается, что b ≥ a.)

54