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

редактировать
Дифференциальное уравнение, содержащее одну или несколько функций одной независимой переменной и ее производных

В математике, обыкновенное дифференциальное уравнение (ODE ) - это дифференциальное уравнение, содержащее одну или несколько функций от одной независимой переменной и производные этих функций. Термин обыкновенный используется в отличие от термина уравнение в частных производных, которое может относиться к более чем одной независимой переменной.

Содержание

  • 1 Дифференциальные уравнения
  • 2 Предпосылки
  • 3 Определения
    • 3.1 Общее определение
    • 3.2 Система ОДУ
    • 3.3 Решения
  • 4 Теории
    • 4.1 Сингулярные решения
    • 4.2 Приведение к квадратурам
    • 4.3 Фуксова теория
    • 4.4 Теория Ли
    • 4.5 Теория Штурма – Лиувилля
  • 5 Существование и единственность решений
    • 5.1 Локальная теорема существования и единственности упрощенная
    • 5.2 Глобальная единственность и максимальная область решения
  • 6 Снижение порядка
    • 6.1 Сведение к система первого порядка
  • 7 Сводка точных решений
  • 8 Метод предположений
  • 9 Программное обеспечение для решения ODE
  • 10 См. также
  • 11 Примечания
  • 12 Ссылки
  • 13 Библиография
  • 14 Внешние ссылки

Дифференциальные уравнения

A линейное дифференциальное уравнение - это дифференциальное уравнение, которое определяется линейным полиномом в неизвестная функция и ее производные, то есть уравнение вида

a 0 (x) y + a 1 (x) y ′ + a 2 (x) y ″ + ⋯ + an ( Икс) Y (N) + б (Икс) знак равно 0, {\ Displaystyle a_ {0} (x) Y + a_ {1} (x) y '+ a_ {2} (x) y' '+ \ cdots + a_ {n} (x) y ^ {(n)} + b (x) = 0,}{\displaystyle a_{0}(x)y+a_{1}(x)y'+a_{2}(x)y''+\cdots +a_{n}(x)y^{(n)}+b(x)=0,}

где a 0 (x) {\ displaystyle a_ {0} (x)}{\ displaystyle a_ {0} (x)} ,..., an (x) {\ displaystyle a_ {n} (x)}{\ displaystyle a_ {n} ( x)} и b (x) {\ displaystyle b (x)}b(x)- произвольные дифференцируемые функции, которые не обязательно должны быть линейными, и y ′,…, y (n) {\ displaystyle y ', \ ldots, y ^ {(n)} }{\displaystyle y',\ldots,y^{(n)}}- последовательные производные неизвестной функции y от переменной x.

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

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

Предпосылки

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

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

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

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

Простым примером является второй закон Ньютона - взаимосвязь между смещением x и время t объекта под действием силы F задаются дифференциальным уравнением

md 2 x (t) dt 2 = F (x (t)) {\ displaystyle m {\ frac {\ mathrm {d } ^ {2} x (t)} {\ mathrm {d} t ^ {2}}} = F (x (t)) \,}{\ displaystyle m {\ frac {\ mathrm {d} ^ {2} x (t)} {\ mathrm {d} t ^ {2}}} = F (x (t)) \,}

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

Определения

Далее пусть y будет зависимая переменная и x независимая переменная, а y = f (x) - неизвестная функция от x. Обозначения для различения различаются в зависимости от автора и от того, какие обозначения наиболее полезны для решения поставленной задачи. В этом контексте нотация Лейбница (dy / dx, dy / dx,..., dy / dx) более полезна для дифференциации и интегрирования, тогда как нотация Лагранжа (y ′, y ′ ′,..., y) более полезно для компактного представления производных любого порядка, а ньютоновская нотация (y ˙, y ¨, y..) {\ displaystyle ({\ dot {y}}, {\ ddot {y}}, {\ overset {...} {y}})}{\ displaystyle ({\ dot {y}}, {\ ddot {y}}, {\ overset {...} {y}})} часто используется в физике для представления производных низкого порядка по времени.

Общее определение

Дано F, функция от x, y и производные от y. Тогда уравнение вида

F (x, y, y ',…, y (n - 1)) = y (n) {\ displaystyle F \ left (x, y, y', \ ldots, y ^ {(n-1)} \ right) = y ^ {(n)}}{\displaystyle F\left(x,y,y',\ldots,y^{(n-1)}\right)=y^{(n)}}

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

В более общем смысле, неявное обыкновенное дифференциальное уравнение порядка n имеет вид:

F (x, y, y ′, y ″,…, y (n)) = 0 {\ displaystyle F \ left (x, y, y ', y' ', \ \ ldots, \ y ^ {(n)} \ right) = 0}{\displaystyle F\left(x,y,y',y'',\ \ldots,\ y^{(n)}\right)=0}

Существуют другие классификации:

Автономный
Дифференциальное уравнение, не зависящее от x, является называется автономным.
Линейным
Дифференциальное уравнение называется линейным, если F может быть записано как линейная комбинация производных y:
y (n) = ∑ я знак равно 0 N - 1 ai (x) y (i) + r (x) {\ displaystyle y ^ {(n)} = \ sum _ {i = 0} ^ {n-1} a_ {i} (x) y ^ {(i)} + r (x)}{\ displaystyle y ^ {(n)} = \ sum _ {i = 0} ^ {n-1} a_ {i} (x) y ^ {(i)} + r (x)}
где a i (x) и r (x) являются непрерывными функциями от x. Функция r (x) называется исходным членом, что приводит к еще двум важным классификациям:
Однородный
Если r (x) = 0, и, следовательно, одно «автоматическое» решение является тривиальным решением, y = 0. Решением линейного однородного уравнения является дополнительная функция, обозначаемая здесь y c.
Неоднородная (или неоднородная)
Если r (x) ≠ 0. Дополнительное решение к дополнительной функции является частный интеграл, обозначенный здесь y p.

Общее решение линейного уравнения может быть записано как y = y c + y p.

Non- linear
Дифференциальное уравнение, которое не может быть записано в форме линейной комбинации.

Система ОДУ

Ряд связанных дифференциальных уравнений образуют систему уравнений. Если y - вектор, элементы которого являются функциями; y (x) = [y 1 (x), y 2 (x),..., y m (x)], и F является векторнозначной функцией от y и его производных, тогда

y (n) = F (x, y, y ′, y ″,…, y (n - 1)) {\ displaystyle \ mathbf {y} ^ {(n)} = \ mathbf {F} \ left (x, \ mathbf {y}, \ mathbf {y } ', \ mathbf {y}' ', \ ldots, \ mathbf {y} ^ {(n-1)} \ right)}{\displaystyle \mathbf {y} ^{(n)}=\mathbf {F} \left(x,\mathbf {y},\mathbf {y} ',\mathbf {y} '',\ldots,\mathbf {y} ^{(n-1)}\right)}

- явная система обыкновенных дифференциальных уравнений порядка n и размерности m. В формате вектора-столбца :

(y 1 (n) y 2 (n) ⋮ ym (n)) = (f 1 (x, y, y ′, y ″,…, y ( n - 1)) f 2 (x, y, y ′, y ″,…, y (n - 1)) ⋮ fm (x, y, y ′, y ″,…, y (n - 1))) {\ displaystyle {\ begin {pmatrix} y_ {1} ^ {(n)} \\ y_ {2} ^ {(n)} \\\ vdots \\ y_ {m} ^ {(n)} \ end { pmatrix}} = {\ begin {pmatrix} f_ {1} \ left (x, \ mathbf {y}, \ mathbf {y} ', \ mathbf {y}' ', \ ldots, \ mathbf {y} ^ { (n-1)} \ right) \\ f_ {2} \ left (x, \ mathbf {y}, \ mathbf {y} ', \ mathbf {y}' ', \ ldots, \ mathbf {y} ^ {(n-1)} \ right) \\\ vdots \\ f_ {m} \ left (x, \ mathbf {y}, \ mathbf {y} ', \ mathbf {y}' ', \ ldots, \ mathbf {y} ^ {(n-1)} \ right) \ end {pmatrix}}}{\displaystyle {\begin{pmatrix}y_{1}^{(n)}\\y_{2}^{(n)}\\\vdots \\y_{m}^{(n)}\end{pmatrix}}={\begin{pmatrix}f_{1}\left(x,\mathbf {y},\mathbf {y} ',\mathbf {y} '',\ldots,\mathbf {y} ^{(n-1)}\right)\\f_{2}\left(x,\mathbf {y},\mathbf {y} ',\mathbf {y} '',\ldots,\mathbf {y} ^{(n-1)}\right)\\\vdots \\f_{m}\left(x,\mathbf {y},\mathbf {y} ',\mathbf {y} '',\ldots,\mathbf {y} ^{(n-1)}\right)\end{pmatrix}}}

Они не обязательно линейны. Неявный аналог:

F (x, y, y ', y ″,…, y (n)) = 0 {\ displaystyle \ mathbf {F} \ left (x, \ mathbf {y}, \ mathbf {y} ', \ mathbf {y}' ', \ ldots, \ mathbf {y} ^ {(n)} \ right) = {\ boldsymbol {0}}}{\displaystyle \mathbf {F} \left(x,\mathbf {y},\mathbf {y} ',\mathbf {y} '',\ldots,\mathbf {y} ^{(n)}\right)={\boldsymbol {0}}}

где 0 = (0, 0,..., 0) - нулевой вектор . В матричной форме

(f 1 (x, y, y ′, y ″,…, y (n)) f 2 (x, y, y ′, y ″,…, y (n)) ⋮ fm ( x, y, y ′, y ″,…, y (n))) = (0 0 ⋮ 0) {\ displaystyle {\ begin {pmatrix} f_ {1} (x, \ mathbf {y}, \ mathbf { y} ', \ mathbf {y}' ', \ ldots, \ mathbf {y} ^ {(n)}) \\ f_ {2} (x, \ mathbf {y}, \ mathbf {y}', \ mathbf {y} '', \ ldots, \ mathbf {y} ^ {(n)}) \\\ vdots \\ f_ {m} (x, \ mathbf {y}, \ mathbf {y} ', \ mathbf {y} '', \ ldots, \ mathbf {y} ^ {(n)}) \ end {pmatrix}} = {\ begin {pmatrix} 0 \\ 0 \\\ vdots \\ 0 \ end {pmatrix} }}{\displaystyle {\begin{pmatrix}f_{1}(x,\mathbf {y},\mathbf {y} ',\mathbf {y} '',\ldots,\mathbf {y} ^{(n)})\\f_{2}(x,\mathbf {y},\mathbf {y} ',\mathbf {y} '',\ldots,\mathbf {y} ^{(n)})\\\vdots \\f_{m}(x,\mathbf {y},\mathbf {y} ',\mathbf {y} '',\ldots,\mathbf {y} ^{(n)})\end{pmatrix}}={\begin{pmatrix}0\\0\\\vdots \\0\end{pmatrix}}}

Для системы вида F (x, y, y ′) = 0 {\ displaystyle \ mathbf {F} \ left (x, \ mathbf {y}, \ mathbf {y} ' \ right) = {\ boldsymbol {0}}}\mathbf {F} \left(x,\mathbf {y},\mathbf {y} '\right)={\boldsymbol {0}}, некоторые источники также требуют, чтобы матрица Якоби ∂ F (x, u, v) ∂ v {\ displaystyle {\ frac {\ partial \ mathbf {F} (x, \ mathbf {u}, \ mathbf {v})} {\ partial \ mathbf {v}}}}{\ frac {\ partial \ mathbf {F} (x, \ mathbf {u}, \ mathbf {v})} {\ partial \ mathbf {v}}} быть не- единственное число, чтобы называть это неявным ODE [системой]; неявная система ОДУ, удовлетворяющая этому условию неособенности Якоби, может быть преобразована в явную систему ОДУ. В тех же источниках неявные системы ОДУ с сингулярным якобианом называются дифференциально-алгебраическими уравнениями (ДАУ). Это различие не просто терминологическое; DAE имеют принципиально разные характеристики и, как правило, более сложные для решения, чем (неособые) системы ODE. Предположительно для дополнительных производных матрица Гессе и т. Д. Также считаются невырожденными в соответствии с этой схемой, хотя обратите внимание, что любое ОДУ порядка выше единицы может быть [и обычно] переписано как система ОДУ первого порядка, что делает критерий якобианской сингулярности достаточным для того, чтобы эта таксономия была всеобъемлющей на всех порядках.

Поведение системы ОДУ можно визуализировать с помощью фазового портрета.

Решения

Заданное дифференциальное уравнение

F (x, y, y ′,…, Y (n)) = 0 {\ displaystyle F \ left (x, y, y ', \ ldots, y ^ {(n)} \ right) = 0}{\displaystyle F\left(x,y,y',\ldots,y^{(n)}\right)=0}

функция u: I ⊂ R→ R, где I - интервал, называется решением или интегральной кривой для F, если u n-кратно дифференцируемо на I, и

F (x, u, u ′,…, u (n)) = 0 x ∈ I. {\ Displaystyle F (x, u, u ', \ \ ldots, \ u ^ {(n)}) = 0 \ quad x \ in I.}{\displaystyle F(x,u,u',\ \ldots,\ u^{(n)})=0\quad x\in I.}

Даны два решения u: J ⊂ R→ Rи v : I ⊂ R→ R, u называется расширением v, если I ⊂ J и

u (x) = v (x) x ∈ I. {\ displaystyle u (x) = v (x) \ quad x \ in I. \,}u (x) = v (x) \ quad x \ in I. \,

Решение, не имеющее расширения, называется максимальным решением. Решение, определенное на всем R, называется глобальным решением.

Общее решение уравнения n-го порядка - это решение, содержащее n произвольных независимых постоянных интегрирования. Конкретное решение выводится из общего решения путем установки констант на конкретные значения, часто выбираемые для выполнения набора «начальных условий или граничных условий ». сингулярное решение - это решение, которое не может быть получено путем присвоения определенных значений произвольным константам в общем решении.

Теории

Особые решения

Теория сингулярных решений обыкновенных и дифференциальных уравнений с частными производными была предметом исследований со времен Лейбница, но только с середины XIX века ей стало уделяться особое внимание. Ценная, но малоизвестная работа по этому поводу - это работа Хаутена (1854 г.). Дарбу (с 1873 г.) был лидером в теории, и в геометрической интерпретации этих решений он открыл область, в которой работали различные писатели, в частности Казорати и Кэли. Последним обязана (1872 г.) теория сингулярных решений дифференциальных уравнений первого порядка, принятая примерно в 1900 г.

Сведение к квадратурам

Первобытная попытка работы с дифференциальными уравнениями имела место в просмотреть сокращение до квадратур. Как алгебраисты восемнадцатого века надеялись найти метод решения общего уравнения n-й степени, так и аналитики надеялись найти общий метод интегрирования любого дифференциального уравнения. Гаусс (1799) показал, однако, что сложные дифференциальные уравнения требуют комплексных чисел. Таким образом, аналитики начали заменять изучение функций, открывая тем самым новое плодородное поле. Коши был первым, кто оценил важность этой точки зрения. После этого реальный вопрос был уже не в том, возможно ли решение с помощью известных функций или их интегралов, а в том, достаточно ли данного дифференциального уравнения для определения функции независимой переменной или переменных, и, если да, то каковы характерные свойства.

Теория Фукса

Два мемуара Фукса вдохновили на новый подход, который впоследствии был разработан Томе и Фробениусом. Колле внес значительный вклад, начиная с 1869 года. Его метод интеграции нелинейной системы был передан Бертрану в 1868 году. Клебш (1873) атаковал теорию по направлениям, параллельным тем, что в его теории Абелевы интегралы. Поскольку последние могут быть классифицированы в соответствии со свойствами фундаментальной кривой, которая остается неизменной при рациональном преобразовании, Клебш предложил классифицировать трансцендентные функции, определяемые дифференциальными уравнениями, в соответствии с инвариантными свойствами соответствующих поверхностей f = 0 как рациональные. -один трансформации.

Теория Ли

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

Теория групп Ли дифференциальных уравнений была сертифицирована, а именно: (1) она объединяет многие специальные методы, известные для решения дифференциальных уравнений, и (2) что она предоставляет новые мощные способы поиска решений. Теория имеет приложения как к обыкновенным, так и к дифференциальным уравнениям с частными производными.

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

Методы симметрии применялись к дифференциальным уравнениям, возникающим в математике, физике, технике и других дисциплинах.

Теория Штурма – Лиувилля

Теория Штурма – Лиувилля - это теория особого типа линейных обыкновенных дифференциальных уравнений второго порядка. Их решения основаны на собственных значениях и соответствующих собственных функциях линейных операторов, определенных через однородные линейные уравнения второго порядка . Проблемы обозначены как проблемы Штурма-Лиувилля (SLP) и названы в честь J.C.F. Штурм и Дж. Лиувилля, изучавший их в середине 1800-х гг. SLP имеют бесконечное количество собственных значений, а соответствующие собственные функции образуют полный ортогональный набор, что делает возможным ортогональное расширение. Это ключевая идея прикладной математики, физики и инженерии. SLP также полезны при анализе некоторых дифференциальных уравнений в частных производных.

Существование и уникальность решений

Существует несколько теорем, которые устанавливают существование и уникальность решений проблем с начальным значением, включающих ОДУ как локально, так и глобально. Две основные теоремы:

ТеоремаПредположениеЗаключение
Теорема существования Пеано F непрерывное только локальное существование
Теорема Пикара – Линделёфа F Липшицева непрерывное локальное существование и единственность

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

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

Локальная теорема существования и единственности упрощена

Теорема может быть сформулирована просто следующим образом. Для уравнения и задачи начального значения:

y ′ = F (x, y), y 0 = y (x 0) {\ displaystyle y '= F (x, y) \,, \ quad y_ {0} = y (x_ {0})}y'=F(x,y)\,,\quad y_{0}=y(x_{0})

, если F и ∂F / ∂y непрерывны в замкнутом прямоугольнике

R = [x 0 - a, x 0 + a] × [y 0 - b, y 0 + b] {\ displaystyle R = [x_ {0} -a, x_ {0} + a] \ times [y_ {0} -b, y_ {0} + b]}R = [x_ {0} -a, x_ {0} + a] \ times [y_ {0} -b, y_ {0} + b]

в плоскости xy, где a и b являются вещественными (символически: a, b ∈ ℝ), а × обозначает декартово произведение, квадратные скобки обозначают отрезки, тогда существует интервал

I = [x 0 - h, x 0 + h] ⊂ [x 0 - a, x 0 + a] {\ displaystyle I = [x_ {0} -h, x_ {0} + h] \ подмножество [x_ {0} -a, x_ {0} + a]}I = [x_ {0} -h, x_ {0} + h] \ subset [x_ {0} -a, x_ {0} + a]

для некоторого h ∈ ℝ, где может быть найдено решение вышеуказанного уравнения и задачи начального значения. То есть решение есть, и оно уникальное. Поскольку нет ограничений на линейность F, это применимо к нелинейным уравнениям, которые принимают форму F (x, y), а также могут применяться к системам уравнений.

Глобальная уникальность и максимальная область решения

Когда выполняются гипотезы теоремы Пикара – Линделёфа, тогда локальное существование и уникальность может быть расширено до глобального результата. Точнее:

Для каждого начального условия (x 0, y 0) существует уникальный максимальный (возможно бесконечный) открытый интервал

I max = ( x -, x +), x ± ∈ R ∈ {± ∞}, x 0 ∈ I max {\ displaystyle I _ {\ max} = (x _ {-}, x _ {+}), x _ {\ pm} \ in \ mathbb {R} \ cup \ {\ pm \ infty \}, x_ {0} \ in I _ {\ max}}{\ displaystyle I _ {\ max} = (x_ { -}, x _ {+}), x _ {\ pm} \ in \ mathbb {R} \ cup \ {\ pm \ infty \}, x_ {0} \ in I _ {\ max}}

такое, что любое решение, удовлетворяющее этому начальному условию, является ограничением решение, которое удовлетворяет этому начальному условию с областью I max {\ displaystyle I _ {\ max}}I _ {\ max} .

В случае, если x ± ≠ ± ∞ {\ displaystyle x _ {\ pm} \ neq \ pm \ infty}{\ displaystyle x _ {\ pm} \ neq \ pm \ infty} , есть ровно две возможности

  • взрыва за конечное время: lim sup x → x ± ‖ y (x) ‖ → ∞ {\ displaystyle \ limsup _ {x \ to x _ {\ pm}} \ | y (x) \ | \ to \ infty}{\ displaystyle \ limsup _ {x \ to x _ {\ pm}} \ | y ( х) \ | \ к \ infty}
  • выходит из области определения: lim x → x ± y (x) ∈ ∂ Ω ¯ {\ displaystyle \ lim _ { x \ to x _ {\ pm}} y (x) \ \ in \ partial {\ bar {\ Omega}}}\ lim _ {x \ to x _ {\ pm}} y (x) \ \ in \ partial {\ bar {\ Omega}}

где Ω - открытое множество, в котором определено F, и ∂ Ω ¯ { \ Displaystyle \ partial {\ bar {\ Omega }}}\ partial {\ bar {\ Omega}} - его граница.

Обратите внимание, что максимальная область решения

  • всегда представляет собой интервал (чтобы иметь уникальность)
  • может быть меньше, чем R {\ displaystyle \ mathbb {R}}\ mathbb {R}
  • может зависеть от конкретного выбора (x 0, y 0).
Пример.
y ′ = y 2 {\ displaystyle y '= y ^ {2}}y'=y^{2}

Это означает что F (x, y) = y, который является C и, следовательно, локально липшицевым, удовлетворяет теореме Пикара – Линделёфа.

Даже в такой простой ситуации максимальная область решения не может быть полностью R {\ displaystyle \ mathbb {R}}\ mathbb {R} , поскольку решение равно

y (x) = y 0 (x 0 - x) y 0 + 1 {\ displaystyle y (x) = {\ гидроразрыв {y_ {0}} {(x_ {0} -x) y_ {0} +1}}}y (x) = {\ frac {y_ {0}} {(x_ {0} -x) y_ {0} +1}}

с максимальным доменом:

{R y 0 = 0 (- ∞, x 0 + 1 y 0) y 0>0 (x 0 + 1 y 0, + ∞) y 0 < 0 {\displaystyle {\begin{cases}\mathbb {R} y_{0}=0\\[4pt]\left(-\infty,x_{0}+{\frac {1}{y_{0}}}\right)y_{0}>0 \\ [4pt] \ left (x_ {0} + {\ frac {1} {y_ {0}}}, + \ infty \ right) y_ {0} <0\end{cases}}}{\displaystyle {\begin{cases}\mathbb {R} y_{0}=0\\[4pt]\left(-\infty,x_{0}+{\frac {1}{y_{0}}}\right)y_{0}>0 \\ [4pt] \ left (x_ {0} + {\ frac {1} {y_ {0}}}, + \ infty \ right) y_ { 0} <0\end{cases}}}

Это ясно показывает, что максимальный интервал может зависеть от начальных условий. Область y может быть принято как R ∖ (x 0 + 1 / y 0), {\ displaystyle \ mathbb {R} \ setminus (x_ {0} + 1 / y_ {0}),}{\ displaystyle \ mathbb {R} \ setminus (x_ {0} + 1 / y_ {0}),} но это привело бы к области, которая не является интервалом, так что сторона, противоположная начальному условию, будет отключена от начального условия и, следовательно, не будет однозначно определена им.

Максимальный домен не R {\ displaystyle \ mathbb {R}}\ mathbb {R} , потому что

lim x → x ± ‖ y (x) ‖ → ∞, {\ displaystyle \ lim _ {x \ to x _ {\ pm}} \ | y (x) \ | \ to \ infty,}{\ displaystyle \ lim _ {x \ to x _ {\ pm}} \ | y (x) \ | \ to \ infty,}

, который является одним из двух возможных случаев согласно приведенной выше теореме.

Понижение порядка

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

Редукция к системе первого порядка

Любое явное дифференциальное уравнение порядка n,

F (x, y, y ', y ″,…, y (n - 1)) = Y (N) {\ Displaystyle F \ left (x, y, y ', y' ', \ \ ldots, \ y ^ {(n-1)} \ right) = y ^ {(n)}}{\displaystyle F\left(x,y,y',y'',\ \ldots,\ y^{(n-1)}\right)=y^{(n)}}

можно записать в виде системы n дифференциальных уравнений первого порядка путем определения нового семейства неизвестных функций

yi = y (i - 1). {\ displaystyle y_ {i} = y ^ {(i-1)}. \!}y_ {i} = y ^ {(i-1)}. \!

для i = 1, 2,..., n. Тогда n-мерная система связанных дифференциальных уравнений первого порядка имеет вид

y 1 ′ = y 2 y 2 ′ = y 3 ⋮ y n - 1 ′ = y n y n ′ = F (x, y 1,…, y n). {\ displaystyle {\ begin {array} {rcl} y_ {1} '= y_ {2} \\ y_ {2}' = y_ {3} \\ \ vdots \\ y_ {n-1} '= y_ {n} \\ y_ {n}' = F (x, y_ {1}, \ ldots, y_ {n}). \ end {array}}}{\displaystyle {\begin{array}{rcl}y_{1}'=y_{2}\\y_{2}'=y_{3}\\\vdots \\y_{n-1}'=y_{n}\\y_{n}'=F(x,y_{1},\ldots,y_{n}).\end{array}}}

более компактно в векторной записи:

y ′ = F (x, y) {\ displaystyle \ mathbf {y} '= \ mathbf {F} (x, \ mathbf {y})}\mathbf {y} '=\mathbf {F} (x,\mathbf {y})

где

y = (y 1, …, Yn), F (x, y 1,…, yn) = (y 2,…, yn, F (x, y 1,…, yn)). {\ displaystyle \ mathbf {y} = (y_ {1}, \ ldots, y_ {n}), \ quad \ mathbf {F} (x, y_ {1}, \ ldots, y_ {n}) = (y_ {2}, \ ldots, y_ {n}, F (x, y_ {1}, \ ldots, y_ {n})).}{\ displaystyle \ mathbf {y} = (y_ {1}, \ ldots, y_ {n}), \ quad \ mathbf {F} (x, y_ {1}, \ ldots, y_ {n}) = (y_ {2}, \ ldots, y_ {n}, F (x, y_ { 1}, \ ldots, y_ {n})).}

Сводка точных решений

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

В таблице ниже P (x), Q (x), P (y), Q (y) и M (x, y), N (x, y) - любые интегрируемые функции от x, y, а также b и c являются действительными заданными константами, а C 1, C 2,... являются произвольными константами (complex в целом). Дифференциальные уравнения представлены в их эквивалентной и альтернативной формах, которые приводят к решению путем интегрирования.

В интегральных решениях λ и ε - фиктивные переменные интегрирования (континуальные аналоги индексов в суммировании ), а обозначение ∫F (λ) dλ просто означает интегрирование F ( λ) по λ, то после интегрирования подставляем λ = x без добавления констант (явно указанных).

ТипДифференциальное уравнениеМетод решенияОбщее решение
РазделимоеПервого порядка, разделимое по x и y (общий случай, специальные случаи см. ниже)

P 1 (x) Q 1 (y) + P 2 (x) Q 2 (y) dydx = 0 P 1 (x) Q 1 (y) dx + P 2 (Икс) Q 2 (Y) dy знак равно 0 {\ Displaystyle {\ begin {align} P_ {1} (x) Q_ {1} (y) + P_ {2} (x) Q_ {2} (y) \, {\ frac {dy} {dx}} = 0 \\ P_ {1} (x) Q_ {1} (y) \, dx + P_ {2} (x) Q_ {2} (y) \, dy = 0 \ end {align}}}{\ displaystyle {\ begin {align} P_ {1} (x) Q_ {1} (y) + P_ {2} (x) Q_ {2} (y) \, {\ frac {dy} {dx}} = 0 \\ P_ {1} (x) Q_ {1} (y) \, dx + P_ {2} (x) Q_ {2} (у) \, dy = 0 \ конец {выровнено}}}

Разделение переменных (разделить на P 2Q1).∫ Икс п 1 (λ) п 2 (λ) d λ + ∫ Y Q 2 (λ) Q 1 (λ) d λ знак равно С {\ Displaystyle \ int ^ {x} {\ гидроразрыва {P_ {1 } (\ lambda)} {P_ {2} (\ lambda)}} \, d \ lambda + \ int ^ {y} {\ frac {Q_ {2} (\ lambda)} {Q_ {1} (\ lambda)}} \, d \ lambda = C \, \!}\ int ^ {x} {\ frac {P_ {1} (\ lambda) } {P_ {2} (\ lambda)}} \, d \ lambda + \ int ^ {y} {\ frac {Q_ {2} (\ lambda)} {Q_ {1} (\ lambda)}} \, d \ лямбда = C \, \!
Первый порядок, разделимый по x

dydx = F (x) dy = F (x) dx {\ displaystyle {\ begin {align} {\ frac {dy} {dx}} = F (x) \\ dy = F (x) \, dx \ end {align}}}{\ displaystyle {\ begin {align} {\ frac {dy} {dx}} = F (x) \\ dy = F (x) \, dx \ end {выровнено}}}

Прямое интегрирование.Y = ∫ Икс F (λ) d λ + C {\ Displaystyle y = \ int ^ {x} F (\ lambda) \, d \ lambda + C \, \!}y = \ int ^ {x} F (\ lambda) \, d \ lambda + C \, \!
Первого порядка, автономный, разделяемый по y

dydx = F (y) dy = F (y) dx {\ displaystyle {\ begin {align} {\ frac {dy} {dx}} = F (y) \ \ dy = F (y) \, dx \ end {align}}}{\ displaystyle {\ begin {align} {\ frac { dy} {dx}} = F (y) \\ dy = F (y) \, dx \ end {align}}}

Разделение переменных (разделить на F).Икс знак равно ∫ ярд λ F (λ) + C {\ displaystyle x = \ int ^ {y} {\ frac {d \ lambda} {F (\ lambda)}} + C \, \!}x = \ дюйм t ^ {y} {\ гидроразрыва {d \ лямбда} {F (\ lambda)}} + C \, \!
Первый порядок, разделяемый по x и y

P (y) dydx + Q (x) = 0 P (y) dy + Q (x) dx = 0 {\ displaystyle {\ begin {align} P (y) {\ frac {dy} {dx}} + Q (x) = 0 \\ P (y) \, dy + Q (x) \, dx = 0 \ end {align}}}{\ displaystyle {\ begin {align} P (y) {\ frac {dy} {dx}} + Q (x) знак равно 0 \\ P (y) \, dy + Q (x) \, dx = 0 \ end {align}}}

Интегрируйте во всем.∫ Y П (λ) d λ + ∫ Икс Q (λ) d λ знак равно С {\ Displaystyle \ int ^ {y} P (\ lambda) \, d \ lambda + \ int ^ {x} Q ( \ lambda) \, d \ lambda = C \, \!}{\ displaystyle \ int ^ {y} P (\ lambda) \, d \ lambda + \ int ^ {x} Q (\ lambda) \, d \ lambda = C \, \!}
Общие первого порядкаОднородные первого порядка

dydx = F (yx) {\ displaystyle {\ frac {dy} {dx}} = F \ left ({\ frac {y} {x}} \ right) \, \!}{\ displaystyle {\ frac {dy} {dx}} = F \ left ({\ frac {y } {x}} \ right) \, \!}

Установите y = ux, затем решите, разделив переменные в u и x.пер ⁡ (С Икс) знак равно ∫ Y / xd λ F (λ) - λ {\ displaystyle \ ln (Cx) = \ int ^ {y / x} {\ frac {d \ lambda} {F (\ lambda) - \ lambda}} \, \!}\ ln (Cx) = \ int ^ {y / x} {\ frac {d \ lambda} {F (\ lambda) - \ lambda}} \, \!
Разделимость первого порядка

y M (xy) + x N (xy) dydx = 0 y M (xy) dx + x N ( ху) dy знак равно 0 {\ displaystyle {\ begin {align} yM (xy) + xN (xy) \, {\ frac {dy} {dx}} = 0 \\ yM (xy) \, dx + xN ( xy) \, dy = 0 \ end {align}}}{\ displaystyle {\ begin {выровнено} yM (xy) + xN (xy) \, {\ frac {dy} {dx}} = 0 \\ yM (xy) \, dx + xN (xy) \, dy = 0 \ end {align}}}

Разделение переменных (разделить на xy).

пер ⁡ (С Икс) знак равно ∫ ху N (λ) d λ λ [N (λ) - M (λ)] {\ Displaystyle \ ln (Cx) = \ int ^ {xy} {\ гидроразрыва {N (\ lambda) \, d \ lambda} {\ lambda [N (\ lambda) -M (\ lambda)]}} \, \!}\ ln (Cx) = \ int ^ {xy} {\ frac {N (\ lambda) \, d \ lambda} {\ лямбда [N (\ lambda) -M (\ lambda)]}} \, \!

Если N = M, решением будет xy = C.

Точный дифференциал, первого порядка

M (x, y) dydx + N (x, y) = 0 M (x, y) dy + N (x, y) dx = 0 { \ Displaystyle {\ begin {align} M (x, y) {\ frac {dy} {dx}} + N (x, y) = 0 \\ M (x, y) \, dy + N (x, y) \, dx = 0 \ end {align}}}{\ displaystyle {\ begin {выровнено} M (x, y) {\ frac {dy} {dx}} + N (x, y) = 0 \\ M (x, y) \, dy + N (x, y) \, dx = 0 \ конец {выровнено}}}

где ∂ M ∂ x = ∂ N ∂ y {\ displaystyle {\ frac {\ partial M} {\ partial x}} = {\ frac {\ partial N} {\ partial y}} \, \!}{\ frac {\ partial M} {\ partial x}} = {\ frac {\ partial N} {\ partial y}} \, \!

Интегрировать всюду.F (x, y) знак равно ∫ Y M (x, λ) d λ + ∫ x N (λ, y) d λ + Y (y) + X (x) = C {\ displaystyle {\ begin { выровнено} F (x, y) = \ int ^ {y} M (x, \ lambda) \, d \ lambda + \ int ^ {x} N (\ lambda, y) \, d \ lambda \\ + Y (y) + X (x) = C \ end {align}}}{\ displaystyle {\ begin {align} F (x, y) = \ int ^ {y} M (x, \ lambda) \, d \ lambda + \ int ^ {x} N (\ lambda, y) \, d \ lambda \\ + Y (y) + X (x) = C \ end {align}}}

где Y (y) и X (x) - это функции от интегралов, а не постоянные значения, которые устанавливаются, чтобы сделать конечную функцию F (x, y) удовлетворяют исходному уравнению.

Неточный дифференциал, первый порядок

M (x, y) dydx + N (x, y) = 0 M (x, y) dy + N (x, y) dx = 0 {\ displaystyle {\ begin {align} M (x, y) {\ frac {dy} {dx}} + N (x, y) = 0 \\ M (x, y) \, dy + N ( x, y) \, dx = 0 \ end {align}}}{\ displaystyle {\ begin {выровнено} M (x, y) {\ frac {dy} {dx}} + N (x, y) = 0 \\ M (x, y) \, dy + N (x, y) \, dx = 0 \ конец {выровнено}}}

где ∂ M ∂ x ≠ ∂ N ∂ y {\ displaystyle {\ frac {\ partial M} {\ partial x}} \ neq {\ frac {\ partial N} {\ partial y}} \, \!}{\ frac {\ partial M} {\ partial x}} \ neq {\ frac {\ partial N} {\ partial y}} \, \!

Коэффициент интегрирования μ (x, y), удовлетворяющий

∂ (μ M) ∂ x = ∂ (μ N) ∂ Y {\ Displaystyle {\ frac {\ partial (\ mu M)} {\ partial x}} = {\ frac {\ partial (\ mu N)} {\ partial y}} \, \! }{\ frac {\ partial (\ mu M)} {\ partial x}} = {\ frac {\ partial (\ mu N)} {\ partial y}} \, \ !

Если можно найти μ (x, y):

F (x, y) = ∫ y μ (x, λ) M (x, λ) d λ + ∫ x μ (λ, y) N (λ, y) d λ + Y (y) + X (x) = C {\ displaystyle {\ begin {align} F (x, y) = \ int ^ {y} \ mu (x, \ lambda) M (x, \ lambda) \, d \ lambda + \ int ^ {x} \ mu (\ lambda, y) N (\ lambda, y) \, d \ lambda \\ + Y (y) + X (x) = C \ end {align}}}{\ displaystyle {\ begin {align} F (x, y) = \ int ^ {y} \ mu (x, \ lambda) M (x, \ lambda) \, d \ lambda + \ int ^ {x} \ mu (\ lambda, y) N (\ lambda, y) \, d \ лямбда \\ + Y (y) + X (x) = C \ end {align}}}

Обычный второй порядокВторой порядок, автономный

d 2 ydx 2 = F (y) {\ displaystyle {\ frac {d ^ {2} y} {dx ^ {2}}} = F (y) \, \!}{\ гидроразрыва {d ^ {2} y} {dx ^ {2}}} = F (y) \, \!

Умножьте обе части уравнения на 2dy / dx, замените 2 dydxd 2 ydx 2 = ddx (dydx) 2 {\ displaystyle 2 {\ frac {dy} {dx}} {\ frac {d ^ {2} y} {dx ^ {2}}} = {\ frac {d} {dx}} \ left ({\ frac {dy} {dx}} \ right) ^ {2} \, \!}2 {\ frac {dy} {dx}} {\ frac {d ^ {2} y} {dx ^ {2}}} = {\ frac {d} {dx}} \ left ({\ frac {dy} {dx}} \ right) ^ {2} \, \! , затем проинтегрируйте дважды.Икс = ± ∫ ярд λ 2 ∫ λ F (ε) d ε + C 1 + C 2 {\ displaystyle x = \ pm \ int ^ {y} {\ frac {d \ lambda} {\ sqrt {2 \ int ^ {\ lambda} F (\ varepsilon) \, d \ varepsilon + C_ {1}}}} + C_ {2} \, \!}{\ displaystyle x = \ pm \ int ^ {y} {\ frac {d \ lambda} {\ sqrt {2 \ int ^ {\ lambda} F (\ varepsilon) \, d \ varepsilon + C_ {1}}}} + C_ {2} \, \!}
Линейно до n-го порядкаПервый- упорядоченные, линейные, неоднородные, функциональные коэффициенты

dydx + P (x) y = Q (x) {\ displaystyle {\ frac {dy} {dx}} + P (x) y = Q (x) \, \!}{\ displaystyle {\ frac {dy} {dx}} + P (x) y = Q (x) \, \!}

Интегрирующий множитель: e ∫ x P (λ) d λ. {\ Displaystyle е ^ {\ int ^ {x} P (\ lambda) \, d \ lambda}.}{\ displaystyle e ^ {\ int ^ {x } P (\ lambda) \, d \ lambda}.} y = e - ∫ x P (λ) d λ [∫ xe ∫ λ P (ε) d ε Q (λ) d λ + C] {\ Displaystyle у = е ^ {- \ int ^ {x} P (\ lambda) \, d \ lambda} \ left [\ int ^ {x} e ^ {\ int ^ {\ lambda} P (\ varepsilon) \, d \ varepsilon} Q (\ lambda) \, d \ lambda + C \ right]}{\ displaystyle y = e ^ {- \ int ^ {x} P (\ lambda) \, d \ lambda} \ left [\ int ^ {x } e ^ {\ int ^ {\ lambda} P (\ varepsilon) \, d \ varepsilon} Q (\ lambda) \, d \ lambda + C \ right]}
Второй порядок, линейный, неоднородный, функциональные коэффициенты

d 2 ydx 2 + 2 p (x) dydx + (p (x) 2 + p '(x)) y = q (x) {\ displaystyle {d ^ {2} y \ over dx ^ {2}} + 2p (x) {dy \ over dx} + (p (x) ^ {2} + p '(x)) y = q (x)}{\displaystyle {d^{2}y \over dx^{2}}+2p(x){dy \over dx}+(p(x)^{2}+p'(x))y=q(x)}

Интегрирующий множитель: e ∫ x P (λ) d λ {\ displaystyle e ^ {\ int ^ {x} P (\ lambda) \, d \ lambda}}е ^ {\ int ^ {x} P (\ lambda) \, d \ lambda} y = e - ∫ x P (λ) d λ [∫ x (∫ ξ e ∫ λ п (ε) d ε Q (λ) d λ) d ξ + C 1 Икс + C 2] {\ Displaystyle y = e ^ {- \ int ^ {x} P (\ lambda) \, d \ lambda} \ left [\ int ^ {x} (\ int ^ {\ xi} e ^ {\ int ^ {\ lambda} P (\ varepsilon) \, d \ varepsilon} Q (\ lambda) \, d \ lambda) d \ xi + C_ {1} x + C_ {2} \ right]}{\ displaystyle y = e ^ {- \ int ^ {x} P ( \ lambda) \, d \ lambda} \ left [\ int ^ {x} (\ int ^ {\ xi} e ^ {\ int ^ {\ lambda} P (\ varepsilon) \, d \ varepsilon} Q (\ lambda) \, d \ lambda) d \ xi + C_ {1} x + C_ {2} \ right]}
Линейные, неоднородные, постоянные коэффициенты второго порядка

d 2 ydx 2 + bdydx + cy = r (x) { \ displaysty le {\ frac {d ^ {2} y} {dx ^ {2}}} + b {\ frac {dy} {dx}} + cy = r (x) \, \!}{\ frac {d ^ {2} y} {dx ^ {2}}} + b {\ frac {dy} {dx}} + cy = r (x) \, \!

Дополнительная функция y c : принять y c = e, подставить и решить многочлен от α, чтобы найти линейно независимые функции e α jx {\ displaystyle e ^ {\ alpha _ {j} x}}е ^ {\ альфа _ {j} x} .

Частный интеграл y p : в общем метод изменения параметров, хотя для очень простой проверки r (x) может работать

y = yc + yp {\ displaystyle y = y_ {c} + y_ {p}}{\ displaystyle y = y_ {c} + y_ {p}}

Если b>4c, то

yc = C 1 e - x 2 (b + b 2 - 4 с) + С 2 е - Икс 2 (b - b 2-4 с) {\ displaystyle y_ {c} = C_ {1} e ^ {- {\ frac {x} {2}} \, \ left ( b + {\ sqrt {b ^ {2} -4c}} \ right)} + C_ {2} e ^ {- {\ frac {x} {2}} \, \ left (b - {\ sqrt {b ^ {2} -4c}} \ right)}}{\ displaystyle y_ {c} = C_ {1} e ^ {- {\ frac {x} {2}} \, \ left (b + {\ sqrt {b ^ {2 } -4c}} \ right)} + C_ {2} e ^ {- {\ frac {x} {2}} \, \ left (b - {\ sqrt {b ^ {2} -4c}} \ right)}}

Если b = 4c, то

yc = (C 1 x + C 2) e - bx 2 {\ displaystyle y_ {c} = (C_ {1 } x + C_ {2}) e ​​^ {- {\ frac {bx} {2}}}}{\ displaystyle y_ {c} = (C_ {1} x + C_ {2}) e ​​^ {- {\ frac {bx} {2}}} }

Если b < 4c, then

yc = e - bx 2 [C 1 sin ⁡ (x 4 c - b 2 2) + C 2 соз ⁡ (Икс 4 C - b 2 2)] {\ Displaystyle Y_ {C} = e ^ {- {\ frac {bx} {2}}} \ left [C_ {1} \ sin \ left (x \, {\ frac {\ sqrt {4c-b ^ {2}}} {2}} \ right) + C_ {2} \ cos \ left (x \, {\ frac {\ sqrt {4c-b ^ {2}}} {2}} \ right) \ right]}{\ displaystyle y_ {c} = e ^ {- {\ frac {bx} {2}}} \ left [C_ {1} \ sin \ left (x \, {\ frac {\ sqrt {4c-b ^ {) 2}}} {2}} \ right) + C_ {2} \ cos \ left (x \, {\ frac {\ sqrt {4c-b ^ {2}}} {2}} \ right) \ right] }
n-го порядка, линейный, неоднородный, постоянные коэффициенты

∑ j = 0 nbjdjydxj = r (x) {\ displaystyle \ sum _ {j = 0} ^ {n} b_ {j} {\ frac {d ^ {j} y} {dx ^ {j}}} = r (x) \, \!}\ sum _ {j = 0} ^ {n} b_ {j} {\ frac {d ^ {j} y} {dx ^ {j}}} = r (x) \, \!

Дополнительная функция y c : принять y c = e, подставить и решить многочлен от α, чтобы найти линейно независимые функции e α jx {\ displaystyle e ^ {\ alpha _ {j} x}}е ^ {\ альфа _ {j} x} .

Частный интеграл y p : в общем случае метод изменения параметров, хотя для очень простой проверки r (x) может работать.

y = yc + yp {\ displaystyle y = y_ {c} + y_ {p}}{\ displaystyle y = y_ {c} + y_ {p}}

Поскольку α j - решения полинома степени n: ∏ j = 1 n (α - α j) = 0 {\ displaystyle \ prod _ {j = 1} ^ {n} (\ alpha - \ alpha _ {j}) = 0 \, \!}{\ displaystyle \ prod _ {j = 1} ^ {n} (\ alpha - \ alpha _ {j}) = 0 \, \!} , тогда:

для α j все разные,

yc = ∑ j = 1 n C je α jx {\ displaystyle y_ {c} = \ sum _ {j = 1} ^ {n} C_ {j} e ^ {\ alpha _ {j} x} \, \!}y_ {c} = \ sum _ {j = 1} ^ {n} C_ {j} e ^ {\ alpha _ {j} x} \, \!

для каждого корня α j re повторяется k j раз,

yc = ∑ j = 1 n (∑ ℓ = 1 kj C j, ℓ x ℓ - 1) e α jx {\ displaystyle y_ {c} = \ sum _ {j = 1} ^ {n} \ left (\ sum _ {\ ell = 1} ^ {k_ {j}} C_ {j, \ ell} x ^ {\ ell -1} \ right) e ^ {\ alpha _ {j} x} \, \!}{\ displaystyle y_ {c} = \ sum _ {j = 1} ^ {n} \ left (\ sum _ {\ ell = 1} ^ {k_ {j}} C_ {j, \ ell} x ^ {\ ell -1} \ right) e ^ {\ alpha _ {j} x} \, \!}

для некоторого комплекса α j, затем установив α = χ j + iγ j, и использование формулы Эйлера, позволяет записать некоторые термины из предыдущих результатов в виде

C je α jx = C je χ jx cos ⁡ (γ jx + φ j) {\ displaystyle C_ {j } e ^ {\ alpha _ {j} x} = C_ {j} e ^ {\ chi _ {j} x} \ cos (\ gamma _ {j} x + \ varphi _ {j}) \, \!}{\ displaystyle C_ {j} e ^ {\ alpha _ {j} x} = C_ {j} e ^ {\ chi _ {j} x} \ соз (\ gamma _ {j} x + \ varphi _ {j}) \, \!}

где ϕ j - произвольная постоянная (фазовый сдвиг).

Метод угадывания

Когда все другие методы решения ODE терпят неудачу, или в тех случаях, когда у нас есть некоторая интуиция относительно того, как может выглядеть решение DE, иногда возможно решить DE просто угадав решение и проверив его правильность. Чтобы использовать этот метод, мы просто угадываем решение дифференциального уравнения, а затем подставляем решение в дифференциальное уравнение, чтобы проверить, удовлетворяет ли оно уравнению. Если это так, то у нас есть конкретное решение для DE, в противном случае мы начинаем заново и пробуем другое предположение. Например, мы могли догадаться, что решение DE имеет форму: y = A e α t {\ displaystyle y = Ae ^ {\ alpha t}}{\ displaystyle y = Ae ^ {\ alpha t}} , поскольку это очень распространенное решение который физически ведет себя синусоидальным образом.

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

общее решение = однородное решение + частное решение {\ displaystyle {\ text {полное решение}} = {\ text {однородное решение }} + {\ text {особое решение}}}{\ displaystyle {\ text {total solution}} = {\ text { однородное решение}} + {\ text {особое решение}}}

Программное обеспечение для решения ODE

  • Maxima, система компьютерной алгебры с открытым исходным кодом.
  • COPASI, бесплатная (Artistic License 2.0 ) программный пакет для интеграции и анализа ODE.
  • MATLAB, приложение для технических вычислений (MATrix LABoratory)
  • GNU Octave, язык высокого уровня, в первую очередь предназначен для численных вычислений.
  • Scilab, приложение с открытым исходным кодом для численных вычислений.
  • Maple, проприетарное приложение для символьных вычислений.
  • Mathematica, проприетарное приложение, в первую очередь предназначенное для символьные вычисления.
  • Julia (язык программирования), язык высокого уровня, в первую очередь предназначенный для числовых вычислений.
  • SageMath, ope Приложение n-source, использующее синтаксис, подобный Python, с широким спектром возможностей, охватывающих несколько областей математики.
  • SciPy, пакет Python, включающий модуль интеграции ODE.
  • Chebfun, пакет с открытым исходным кодом, написанный на MATLAB, для вычислений с функциями с точностью до 15 цифр.
  • GNU R, вычислительная среда с открытым исходным кодом, в первую очередь предназначенная для статистики, которая включает пакеты для решения ODE.

См. Также

Примечания

Ссылки

  • Холлидей, Дэвид ; Резник, Роберт (1977), Physics (3-е изд.), Нью-Йорк: Wiley, ISBN 0-471-71716-9
  • Харпер, Чарли (1976), Введение в математическую физику, Нью-Джерси: Прентис-Холл, ISBN 0-13-487538-9
  • Крейсциг, Эрвин (1972), Advanced Engineering Mathematics (3-е изд.), Нью-Йорк: Wiley, ISBN 0-471-50728-8.
  • Полянин, А.Д. и В.Ф. Зайцев, Справочник по точным решениям обыкновенных дифференциальных уравнений (2-е издание) », Chapman Hall / CRC Press, Boca Raton, 2003. ISBN 1- 58488-297-2
  • Симмонс, Джордж Ф. (1972), Дифференциальные уравнения с приложениями и историческими примечаниями, Нью-Йорк: McGraw-Hill, LCCN 75173716
  • Типлер, Пол А. (1991), Физика для ученых и инженеров: Расширенная версия (3-е изд.), Нью-Йорк: Worth Publishers, ISBN 0-87901-432-6
  • Boscain, Ugo; Chitour, Yacine (2011), Introduction à l'automatique (PD F) (на французском языке)
  • Дреснер, Лоуренс (1999), Приложения теории обыкновенных и дифференциальных уравнений с частными производными Ли, Бристоль и Филадельфия: Издательский институт физики, ISBN 978-0750305303

Библиография

Внешние ссылки

В Wikibooks есть книга по теме: Исчисление / обыкновенные дифференциальные уравнения
Викимедиа У Commons есть материалы, связанные с обыкновенными дифференциальными уравнениями.
Последняя правка сделана 2021-06-01 14:14:59
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте