Обыкновенные дифференциальные уравнения

Обыкновенное дифференциальное уравнение (ОДУ) — это уравнение, содержащее одну независимую переменную, неизвестную функцию этой переменной и её производные различных порядков. Порядок уравнения определяется наивысшим порядком входящей в него производной. ОДУ являются основным математическим инструментом для описания эволюционных процессов в физике, химии, биологии, экономике и технике, где состояние системы изменяется во времени или пространстве.

История

Первые дифференциальные уравнения возникли в конце XVII века в работах Исаака Ньютона и Готфрида Вильгельма Лейбница в связи с задачами механики (законы движения планет, колебания маятника). Ньютон в «Математических началах натуральной философии» (1687) сформулировал дифференциальные уравнения движения, хотя и в геометрической форме. Лейбниц ввёл современные обозначения производных и интегралов.

В XVIII веке математики, такие как Леонард Эйлер, Даниил Бернулли и Жан Лерон Д’Аламбер, разработали методы решения многих типов ОДУ. Эйлер систематизировал теорию линейных уравнений с постоянными коэффициентами, ввёл понятие характеристического уравнения. В XIX веке Огюстен Луи Коши заложил основы теории существования и единственности решений, а Софус Ли разработал теорию непрерывных групп симметрий, позволившую классифицировать интегрируемые уравнения. К концу XIX века сформировалась качественная теория дифференциальных уравнений (Анри Пуанкаре, Александр Ляпунов), изучающая поведение решений без их явного нахождения.

Классификация

ОДУ классифицируются по нескольким признакам.

По порядку

Порядок уравнения равен наивысшей производной неизвестной функции. Например:

• Уравнение первого порядка: \( y' = f(x, y) \).
• Уравнение второго порядка: \( y'' + p(x)y' + q(x)y = r(x) \).

По линейности

• Линейные ОДУ: функция и её производные входят в уравнение линейно (не умножаются друг на друга, не возводятся в степень). Общий вид: \( a_n(x) y^{(n)} + a_{n-1}(x) y^{(n-1)} + \dots + a_0(x) y = f(x) \). Если \( f(x) = 0 \), уравнение называется однородным, иначе — неоднородным.
• Нелинейные ОДУ: содержат нелинейные комбинации функции и её производных (например, \( y y' = 1 \), \( (y')^2 + y = 0 \)).

По типу коэффициентов

• С постоянными коэффициентами (например, \( y'' + 3y' + 2y = 0 \)).
• С переменными коэффициентами (например, уравнение Эйлера: \( x^2 y'' + a x y' + b y = 0 \)).

Автономные и неавтономные

• Автономные: правая часть явно не зависит от независимой переменной (например, \( y' = f(y) \)). Такие уравнения описывают системы, свойства которых не меняются со временем.
• Неавтономные: правая часть зависит от \( x \) (например, \( y' = f(x, y) \)).

Основные понятия

Решение и интегральная кривая

Решением ОДУ на интервале \( I \) называется функция \( y = \varphi(x) \), определённая на \( I \) и обращающая уравнение в тождество при подстановке. График решения называется интегральной кривой. Общее решение содержит произвольные постоянные (число которых равно порядку уравнения). Частное решение получается из общего фиксацией постоянных с помощью начальных или граничных условий.

Задача Коши

Задача Коши (начальная задача) для ОДУ \( n \)-го порядка заключается в нахождении решения, удовлетворяющего \( n \) начальным условиям: \( y(x_0) = y_0 \), \( y'(x_0) = y_1 \), …, \( y^{(n-1)}(x_0) = y_{n-1} \). Теорема Коши (Пикара — Линделёфа) утверждает, что если функция \( f(x, y, y', \dots, y^{(n-1)}) \) непрерывна и удовлетворяет условию Липшица по переменным \( y, y', \dots, y^{(n-1)} \), то задача Коши имеет единственное решение в некоторой окрестности точки \( x_0 \).

Краевая задача

Краевая задача задаётся условиями на значения функции и/или её производных на границах интервала. Например, для уравнения \( y'' + \lambda y = 0 \) на отрезке \( [0, L] \) с условиями \( y(0) = 0 \), \( y(L) = 0 \) — это задача Штурма — Лиувилля. Решения существуют лишь при определённых значениях параметра \( \lambda \) (собственные значения).

Методы решения

Точные методы

Для некоторых типов ОДУ существуют аналитические методы:

• Уравнения с разделяющимися переменными: \( y' = f(x)g(y) \) — интегрирование после разделения.
• Линейные уравнения первого порядка: \( y' + p(x)y = q(x) \) — метод вариации постоянной или интегрирующий множитель.
• Уравнения в полных дифференциалах: \( M(x,y)dx + N(x,y)dy = 0 \) при условии \( \frac{\partial M}{\partial y} = \frac{\partial N}{\partial x} \).
• Линейные уравнения с постоянными коэффициентами: решение ищется в виде \( y = e^{rx} \) (характеристическое уравнение).
• Уравнения Эйлера: сводятся к линейным с постоянными коэффициентами заменой \( x = e^t \).
• Уравнения Бернулли: \( y' + p(x)y = q(x) y^n \) — сводятся к линейным подстановкой \( z = y^{1-n} \).

Численные методы

Когда точное решение получить невозможно (большинство нелинейных уравнений), применяются численные методы:

• Метод Эйлера: простейший одношаговый метод, основанный на аппроксимации производной конечной разностью.
• Методы Рунге — Кутты: семейство методов различного порядка точности (наиболее популярен метод Рунге — Кутты 4-го порядка).
• Многошаговые методы (методы Адамса): используют информацию о нескольких предыдущих точках.
• Неявные методы (например, метод трапеций): более устойчивы для жёстких систем.

Качественные методы

Качественная теория изучает поведение решений без их явного нахождения:

• Фазовый портрет: для автономных систем второго порядка изображается на плоскости (\( y, y' \)).
• Особые точки: точки, где правая часть обращается в нуль (например, устойчивый узел, седло, центр).
• Предельные циклы: замкнутые траектории, к которым стремятся соседние решения (характерны для автоколебательных систем, например, уравнение Ван дер Поля).

Применение

ОДУ используются для моделирования в различных областях.

Физика и механика

• Второй закон Ньютона: \( m \frac{d^2 x}{dt^2} = F(x, v, t) \) — основа классической механики.
• Уравнения гармонического осциллятора: \( \ddot{x} + \omega^2 x = 0 \) (пружинный маятник, колебательный контур).
• Уравнения радиоактивного распада: \( \frac{dN}{dt} = -\lambda N \).
• Уравнения движения планет (задача двух тел).

Биология и химия

• Модель Мальтуса (экспоненциальный рост популяции): \( \frac{dP}{dt} = rP \).
• Логистическое уравнение (Ферхюльста): \( \frac{dP}{dt} = rP(1 - \frac{P}{K}) \).
• Модель Лотки — Вольтерры (хищник — жертва): система двух ОДУ.
• Кинетика химических реакций: скорости реакций описываются ОДУ на концентрации (например, реакция первого порядка: \( \frac{d[A]}{dt} = -k[A] \)).

Экономика

• Модель Солоу (экономический рост): \( \dot{k} = s f(k) - (n + \delta) k \).
• Модель ценообразования (динамика спроса и предложения).

Техника

• Уравнения электрических цепей (RL, RC, RLC-цепи).
• Уравнения теплопроводности (после дискретизации по пространству — система ОДУ).
• Уравнения динамики механических систем (робототехника, вибрации).

Примеры

1. Уравнение роста численности населения (Мальтус): \( \frac{dP}{dt} = 0.02 P \), где \( P(t) \) — численность в момент \( t \). Решение: \( P(t) = P_0 e^{0.02 t} \).
2. Уравнение затухающих колебаний: \( \ddot{x} + 2\beta \dot{x} + \omega_0^2 x = 0 \). При \( \beta < \omega_0 \) — затухающие колебания, при \( \beta \ge \omega_0 \) — апериодический режим.
3. Уравнение Бернулли: \( y' + \frac{2}{x} y = x^2 y^3 \). Подстановка \( z = y^{-2} \) сводит к линейному уравнению.

Интересные факты

• Уравнение \( y' = y^2 + 1 \) имеет решение \( y = \tan(x + C) \), которое обращается в бесконечность за конечное время (явление «blow-up»).
• Система Лоренца (три ОДУ) демонстрирует хаотическое поведение при определённых параметрах, что стало основой теории хаоса.
• Для многих нелинейных ОДУ (например, \( y'' + \sin y = 0 \) — уравнение математического маятника) решения выражаются через эллиптические функции, которые не являются элементарными.

Источники

• Арнольд В. И. «Обыкновенные дифференциальные уравнения». — М.: МЦНМО, 2012.
• Петровский И. Г. «Лекции по теории обыкновенных дифференциальных уравнений». — М.: Наука, 1970.
• Понтрягин Л. С. «Обыкновенные дифференциальные уравнения». — М.: Наука, 1974.
• Тихонов А. Н., Васильева А. Б., Свешников А. Г. «Дифференциальные уравнения». — М.: Наука, 1985.
• Эльсгольц Л. Э. «Дифференциальные уравнения и вариационное исчисление». — М.: Наука, 1969.