Простейшие дифференциальные уравнения

Простейшие дифференциальные уравнения

Содержание

Введение
. Движение тела под действием постоянной силы

. Уравнение гармонического осциллятора
. Математический маятник
. Движение планет вокруг Солнца
Литература
Введение

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

Иссаку Ньютону принадлежит огромное число фундаментальных научных результатов. В те времена, недостаточно проверенные открытия, было модно в целях закрепления приоритета, публиковать в виде некоторых анаграмм (в формуле открытия переставляются буквы способом, известным лишь его автору). Заметим, что многие ученые увлекались расшифровками (обычно неудачными) чужих анаграмм. В их числе был и Иоганн Кеплер.
Из всех своих открытий И. Ньютон опубликовал в виде анаграммы только одно. В переводе на современный язык, оно звучит следующим образом: «Полезно решать дифференциальные уравнения».
Определение. Обыкновенным дифференциальным уравнением называется уравнение, содержащее неизвестную функцию и некоторое число ее производных. Решением дифференциального уравнения называется функция, при подстановке которого в уравнение, оно обращает в тождество.
Таким образом, в отличие от алгебраического уравнения, где неизвестными являются его корни (т.е. некоторая совокупность чисел), в дифференциальном уравнении неизвестной является функция. Например, уравнение математического маятника, которое описывает изменение со временем отклонения его подвеса от вертикального положения (вывод этого уравнения см. далее), имеет вид .
Оно представляет собой дифференциальное уравнение второго порядка относительно неизвестного .
Порядком ОДУ называется порядок старшей производной от неизвестной функции, входящей в рассматриваемое уравнение.
В дифференциальное уравнение могут входить несколько разных производных. Например, при учете влияния сопротивления воздуха или жидкости, уравнение, описывающие колебания физического маятника имеет вид .
В отличие от уравнения (3), в этом уравнении появился член , пропорциональный угловой скорости движения.
Проиллюстрируем появление ОДУ при решении простейших задач механики.

1. Движение тела под действием постоянной силы

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

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

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


,

где - некоторая произвольная постоянная.

Беря неопределенный интеграл от обеих частей последнего уравнения, находим


,

где - еще одна произвольная постоянная, которая возникает за счет второго интегрирования.

Формула приобретает совершенно иной вид, если ввести более привычные обозначения для входящих в нее постоянных величин. Пусть есть постоянное во времени ускорение. Константа имеет смысл начальной координаты тела (поскольку из (1) следует, что ) и ее принято обозначать символом . Произвольная постоянная имеет смысл начальной скорости, поскольку из формулы (6) ясно, что . В силу этого смысла, принято обозначать символом . Тогда наше дифференциальное уравнение принимает хорошо знакомый из школьного курса физики вид


.

Это выражение есть формула для перемещения материальной точки при равноускоренном движении по прямой.

. Уравнение гармонического осциллятора

Рис. 1

Рассмотрим колебания грузика массой m на пружинке с коэффициентом жесткости k, который лежит на плоском горизонтальном столе, предполагая, что трение грузика об поверхности стола отсутствует. Если грузик вывести из положения равновесия, он будет совершать колебания относительно этого положения. Эти колебания мы будем описываем зависящей от времени функцией , считая, что она определяет отклонение грузика из своего положения равновесия в момент времени t.

В горизонтальном направлении на грузик действует только одна сила - сила упругости пружинки, , определенная известным законом Гука


.

Деформация пружины является функцией времени, в силу чего, также является переменной.

Из второго закона Ньютона имеем


,

поскольку ускорение является второй производной от смещения : .

Уравнение (9) можно переписать в форме


,

где . Это уравнение получило название уравнение гармонического осциллятора.

Замечание. В математической литературе, при написании дифференциального уравнения обычно не указывают аргумент (t) около всех, зависящих от него функций. Такая зависимость предполагается по умолчанию. При использовании же математического пакета Maple в (10) необходимо указывать явную зависимость функции .
В отличие от предыдущего примера движения тела под действием постоянной силы в нашем случае сила изменяется с течением времени, и уравнение (10) уже нельзя решить с помощью обычной процедуры интегрирования. Попытаемся угадать решение этого уравнения, зная, что оно описывает некоторый колебательный процесс. В качестве одного из возможных решений уравнения (10) можно выбрать следующую функцию:


.

Дифференцируя функцию (11), имеем

Подставляя выражение (12) в уравнение (10), убеждаемся, что оно удовлетворяется тождественно при любом значении t.

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

с постоянными коэффициентами A и B также является решениеv уравнения гармонического осциллятора.

Можно доказать, что зависящее от двух постоянных решение (13) является общим решением уравнения гармонического осциллятора (10). Это означает, что формула (13) исчерпывает все возможные решения этого уравнения. Иными словами, других частных решений, кроме тех, которые получаются из формулы (13) фиксацией произвольных постоянных А и В, уравнение гармонического осциллятора не имеет.
Заметим, что в физике наиболее часто приходится искать именно некоторые частные решения отдельных ОДУ или их систем. Рассмотрим этот вопрос более подробно.
Возбудить колебания в рассматриваемой нами системе грузика на пружинке можно разными способами. Пусть мы задали следующие начальные условия

Полагая в этом выражении t=0 и учитывая второе начальное условие из (14), получим , отсюда следует, что A=0 и, таким образом, исходное частное решение имеет вид

.

Оно описывает колебательный режим рассматриваемой механической системы, который определяется условиями начального возбуждения (14).

Из школьного курса физики известно, что в формуле (16) a является амплитудой колебаний (она задает максимальную величину отклонения грузика от своего положения равновесия), является циклической частотой, а - фазой колебаний (начальная фаза оказывается при этом равной нулю).
Уравнение гармонического осциллятора (10) является примером линейного ОДУ. Это значит, что неизвестная функция и все ее производные входят в каждый член уравнения в первой степени. Линейные дифференциальные уравнения обладают чрезвычайно важным отличительным свойством: они удовлетворяют принципу суперпозиции. Это значит, что любая линейная комбинация двух каких либо решений линейного ОДУ также является его решением.
В рассматриваемом нами примере уравнения гармонического осциллятора, произвольная линейная комбинация двух частных решений и является не просто каким-то новым решением, но общим решением этого уравнения (оно исчерпывает все возможные его решения).
В общем случае, это не так. Например, если бы мы имели дело с линейным дифференциальным уравнением третьего порядка, (т.е. если бы в уравнение входила бы третья производная ), то линейная комбинация каких-либо двух его частных решений также была бы решением этого уравнения, но не представляла бы собой его общее решение.
В курсе дифференциальных уравнений доказывается теорема о том, что общее решение ОДУ N-ого порядка (линейного или нелинейного) зависит от N произвольных постоянных. В случае нелинейного уравнения эти произвольные постоянные могут входить в общее решение (в отличие от (13)), нелинейным образом.
Принцип суперпозиции играет в теории ОДУ исключительно важную роль, поскольку с его помощью можно построить общее решение дифференциального уравнения в виде суперпозиции его частных решений. Например, для случая линейных ОДУ с постоянными коэффициентами и их систем (уравнение гармонического осциллятора относится именно к этому типу уравнений) в теории дифференциальных уравнений разработан общий метод решения. Суть его заключается в следующем. Ищется частное решение в виде . В результате его подстановки в исходное уравнение, все зависящие от времени множители сокращаются и мы приходим к некоторому характеристическому уравнению, которое для ОДУ N-ого порядка представляет собой алгебраическое уравнение N-ой степени. Решая его, мы находим, тем самым, все возможные частные решения, произвольная линейная комбинация которых и дает общее решение исходного ОДУ. Мы не будем далее останавливаться на этом вопросе, отсылая читателя к соответствующим учебникам по теории дифференциальным уравнениям, в которых можно найти дальнейшие детали, в частности, рассмотрение случая, когда характеристическое уравнение содержит кратные корни.
Если рассматривается линейное ОДУ с переменными коэффициентами, (его коэффициенты зависят от времени), то принцип суперпозиции также справедлив, но построить в явном виде общее решение этого уравнение каким-либо стандартным методом, уже не представляется возможным. Мы вернемся к этому вопросу далее, обсуждая явление параметрического резонанса и связанным с его исследованием уравненем Матье.

3. Математический маятник

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

Рис. 2

В силу нерастяжимости нити мы заведомо знаем, что маятник, точнее соответствующая ему материальная точка с массой m, движется вдоль дуги окружности радиуса l (см. рис.2). Отсюда следует, что движение маятника является одномерным и мы можем в качестве координаты, описывающей это движение, выбрать длину дуги окружности s. Отсчет этой криволинейной координаты начинается в точке равновесия маятника A и ведется в направлении положительных углов отклонения подвеса маятника относительно вертикального положения (т.е. против часовой стрелки).

Обыкновенное дифференциальное уравнение, описывающее движение маятника в момент времени t, который отвечает отклонению его подвеса на угол (см. рис 2.), получим с помощью второго закона Ньютона. На материальную точку массой m действуют две силы - сила притяжения к Земле mg, направленная вертикально вниз, и сила натяжения нити, направленная по радиусу окружности, соответствующему выбранному нами отклонению на угол (t) от вертикального положения. Чтобы применить второй закон Ньютона для описания движения материальной точки m по дуге окружности s, необходимо найти проекцию равнодействующей силы на направление касательной к окружности в той ее точке, где в момент времени t находится маятник. Для этого достаточно найти сумму проекций на касательную в точке B вышеуказанных двух сил. Поскольку сила натяжения нити перпендикулярна касательной, достаточно учесть только проекцию F силы тяжести mg. Из рис.2 имеем F=-mg sin( ) (угол отклонения маятника равен углу обозначенному той же буквой в силовом треугольнике, как углы со взаимно параллельными сторонами). С учетом вышесказанного, из второго закона Ньютона получим


,

здесь есть вторая производная от криволинейного пути перемещения маятника s(t), а угол измеряется в радианах. Принципиально важным является правильность учета знаков проекций силы mg: из рисунка видно, что при отклонении маятника от положения равновесия вправо, эта проекция направлена влево и наоборот (таким образом, она является возвращающей силой).

Литература

1. Кунин С. Вычислительная физика. М.: Мир, 2002.

. Федоренко Р.П. Введение в вычислительную физику. М.: Изд-во МФТИ, 2004.
. Валуев А.А., Норман Г.Э., Подлипчук В.Ю. - В сб.: Математическое моделирование. М.: Наука, 2009, с. 5-40.