- •Е.А.Коган
- •и естественно-научным дисциплинам
- •Москва 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).
- •Решение методом Бубнова
- •Как видим, ошибка, полученная при решении методом Бубнова в двучленном приближении, составляет сотые доли процента.
- •Решение методом наименьших квадратов
- •Решение методом коллокаций
- •В результате находим
13
F(x) = ∫M (x)dx, Ф(y) = ∫N (y)dy .
Конечное (не дифференциальное) соотношение (1.8) и является общим интегралом уравнения (1.7).
Пример. Решить уравнение ex dx + ln ydy = 0.
Очевидно, это уравнение с разделëнными переменными. Интегрируя его, получим
∫ex dx + ∫ln ydy = C.
Следовательно, общий интеграл уравнения будет ex + y(ln y −1) = C.
Дифференциальное уравнение вида
M1 (x)N1 (y)dx + M 2 (x)N2 (y)dy = 0, |
(1.9) |
в котором коэффициенты при дифференциалах распадаются на произведение сомножителей, каждый из которых зависит только от одной переменной, называется дифференциальным уравне-
нием с разделяющимися переменными.
Уравнение (1.9) делением обеих частей на произведение функций N1 (y)M 2 (x) ≠ 0 приводится к уравнению с разделëнны-
ми переменными |
|
|
|
||
|
M1 (x) dx + |
N2 (y)dy = 0, |
|||
|
M 2 (x) |
N1 (y) |
|||
общий интеграл которого |
|
|
|
||
∫ |
M1 (x) |
dx + |
∫ |
N2 (y) |
dy = C. |
|
|||||
|
|
||||
|
M 2 (x) |
|
N1 (y) |
Пример. Решить равнение (1 + y2 )xdx + (1 + x2 )ydy = 0. Разделяем переменные делением на выражение (1+ y2)(1+ x2) ≠ 0:
x |
dx + |
y |
dy = 0. |
|
1+ x2 |
1+ y2 |
|||
|
|
Интегрируем полученное уравнение с разделëнными переменн ы- ми
∫1 +xx2 dx + ∫1 +yy2 dy = C1.
Тогда
14
12ln(1+ x2)+ 12ln(1+ y2) = C1
Так как C1 - произвольная постоянная, принимая еë для упрощ е- ния полученного выражения в виде C1 = 12lnC, представим общий интеграл уравнения в виде (1 + x2 )(1 + y2 ) = C.
1.4. Однородные дифференциальные уравнения
Предварительно введем понятие однородной функции. Функция двух переменных f (x, y) называется однородной
n –го измерения, если при любом k > 0 справедливо равенство
f(kx,ky) = k n f (x, y).
Вчастности, если при изменении аргументов x и y в k раз
вид функции не меняется, то есть f (kx,ky) = f (x, y), то функция f (x, y) называется однородной нулевого измерения.
Соответственно дифференциальное уравнение первого порядка, разрешëнное относительно производной y′ = f (x, y), на-
зывается однородным, если правая часть f (x, y) есть однород-
ная функция нулевого измерения.
Если уравнение первого порядка записано в форме, содержащей дифференциалы M (x,y)dx + N (x,y)dy = 0, то оно будет
однородным, если M (x, y) и N (x, y) - однородные функции од-
ного и того же измерения.
Однородное дифференциальное уравнение подстановкой
|
|
|
|
|
|
|
t = |
y |
|
|
|
||
|
|
|
|
|
|
|
x |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
||
приводится к уравнению с разделëнными переменными. |
|
||||||||||||
Действительно, пусть уравнение |
y′ = f (x, y) |
однородное. |
Тогда |
||||||||||
f (kx,ky) = f (x, y). Полагая |
k =1/ x , получим |
f (kx,ky) = f (x, y) = |
|||||||||||
= f (1, y / x). |
Введем теперь |
новую функцию t = y / x. |
Тогда |
||||||||||
y = tx , |
y |
′ |
= t x +t, |
и уравнение примет вид |
t x +t = f (1,t) =ϕ(t) |
||||||||
|
|
|
′ |
|
|
|
|
|
|
′ |
|
||
или dt |
x =ϕ(t)−t. |
|
dt |
|
= dx - уравнение с разделëнными |
||||||||
ϕ(t) |
−t |
||||||||||||
dx |
|
|
|
|
|
|
|
x |
|
|
переменными.
15
Пример. Решить уравнение
(x + y)dx + (y − x)dy = 0.
Преобразуем уравнение к виду |
|
|
|
x + y |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
y′ = |
|
. |
|
|
|
|
|||||||
|
kx +ky |
|
x + y |
|
|
|
x − y |
|
|
|
|
|||||||||
Так как |
= |
|
, то исходное уравнение однородное. |
|||||||||||||||||
kx −ky |
x − y |
|||||||||||||||||||
|
|
y = tx |
и |
|
y |
′ |
= t x +t. |
|
||||||||||||
Полагаем t = y / x, |
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
|
|
||
Тогда уравнение примет вид |
|
|
|
|
1+t |
|
|
1+t2 |
1 + t 2 |
|||||||||||
′ |
|
x +tx 1+t |
|
dt |
|
|
|
|||||||||||||
|
|
= 1−t |
или x dx |
= |
1−t |
−t = |
1−t , |
xdt = 1 −t dx. |
||||||||||||
t x +t = x −tx |
||||||||||||||||||||
Разделив обе части уравнения на |
x |
1+t2 |
≠ |
0, приходим к уравне- |
||||||||||||||||
|
1−t |
|||||||||||||||||||
нию с разделëнными переменными |
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
1 −t |
|
dt = dx . |
|
|
|
||||||||
|
|
|
|
|
|
1 + t 2 |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
Интегрируя его, находим
arctg t − 12ln (1 + t 2 ) = ln x + C
или
arctg t −ln (x1+t2 ) = C.
Возвращаясь к старой переменной, получим общий интеграл исходного уравнения в виде
C= arctg xy − ln x2 + y2 .
1.5.Дифференциальные уравнения, приводящиеся к однородным
Дифференциальное уравнение вида
|
′ |
|
|
a1 x + b1 y + c1 |
|
||||
y |
= |
|
|
||||||
|
|||||||||
|
f a |
2 |
x + b y + c |
2 |
|
||||
|
|
|
|
|
2 |
|
называется приводящимся к однородному. В частности, к этому классу относится уравнение вида
16
y′ = a1x +b1 y +c1 . a2x +b2 y +c2
Некоторые из коэффициентов (но не одновременно
быть равны нулю.
Следует различать два случая:
(1.10)
c1 и c2) могут
1). Если определитель δ = |
|
a1 |
b1 |
|
≠ 0, |
то уравнение (1.10) |
|
|
|
||||||
|
|
a |
2 |
b |
|
|
|
|
|
|
2 |
|
|
|
|
приводится к однородному подстановкой |
|
||||||
x = u +α, y = v + β , |
(1.11) |
где постоянные a и β определяются из системы уравнений:
a1α + b1β + c1 = 0,
a2α + b2β + c2 = 0.
Заметим, что эту систему уравнений можно записать непосредственно по виду правой части уравнения (1.10), если заменить в ней x на α , y на β и приравнять числитель и знаменатель дроби ну-
лю.
Учитывая, что du = dx, dv = dy, следовательно, dydx = dudv , и
подставляя (1.11) в (1.10), получим однородное уравнение относительно новой функции v(u):
dv = |
a1 (u +α)+b1 (v + β)+c1 |
= |
a1u +b1v + a1α +b1β +c1 |
= |
a1u +b1v |
. |
||||||||||
du |
a |
(u +α)+b |
(v + β)+c |
2 |
|
a |
u +b v + a α +b |
β +c |
2 |
|
a |
u +b v |
|
|||
2 |
2 |
|
|
2 |
|
2 |
2 |
2 |
|
2 |
2 |
|
||||
Полагая далее |
|
|
t = v |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
, |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
u |
|
|
|
|
|
|
|
|
|
приводим последнее уравнение к уравнению с разделяющимися переменными.
2) Если определитель
δ = |
|
a1 |
b1 |
|
|
|
= 0, |
||||
|
|
a |
2 |
b |
|
|
|
|
2 |
|
то уравнение (1.10) сразу приводится к уравнению с разделëн ными переменными заменой u = a1x +b1 y.