- •Е.А.Коган
- •и естественно-научным дисциплинам
- •Москва 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).
- •Решение методом Бубнова
- •Как видим, ошибка, полученная при решении методом Бубнова в двучленном приближении, составляет сотые доли процента.
- •Решение методом наименьших квадратов
- •Решение методом коллокаций
- •В результате находим
|
|
|
|
|
|
|
|
90 |
|
|
|
|
|
|
|
|
где |
|
|
|
|
b3 |
|
|
|
|
|
= b4 |
|
|
|
|
|
Ω2 = |
(b ω) |
|
|
, |
a |
0 |
, |
|
a = Ω2 |
(b −ω2) |
||||||
|
|
|
2 |
+(b −ω2)2 |
|
|
|
b |
|
|
1 |
2 |
|
|||
|
|
|
1 |
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
a |
2 |
= b ωΩ2 |
, |
a = −(a |
0 |
+ a ), |
|
a |
4 |
= −(a ω + a α)β−1. |
||||||
|
|
1 |
|
3 |
|
|
1 |
|
|
|
2 |
3 |
|
|||
График зависимости yк |
от времени показан на рис. I.13. |
|
||||||||||||||
8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
0,1 |
|
0,2 |
|
|
|
|
|
0,3 |
|
|
0,4 |
t |
0,5 |
Рис. I.13. |
|
График колебаний центра масс колеса |
|
2.8. Задачи на собственные значения
Рассмотрим линейное однородное дифференциальное уравнение второго порядка
y′′+λy = 0 |
(2.90) |
|
c однородными граничными условиями |
|
|
y(0) = 0, |
y(l) = 0. |
(2.91) |
Здесь предполагается, что |
0≤ x ≤ L, |
λ −некоторый пара- |
метр, имеющий определëнный физический смысл.
Задача (2.90), (2.91) является однородной краевой задачей. Особенность еë в том, что она имеет нулевое ("тривиальное") решение y = 0, не представляющее физического интереса. Но
кроме того, имеются еще определëнные значения параметра λ, при которых задача имеет не равные тождественно нулю (нетривиальные) решения. Такие значения параметра λ называются
91
собственными значениями, а соответствующие им ненулевые решения y(x) называются собственными функциями. Сама же
задача отыскания собственных значений и собственных функций называется задачей на собственные значения или задачей Штурма-Лиувилля14.15. Эта задача представляет большой интерес для физики и для технических приложений. К ней, например, сводятся многие задачи устойчивости и колебаний упругих систем.
Покажем, прежде всего, что ненулевые решения задачи (2.90), (2.91) существуют при λ > 0. Доказательство от противного.
1) Пусть λ = 0. Тогда уравнение (2.90) примет вид y′′ = 0. Его общее решение, очевидно, будет y = C1x +C2. Определяя да-
лее произвольные постоянные из граничных условий (2.91), находим C1 = C2 = 0, y ≡ 0. Поэтому λ = 0 не является собствен-
ным значением.
2) Пусть λ < 0. Характеристическое уравнение, соответствующее уравнению (2.90) k 2 +λ = 0, имеет действительные и различные корни: k1,2 = ±−λ. Поэтому общее решение уравне-
ния (2.90) запишется в виде y = C1e−λx +C2e−−λx . Подставляя его
в граничные условия (2.91), получим систему однородных линейных алгебраических уравнений
C1 +C2 = 0,
e−λxC1 +e−−λxC2 = 0.
Определитель этой системы
1 e
1
−λx e−−λx
≠ 0, следователь-
но, C1 = C2 = 0. Поэтому задача (2.90), (2.91) имеет только нулевое решение y ≡ 0.
14Ш т у р м Жак Шарль Франсуа (29.09.1803 – 18.12.1855) –
французский математик,
15Л и у в и л л ь Жозеф (24.03.1809 – 08.09.1882) – французский математик.
92
3) Пусть λ > 0. Тогда корни характеристического уравнения будут мнимыми: k1,2 = ±−λ = ±λi. Поэтому частные решения
будут y1 = cos |
|
x, y1 = sin |
|
x, и общее |
решение уравнения |
||||
λ |
λ |
||||||||
(2.90) запишется в виде |
|
|
|
||||||
|
|
y = C1 cos |
|
x +C2 sin |
|
x. |
(2.92) |
||
|
|
λ |
λ |
Из первого граничного условия (при x = 0) получаем C1 = 0, а из второго граничного условия (при x = l) следует C2 sin λl = 0.
Если C2 = 0, то опять получим нулевое решение y = 0. По-
этому, чтобы существовало нетривиальное решение задачи, необходимо принять
sin λl = 0.
Из этого тригонометрического уравнения следует λl = nπ, где n – произвольное целое число. Поэтому собственные значения параметра для задачи (2.90), (2.91) будут
λ |
|
nπ |
2 |
(2.93) |
||
n |
= |
l |
|
(n = ±1,±2,...) |
||
|
|
|
|
|
(так как собственные значения будут различными для разных n, то им приписывается соответствующий индекс).
Соответствующие им собственные функции с точностью до
постоянного множителя C2 определяются по формуле |
|
||
y = sin |
nπx |
. |
(2.94) |
|
|||
|
l |
|
Рассмотрим в качестве примера классическую задачу об устойчивости прямолинейного упругого шарнирноëртопго стержня длиной l , нагруженного осевой сжимающей силой P, линия действия которой совпадает с продольной осью стержня. Постановка и решение этой задачи впервые были даны Л.Эйлером в 1765 году.
Пока сжимающая нагрузка мала, стержень сохраняет прямолинейную форму равновесия. Если увеличивать нагрузку, то при некотором значении сжимающей силы прямолинейная форма равновесия перестаëт быть устойчивой, и возникает близкая к ней искривлëнная форма равновесия (см. рис. I.14).
93
Уравнение равновесия отсечëнной части стержня длиной x может быть записано в виде [20]:
EIy′′+ Py = 0, |
(2.95) |
где y(x) - прогиб стержня, y′′(x)- приближëнное выражение для
кривизны изогнутой оси стержня, EI - изгибная жесткость стержня, Py - изгибающий момент в произвольном сечении стержня.
Полагая
λ = |
P |
, |
(2.96) |
|
EI |
||||
|
|
|
Рис. I.14
приведëм уравнение (2.95) к виду (2.90): y′′+λy = 0.
Так как прогибы на концах шарнирно опëртого стержня равны нулю, то граничные условия запишутся в виде (2.91): y(0) = 0,
y(l) = 0.
Поэтому с математической точки зрения задача об устойчивости упругого шарнирно опëртого стержня сведена к рассмотренной выше задаче на собственные значения (2.90), (2.91).
Тривиальное решение этой задачи (y ≡ 0) соответствует ис-
ходной прямолинейной форме равновесия. Нетривиальное решение, соответствующее искривлëнно й форме равновесия, существует при собственных значениях параметра λ, определяемых со-
гласно (2.93).
Приравнивая выражения (2.93) и (2.96), получим