Русская Википедия:Обыкновенное дифференциальное уравнение
Обыкновенное дифференциальное уравне́ние (ОДУ) — дифференциальное уравнение для функции от одной переменной. (Этим оно отличается от уравнения в частных производных, где неизвестная — функция нескольких переменных.) Таким образом, ОДУ — уравнения вида
- <math>F(x,y,y',y,...,y^{(n)})=0, \qquad (1)</math>
где <math>y(x)</math> — неизвестная функция (возможно, вектор-функция, тогда <math>F</math>, как правило, тоже вектор-функция со значениями в пространстве той же размерности; в этом случае говорят о системе дифференциальных уравнений), зависящая от независимой переменной <math>x</math>, штрих означает дифференцирование по <math>x</math>. Число <math>n</math> (порядок старшей производной, входящей в данное уравнение) называется порядком дифференциального уравнения (1).
Независимая переменная <math>x</math> часто интерпретируется (особенно в дифференциальных уравнениях, возникающих в физических и других естественно-научных задачах) как время, поэтому её часто обозначают буквой <math>t</math>. Переменная <math>y</math> — некоторая величина (или совокупность величин, если <math>y</math> является вектор-функцией), изменяющаяся со временем. Например, <math>y</math> может означать набор координат точки в пространстве; в этом случае уравнение (1) описывает движение точки в пространстве, то есть изменение её координат с течением времени. Независимая переменная <math>x</math> обычно принимает вещественные значения, однако рассматриваются и дифференциальные уравнения, в которых переменная <math>x</math> комплексная (так называемые уравнения с комплексным временем).
Наиболее часто встречаются дифференциальные уравнения вида
- <math>y^{(n)} = f(x,y,y',y,...,y^{(n-1)}), \qquad (2)</math>
в которых старшая производная <math>y^{(n)}</math> выражается в виде функции от переменных <math>x,</math> <math>y</math> и производных <math>y^{(i)}</math> порядков меньше <math>n.</math> Такие дифференциальные уравнения называются нормальными или разрешёнными относительно производной.
В противоположность уравнениям вида (2), дифференциальные уравнения вида (1) называются уравнениями, не разрешёнными относительно производной или неявными дифференциальными уравнениями.
Классическим решением дифференциального уравнения (2) называется <math>n</math> раз дифференцируемая функция <math>y(x)</math>, удовлетворяющая уравнению во всех точках своей области определения. Обычно существует целое множество таких функций, и для выбора одного из них требуется наложить на него дополнительное условие. Начальным условием для уравнения (2) называется условие
- <math>y(x_0)=y_0, \ y'(x_0)=y_0^{(1)}, y(x_0)=y_0^{(2)}, \, \ldots, \, y^{(n-1)}(x_0)=y_0^{(n-1)}, \qquad (3)</math>
где <math>x_0</math> — некоторое фиксированное значение независимой переменной (фиксированный момент времени), а <math>y_0</math> и <math>y_0^{(i)}</math> — соответственно, фиксированные значения функции <math>y</math> и всех её производных до порядка <math>n-1</math> включительно. Дифференциальное уравнение (2) вместе с начальным условием (3) называется начальной задачей или задачей Коши:
- <math>\left\{\begin{array}{lcl}y^{(n)} = f(x,y,y',y,...,y^{(n-1)}), \\ {} \\ y(x_0)=y_0, \ y'(x_0)=y_0^{(1)}, y(x_0)=y_0^{(2)}, \, \ldots, \, y^{(n-1)}(x_0)=y_0^{(n-1)}. \end{array}\right.</math>
Теорема существования и единственности решения обыкновенного дифференциального уравнения описывает совокупность всех решений обыкновенного дифференциального уравнения. Является основным теоретическим положением при изучении обыкновенных дифференциальных уравнений.[1]
Теорема Пикара утверждает, что при достаточно общих ограничениях на функцию <math>f</math>, стоящую в правой части уравнения (2), задача Коши для этого уравнения имеет единственное решение, определённое на некотором интервале оси времени <math>x</math>, содержащем начальное значение <math>x_0</math> (этот интервал, вообще говоря, может не совпадать со всей осью). Основные задачи и результаты теории дифференциальных уравнений: существование и единственность решения различных задач для ОДУ, методы решения простейших ОДУ, качественное исследование решений ОДУ без нахождения их явного вида.
Теорема о зависимости решений обыкновенных дифференциальных уравнений от параметров формулирует зависимость свойств решений обыкновенных дифференциальных уравнений от параметров этих уравнений.
История
Дифференциальные уравнения встречались уже в работах И. Ньютона и Г. Лейбница; термин «дифференциальные уравнения» принадлежит Лейбницу. Ньютон при создании исчисления «флюксий» и «флюент» ставил две задачи: по данному соотношению между флюентами определить соотношение между флюксиями; по данному уравнению, содержащему флюксии, найти соотношение между флюентами. С современной точки зрения, первая из этих задач (вычисление по функциям их производных) относится к дифференциальному исчислению, а вторая составляет содержание теории обыкновенных дифференциальных уравнений. Задачу нахождения неопределённого интеграла F(x) функции f(x) Ньютон рассматривал просто как частный случай его второй задачи. Такой подход был для Ньютона как создателя основ математического естествознания вполне оправданным: в очень большом числе случаев законы природы, управляющие теми или иными процессами, выражаются в форме дифференциальных уравнений, а расчёт течения этих процессов сводится к решению дифференциального уравнения.[2]
Основное открытие Ньютона, то, которое он счёл нужным засекретить и опубликовал лишь в виде анаграммы, состоит в следующем: «Data aequatione quotcunque fluentes quantitae involvente fluxiones invenire et vice versa». В переводе на современный математический язык это означает: «Полезно решать дифференциальные уравнения». В настоящее время теория дифференциальных уравнений представляет собой трудно обозримый конгломерат большого количества разнообразных идей и методов, в высшей степени полезный для всевозможных приложений и постоянно стимулирующий теоретические исследования во всех отделах математики.[3] [4]
Примеры
- Одно из простейших применений дифференциальных уравнений — решение нетривиальной задачи нахождения траектории тела по известным проекциям ускорения. Например, в соответствии со вторым законом Ньютона, ускорение тела пропорционально сумме действующих сил; соответствующее дифференциальное уравнение имеет вид <math>m \ddot{x}= F(x,t)</math>. Зная действующие силы (правая часть), можно решить это уравнение и, учитывая начальные условия (координаты и скорость в начальный момент времени), найти траекторию движения точки.
- Дифференциальное уравнение <math>y'=y</math>, вместе с начальным условием <math>y(0)=1</math>, задаёт экспоненту: <math>y(x)=e^x</math>. Если <math>x</math> обозначает время, то эта функция описывает, например, рост популяции в условиях неограниченности ресурсов, а также и многое другое.
- Решением дифференциального уравнения <math>y'=f(x)</math>, правая часть которого не зависит от неизвестной функции, является неопределённый интеграл
- <math>y(x)=\int\! f(x)\,dx+C,</math>
где <math>C</math> — произвольная константа.
Дифференциальные уравнения первого порядка
Уравнения с разделяющимися переменными
Дифференциальное уравнение <math>\dot{y} = f(x,y)</math> называется уравнением с разделяющимися (отделяющимися) переменными, если его правая часть представима в виде <math>y' = f_1(x)f_2(y)</math>. Тогда, в случае <math>f_2(y)\neq 0</math>, общим решением уравнения является <math>\int\!\frac{dy}{f_2(y)} = \int\! f_1(x)\,dx</math>.
Примеры физических задач, приводящих к уравнениям с разделяющимися переменными
Охлаждение тела
Пусть <math>T</math> — температура тела, <math>T_0</math> — температура окружающей среды (<math>T > T_0</math>). Пусть <math>Q</math> — количество теплоты, <math>c</math> — удельная теплоёмкость. Тогда количество теплоты, передаваемое окружающей среде до выравнивания температур, выражается формулой <math>Q = mc(T-T_0)</math>, или, в дифференциальной форме, <math>dQ = mc\,dT</math>. С другой стороны, скорость отдачи тепла можно выразить в виде <math>dQ = -k(T-T_0)\,dt</math>, где <math>k</math> — некий коэффициент пропорциональности. Исключая из этих двух уравнений <math>dQ</math>, получаем уравнение с разделяющимися переменными:
- <math>mc\,dT = -k(T-T_0)\,dt</math>.
Общим решением этого уравнения является семейство функций <math>T = T_0+C e^{-\frac{kt}{mc}}</math>.
Однородные уравнения
Шаблон:Main Дифференциальное уравнение <math>\dot{y} = f(x,y)</math> называется однородным, если <math>f(x, y)</math> — однородная функция нулевой степени. Функция <math>f(x, y)</math> называется однородной степени <math>k</math>, если для любого <math>\lambda > 0</math> выполняется равенство <math>f(\lambda x, \lambda y) = \lambda ^ k f(x,y)</math>.
Замена <math>y(x) = xz(x)</math> приводит при <math>x > 0</math> однородное уравнение к уравнению с разделяющимися переменными:
- <math>f(x, xz) = x^{0} f(1, z) = f(1, z)</math>
- <math>\dot{y} = x \dot{z} + z</math>
Подставив в исходное уравнение, получаем:
- <math>\dot{z} = \frac{1}{x} (f(1, z) - z)</math>,
что является уравнением с разделяющимися переменными.
Квазиоднородные уравнения
Дифференциальное уравнение <math>\dot{y} = f(x,y)</math> называется квазиоднородным, если для любого <math>\lambda > 0</math> выполняется соотношение <math>f\left(\lambda ^ \alpha x, \lambda ^ \beta y\right) = \lambda ^ {\beta - \alpha} f(x, y)</math>.
Данное уравнение решается заменой <math>y = z ^ {\frac{\beta}{\alpha}}</math>:
- <math>\dot{z} = \frac{\alpha}{\beta} \left(z ^ {- \frac{1}{\alpha}}\right) ^ {\beta - \alpha} f\left(x, z ^ \frac{\beta}{\alpha}\right)</math>
В силу квазиоднородности, положив <math>\lambda = z ^ {- \frac{1}{\alpha}}</math>, получаем:
- <math>\left(z ^ {- \frac{1}{\alpha}}\right) ^ {\beta - \alpha} f \left(x, z ^ \frac{\beta}{\alpha}\right) = f\left(\frac{x}{z}, 1\right)</math>
- <math>\dot{z} = \frac{\alpha}{\beta} f\left(\frac{x}{z}, 1\right)</math>,
что, очевидно, является однородным уравнением.
Линейные уравнения
Шаблон:Main Дифференциальное уравнение <math>y' + a(x)y = b(x)</math> называется линейным и может быть решено тремя методами: методом интегрирующего множителя, методом вариации постоянной или методом Бернулли.
Метод интегрирующего множителя
Пусть задана функция <math>\mu(x)</math> — интегрирующий множитель, в виде:
- <math>\mu(x) = e ^ {\int\! a(x)\,dx}</math>
Умножим обе части исходного уравнения на <math>\mu(x)</math>, получим:
- <math>\dot{y}e ^ {\int\!a(x)\,dx} + ya(x) e ^ {\int\! a(x)\,dx} = b(x) \mu (x)</math>
Легко заметить, что левая часть является производной функции <math>\mu(x) y(x)</math> по <math>x</math>. Поэтому уравнение можно переписать:
- <math>(\mu (x) y(x))' = b(x) \mu (x)</math>
Проинтегрируем:
- <math>y(x) \mu(x) = \int\!b(x) \mu (x)\,dx + C</math>
Таким образом, решение линейного уравнения будет:
- <math>y(x) = e ^ {- \int\! a(x)\,dx}\left(\int\! b(x) \mu (x)\,dx + C\right)</math>
Метод вариации постоянной (метод Лагранжа)
Рассмотрим однородное уравнение <math>\dot{y} + a(x)y = 0</math>. Очевидно, это уравнение с разделяющимися переменными, его решение:
- <math>y(x) = ce ^ {-\int\!a(x)\,dx}</math>
Решения исходного уравнения будем искать в виде:
- <math>y(x) = c(x)e ^ {-\int\!a(x)\,dx}</math>
Подставив полученное решение в исходное уравнение:
- <math>\dot{c} = b(x) e ^ {\int\!a(x)\,dx}</math>,
получаем:
- <math>c(x) = c_1 + \int\!b(x) e ^ {\int\!a(x)\,dx}\,dx</math>,
где <math>c_1</math> — произвольная константа.
Таким образом, решение исходного уравнения можно получить путём подстановки <math>c(x)</math> в решение однородного уравнения:
- <math>y(x) = e ^ {-\int\! a(x)\,dx}\left(c_1 + \int\! b(x)e ^ {\int\! a(x)\,dx}\,dx\right)</math>
Уравнение Бернулли
Шаблон:Main Дифференциальное уравнение <math>\dot{y} + a(x) y = b(x) y^n</math> называется уравнением Бернулли (при <math>n=0</math> или <math>n=1</math> получаем неоднородное или однородное линейное уравнение). При <math>n=2</math> является частным случаем уравнения Риккати. Названо в честь Якоба Бернулли, опубликовавшего это уравнение в 1695 году. Метод решения с помощью замены, сводящей это уравнение к линейному, нашёл его брат Иоганн Бернулли в 1697 году.
Биномиальное дифференциальное уравнение
Это уравнение вида
- <math>\left( y' \right)^m = f(x,y),</math> где <math>m</math> — натуральное число, а <math>f(x,y)</math> — многочлен от двух переменных[5].
Литература
Учебники
- Арнольд В. И. Обыкновенные дифференциальные уравнения, — Любое издание.
- Арнольд В. И. Дополнительные главы теории обыкновенных дифференциальных уравнений, — Любое
издание.
- Арнольд В. И. Геометрические методы в теории обыкновенных дифференциальных уравнений, — Любое издание.
- Арнольд В. И., Ильяшенко Ю. С. Обыкновенные дифференциальные уравнения, — Итоги науки и техн. Сер. Совр. пробл. мат. Фундам. направ., 1985, том 1.
- Петровский И. Г. Лекции по теории обыкновенных дифференциальных уравнений, — Любое издание.
- Понтрягин Л. С. Обыкновенные дифференциальные уравнения, — Любое издание.
- Степанов В. В. Курс дифференциальных уравнений, — Любое издание.
- Трикоми Ф. Дифференциальные уравнения, — Любое издание.
- Филиппов А. Ф. Введение в теорию дифференциальных уравнений, — Любое издание.
- Филипс Г. Дифференциальные уравнения, — Любое издание.
- Хартман Ф. Обыкновенные дифференциальные уравнения, — Любое издание.
- Эльсгольц Л. Э. Дифференциальные уравнения и вариационное исчисление, — Любое издание.
- Шаблон:Книга
- Шаблон:Книга
Задачники
- Филиппов А. Ф. Сборник задач по дифференциальным уравнениям, — Любое издание.
Справочники
- Камке Э. Справочник по обыкновенным дифференциальным уравнениям, — Любое издание.
- Зайцев В. Ф., Полянин А. Д. Справочник по обыкновенным дифференциальным уравнениям, — Любое издание.
Примечания
- ↑ Л.С. Понтрягин Дифференциальные уравнения и их приложения. - Шаблон:М., Наука, 1988. - c. 15
- ↑ БСЭ. Дифференциальные уравнения.
- ↑ Арнольд В. И. Дополнительные главы теории обыкновенных дифференциальных уравнений.
- ↑ Арнольд В. И. Геометрические методы в теории обыкновенных дифференциальных уравнений.
- ↑ Шаблон:Книга
Шаблон:Выбор языка Шаблон:Математическая физика