- •Е.А.Коган
- •и естественно-научным дисциплинам
- •Москва 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).
- •Решение методом Бубнова
- •Как видим, ошибка, полученная при решении методом Бубнова в двучленном приближении, составляет сотые доли процента.
- •Решение методом наименьших квадратов
- •Решение методом коллокаций
- •В результате находим
40
Таким образом, с математической точки зрения задача об изгибе консольно защемлëнного стержня сводится к задаче Коши (2.2), (2.3), а задача об изгибе шарнирноëопртого стержня – к краевой задаче (2.2), (2.4).
2.1.Интегрирование дифференциальных уравнений n – го порядка методом понижения порядка
Если правая часть уравнения (2.1) является известной непрерывной функцией f (x) только одной переменной x или не со-
держит искомую функцию y: f (x,y′,y′′, ,y(n−1) ), или не содер-
жит явно независимую переменную x: f (y,y′,y′′, ,y(n−1) ), то для
решения уравнения (2.1) может быть применëн метод понижения порядка.
1. y(n) (x) = f (x). (2.5)
Это простейшее уравнение n -го порядка, общее решение которого получается в квадратурах последовательным интегрированием n раз. При каждом интегрировании порядок уравнения понижается на единицу, и появляется произвольная постоянная:
y(n−1) = ∫ f (x)dx +C1,
y(n−2) = ∫dx∫ f (x)dx +C1x +C2,
………………………………………….
В результате общее решение уравнения будет иметь вид
y = ∫∫ ∫ f (x)dx + C1 xn−1 + C2xn−2 + + Cn−1 x + Cn , |
(2.6) |
где в правой части - n-кратный интеграл от функции f (x) |
и мно- |
гочлен (n-1)-ой степени от x, коэффициентами которого являются
n произвольных постоянных. |
|
|
|
|
|
|
|
|||||||
|
Если |
для уравнения (2.5) решается задача Коши с началь- |
||||||||||||
ными условиями |
y(x0 ) = y0, |
y′(x0 ) = y0′,…, y(n−1) (x0) = y0 |
(n−1), то |
|||||||||||
частное решение уравнения может быть записано в виде |
|
|||||||||||||
|
x |
|
x |
|
x |
|
y0(n−1) |
|
(x − x )n−1 + |
y0(n−2) |
|
(x − x )n−2 |
|
|
y = |
∫ |
dx |
∫ |
dx.... |
∫ |
f (x)dx + |
|
|
+...+ |
|||||
|
|
(n −2)! |
||||||||||||
|
|
|
|
(n −1)! |
0 |
0 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
x0 |
|
x0 |
|
x0 |
|
|
|
|
|
|
|
|
|
+ y0′(x − x0 )+ y0. |
|
|
|
|
|
|
|
|
41
Пример. Найти максимальный прогиб шарнирно опëртого стержня при действии на него равномерно распределëнной поп е- речной нагрузки интенсивности q(x) = q0 (см. рис. 5б).
Изгибающий момент в произвольном поперечном сечении стержня с абсциссой x. равен алгебраической сумме моментов сил, действующих по одну сторону от сечения:
M (x) = q20l x − q02x2 ,
где первое слагаемое – момент от опорной реакции, а второе – момент от равномерно распределëнной поперечной нагрузки. Поэтому уравнение изгиба стержня (2.2) примет вид
|
1 |
|
|
|
2 |
|
y′′ = |
q0l x − q0x |
. |
||||
|
||||||
|
|
|
2 |
2 |
|
|
|
EI |
|
Это уравнение относится к типу уравнений (2.5). Интегрируя его дважды, получим общее решение в виде
|
q0 |
|
|
3 |
|
x |
4 |
|
|
y = |
lx |
|
− |
|
+C x +C |
. |
|||
|
|
|
|
||||||
|
|
|
6 |
12 |
1 |
2 |
|||
|
2EI |
|
|
Произвольные постоянные находим из граничных условий (2.4):
|
C |
= − |
|
l3 |
, |
|
C |
|
= 0. |
|
|
|
||||
|
12 |
|
|
|
|
|
||||||||||
|
1 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
||
Подставляя их в общее решение, получим |
|
|
|
|||||||||||||
y = |
q |
l4 |
x 3 |
|
x |
4 |
x |
(2.7) |
||||||||
0 |
|
|
2 |
|
|
− |
|
|
|
− |
|
. |
||||
|
|
|
|
|||||||||||||
|
24EI |
l |
|
|
l |
|
l |
|
Максимальный прогиб стержня при действии равномерно
распределëнной нагрузки при |
x = l / 2 |
|
будет равен |
|
||||||||||||
|
|
|
|
y |
|
= − |
5q |
l4 |
|
. |
|
|
|
|||
|
|
|
|
max |
|
|
0 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
384EI |
|
|
|
||||||
2. |
y |
(n) |
= |
|
|
|
′ |
|
′′ |
, |
,y |
(n−1) |
). |
(2.8) |
||
|
f (x,y ,y |
|
|
Уравнение (2.8) не содержит искомой функции y(x).
Рассмотрим процедуру интегрирования уравнения данного типа на примере уравнения второго порядка
y |
′′ |
= f (x,y ). |
(2.9) |
|
′ |
|
|
42 |
|
Понижение порядка достигается подстановкой |
(2.10) |
|
y (x) = z(x). |
||
′ |
|
|
Тогда |
|
(2.11) |
y (x) = z (x), |
||
′′ |
′ |
|
и уравнение (2.9) приводится к уравнению первого |
порядка от- |
|
носительно функции z(x): |
|
|
′ |
|
(2.12) |
z (x) = f (x,z). |
Интегрируя уравнение (2.12), находим его общий интеграл в
виде
z(x) = ϕ(x,C1 ), |
(2.13) |
где С1 - произвольная постоянная. Далее в (2.13) заменяем левую
часть согласно (2.10). и вновь получаем уравнение первого порядка относительно искомой функции y:
y (x) = ϕ(x,C1 ), |
(2.14) |
′ |
|
Интегрируя уравнение (2.14), находим общее решение исходного уравнения (2.9) в виде
y(x) = ξ(x,C1 ,C2 ). |
(2.15) |
Замечание. Если уравнение (2.9) не содержит ни искомой функции y, ни её производных до ( k- 1) - го порядка включительно, то есть имеет вид
y(n) = f (x, y(k ), y(k +1), , y(n−1) ),
то его порядок может быть понижен сразу на k единиц подстановкой
y(k ) = z(x).
Пример. Решить уравнение x3 yI V + x2 y′′′ =1. |
|
||||
Это уравнение не |
содержит явно искомой функции y(x) и |
||||
её первых производных |
y , |
y |
′′ |
и относится ко второму из рас- |
|
|
′ |
|
|
|
|
смотренных нами типов. Применяя подстановку |
(2.16) |
||||
|
y |
|
(x) = z(x), |
||
|
|
′′′ |
|
|
получаем линейное неоднородное уравнение первого порядка
43
x3 z′+ x2z =1
или |
1 |
1 |
|
|
z′+ |
|
|
||
x z = |
|
. |
(2.17) |
|
x3 |
Интегрируем его методом вариации произвольной постоянной. Сначала решаем соответствующее однородное уравнение
z′+ |
1 |
|
x z = 0. |
(2.18) |
Разделяем в нëм переменные:
dzz = −dxx .
После интегрирования получим
ln z = −ln x + lnC1 = ln Cx1 .
Следовательно, общее решение уравнения (2.18) будет
z = |
C1 . |
(2.19) |
|
x |
|
Далее ищем решение уравнения (2.17) в форме, аналогичной по структуре выражению (2.19), но произвольную постоянную в (2.19) заменяем неизвестной функцией v(x):
|
|
|
|
|
z = |
v(x) |
. |
|
|
(2.20) |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
x |
|
|
|
||
Подставляя (2.20) в (2.17), получим |
|
|
|
||||||||
|
|
′ |
|
|
|
|
|
1 |
|
||
|
|
v (x)x −v(x) + v(x) |
= |
, |
|||||||
|
dv(x) |
|
1 |
x2 |
|
x2 |
|
x3 |
|
||
откуда следует |
= |
и, |
после разделения переменных, |
||||||||
dx |
x2 |
||||||||||
|
|
|
|
|
|
|
|
|
dv(x) = dxx2 .
Интегрируя это уравнение, находим v(x) = −1/ x +C1. Поэтому согласно (2.20) имеем
z = − |
1 |
+ C1 . |
(2.21) |
|
x2 |
||||
|
x |
|
44
Заменяя в выражении (2.21) z по формуле (2.16), приходим к уравнению третьего порядка относительно искомой функции y(x):
y |
′′′ |
1 |
C1 |
|
||
|
|
+ x . |
|
|||
= − x2 |
(2.22) |
|||||
|
||||||
Уравнение (2.22) содержит |
в правой |
части известную |
функцию от x и относится к первому из рассмотренных нами типов (см. (2.5)). Интегрируя его последовательно три раза, окончательно получим общее решение исходного уравнения, содержащее 4 произвольных постоянных :
|
|
|
|
y′′ = 1 + C1 ln x + C2, |
|
|
|
|
|||||||||
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
y′ = ln x +C1x(lnx −1)+C2x +C3 , |
|
|
||||||||||||||
|
|
|
|
|
x2 |
|
3 |
|
|
|
|
x2 |
|
|
|
||
|
y = x(lnx −1) + C |
|
|
ln x − |
|
|
+ C |
|
|
+ C |
x + C |
. |
|||||
|
|
|
2 |
2 2 |
|||||||||||||
|
|
|
|
1 2 |
|
|
|
|
|
|
3 |
4 |
|
||||
3. |
y |
(n) |
= |
f (y, |
|
′ |
′′ |
|
y |
(n−1) |
). |
|
|
|
(2.23) |
||
|
y , y |
, , |
|
|
|
|
Это уравнение не содержит явно независимой переменной x. В частном случае уравнение второго порядка данного типа будет
y |
′′ |
= f (y,y ). |
(2.24) |
|
′ |
|
|
Понижение порядка достигается подстановкой |
(2.25) |
||
y (x) = z(y). |
|||
|
′ |
|
|
Тогда по правилу вычисления производной от сложной функции
|
′′ |
dz dy |
dz |
|
y |
(x) = dy dx |
= dy z . |
(2.26) |
Поэтому уравнение (2.24) приводится к уравнению первого порядка относительно функции z(y):
z dz |
= f (y,z). |
(2.27) |
dy |
|
|
Интегрируя уравнение (2.27), находим его общее решение в виде
z(y) =ϕ(y,C1 ), |
(2.28) |
где C1 - произвольная постоянная. Далее в (2.28) заменяем левую часть согласно (2.25) и вновь приходим к уравнению первого по-
45
рядка с разделяющимися переменными относительно искомой
функции y : y (x) =ϕ (y,C1 ), следовательно, |
|
|||
′ |
|
|
||
|
dy |
= dx. |
(2.29) |
|
ϕ(y,C ) |
||||
|
|
|||
1 |
|
|
Интегрируя уравнение (2.29), окончательно получим общий интеграл исходного уравнения (2.24) в виде
|
|
|
|
|
y(x) = ξ(x,C1 ,C2 ). |
|
(2.30) |
|||||||||
Пример. Решить задачу Коши: |
|
|
|
|
|
|
|
|
|
|||||||
16y |
3 |
y |
′′ |
= 4 y |
4 |
−1, |
y(0) = |
|
2 |
|
|
′ |
|
2 |
|
|
|
|
2 |
, |
2 . |
||||||||||||
|
|
|
y (0) = |
Уравнение не содержит явно независимой переменной x и потому относится к третьему из перечисленных типов. Принимаем y′(x) = z(y), тогда с учетом (2.26) получим уравнение первого
порядка с разделяющимися переменными
16y3 z dydz = 4 y4 −1.
Приводим его к виду zdz = 4 y4 −3 1 dy. Интегрируя, получим
16y
|
z2 |
|
y2 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
y2 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||
|
|
= |
|
+ |
|
|
|
+C |
|
, |
z = ± |
|
|
|
|
|
+ |
|
|
|
|
|
|
+C |
, где C |
= 2C |
. |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
2 |
|
8 |
|
16y2 |
1 |
|
|
|
|
|
|
|
4 |
|
|
|
|
16y2 |
|
|
|
1 |
|
1 |
1 |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
Заменяя z |
|
на y′ |
согласно (2.25), |
|
|
|
вновь получаем уравнение |
||||||||||||||||||||||||||||||
первого порядка |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
y′ = ± |
4 y4 |
+16C y2 + |
1 |
|
. |
|
|
|
(2.31) |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
Определим произвольную постоянную C1 |
из начальных ус- |
||||||||||||||||||||||||||||||||||
ловий. При x = 0 равенство (2.31) примет вид |
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
1 |
+16 |
1 C +1 |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
4 |
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
= ± |
|
|
|
|
|
|
|
|
|
2 |
1 |
|
|
|
|
, |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
2 |
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
откуда получим 2 = ±2+8C1 , C1 = 14 . Подставляя C1 в (2.31), приходим к уравнению
|
|
|
|
|
|
|
|
|
46 |
|
|
|
|
|
|
|
2y2 +1 |
|
|
y |
′ |
|
|
4 y4 |
+ 4 y2 +1 |
|
|
|
|
= |
|
|
|
|
|
= 4 y (знак минус отброшен, так как при |
|||
|
|
4 y |
|||||||
|
|
|
этом не удовлетворяются начальные условия).
Разделяя переменные, имеем |
4 ydy |
|
= dx , следовательно, после |
|
2y2 +1 |
||||
|
|
интегрирования получим ln(2y2 +1) = x +ln C2, 2y2 +1 = C2ex . Из первого начального условия находим C2 = 2.
В результате решение задачи Коши запишется в виде y = ex − 12.
Пример. Решить задачу Коши: y |
′′ |
= 2y |
3 |
, |
y(0) =1, |
′ |
||||||
|
|
|
y (0) =1. |
|||||||||
Полагая, как и выше, |
′ |
|
|
y |
′′ |
dz |
получим |
|||||
y (x) = z(y), |
|
(x) = dy z , |
||||||||||
уравнение первого порядка |
z dz = 2y3 |
. |
Разделяя переменные и |
|||||||||
|
dy |
|
|
|
|
|
|
|
|
|
|
|
интегрируя, получим z = ± |
|
|
|
|
|
|
y′ = |
|
(знак |
|||
y4 +C1 |
. |
|
|
|
|
|
y4 +C1 |
минус отброшен, так как при этом не удовлетворяются начальные условия).
Используя начальные условия, находим C1 = 0. Поэтому y′ = y2 , и уравнение легко интегрируется: − 1y = x +C2. Далее из первого
начального условия находим C2 = −1. Поэтому решение задачи Коши запишется в виде
y =1−1 x .
Замечание. При решении задачи Коши часто целесообразно (как в приведëнных примерах) определять значения произвол ь- ных постоянных в процессе решения, а не после нахождения общего решения исходного уравнения, так как при произвольных значениях C1,C2,...интегрирование может оказаться весьма слож-
ным, а то и вообще невозможным. Действительно, например, в последней задаче решение уравнения y′ = y4 +C1 приводит по-