Решение дифференциальных уравнений 1 порядка методом лагранжа

Метод Лагранжа (вариации постоянной). Линейные дифференциальные уравнения первого порядка.

Рассмотрим линейное неоднородное дифференциальное уравнение первого порядка:
(1) .
Существует три способа решения этого уравнения:

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

Метод вариации постоянной (Лагранжа)

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

Шаг 1 Решение однородного уравнения

Ищем решение однородного уравнения:

Это уравнение с разделяющимися переменными

Разделяем переменные — умножаем на dx , делим на y :

Интегрируем:

Интеграл по y — табличный:

Тогда

Потенцируем:

Заменим постоянную e C на C и уберем знак модуля, что сводится к умножению на постоянную ±1 , которую включим в C :

Шаг 2 Заменим постоянную C на функцию

Теперь заменим постоянную C на функцию от x :
C → u ( x )
То есть, будем искать решение исходного уравнения (1) в виде:
(2)
Находим производную.

По правилу дифференцирования сложной функции:
.
По правилу дифференцирования произведения:

.
Подставляем в исходное уравнение (1):
(1) ;

.
Два члена сокращаются:
;
.
Интегрируем:
.
Подставляем в (2):
.
В результате получаем общее решение линейного дифференциального уравнения первого порядка:
.

Пример решения линейного дифференциального уравнения первого порядка методом Лагранжа

Решаем однородное уравнение:

Разделяем переменные:

Умножим на :

Интегрируем:

Интегралы табличные:

Потенцируем:

Заменим постоянную e C на C и убираем знаки модуля:

Отсюда:

Заменим постоянную C на функцию от x :
C → u ( x )

Находим производную:
.
Подставляем в исходное уравнение:
;
;
Или:
;
.
Интегрируем:
;
Решение уравнения:
.

Общее решение уравнения:
.

Автор: Олег Одинцов . Опубликовано: 27-07-2012 Изменено: 01-03-2015

ЛДУ с переменными коэффициентами. Метод Лагранжа

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

Если известно частное решение уравнения

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

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

Общее решение уравнения

есть сумма какого-нибудь его частного решения и общего решения соответствующего однородного уравнения (32).

Если известна фундаментальная система соответствующего однородного уравнения (32), то общее решение неоднородного уравнения (33) может быть найдено методом вариации постоянных ( метод Лагранжа ).

Общее решение уравнения (32) имеет вид

где — произвольные постоянные.

Будем искать решение уравнения (33) в виде

где — некоторые пока неизвестные функции от . Для их определения получаем систему

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

где — произвольные постоянные. Внося найденные значения в (34), получаем общее решения уравнения (33).

В частности, для уравнения второго порядка

Решая (36) относительно и , получаем

где и — постоянные интегрирования.

Замечание. Для уравнения , где , система (36) будет выглядеть так:

Пример 1. Найти общее решение уравнения , если есть его частное решение.

Решение. Положим , где — новая неизвестная функция от , тогда

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

Но так как есть частное решение данного уравнения, то , поэтому имеем

Но , а значит , и уравнение (37) примет вид

Перепишем его в виде . Отсюда имеем , откуда

Интегрируя это уравнение, найдем и, следовательно, общее решение данного уравнения будет

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

Решение. Общее решение соответствующего однородного уравнения имеет вид (см. пример 1)

и следовательно, его фундаментальная система решений будет

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

где — постоянные неизвестные функции от , подлежащие определению. Для их нахождения составим следующую систему:

Отсюда находим: . Интегрируя, получаем

Подставляя эти значения и в выражение для , найдем общее решение данного уравнения

Пример 3. Решить уравнение .

Решение. Соответствующее однородное уравнение будет . Его характеристическое уравнение имеет мнимые корни , и общее решение однородного уравнения имеет вид

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

где и — неизвестные функции от . Для их нахождения составим систему

Разрешаем эту систему относительно и :

Подставляя выражения и в (38), получаем общее решение данного уравнения

Здесь есть частное решение исходного неоднородного уравнения.

Пример 4. Зная фундаментальную систему решений соответствующего однородного уравнения, найти частное решение уравнения

Решение. Применяя метод вариации постоянных, находим общее решение уравнения (39):

При первые два слагаемых правой части (40) стремятся к бесконечности, причем при любых , неравных нулю одновременно, функция есть бесконечно большая функция при . Третье слагаемое правой части (40) имеет пределом ноль при , что легко установить с помощью правила Лопиталя. Таким образом, функция , которая получается из (40) при и , будет решением уравнения (39), удовлетворяющим условию .

Составление дифференциального уравнения по заданной фундаментальной системе решений

Рассмотрим линейно независимую на отрезке систему функций

имеющих все производные до n-го порядка включительно. Тогда уравнение

где — неизвестная функция, будет линейным дифференциальным уравнением, для которого, как нетрудно видеть, функции составляют фундаментальную систему решений. Коэффициент при в (42) есть определитель Вронского системы (41). Те точки, в которых этот определитель обращается в ноль, будут особыми точками построенного уравнения — в этих точках обращается в ноль коэффициент при старшей производной .

Пример 1. Составить дифференциальное уравнение, для которого образуют фундаментальную систему решений.

Решение. Применяя формулу (42), получаем

Раскрывая определитель в левой части (43) по элементам третьего столбца, будем иметь . Это и есть искомое дифференциальное уравнение.

Пример 2. Составить дифференциальное уравнение, для которого функции фундаментальную систему решений образуют функции .

Решение. Составим уравнение вида (42):

Раскрывая последний определитель по элементам 3-го столбца, будем иметь

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

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

видим, что коэффициент при терпит разрыв при , так что в точке непрерывность коэффициентов уравнения (45) нарушается.

Разные задачи

Пусть — фундаментальная система линейного однородного уравнения

Тогда имеет место формула Остроградского–Лиувилля

где — определитель Вронского, а — любое значение из отрезка , на котором непрерывны коэффициенты уравнения.

Пример 1. Показать, что линейное дифференциальное уравнение имеет решение вида , где — некоторый многочлен. Показать, что второе решение этого уравнения имеет вид , где — также многочлен.

Решение. Будем искать решение в виде многочлена, например, первой степени: . Подставляя в уравнение, найдем, что . Пусть , тогда ;. таким образом, многочлен будет решением данного уравнения. Перепишем данное уравнение в виде

Пусть — второе частное решение данного уравнения, линейно независимое с первым. Находим определитель Вронского системы решений

здесь . Применяя формулу Остроградского–Лиувилля, будем иметь

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

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

В данной теме поговорим о способах решения линейных неоднородных дифференциальных уравнений вида y ‘ = P ( x ) · y = Q ( x ) . Начнем с метода вариации произвольной постоянной и покажем способ применения этого метода для решения задачи Коши. Продолжим рассмотрением метода, который предполагает представление произвольной постоянной у как произведения двух функций u ( x ) и v ( x ) . В разделе мы приводим большое количество задач по теме с детальным разбором решения.

На тот случай, если применяемые при разборе темы термины и понятия окажутся незнакомыми для вас, мы рекомендуем заглядывать в раздел «Основные термины и определения теории дифференциальных уравнений».

Метод вариации произвольной постоянной для решения ЛНДУ первого порядка

Для краткости будет обозначать линейное неоднородное дифференциальное уравнение аббревиатурой ЛНДУ, а линейное однородное дифференциальное уравнение (ЛОДУ).

ЛНДУ вида y ‘ = P ( x ) · y = Q ( x ) соответствует ЛОДУ вида y ‘ = P ( x ) · y = 0 , при Q ( x ) = 0 . Если посмотреть на дифференциальное уравнение y ‘ = P ( x ) · y = 0 , становится понятно, что мы имеем дело с уравнением с разделяющимися переменными. Мы можем его проинтегрировать: y ‘ = P ( x ) · y = 0 ⇔ d y y = — P ( x ) d x , y ≠ 0 ∫ d y y = — ∫ P ( x ) d x ⇔ ln y + C 1 = — ∫ P ( x ) d x ⇔ ln y = ln C — ∫ P ( x ) d x , ln C = — C 1 , C ≠ 0 ⇔ e ln y = e ln C — ∫ P ( x ) d x ⇔ y = C · e — ∫ P ( x ) d x

Мы можем утверждать, что значение переменной y = 0 тоже является решением, так как при этом значении переменной уравнение y ‘ = P ( x ) · y = 0 обращается в тождество. Этому случаю соответствует решение y = C · e — ∫ P ( x ) d x при значении C = 0 .

Получается, что y = C · e — ∫ P ( x ) d x — общее решение ЛОДУ, где С – произвольная постоянная.

y = C · e — ∫ P ( x ) d x — это решение ЛОДУ y ‘ = P ( x ) · y = 0 .

Для того, чтобы найти общее решение неоднородного уравнения y ‘ = P ( x ) · y = Q ( x ) , будем считать С не константой, а функцией аргумента х . Фактически, мы примем y = C ( x ) · e — ∫ P ( x ) d x общим решением ЛНДУ.

Подставим y = C ( x ) · e — ∫ P ( x ) d x в дифференциальное уравнение y ‘ = P ( x ) · y = Q ( x ) . Оно при этом обращается в тождество:

y ‘ = P ( x ) · y = Q ( x ) C x · e — ∫ P ( x ) d x + P ( x ) · C ( x ) · e — ∫ P ( x ) d x = Q ( x )

Теперь обратимся к правилу дифференцирования произведения. Получаем:

C ‘ ( x ) · e — ∫ P ( x ) d x + C ( x ) · e — ∫ P ( x ) d x + P ( x ) · C ( x ) · e — ∫ P ( x ) d x = Q ( x )

Производная сложной функции e — ∫ P ( x ) d x ‘ равна e — ∫ P ( x ) d x · — ∫ P ( x ) d x ‘ .

Теперь вспомним свойства неопределенного интеграла. Получаем:

e — ∫ P ( x ) d x · — ∫ P ( x ) d x ‘ = — e — ∫ P ( x ) d x · P ( x )

Теперь выполним переход:

C ‘ ( x ) · e — ∫ P ( x ) d x + C ( x ) · e — ∫ P ( x ) d x ‘ + P ( x ) · C ( x ) · e — ∫ P ( x ) d x = Q ( x ) C ‘ ( x ) · e — ∫ P ( x ) d x — P ( x ) · C ( x ) · e — ∫ P ( x ) d x + P ( x ) · C ( x ) · e — ∫ P ( x ) d x = Q ( x ) C ‘ ( x ) · e — ∫ P ( x ) d x = Q ( x )

Так мы пришли к простейшему дифференциальному уравнению первого порядка. В ходе решения этого уравнения мы определим функцию C ( x ) . Это позволит нам записать решение исходного ЛНДУ первого порядка следующим образом:

y = C ( x ) · e — ∫ P ( x ) d x

Подведем итог

Метод вариации произвольной постоянной при решении ЛНДУ предполагает проведение трех этапов:

• нахождение общего решения соответствующего ЛОДУ y ‘ + P ( x ) · y = 0 в виде y = C · e — ∫ P ( x ) d x ;
• варьирование произвольной постоянной С , что заключается в замене ее функцией С ( x ) ;
• подстановка функции y = C ( x ) · e — ∫ P ( x ) d x в исходное дифференциальное уравнение, откуда мы можем вычислить C ( x ) и записать ответ.

Теперь применим этот алгоритм к решению задачи.

Найдите решение задачи Коши y ‘ — 2 x y 1 + x 2 = 1 + x 2 , y ( 1 ) = 3 .

Нам нужно отыскать частное решение ЛНДУ y ‘ — 2 x y 1 + x 2 = 1 + x 2 при начальном условии y ( 1 ) = 3 .

В нашем примере P ( x ) = — 2 x 1 + x 2 и Q ( x ) = x 2 + 1 . Начнем с того, что найдем общее решение ЛОДУ. После этого применим метод вариации произвольной постоянной и определим общее решение ЛНДУ. Это позволит нам найти искомое частное решение.

Общим решением соответствующего ЛОДУ y ‘ — 2 x y 1 + x 2 = 0 будет семейство функций y = C · ( x 2 + 1 ) , где С – произвольная постоянная.

Варьируем произвольную постоянную y = C ( x ) · ( x 2 + 1 ) и подставляем эту функцию в исходное уравнение:
y ‘ — 2 x y 1 + x 2 = 1 + x 2 C x · ( x 2 + 1 ‘ — 2 x · C ( x ) · ( x 2 + 1 ) 1 + x 2 = 1 + x 2 C ‘ ( x ) · ( x 2 + 1 ) + C ( x ) · 2 x — 2 x · C ( x ) = 1 + x 2 C ‘ ( x ) = 1 ,

откуда C ( x ) = ∫ d x = x + C 1 , где C 1 – произвольная постоянная.

Это значит, что y = C ( x ) · ( x 2 + 1 ) = ( x + C 1 ) · ( x 2 + 1 ) — общее решение неоднородного уравнения.

Теперь приступим к отысканию частного решения, которое будет удовлетворять начальному условию y ( 1 ) = 3 .

Так как y = ( x + C 1 ) · ( x 2 + 1 ) , то y ( 1 ) = ( 1 + C 1 ) · ( 1 2 + 1 ) = 2 · ( 1 + C 1 ) . Обратившись к начальному условию, получаем уравнение 2 · ( 1 + C 1 ) = 3 , откуда C 1 = 1 2 . Следовательно, искомое решение задачи Коши имеет вид y = x + 1 2 · ( x 2 + 1 )

Теперь рассмотрим еще один метод решения линейных неоднородных дифференциальных уравнений y ‘ + P ( x ) · y = Q ( x ) .

Еще один метод решения ЛНДУ первого порядка

Мы можем представить неизвестную функцию как произведение y = u ⋅ v , где u и v – функции аргумента x .

Мы можем подставить эту функцию в ЛНДУ первого порядка. Имеем:

y ‘ + P ( x ) · y = Q ( x ) ( u · v ) ‘ + P ( x ) · u · v = Q ( x ) u ‘ · v + u · v ‘ + P ( x ) · u · v = Q ( x ) u ‘ · v + u · ( v ‘ + P ( x ) · v ) = Q ( x )

Если найти такое v , чтобы оно было ненулевым частным решением дифференциального уравнения v ‘ + P ( x ) · v = 0 , то u можно будет определить из уравнения с разделяющимися переменными u ‘ · v = Q ( x ) .

Рассмотрим этот алгоритм решения на предыдущем примере. Это позволит нам сосредоточиться на главном, не отвлекаясь на второстепенные детали.

Найдите общее решение линейного неоднородного дифференциального уравнения y ‘ — 2 x y 1 + x 2 = 1 + x 2 .

Пусть y = u ⋅ v , тогда
y ‘ — 2 x y x 2 + 1 = x 2 + 1 ⇔ ( u · v ) — 2 x · u · v x 2 + 1 = x 2 + 1 u ‘ · v + u · v ‘ — 2 x · u · v x 2 + 1 = x 2 + 1 u ‘ · v + u · v ‘ — 2 x · v x 2 + 1 = x 2 + 1

Находим такое v , отличное от нуля, чтобы выражение в скобках обращалось в ноль. Иными словами, находим частное решение дифференциального уравнения v ‘ — 2 x · v x 2 + 1 = 0 .
v ‘ — 2 x · v x 2 + 1 = 0 ⇔ d v d x = 2 x · v x 2 + 1 ⇒ d v v = 2 x d x x 2 + 1 ⇔ d v v = d ( x 2 + 1 ) x 2 + 1 ∫ d v v = ∫ d ( x 2 + 1 ) x 2 + 1 ln v + C 1 = ln ( x 2 + 1 ) + C 2

Возьмем частное решение v = x 2 + 1 , соответствующее C 2 – С 1 = 0 .

Для этого частного решения имеем
u ‘ · v + u · v ‘ — 2 x · v x 2 + 1 = x 2 + 1 ⇔ u ‘ · ( x 2 + 1 ) + u · 0 = x 2 + 1 ⇔ u ‘ = 1 ⇔ u = x + C

Следовательно, общее решение исходного линейного неоднородного дифференциального уравнения есть y = u · v = ( x + C ) · ( x 2 + 1 )

Ответы в обоих случаях совпадают. Это значит, что оба метода решения, которые мы привели в статье, равнозначны. Выбирать, какой из них применить для решения задачи, вам.