В математике обыкновенное дифференциальное уравнение ( ОДУ ) — это дифференциальное уравнение (ДУ), зависящее только от одной независимой переменной . Как и в случае с другими ДУ, его неизвестные состоят из одной (или нескольких) функций и включают производные этих функций. [1] Термин «обыкновенный» используется в отличие от уравнений в частных производных (УЧП), которые могут быть относительно более чем одной независимой переменной, [2] и, реже, в отличие от стохастических дифференциальных уравнений (СДУ), где прогрессия случайна. [3]
Линейное дифференциальное уравнение — это дифференциальное уравнение, которое определяется линейным полиномом относительно неизвестной функции и ее производных, то есть уравнением вида
где , ..., и — произвольные дифференцируемые функции , которые не обязательно должны быть линейными, а — последовательные производные неизвестной функции y от переменной x . [4]
Среди обыкновенных дифференциальных уравнений линейные дифференциальные уравнения играют заметную роль по нескольким причинам. Большинство элементарных и специальных функций, которые встречаются в физике и прикладной математике, являются решениями линейных дифференциальных уравнений (см. Голономная функция ). Когда физические явления моделируются нелинейными уравнениями, они обычно аппроксимируются линейными дифференциальными уравнениями для более простого решения. Несколько нелинейных ОДУ, которые могут быть решены явно, обычно решаются путем преобразования уравнения в эквивалентное линейное ОДУ (см., например, уравнение Риккати ). [5]
Некоторые ОДУ могут быть решены явно в терминах известных функций и интегралов . Когда это невозможно, уравнение для вычисления ряда Тейлора решений может быть полезным. Для прикладных задач численные методы для обыкновенных дифференциальных уравнений могут предоставить приближение решения.
Обыкновенные дифференциальные уравнения (ОДУ) возникают во многих контекстах математики , социальных и естественных наук . Математические описания изменений используют дифференциалы и производные. Различные дифференциалы, производные и функции становятся связанными посредством уравнений, так что дифференциальное уравнение является результатом, который описывает динамически изменяющиеся явления, эволюцию и вариацию. Часто величины определяются как скорость изменения других величин (например, производные смещения по времени) или градиенты величин, и именно так они входят в дифференциальные уравнения. [6]
Конкретные математические области включают геометрию и аналитическую механику . Научные области включают большую часть физики и астрономии (небесная механика), метеорологию (моделирование погоды), химию (скорости реакций), [7] биологию (инфекционные заболевания, генетическая изменчивость), экологию и моделирование популяций (конкуренция популяций), экономику (тенденции акций, процентные ставки и рыночное равновесное изменение цен).
Многие математики изучали дифференциальные уравнения и внесли свой вклад в эту область, в том числе Ньютон , Лейбниц , семья Бернулли , Риккати , Клеро , Д'Аламбер и Эйлер .
Простым примером является второй закон движения Ньютона — связь между смещением x и временем t объекта под действием силы F задается дифференциальным уравнением
которая ограничивает движение частицы постоянной массы m . В общем случае F является функцией положения x ( t ) частицы в момент времени t . Неизвестная функция x ( t ) появляется в обеих частях дифференциального уравнения и указывается в обозначении F ( x ( t )). [8] [9] [10] [11]
В дальнейшем y — зависимая переменная, представляющая неизвестную функцию y = f ( x ) независимой переменной x . Обозначение для дифференциации различается в зависимости от автора и от того, какое обозначение наиболее полезно для поставленной задачи. В этом контексте обозначение Лейбница ( dy/дх , д 2 г/дх 2 , …, д н г/дх н ) более полезна для дифференцирования и интегрирования , тогда как обозначение Лагранжа ( y ′, y ′′, …, y ( n ) ) более полезно для компактного представления производных более высокого порядка , а обозначение Ньютона часто используется в физике для представления производных более низкого порядка по времени.
Дано F , функция x , y и производные y . Тогда уравнение вида
называется явным обыкновенным дифференциальным уравнением порядка n . [12] [13]
В более общем случае неявное обыкновенное дифференциальное уравнение порядка n имеет вид: [14]
Существуют и другие классификации:
Ряд связанных дифференциальных уравнений образуют систему уравнений. Если y — вектор, элементами которого являются функции; y ( x ) = [ y 1 ( x ), y 2 ( x ),..., y m ( x )], а F — векторнозначная функция y и ее производных, то
является явной системой обыкновенных дифференциальных уравнений порядка n и размерности m . В форме вектор - столбца :
Они не обязательно линейны. Неявный аналог:
где 0 = (0, 0, ..., 0) — нулевой вектор . В матричной форме
Для системы вида некоторые источники также требуют, чтобы матрица Якоби была невырожденной , чтобы назвать это неявным ОДУ [системой]; неявная система ОДУ, удовлетворяющая этому условию невырожденности Якоби, может быть преобразована в явную систему ОДУ. В тех же источниках неявные системы ОДУ с вырожденным якобианом называются дифференциально-алгебраическими уравнениями (ДАУ). Это различие не просто терминологическое; ДАУ имеют принципиально разные характеристики и, как правило, более сложны для решения, чем (невырожденные) системы ОДУ. [18] [19] [20] Предположительно для дополнительных производных матрица Гессе и т. д. также предполагаются невырожденными в соответствии с этой схемой, [ требуется ссылка ], хотя следует отметить, что любое ОДУ порядка больше единицы может быть (и обычно так и есть) переписано как система ОДУ первого порядка, [21] что делает критерий вырожденности Якоби достаточной для того, чтобы эта таксономия была всеобъемлющей во всех порядках.
Поведение системы ОДУ можно визуализировать с помощью фазового портрета .
Дано дифференциальное уравнение
Функция u : I ⊂ R → R , где I — интервал, называется решением или интегральной кривой для F , если u n - кратно дифференцируема на I , и
Если даны два решения u : J ⊂ R → R и v : I ⊂ R → R , u называется расширением v , если I ⊂ J и
Решение, не имеющее расширения, называется максимальным решением . Решение, определенное на всем R, называется глобальным решением .
Общее решение уравнения n-го порядка — это решение, содержащее n произвольных независимых констант интегрирования . Частное решение выводится из общего решения путем задания константам конкретных значений, часто выбираемых для выполнения набора « начальных условий или граничных условий ». [22] Особое решение — это решение, которое не может быть получено путем присвоения определенных значений произвольным константам в общем решении. [23]
В контексте линейного ОДУ термин частное решение может также относиться к любому решению ОДУ (не обязательно удовлетворяющему начальным условиям), которое затем добавляется к однородному решению (общему решению однородного ОДУ), которое затем образует общее решение исходного ОДУ. Это терминология, используемая в разделе метода угадывания в этой статье, и часто используется при обсуждении метода неопределенных коэффициентов и вариации параметров .
Для нелинейных автономных ОДУ при некоторых условиях возможно разработать решения конечной продолжительности, [24] имея в виду, что из своей собственной динамики система достигнет нулевого значения в конечное время и останется там в нуле навсегда. Эти решения конечной продолжительности не могут быть аналитическими функциями на всей действительной прямой, и поскольку они будут нелипшицевыми функциями в конечное время, они не включены в теорему единственности решений дифференциальных уравнений Липшица.
Например, уравнение:
Допускает решение с конечной продолжительностью:
Теория сингулярных решений обыкновенных и частных дифференциальных уравнений была предметом исследований со времен Лейбница, но только с середины девятнадцатого века ей стало уделяться особое внимание. Ценная, но малоизвестная работа по этой теме — работа Хоутэна (1854). Дарбу (с 1873) был лидером в этой теории, и в геометрической интерпретации этих решений он открыл область, над которой работали различные авторы, в частности Казорати и Кэли . Последнему принадлежит (1872) теория сингулярных решений дифференциальных уравнений первого порядка, принятая около 1900 года.
Первобытная попытка иметь дело с дифференциальными уравнениями имела в виду сведение к квадратурам . Как надежда алгебраистов восемнадцатого века была найти метод решения общего уравнения n -й степени, так и надежда аналитиков была найти общий метод интегрирования любого дифференциального уравнения. Гаусс (1799) показал, однако, что сложные дифференциальные уравнения требуют комплексных чисел . Поэтому аналитики начали заменять изучение функций, тем самым открывая новое и плодородное поле. Коши был первым, кто оценил важность этой точки зрения. После этого реальный вопрос больше не заключался в том, возможно ли решение с помощью известных функций или их интегралов, а в том, достаточно ли данного дифференциального уравнения для определения функции независимой переменной или переменных, и, если да, каковы характеристические свойства.
Два мемуара Фукса [25] вдохновили на новый подход, впоследствии разработанный Томе и Фробениусом . Колле был видным деятелем, начиная с 1869 года. Его метод интегрирования нелинейной системы был передан Бертрану в 1868 году. Клебш (1873) атаковал теорию по линиям, параллельным тем, что были в его теории абелевых интегралов . Поскольку последняя может быть классифицирована в соответствии со свойствами фундаментальной кривой, которая остается неизменной при рациональном преобразовании, Клебш предложил классифицировать трансцендентные функции, определяемые дифференциальными уравнениями, в соответствии с инвариантными свойствами соответствующих поверхностей f = 0 при рациональных взаимно-однозначных преобразованиях.
Начиная с 1870 года, работа Софуса Ли поставила теорию дифференциальных уравнений на более прочную основу. Он показал, что теории интегрирования математиков старшего поколения могут, используя группы Ли , быть отнесены к общему источнику, и что обыкновенные дифференциальные уравнения, которые допускают те же самые бесконечно малые преобразования, представляют сопоставимые трудности интегрирования. Он также подчеркнул тему преобразований контакта .
Теория групп Ли дифференциальных уравнений была сертифицирована, а именно: (1) что она объединяет множество известных методов ad hoc для решения дифференциальных уравнений, и (2) что она предоставляет мощные новые способы поиска решений. Теория применима как к обычным, так и к частным дифференциальным уравнениям. [26]
Общий подход к решению использует свойство симметрии дифференциальных уравнений, непрерывные бесконечно малые преобразования решений в решения ( теория Ли ). Непрерывная теория групп , алгебры Ли и дифференциальная геометрия используются для понимания структуры линейных и нелинейных (частных) дифференциальных уравнений для генерации интегрируемых уравнений, для нахождения их пар Лакса , операторов рекурсии, преобразования Бэклунда и, наконец, для нахождения точных аналитических решений ДУ.
Методы симметрии применяются к дифференциальным уравнениям, возникающим в математике, физике, технике и других дисциплинах.
Теория Штурма–Лиувилля — это теория специального типа линейных обыкновенных дифференциальных уравнений второго порядка. Их решения основаны на собственных значениях и соответствующих собственных функциях линейных операторов, определенных с помощью однородных линейных уравнений второго порядка . Задачи определяются как задачи Штурма–Лиувилля (SLP) и названы в честь J. C. F. Sturm и J. Liouville , которые изучали их в середине 1800-х годов. SLP имеют бесконечное число собственных значений, а соответствующие собственные функции образуют полный ортогональный набор, что делает возможными ортогональные разложения. Это ключевая идея в прикладной математике, физике и технике. [27] SLP также полезны при анализе некоторых уравнений с частными производными.
Существует несколько теорем, которые устанавливают существование и единственность решений задач начального значения, включающих ОДУ как локально, так и глобально. Две основные теоремы:
В своей базовой форме обе эти теоремы гарантируют только локальные результаты, хотя последняя может быть расширена для получения глобального результата, например, если выполнены условия неравенства Грёнвалля .
Кроме того, теоремы единственности, подобные теореме Липшица, приведенной выше, не применимы к системам ДАУ , которые могут иметь несколько решений, вытекающих только из их (нелинейной) алгебраической части. [28]
Теорему можно сформулировать просто следующим образом. [29] Для уравнения и задачи начального значения: если F и ∂ F /∂ y непрерывны в замкнутом прямоугольнике в плоскости xy , где a и b являются действительными (символически: a , b ∈ R ) и × обозначает декартово произведение , квадратные скобки обозначают замкнутые интервалы , то существует интервал для некоторого h ∈ R , где может быть найдено решение вышеуказанного уравнения и задачи начального значения. То есть, существует решение, и оно единственно. Поскольку нет ограничения на то, что F должна быть линейной, это применимо к нелинейным уравнениям, которые принимают форму F ( x , y ), и это также может быть применено к системам уравнений.
Когда гипотезы теоремы Пикара–Линделёфа удовлетворены, тогда локальное существование и единственность могут быть расширены до глобального результата. Точнее: [30]
Для каждого начального условия ( x 0 , y 0 ) существует уникальный максимальный (возможно бесконечный) открытый интервал
таким образом, что любое решение, удовлетворяющее этому начальному условию, является ограничением решения, удовлетворяющего этому начальному условию с областью определения .
В случае, если , есть ровно две возможности
где Ω — открытое множество, в котором определено F , а — его граница.
Обратите внимание, что максимальная область решения
Это означает, что F ( x, y ) = y 2 , что соответствует C 1 и, следовательно, локально Липшицеву функтору, удовлетворяющему теореме Пикара–Линделёфа.
Даже в такой простой ситуации максимальная область решения не может быть всей, поскольку решение
который имеет максимальный домен:
Это ясно показывает, что максимальный интервал может зависеть от начальных условий. Область y можно было бы взять как есть, но это привело бы к области, которая не является интервалом, так что сторона, противоположная начальному условию, была бы отсоединена от начального условия и, следовательно, не определялась бы им однозначно.
Максимальный домен не потому, что
что является одним из двух возможных случаев согласно приведенной выше теореме.
Дифференциальные уравнения обычно легче решать, если можно понизить порядок уравнения.
Любое явное дифференциальное уравнение порядка n ,
можно записать как систему n дифференциальных уравнений первого порядка, определив новое семейство неизвестных функций
для i = 1, 2, ..., n . Тогда n- мерная система связанных дифференциальных уравнений первого порядка имеет вид
более компактно в векторной записи:
где
Некоторые дифференциальные уравнения имеют решения, которые можно записать в точной и замкнутой форме. Здесь приведены несколько важных классов.
В таблице ниже P ( x ) , Q ( x ) , P ( y ) , Q ( y ) и M ( x , y ) , N ( x , y ) — любые интегрируемые функции x , y ; b и c — действительные заданные константы; C 1 , C 2 , ... — произвольные константы ( в общем случае комплексные ). Дифференциальные уравнения находятся в их эквивалентных и альтернативных формах, которые приводят к решению посредством интегрирования.
В интегральных решениях λ и ε являются фиктивными переменными интегрирования (континуальными аналогами индексов при суммировании ), а запись ∫ x F ( λ ) dλ просто означает интегрирование F ( λ ) по λ , а затем после интегрирования подстановку λ = x , без добавления констант (что явно указано).
Когда все другие методы решения ОДУ терпят неудачу или в случаях, когда у нас есть некоторая интуиция относительно того, как может выглядеть решение ДУ, иногда можно решить ДУ, просто угадывая решение и проверяя его правильность. Чтобы использовать этот метод, мы просто угадываем решение дифференциального уравнения, а затем подставляем решение в дифференциальное уравнение, чтобы проверить, удовлетворяет ли оно уравнению. Если это так, то у нас есть частное решение ДУ, в противном случае мы начинаем снова и пробуем еще одну догадку. Например, мы могли бы угадать, что решение ДУ имеет вид: поскольку это очень распространенное решение, которое физически ведет себя синусоидальным образом.
В случае неоднородного ОДУ первого порядка нам нужно сначала найти решение однородной части ОДУ, иначе известной как связанное однородное уравнение, а затем найти решение всего неоднородного уравнения методом угадывания. Наконец, мы складываем оба этих решения вместе, чтобы получить общее решение ОДУ, то есть: