- •Е.А.Коган
- •и естественно-научным дисциплинам
- •Москва 2010
- •Настоящее пособие посвящено двум важным и тесно связанным между собой разделам математики: обыкновенным дифференциальным уравнениям и вариационному исчислению.
- •Подробное изложение теории дифференциальных уравнений и вариационного исчисления содержится во многих учебниках и пособиях, некоторые из которых, рекомендуемые студентам втузов для самостоятельного изучения, приведены в списке литературы.
- •1.1. Основные понятия
- •Обыкновенное дифференциальное уравнение n-го порядка в самом общем виде записывается так:
- •Для дифференциальных уравнений первого порядка различают общее, частное и особое решения, а также общий, частный и особый интегралы.
- •Поэтому общее решение дифференциального уравнения можно определить как множество всех частных решений уравнения.
- •Особым решением дифференциального уравнения называется решение, которое не может быть получено из общего решения ни
- •при одном частном значении произвольной постоянной.
- •Часто при интегрировании уравнения первого порядка не удается найти его общее решение в явном виде, а получается конечное (не дифференциальное) соотношение вида
- •Знание изоклин позволяет во многих случаях даже для не интегрируемых явно дифференциальных уравнений получить графическое решение задачи Коши и выявить характер интегральных кривых.
- •Пример. Построить методом изоклин интегральную кривую уравнения
- •Таким образом, изоклинами данного уравнения являются концентрические окружности с центром в начале координат, причëм угловые коэффициенты касательных к искомым интегральным кривым равны радиусам этих окружностей С.
- •Очевидно, это уравнение с разделëнными переменными. Интегрируя его, получим
- •Следовательно, общий интеграл уравнения будет
- •Тогда
- •Тогда уравнение примет вид
- •Дифференциальное уравнение вида
- •называется приводящимся к однородному. В частности, к этому классу относится уравнение вида
- •В результате уравнение приводится к однородному
- •Возвращаясь к старой переменной, получим
- •Общее решение уравнения примет вид
- •После разделения переменных получим
- •Тогда уравнение (1.18) примет вид
- •Следовательно,
- •Подставляя (1.19) и (1.20) в (1.16), получим общее решение уравнения (1.12) в виде (1.15).
- •Уравнение вида
- •Следовательно, его общий интеграл, а значит, и общий интеграл уравнения (1.24) имеет вид
- •В результате общее решение уравнения будет иметь вид
- •После интегрирования получим
- •Следовательно, общее решение уравнения (2.18) будет
- •Далее ищем решение уравнения (2.17) в форме, аналогичной по структуре выражению (2.19), но произвольную постоянную в (2.19) заменяем неизвестной функцией
- •Подставляя (2.20) в (2.17), получим
- •В качестве практического примера применения метода понижения порядка рассмотрим задачу об осесимметричном изгибе упругих круговых пластин.
- •Определителем Вронского или вронскианом называется функциональный определитель вида:
- •Поэтому теорема о структуре общего решения линейного однородного дифференциального уравнения n-го порядка может быть переформулирована так:
- •а общее решение запишется в виде
- •Если известно общее решение однородного линейного дифференциального уравнения, соответствующего заданному неоднородному, то его частное решение может быть найдено методом вариации произвольных постоянных (методом Лагранжа).
- •Рассмотрим его реализацию для линейных дифференциальных уравнений второго порядка:
- •Интегрируя дифференциальные уравнения первого порядка (2.66), находим
- •Пример. Решить уравнение
- •Определитель этой системы
- •Интегрируя уравнения (2.71), находим
- •Подставляя (2.72) в (2.70), получим общее решение уравнения в виде
- •Рассмотрим линейное однородное дифференциальное уравнение второго порядка
- •Пример. Найти собственные значения и собственные функции однородной краевой задачи для уравнения (2.90) при граничных условиях
- •Общее решение уравнения согласно (2.92) имеет вид
- •Вычисляем
- •Раскрывая определитель, получим уравнение относительно параметра
- •корни которого являются собственными значениями задачи:
- •Уравнением Эйлера называется линейное уравнение вида
- •Пример. Решить уравнение
- •Уравнение (2.103) есть уравнение Эйлера второго порядка. Применим замену независимой переменной и производных по формулам (2.101), (2.102). Тогда уравнение примет вид
- •Его общее решение
- •Эти решения находятся методом подбора (см. выше) и имеют вид
- •Поэтому общее решение уравнения (2.104) будет
- •Пример. Решить задачу Коши для уравнения
- •Общее решение уравнения имеет вид
- •Убедимся в этом на примере системы двух уравнений относительно двух неизвестных функций
- •Будем искать частное решение однородной системы в виде
- •Пример. Решить методом Эйлера систему
- •Характеристическое уравнение системы
- •Общее решение однородной системы (3.18) запишется в виде
- •Этот метод применим к решению систем неоднородных линейных уравнений n-го порядка. Ограничимся для простоты нормальной системой двух линейных уравнений с постоянными коэффициентами
- •Пусть общее решение однородной системы уравнений известно:
- •Пример. Решить систему
- •Принимаем частное решение системы (3.31) в виде
- •Интегрируя эти уравнения, получим
- •Общее решение системы запишется в виде:
- •Пусть требуется решить краевую задачу для уравнения
- •Будем искать решение задачи в виде
- •Задача сводится к решению дифференциального уравнения
- •Выбранные функции удовлетворяют всем перечисленным выше требованиям. Они линейно независимы, непрерывно дифференцируемы и удовлетворяют граничным условиям (4.12).
- •Решение методом Бубнова
- •Как видим, ошибка, полученная при решении методом Бубнова в двучленном приближении, составляет сотые доли процента.
- •Решение методом наименьших квадратов
- •Решение методом коллокаций
- •В результате находим
99
2.9.2. Решение задачи Коши методом степенных рядов
Пусть дано линейное дифференциальное уравнение n-го порядка
y(n) + a1 (x) y(n−1) + a2(x) y(n−2) + an−1 (x) y′+ an (x) y = 0,
коэффициенты которого являются аналитическими функциями (то есть представимы в виде степенных рядов), и пусть начальные условия имеют вид
y(x0) = y0, |
y′(x0) = y0′, ..., y(n−1) (x0) = y0(n−1), |
то есть сформулирована задача Коши.
Будем искать решение уравнения в виде бесконечного степенного ряда с неопределёнными коэффициентами
∞
y(x) = ∑cn xn = c0 +c1x +c2x2 +...+cn xn +...
n=0
(здесь для простоты принято x0 = 0).
Формальный алгоритм определения коэффициентов ряда следующий: ряд подставляется в уравнение. В получающемся тождественном равенстве коэффициенты при различных степенях x приравниваются нулю. В результате получается система уравнений для определения коэффициентов cn . Полученное ре-
шение исследуют на сходимость и на возможность почленного дифференцирования. В области, в которой ряд сходится и допускает n – кратное дифференцирование, он и является искомым решением.
Пример. Решить задачу Коши для уравнения
y′′− 2xy′− 4 y = 0 |
(2.106) |
при начальных условиях |
(2.107) |
y(0) = 0, y (0) =1. |
|
′ |
|
Представим решение задачи (2.106), (2.107) степенным рядом с неопределёнными коэффициентами:
∞ |
|
y = ∑an xn =a0 + a1 x + a2x2 + + an xn + |
(2.108) |
n=0
Подставляем искомое решение (2.108) в уравнение (2. 106), учитывая соотношения
y′ = a1 + 2a2x +3a3 x2 + + nan xn−1 +
100
y′′ = 2a2 +3 2a3 x + 4 3a4 x2 + + n(n −1)an xn−2 +
Тогда имеем:
[2a2 +3 2a3 x + 4 3a4 x2 + + n(n −1)an xn−2 + ] −2x(a1 +
+2a2x +3a3 x2 + + nan xn−1 + )−4(a0 + a1x + a2x2 + a3 x3 +
++ an xn + ) = 0.
Группируя подобные члены при одинаковых степенях x, получим тождественное равенство:
2a2 − 4a0 + (3 2a3 − 2a1 − 4a1 )x + (4 3a4 − 2 2a2 − 4a2 )x2 + (5 4a5 −
−3 2a3 − 4a3 )x3 +...+[n(n −1)an − 2(n − 2)an−2 − 4an−2 ]xn−2 +... ≡ 0.
В левой части этого тождества имеем степенной ряд. Ряд будет тождественно равен нулю, если его коэффициенты будут равны нулю. Приравнивая коэффициенты при различных степенях x нулю, получаем систему уравнений:
2a2 − 4a0 = 0, 3 2a3 − 6a1 = 0,
4 3a4 −8a2 = 0, 5 4a5 −10a3 = 0,
……………………….
n(n −1)an − 2nan−2 = 0.
Из этой системы следуют рекуррентные формулы:
a2 = 2a0, a3 = a1 ,
………….. |
|
|||
an = |
2an−2 |
, |
(2.109) |
|
n −1 |
||||
|
|
|
позволяющие выразить последующие коэффициенты ряда через предыдущие.
Подчиняя искомое решение (2.108) начальным условиям, получаем a0 = 0, a1 =1. Остальные коэффициенты ряда опреде-
ляются по рекуррентным формулам (2.109):
a2 = a0 = 0, a3 = a1 =1, a4 = 0, a5 =1/ 2!, a6 = 0, a7 =1/ 3! и т. д.
101
После подстановки коэффициентов в (2.108) находим окончательное выражение
y = x + |
x3 |
+ |
x5 |
+ |
x7 |
+ + |
x2 n+1 |
+ |
(2.110) |
|
1! |
2! |
3! |
n! |
|||||||
|
|
|
|
|
|
Применяя признак Даламбера17, легко установить, что ряд
(2.110) сходится на всей числовой оси. Следовательно, выражение (2.110) представляет решение исходной задачи при
x(−∞;∞).
2.9.3.Построение общего решения линейного неоднородного уравнения методом степенных рядов
Пусть требуется найти общее решение линейного неоднородного дифференциального уравнения
y′′ − xy′ − y = 2. |
(2.111) |
Общее решение уравнения имеет вид |
|
y = C1 y1 +C2 y2 + y , |
(2.112) |
где y - частное решение неоднородного уравнения, |
C1 и C2 - |
произвольные постоянные, y1 и y2 - частные линейно независимые решения однородного уравнения
y′′ − xy′ − y = 0. |
(2.113) |
Чтобы построить общее решение уравнения (2.111), |
будем р е- |
шать задачу Коши для этого уравнения при начальных условиях:
|
|
y(0) = C1 , |
y (0) = C2 . |
(2.114) |
|||||
|
|
|
|
|
|
′ |
|
|
|
Решение задачи (2.113), |
(2.114) |
представим рядом |
|
||||||
y = a |
0 |
+ a x + a |
x2 |
+ a x3 + + a |
xn + |
(2.115) |
|||
|
1 |
2 |
|
3 |
|
n |
|
|
|
Легко видеть, |
что a0 |
= y(0), |
a1 |
= y1′(0), то есть первые два |
коэффициента разложения решения в ряд (2.115) представляют собой соответственно начальное значение функции и её первой производной. Поэтому в соответствии с начальными условиями (2.114) можно принять
17 Д а л а м б е р Жан Лоран (16.11.1717 – 29.10.1783) – французский математик, механик, физик.
102
a0 = C1 , a1 = C2 . |
(2.116) |
Для определения последующих коэффициентов ряда подставляем его в уравнение (2.111)
[2a2 +3 2a3 x + + n (n −1)an xn−2 + ] − x[a1 + 2a2x +3a3 x2 +
+ + nan xn−1 + ] −[a0 + a1x + a2x2 + a3 x3 + + an xn + ] = 2,
следовательно,
2a2 −a0 +(3 2a3 −2a1 )x +(4 3a4 −3a2)x2 +(5 4a5 −4a3 )x3 +...+ +[n(n −1)an −(n −1)an−2] xn−2 + = 2.
Приравнивая коэффициенты при одинаковых степенях x в левой и правой частях полученного равенства, приходим к соотношениям вида
2a2 − a0 = 2, 3 2a3 − 2a1 = 0,
4 3a4 −3a2 = 0,
…………………………
n(n −1)an −(n −1)an−2 = 0,
В результате рекуррентные формулы для коэффициентов примут вид:
a2 |
= a0 / 2+1, |
|
a3 |
= a1 / 3, |
|
a4 |
= a2 / 4, |
|
……………. |
|
|
an |
= an−2 / n. |
(2.117) |
По этим формулам с учётом (2.116) получаем:
a2 = С1 / 2+1, a3 = С2 / 3,
a4 = С1 /(2 4) +1/ 4, a5 = C2 /(3 5),
a6 = C1 /(2 4 6) +1/(4 6)
и так далее. Подставляя найденные коэффициенты в искомое решение (2.115) и группируя подобные члены, получим общее ре-
103
шение уравнения (2.111) в форме: |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
y = C |
|
|
|
|
|
|
+ x |
|
|
+ |
|
x |
|
|
|
+C |
|
|
|
+ |
x |
|
+ |
|
x |
|
|
||
1+ x |
2 |
4 |
|
6 |
|
+ |
x + x |
3 |
5 |
|
7 |
+ + |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 2 4 |
|
2 4 6 |
3 |
|
3 5 3 |
5 7 |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
+ x2 + |
x4 |
+ |
|
|
x6 |
+ |
|
|
x8 |
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
(2.118) |
|||||
|
|
|
|
|
4 6 |
8 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
4 |
|
|
|
4 6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Первые два ряда в (2.118) сходятся на всей числовой оси и определяют функции y1 и y2, линейно независимые в окрестно-
сти точки x=0. Ряд
y = x2 + |
x4 |
+ |
x6 |
+ |
x8 |
|
+ |
|
4 6 |
4 6 |
8 |
||||
4 |
|
|
|
также сходится на всей числовой оси и удовлетворяет уравнению (2.111). Поэтому функция y является разложением в сте-
пенной ряд частного решения исходного линейного неоднородного уравнения. Заметим, что в данном примере частное решение y можно было проще найти методом подбора. Оно, очевидно,
равно y = −2. Тогда функции y1 и y2 следует находить как ча-
стные решения однородного уравнения (2.113), используя рекуррентные соотношения (2.117) и формулы (2.116).
2.9.4. Разложение решения задачи Коши в ряд Тейлора
Пусть требуется найти решение дифференциального уравнения n-го порядка, разрешëнного относительно старшей производной
y |
(n) |
′ |
′′ |
(n−1) |
), |
(2.119) |
|
= f (x,y,y ,y |
,...,y |
|
удовлетворяющее начальным условиям
y(x0 ) = y0 , y′(x0 ) = y0′, ...y(n−1) (x0 ) = y0(n−1). (2.120)
Будем искать решение задачи в виде разложения в ряд Тейлора
|
∞ |
(n) |
(x0) |
|
|
|
y′(x0) |
|
|
|||
y(x) = ∑ |
y |
|
(x − x0)n = y(x0)+ |
(x − x0)+ |
|
|||||||
|
n! |
|
|
|||||||||
|
n=0 |
|
|
|
1! |
|
|
|
||||
+ |
y′′(x0) |
(x − x )2 |
+...+ |
y(n) (x0) |
(x − x )n |
+... |
(2.121) |
|||||
|
|
|||||||||||
2! |
|
|
|
0 |
|
n! |
0 |
|
|
|
||
|
|
|
|
|
|
|
|
|
104
В частности, при x0 = 0 получаем разложение в ряд по степеням x (ряд Маклорена18):
|
∞ |
|
y |
(n) |
(0) |
|
|
y′(0) |
|
y′′(0) |
|
|
||
y(x) = ∑ |
|
xn =y(0) |
+ |
x + |
x2 |
+...+ |
||||||||
|
|
|
|
|
|
|||||||||
|
n=0 |
|
n! |
1! |
2! |
|
(2.122) |
|||||||
|
y(n) (0) |
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
x |
n |
+... |
|
|
|
|
|
|
|
||||
n! |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Задача сводится к отысканию коэффициентов ряда (2.121). Как видно, первые n коэффициентов ряда могут быть найдены непосредственно из начальных условий (2.120); (n+1)-ый коэф-
фициент разложения решения в ряд: y(n) (x0 ) определяется из
дифференциального уравнения (2.119) после подстановки в правую часть начальных условий (2.120). Последующие коэффициенты ряда (2.121) находятся дифференцированием правой части уравнения (2.119) (обычно с использованием правила дифференцирования сложной функции) и начальных условий (2.120). Если полученный ряд сходится, то в интервале сходимости он является решением задачи ( 2.119), (2.120).
Заметим, что этот метод является достаточно универсальным и применим к нелинейным уравнениям с переменными коэффициентами, но вычисление производных высокого порядка может оказаться трудоёмким из-за необходимости последовательного дифференцирования сложных функций.
Пример. Решить задачу Коши для уравнения
y′′ = y′2 − xey |
(2.123) |
при начальных условиях y(0) =1, y′(0) =1.
Принимаем решение в виде (2.122). В рассматриваемой задаче первые два коэффициента ряда определяются из начальных условий. Третий коэффициент получается непосредственно из дифференциального уравнения (2.123) при x = 0 и равен y′′(0) =1.Дифференцируя последовательно уравнение (2.123),
получим
18 М а к л о р е н Колин (1698 – 14.06.1746) – шотландский математик.
|
|
|
|
|
|
|
105 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
′′′ |
|
|
′ ′′ |
− e |
y |
− xe |
y |
|
′ |
|
|
|
|
|
|||
|
|
(x) = 2y y |
|
|
y , |
|
|
|
|
|
||||||||||
y |
IV |
(x) = 2y |
′′2 |
|
′ |
′′′ |
− |
2e |
y |
y |
′ |
− xe |
y |
y |
′2 |
− xe |
y |
′′ |
||
|
|
+ 2y y |
|
|
|
|
|
|
y . |
|||||||||||
Поэтому при x=0 находим |
′′′ |
|
|
|
|
|
|
|
|
|
|
|
|
y |
IV |
(0) = 6−4e. |
||||
y (0) = 2 1 1−e = 2−e, |
|
Ограничиваясь первыми пятью членами ряда (2.122), получим решение задачи в виде
y(x) =1 + x + 21! x2 + 23−! e x3 + 6−4!4e x4 +...
Полученное выражение является решением при условии сходимости ряда и даёт достаточно точное значение функции y в малой окрестности точки x = 0.
106
3. СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
Решение многих прикладных задач, как правило, приводит к системам дифференциальных уравнений. Кроме того, при численном решении дифференциальных уравнений высокого порядка обычно применяется процедура сведения дифференциального уравнения n-го порядка относительно одной неизвестной функции к системе дифференциальных уравнений относительно n неизвестных функций.
Она базируется на следующей теореме: Дифференциальное уравнение n-го порядка y(n) = f (x, y, y′, y′′,..., y(n−1) ), разрешëнное
относительно старшей производной, может быть сведено к системе n дифференциальных уравнений первого порядка. Это достигается следующей процедурой замены функции и ëе произво д- ных разного порядка новыми функциями:
y = y1,
y′ = y1′ = y2,
y′′ = y1′′= y′2 = y3 ,
..............................
y(n−1) =.....= y′n-1 = yn ,
y(n) = y′n = f (x, y1, y2,...yn ).
В результате, исходное дифференциальное уравнение оказывется эквивалентным следующей системе дифференциальных уравнений первого порядка относительно n неизвестных функций
y1 (x), y2(x), , yn (x).
y′ |
= y |
|
, |
1 |
|
2 |
|
y2′ = y3 , |
|||
.................... |
|||
|
|
|
|
|
|
|
|
yn′−1 = yn , |
|||
′ |
= f (x, y1, y2,..., yn ). |
||
yn |
107
Система дифференциальных уравнений первого порядка, разрешëнных относительно производных от неизвестных функций
yi′= fi (x, y1 , y2, , yn ), (i =1,2,...,n) |
(3.1) |
называется системой, записанной в нормальной форме Коши или нормальной системой.
В частности, нормальная система линейных дифференциальных уравнений с постоянными коэффициентами записывается в форме
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
yi′ = ∑aik yk + |
fi (x) |
|
(i=1,2,…,n) |
|||||||||||||
или |
|
|
|
k =1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
= a y + a y + + a y + f (x), |
|||||||||||||||||
|
dx |
11 1 |
12 2 |
|
|
|
1n n |
|
1 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy2 |
= a21y1 + a22 y2 + + a2n yn + f2 (x), |
||||||||||||||||
|
|
|
|
|||||||||||||||
|
|
|
||||||||||||||||
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
...................................................... |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dyn |
= a |
y |
+ a |
n2 |
y |
2 |
+ |
+ a |
n n |
y |
n |
+ f |
n |
(x). |
|||
|
|
|
|
n1 1 |
|
|
|
|
|
|
|
|||||||
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Если |
fi (x) = 0 |
(i =1,2,...,n), |
|
то |
система |
называется одно- |
||||||||||||
родной. Если хотя бы одна из функций |
fi (x) ≠ 0, то система на- |
|||||||||||||||||
зывается неоднородной. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Общим решением системы (3.1) на интервале (a, b) изме- |
||||||||||||||||||
нения аргумента |
x |
называется |
|
всякая |
совокупность функций |
|||||||||||||
yi (x,C1,C2, |
...,Cn ) |
(i =1,2,...,n), дифференцируемых на интервале |
(a, b), и обращающих каждое из уравнений системы (3.1) в тождество.
Число произвольных постоянных, входящих в общее решение нормальной системы уравнений, равно числу неизвестных функций. Решение, получающееся из общего при конкретных значениях произвольных постоянных, называется частным решением. Для получения частного (единственного) решения системы произвольные постоянные определяются из начальных или граничных условий.