Дифференциальные уравнения: интуитивное объяснение через поле направлений

В обычном уравнении нужно найти число. Оно даёт равенство вроде 2x + 3 = 11 и спрашивает, при каком значении это равенство верно. В дифференциальном уравнении нужно найти кривую. Оно задаёт локальное правило изменения и спрашивает, какие траектории могут следовать ему на всём интервале.

Именно поэтому предмет поначалу кажется непривычным. Неизвестной является функция, ответ часто содержит произвольную постоянную, а одной точки на кривой недостаточно. Уравнение описывает, как кривая должна двигаться на всём интервале.

Явные дифференциальные уравнения первого порядка проще понять через поле направлений, составленное из коротких отрезков. Для уравнения вида y' = F(x, y) каждой точке, где определена F, соответствует отрезок, показывающий, в какую сторону должна идти кривая. Решение представляет собой траекторию, которая следует этим направлениям. Когда вы видите такое поле, семейства решений, начальные условия и экспоненциальный рост перестают казаться отдельными приёмами.

Дифференциальное уравнение задаёт правило изменения

Производная показывает, насколько быстро изменяется величина. Дифференциальное уравнение разворачивает эту идею в обратную сторону. Оно задаёт связь с производной и просит восстановить саму величину.

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

$$y' = 2y$$

Штрих означает производную y по x. Уравнение говорит, что скорость изменения всегда вдвое больше текущего значения. Если y равно 1, наклон должен быть равен 2. Если y равно 3, наклон должен быть равен 6. Если y отрицательно, наклон тоже должен быть отрицательным. Правило не сообщает напрямую высоту кривой при конкретном x. Оно описывает, как значение функции и её наклон должны изменяться вместе.

В этом и состоит главное изменение постановки:

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

Уравнение y' = 2y имеет первый порядок, потому что его старшая производная является первой производной. Уравнение с y'' имеет второй порядок. Порядок показывает, производная какого наивысшего порядка входит в уравнение, а не насколько сложным оно выглядит.

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

Решение является функцией, а не числом

Чтобы решить y' = 2y, нужно найти функцию, производная которой вдвое больше самой функции. Именно так ведут себя экспоненциальные функции. Попробуем

$$y = Ce^{2x}$$

где C является любой постоянной. Возьмём производную:

$$y' = 2Ce^{2x}$$

Поскольку сама функция y равна $Ce^{2x}$, получаем y' = 2y. Предложенная функция удовлетворяет уравнению при любом x, а значит, является решением.

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

Постоянная C не означает, что решение не доведено до конца. Она является частью ответа. Уравнение описывает семейство кривых:

• $C = 1$ даёт $y = e^{2x}$.
• $C = 3$ даёт $y = 3e^{2x}$.
• $C = -2$ даёт $y = -2e^{2x}$.
• C = 0 даёт постоянное решение y = 0.

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

По той же причине при интегрировании появляется аддитивная постоянная. Дифференцирование стирает постоянные вертикальные сдвиги, поскольку их производная равна нулю. Поэтому обратный процесс должен сохранять все возможные сдвиги. Последующие алгебраические преобразования могут изменить вид параметра. В нашем примере после возведения в степень $e^{C}$ превращается в множитель перед $e^{2x}$.

Поле направлений превращает правило в ландшафт

Явное уравнение первого порядка y' = F(x, y) можно понять ещё до того, как вы найдёте его символическое решение. Возьмите сетку точек, в которых определена F. В каждой точке вычислите требуемый наклон и нарисуйте короткий отрезок с таким наклоном. Полученная картина называется полем наклонов, или полем направлений.

Для y' = 2y наклон зависит только от высоты:

• Вдоль y = 0 все отрезки горизонтальны.
• Вдоль y = 1 каждый отрезок имеет наклон 2.
• Вдоль y = -1 каждый отрезок имеет наклон -2.
• Чем выше ряд с положительным y, тем круче отрезки направлены вверх.
• Чем ниже ряд с отрицательным y, тем круче отрезки направлены вниз.

Поле выглядит как горизонтальные полосы с одинаковым наклоном внутри каждой полосы. Кривая решения может начинаться в любой точке, но затем обязана следовать локальным отрезкам. Положительные решения изгибаются вверх и растут. Отрицательные решения уходят вниз от нуля. Нулевое решение остаётся горизонтальным, потому что все стрелки на этой линии горизонтальны.

Эта картина отделяет уравнение от любого отдельного решения. Дифференциальное уравнение создаёт всё поле. Каждое решение представляет собой одну нить, проходящую через него.

Поле также позволяет предсказывать поведение без явной формулы. Если в некоторой области отрезки направлены вверх, решения там возрастают. Если вдоль y = c они горизонтальны при всех нужных x, то y = c является постоянным стационарным решением. Если соседние отрезки направлены к этой линии, близкие решения возвращаются к ней. Если они направлены от линии, малые возмущения растут. Значительная часть прикладной теории дифференциальных уравнений строится на таком качественном анализе: сначала понять, как могут вести себя траектории, и только затем выводить точные формулы.

В уравнении y' = 2y правая часть гладко зависит от y, поэтому стандартные условия существования и единственности гарантируют ровно одно решение, проходящее через каждую точку. Поле направлений иллюстрирует этот результат, но не доказывает его. Общей точки и одинакового наклона недостаточно, чтобы гарантировать совпадение решений для любого уравнения: менее регулярные правила могут позволить им встретиться, а затем разойтись.

Начальные условия выбирают одну траекторию

Добавим к тому же уравнению условие y(0) = 3. Это начальное условие: при x = 0 значение функции равно 3. Начнём с семейства решений

$$y = Ce^{2x}$$

и подставим известную точку:

$$3 = Ce^{0} = C$$

Значит, частное решение имеет вид

$$y = 3e^{2x}$$

Дифференциальное уравнение задало форму семейства. Начальное условие выбрало один его элемент. На поле направлений это всё равно что поставить палец в точку (0, 3): остаётся только проходящая через неё траектория.

Для регулярных задач Коши необходимое количество условий обычно соответствует порядку уравнения. Уравнение первого порядка y' = v(x, y), в котором скорость уже задана, обычно требует одного начального положения. Для уравнения второго порядка, задающего ускорение, как правило, нужны и начальное положение, и начальная скорость.

Например, нормированное уравнение незатухающего пружинного осциллятора

$$y'' + y = 0$$

говорит, что ускорение направлено к положению равновесия. Физическая форма имеет вид m y'' + k y = 0. Если выбрать единицы так, чтобы k/m = 1, получится приведённое выше нормированное уравнение. Его решения представляют собой комбинации синуса и косинуса. Одного начального положения недостаточно, чтобы понять, отпустили ли пружину из состояния покоя или она уже двигалась. Добавьте начальную скорость, и движение станет определённым.

Разделение переменных разносит x и y по разным сторонам

Уравнение y' = 2y можно решить систематически, а не только угадать ответ. Запишем производную в виде $\frac{dy}{dx}$:

$$\frac{dy}{dx} = 2y$$

Перенесём часть с y к dy, а часть с x к dx:

$$\frac{dy}{y} = 2 dx$$

Теперь проинтегрируем обе стороны:

$$\ln|y| = 2x + C$$

Возводя e в степень обеих частей, получаем

$$y = Ae^{2x}$$

На этом этапе A не равно нулю, поскольку при делении на y мы исключили y = 0. Проверим исключённую функцию в исходном уравнении: y = 0 действительно подходит. Поэтому полное семейство можно записать как y = $Ae^{2x}$ для любой вещественной постоянной A, включая ноль.

Метод называется разделением переменных, потому что перед интегрированием все множители с y собирают рядом с dy, а все множители с x рядом с dx. Он работает для уравнений, которые можно привести к виду

$$\frac{dy}{dx} = f(x)g(y)$$

но в алгебре скрыта ловушка. Деление на g(y) может отбросить решения, для которых g(y) = 0. В нашем примере деление на y временно исключило стационарное решение y = 0.

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

Рост, убывание и колебания задаются разными правилами

Небольшие изменения в дифференциальном уравнении приводят к совершенно разным типам движения.

Модель

$$y' = ky$$

описывает экспоненциальный рост положительных величин при положительном k и экспоненциальное убывание к нулю при отрицательном k. В более общем случае модуль любого ненулевого решения растёт при положительном k и уменьшается при отрицательном k. Модель возникает, когда скорость пропорциональна текущему количеству. Знакомые примеры включают идеализированный рост популяции, непрерывное начисление процентов и радиоактивный распад. Одна и та же формула не означает, что любая реальная популяция или инвестиция всегда меняется экспоненциально. Она говорит лишь о том, что экспоненциальное поведение следует из предположения о постоянной относительной скорости, пока это предположение остаётся верным.

В отличие от него, уравнение пружины y'' + y = 0 описывает колебания. Его решения строятся из синуса и косинуса, потому что после двойного дифференцирования функция меняет знак. Сравните это уравнение с

$$y'' - y = 0$$

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

В реальные модели добавляют члены, отвечающие за внешние силы, трение, ограниченные ресурсы или взаимодействие нескольких величин. В логистической модели популяции рост замедляется возле ёмкости среды. У затухающей пружины амплитуда со временем уменьшается. Система дифференциальных уравнений может одновременно отслеживать несколько изменяющихся величин. Каждый дополнительный член меняет локальное правило, а вместе с ним и траектории, по которым способна двигаться система.

Откуда берутся ошибки

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

Вторая ошибка состоит в проверке только начального условия. Кривая не становится решением y' = 2y только потому, что проходит через (0, 3). Через эту точку проходит бесконечно много кривых. Кривая должна также удовлетворять дифференциальному уравнению на всём интервале.

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

Четвёртая ошибка состоит в потере стационарного решения при делении. Если вы делите на y, y - 1 или любое выражение с неизвестной функцией, остановитесь и проверьте в исходном уравнении все случаи, когда это выражение равно нулю.

Пятая ошибка состоит в том, что модель принимают за сам объект. Уравнение y' = ky приводит к экспоненциальному росту или убыванию, потому что предполагает постоянную относительную скорость. Реальная система следует этой кривой лишь до тех пор, пока предположение остаётся полезным приближением. Правильное решение уравнения не доказывает, что оно правильно моделирует систему.

Как здесь помогает Math Zen

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

Math Zen включает задачи на уравнения с разделяющимися переменными, линейные уравнения первого порядка, уравнения второго порядка, однородные уравнения, уравнения в полных дифференциалах, уравнения Бернулли и системы обыкновенных дифференциальных уравнений. Перед вычислениями задайте те же три вопроса, что и в этой статье: какого изменения требует уравнение, какое семейство траекторий может ему подчиняться и какие условия выбирают нужную траекторию?

Полезный цикл практики прост: классифицируйте уравнение, предскажите направление его траекторий, решите или проверьте его, а затем проверьте начальное условие и все исключённые стационарные решения. Вернитесь к этой теме после перерыва, используя тот же подход с интервальным повторением, что и в остальных разделах Math Zen. Цель состоит не в том, чтобы запомнить набор рецептов, а в том, чтобы сразу видеть правило изменения перед собой.

Главное

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

Для явных уравнений первого порядка всю тему связывает воедино поле направлений. Локальное правило заполняет область определения допустимыми направлениями, а решение следует им. Разделение переменных, задачи Коши, рост, убывание и колебания дают разные способы читать этот ландшафт и прокладывать по нему траектории. Дифференциальное уравнение превращает локальное правило изменения в глобальную кривую.