Решить дифференциальное уравнение методом вариации произвольных постоянных

Метод вариации произвольной постоянной решения линейных неоднородных уравнений

Пример №1 . Найдём общее решение уравнения y» + 4y’ + 3y = 9e -3 x . Рассмотрим соответствующее однородное уравнение y» + 4y’ + 3y = 0. Корни его характеристического уравнения r 2 + 4r + 3 = 0 равны -1 и -3. Поэтому фундаментальная система решений однородного уравнения состоит из функций y₁ = e — x и y₂ = e -3 x . Решение неоднородного уравнения ищем в виде y = C₁(x)e — x + C₂(x)e -3 x . Для нахождения производных C’₁, C’₂ составляем систему уравнений (8)
C′₁·e -x +C′₂·e -3x =0
-C′₁·e -x -3C′₂·e -3x =9e -3x
решая которую, находим , Интегрируя полученные функции, имеем
Окончательно получим

Пример №2 . Решить линейные дифференциальные уравнения второго порядка с постоянными коэффициентами методом вариации произвольных постоянных:

y(0) =1 + 3ln3
y’(0) = 10ln3

Решение:
Данное дифференциальное уравнение относится к линейным дифференциальным уравнениям с постоянными коэффициентами.
Решение уравнения будем искать в виде y = e rx . Для этого составляем характеристическое уравнение линейного однородного дифференциального уравнения с постоянными коэффициентами:
r 2 -6 r + 8 = 0
D = (-6) 2 — 4·1·8 = 4

Корни характеристического уравнения: r₁ = 4, r₂ = 2
Следовательно, фундаментальную систему решений составляют функции: y₁=e 4x , y₂=e 2x
Общее решение однородного уравнения имеет вид: y =C₁·e 4x +C₂·e 2x
Поиск частного решения методом вариации произвольной постоянной.
Для нахождения производных C’_i составляем систему уравнений:
C′₁·e 4x +C′₂·e 2x =0
C′₁(4e 4x ) + C′₂(2e 2x ) = 4/(2+e -2x )
Выразим C’₁ из первого уравнения:
C’₁ = -c₂e -2x
и подставим во второе. В итоге получаем:
C’₁ = 2/(e 2x +2e 4x )
C’₂ = -2e 2x /(e 2x +2e 4x )
Интегрируем полученные функции C’_i:
C₁ = 2ln(e -2x +2) — e -2x + C * ₁
C₂ = ln(2e 2x +1) – 2x+ C * ₂

Поскольку y =C₁·e 4x +C₂·e 2x , то записываем полученные выражения в виде:
C₁ = (2ln(e -2x +2) — e -2x + C * ₁) e 4x = 2 e 4x ln(e -2x +2) — e 2x + C * ₁ e 4x
C₂ = (ln(2e 2x +1) – 2x+ C * ₂)e 2x = e 2x ln(2e 2x +1) – 2x e 2x + C * ₂ e 2x
Таким образом, общее решение дифференциального уравнения имеет вид:
y = 2 e 4x ln(e -2x +2) — e 2x + C * ₁ e 4x + e 2x ln(2e 2x +1) – 2x e 2x + C * ₂ e 2x
или
y = 2 e 4x ln(e -2x +2) — e 2x + e 2x ln(2e 2x +1) – 2x e 2x + C * ₁ e 4x + C * ₂ e 2x

Найдем частное решение при условии:
y(0) =1 + 3ln3
y’(0) = 10ln3

Подставляя x = 0, в найденное уравнение, получим:
y(0) = 2 ln(3) — 1 + ln(3) + C * ₁ + C * ₂ = 3 ln(3) — 1 + C * ₁ + C * ₂ = 1 + 3ln3
Находим первую производную от полученного общего решения:
y’ = 2e 2x (2C₁ e 2x + C₂ -2x +4 e 2x ln(e -2x +2)+ ln(2e 2x +1)-2)
Подставляя x = 0, получим:
y’(0) = 2(2C₁ + C₂ +4 ln(3)+ ln(3)-2) = 4C₁ + 2C₂ +10 ln(3) -4 = 10ln3

Получаем систему из двух уравнений:
3 ln(3) — 1 + C * ₁ + C * ₂ = 1 + 3ln3
4C₁ + 2C₂ +10 ln(3) -4 = 10ln3
или
C * ₁ + C * ₂ = 2
4C₁ + 2C₂ = 4
или
C * ₁ + C * ₂ = 2
2C₁ + C₂ = 2
Откуда: C₁ = 0, C * ₂ = 2
Частное решение запишется как:
y = 2e 4x ·ln(e -2x +2) — e 2x + e 2x ·ln(2e 2x +1) – 2x·e 2x + 2·e 2x

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

Здесь мы применим метод вариации постоянных Лагранжа для решения линейных неоднородных дифференциальных уравнений второго порядка. Подробное описание этого метода для решения уравнений произвольного порядка изложено на странице
Решение линейных неоднородных дифференциальных уравнений высших порядков методом Лагранжа >>> .

Пример 1

Решить дифференциальное уравнение второго порядка с постоянными коэффициентами методом вариации постоянных Лагранжа:
(1)

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

Вначале мы решаем однородное дифференциальное уравнение:
(2)
Ищем решение в виде . Составляем характеристическое уравнение:

Это уравнение второго порядка.

Решаем квадратное уравнение:
.
Корни кратные: . Фундаментальная система решений уравнения (2) имеет вид:
(3) .
Отсюда получаем общее решение однородного уравнения (2):
(4) .

Шаг 2. Вариация постоянных – замена постоянных функциями

Варьируем постоянные C ₁ и C ₂ . То есть заменим в (4) постоянные и на функции:
.
Ищем решение исходного уравнения (1) в виде:
(5) .

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

.
Поскольку и удовлетворяют однородному уравнению (2), то сумма членов в каждом столбце последних трех строк дает нуль и предыдущее уравнение приобретает вид:
(7) .
Здесь .

Вместе с уравнением (6) мы получаем систему уравнений для определения функций и :
(6) :
(7) .

Решение системы уравнений

Решаем систему уравнений (6-7). Выпишем выражения для функций и :
.
Находим их производные:
;
.

Решаем систему уравнений (6-7) методом Крамера. Вычисляем определитель матрицы системы:

.
По формулам Крамера находим:
;
.

Итак, мы нашли производные функций:
;
.
Интегрируем (см. Методы интегрирования корней). Делаем подстановку
; ; ; .

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

;
.

Пример 2

Решить дифференциальное уравнение методом вариации постоянных Лагранжа:
(8)

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

Решаем однородное дифференциальное уравнение:

(9)
Ищем решение в виде . Составляем характеристическое уравнение:

Это уравнение имеет комплексные корни:
.
Фундаментальная система решений, соответствующая этим корням, имеет вид:
(10) .
Общее решение однородного уравнения (9):
(11) .

Шаг 2. Вариация постоянных – замена постоянных функциями

Теперь варьируем постоянные C ₁ и C ₂ . То есть заменим в (11) постоянные на функции:
.
Ищем решение исходного уравнения (8) в виде:
(12) .

Далее ход решения получается таким же, как в примере 1. Мы приходим к следующей системе уравнений для определения функций и :
(13) :
(14) .
Здесь .

Решение системы уравнений

Решаем эту систему. Выпишем выражения функций и :
.
Из таблицы производных находим:
;
.

Решаем систему уравнений (13-14) методом Крамера. Определитель матрицы системы:

.
По формулам Крамера находим:
;
.

Первый интеграл немного сложней (см. Интегрирование тригонометрических рациональных функций). Делаем подстановку :

.
Поскольку , то знак модуля под знаком логарифма можно опустить. Умножим числитель и знаменатель на :
.
Тогда
.

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

.

Автор: Олег Одинцов . Опубликовано: 05-08-2013 Изменено: 19-06-2017

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

В данной теме поговорим о способах решения линейных неоднородных дифференциальных уравнений вида 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 )

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