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

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

a ₀(x) y ¢ + a ₁(x) y = B (x). (1)

При a ₀ ¹ 0 его можно представить в виде:

y ¢ + a (x) y = b (x), (2)

где a (x) = a ₁(x)/ a ₀(x) и b (x) = B (x)/ a ₀(x).

Если правые части (1) и (2) равны нулю, то эти уравнения называются однородными, в противном случае – неоднородными.

Если в уравнении (1) a ₀(x) = a ₀ и a ₁(x) = a ₁, то есть эти функции являются константами, то уравнение (1) называется линейным дифференциальным уравнением первого порядка с постоянными коэффициентами.

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

y ¢ + ay = 0. (3)

Перепишем его в виде: или . Последнюю формулу можно рассматривать как равенство дифференциалов функций одного и того же аргумента x. Интегрируя это равенство, получаем ln y = – ax + C, или y = e^{–ax + C}, где C ‑ произвольная константа. Если теперь ввести обозначение e^C = A, то можно представить так называемое общее решение уравнения (3) в виде:

y = A e^–ax. (4)

Это решение зависит от неопределенной константы A, придавая которой различные значения, можно получить все множество интегральных кривых уравнения (3). Если мы хотим найти интегральную кривую, проходящую через точку (x ₁, y ₁), то нужно подставить координаты точки в формулу (4) и определить значение константы A. С этим значением константы A формула (4) будет определять лишь одну интегральную кривую или так называемое частное решение уравнения (3).

Как правило, задача ставится так: найти решение уравнения (3) при условии

y (0) = y ₀. (5)

Последняя формула называется начальным условием для уравнения (3).

Дифференциальное уравнение (3) при начальном условии (5) имеет единственное решение, которое определяется формулой

y (x) = y ₀ e^–ax. (6)

Заметим, что для задания начального условия, вообще говоря, не обязательно выбирать значение аргумента x, равное нулю. Как сказано выше, выделить единственное решение из множества, задаваемого формулой (4) (то есть определить константу А), можно с помощью любого соотношения y (x ₁) = y ₁, считая его начальным условием.

Если в уравнении (3) a = 0, то интегрирование приводит к решению y (x) = C, то есть к константе, которая при начальном условии (5) равна y ₀. Таким образом решение y (x) сохраняет начальное значение y ₀ при изменении x.

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

y¢ + ay = b, (b = cost) (7)

с начальным условием y (0) = y ₀.

Введем новую неизвестную (считаем, что a ¹ 0).

Теперь уравнение (7) примет вид или z¢ + az = 0. Как было показано выше, решением последнего уравнения является функция z = z ₀ e^–ax, где . Возвращаясь к изначальной неизвестной, получаем решение уравнения (7) при заданном начальном условии:

. (8)

Если в уравнении (7) a = 0, то его решением при заданном начальном условии будет функция y (x) = bx + y ₀.

Заметим, что решение (8) состоит из двух частей: y_h = Ae^–ax ‑ решения однородного уравнения y¢ + ay = 0 и y ₀(x) = b / a ‑ решения, которое назовем равновесным и которое получается, если в уравнении (7) положить y¢ = 0. Такое представление позволяет рассматривать решение (8) уравнения (7) как сумму равновесного или фиксированного значения y_e и отклонения или девиации y_h тра­ектории y (x) от равновесного значения. Это отклонение возрастает экспоненциально с ростом x при a < 0 и стремится к нулю при a > 0. В первом случае (a < 0) решение называется неустойчивым, а во втором – устойчивым (асимптотически устойчивым).

Как показано на рисунках 1 и 2, отклонение y_h = (y ₀ – y_e) e^–ax от уровня равновесия уменьшается с ростом x при a > 0 и увеличивается с ростом x при a < 0.

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

у ¢ + a (x) y = b (x). (9)

Здесь a (x) ‑ некоторая функция аргумента x. Как мы это делали раньше, вначале будем искать решение однородного уравнения, положив функцию b (x) в правой части (9) равной нулю. Представив уравнение у ¢ + a (x) y = 0 в виде

,

после интегрирования получаем

или

. (10)

Здесь A ‑ неопределенная константа, которую можно найти из начального условия y (0) = 0.

Пример 4.

Решить уравнение y’ + 2 xy = 0 при начальном условии y (0) = 3.

Решение.

В этом случае

a (x) = 2 x,

и начальное условие определяет A = 3. Искомое решение имеет вид

.

Перейдем к решению неоднородного линейного дифференциального уравнения первого порядка с переменными коэффициентами. Положим в формуле (10) A = A (x), то есть будем считать множитель A некоторой функцией от x. Этот метод называется методом вариации произвольной постоянной, и с его помощью мы попытаемся решить уравнение (9) при условии, что b (x) есть некоторая функция, не равная тождественно нулю. Из формулы (10) получаем:

; .

После подстановки этих выражений уравнение (9) принимает вид

,

откуда следует уравнение относительно функции A(x):

,

с решением

.

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

. (11)

Пример 5.

Решить уравнение при начальном условии y (1) = 2.

(Заметим, что в данном случае нельзя задавать начальное условие при x = 0, так как это значение не принадлежит области B определения функции F).

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

В нашем уравнении . Решение однородного уравнения получается из формулы (10):

. (12)

Реализуем теперь вариацию произвольной константы A, считая, что A = A (x) есть некоторая функция аргумента x. Тогда , и подставив это выражение вместе с приведенным выше выражением для y в исходное уравнение, получим:

,

откуда следует, что A ¢(x) = x ² или . Если теперь подставить это в формулу (12), то получится общее решение исходного уравнения: . С помощью начального условия найдем значение неопределенной константы C и выпишем решение поставленной задачи: .

Пример 6.

Решение предыдущего примера 5. с использованием MATHCAD

Найти частное решение дифференциального уравнения, удовлетворяющего начальным условиям

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

Разделим переменные

Проинтегрируем последнее выражение

Приравнивая полученные результаты интегрирования, и разрешая уравнение относительно y, получим

.

2) Используем метод вариации для нахождения постоянной С:

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

,

считая С – переменной величиной.

Дифференцируем y(x)

Подставляем y(x) и в исходное дифференциальное уравнение

После сокращения, получаем дифференциальное уравнение относительно С(х)

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

Интегрируем последнее выражение

Записываем общий вид частного решения

.

Тогда, общее решение дифференциального неоднородного уравнения примет вид:

.

Подставляя начальные условия, находим значение постоянной величины С₁.

.

Следовательно, частное решение уравнения при заданных начальных условиях имеет вид

.

Поиск по сайту: