книги / Обыкновенные дифференциальные уравнения
..pdfоткуда, |
допуская, что значения yn −1 D при |
|
x − x0 |
|
≤ h , получим |
|||
|
|
|||||||
|
yn − y0 |
|
≤ |
Mh ≤ b , т.е. n-е приближение также не выходит из той же области. |
||||
|
|
Таким образом, методом полной индукции доказано, что ни одно из последовательных приближений не выйдет из области D, если
x0 − h ≤ x ≤ x0 + h .
Теперь нужно показать, что существует предел последовательности yn
,
Y ( x) = lim yn ( x) |
, |
(1.47) |
n → ∞ |
и что предельная функция удовлетворяет уравнению (1.42) и начальным условиям. Для того чтобы показать существование предела (1.47) последовательности {yn}, достаточно доказать сходимость ряда:
y0 + (y1 − y0 ) + (y2 − y1 ) + ... + (yn − yn−1 ) + ... ; |
(1.48) |
действительно, n-я частная сумма этого ряда Sn = yn , и если мы докажем, что ряд (1.48) сходится к функции Y (x) , то этим самым и будет доказано соотношение (1.47).
Оценим абсолютные величины членов ряда (1.18):
|
|
|
|
|
|
|
|
|
|
|
|
y1 − y0 |
|
= |
|
∫x f (x, y0 )dx |
|
≤ M |
|
x1 − x0 |
|
; |
|
|
|
|
|
|
|
|
(1.49, a) |
|||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
далее |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
y 2 − y1 |
|
= |
|
|
∫ {f (x , y1 ) − f (x , y 0 )}dx |
|
|
≤ M |
∫ |
|
|
|
|
f |
(x , y1 ) − f (x , y 0 ) |
|
dx |
. |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
x 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
На основании условия Липшица (1.44) подынтегральная функция |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
удовлетворяет неравенству |
|
|
|
|
|
f (x, y1 ) − f (x, y0 ) |
|
≤ N |
|
y1 − y0 |
|
|
; подставляя сюда |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
только что найденную оценку для |
|
|
|
|
y1 − y 0 |
|
|
|
|
, находим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
MN |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
y2 − y1 |
|
|
≤ N |
∫ M |
|
x − x0 |
|
dx |
= |
|
|
x − x0 |
|
|
|
|
|
|
|
(1.49, б) |
|||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
1 2 |
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
Аналогично получим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
x |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
y3 − y2 |
|
≤ |
∫ |
|
f (x, y2 ) − f |
(x, y1 ) |
|
dx |
≤ |
∫ N |
|
|
y2 − y1 |
|
dx |
≤ |
MN |
|
|
|
|
∫ (x − x0 )2 dx |
, |
||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 2 |
|
x0 |
|
|||||||||||||||||||
и, наконец, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y3 − y2 |
|
≤ |
MN 3 |
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(1,49, в) |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x − x0 |
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Докажем, что для любого целого положительного n справедливо неравенство
31
yn − yn −1 |
|
≤ |
MN n −1 |
|
|
x − x0 |
|
n |
|
|
|
|
|
||||||
|
n! |
|
|
|
. |
(1,49, n) |
|||
|
|
|
Для этого применим метод полной индукции: допустив, что неравенство (1.49, n) справедливо для n, докажем его справедливость для n + 1. Мы имеем:
|
|
|
x |
|
x |
||||||||
yn +1 − yn |
|
≤ |
∫ |
|
f (x, yn ) − f (x, yn −1 ) |
|
dx |
≤ N |
∫ |
|
yn − yn −1 |
|
dx |
|
|
|
|
|
|||||||||
|
|
|
x0 |
|
x0 |
(последнее неравенство получаем на основании условия Липшица).
Подставив в последний интеграл вместо |
|
y n − y n −1 |
|
выражение (1.49n) (отчего |
||||||||||||||||||||
|
|
|||||||||||||||||||||||
неравенство усилится), получаем: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
≤ N |
MN |
n −1 |
|
x |
|
|
|
n dx |
|
MN |
n |
|
n +1 , |
||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
yn +1 − yn |
|
|
|
∫ |
|
x − x0 |
|
= |
|
|
|
x − x0 |
|
||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||
n! |
(n + |
1)! |
||||||||||||||||||||||
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
т.е. неравенство (1.49n) справедливо и для n + 1. А поскольку данное неравенство оправдывается для n = 2 и n = 3, то оно верно для каждого натурального числа n. Заменяя еще x − x 0 его наибольшим допустимым
значением h, мы приходим к заключению, что каждый член ряда (1.48) не превосходит соответствующего члена числового ряда с положительными членами:
|
|
|
|
|
|
|
|
|
Nh 2 |
|
|
N |
2 h3 |
|
|
N n −1hn |
|
|
||||||
|
y0 |
|
|
+ Mh + |
M |
|
|
|
+ M |
|
|
+ ... + M |
|
|
|
|
|
+ ... . |
(1.50) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
2! |
|
|
|
3! |
|
|
|
|
n! |
|||||||||
Применяя к последнему ряду признак Даламбера, находим: |
|
|||||||||||||||||||||||
|
lim |
|
u |
n+1 |
= lim |
|
MN nhn+1 n! |
= lim |
Nh |
|
= 0 |
< 1. |
|
|||||||||||
|
|
un |
|
|
|
|
|
|
n + 1 |
|
||||||||||||||
n→ ∞ |
|
|
n→ ∞ (n + 1)!MN n−1hn |
n→ |
∞ |
|
|
|
|
Итак, ряд (1.50) сходится. Поскольку все члены ряда (1.48) по абсолютной величине меньше членов числового ряда (1.50), то ряд (1.48) не только сходится, но и, в силу критерия Вейерштрасса, сходится равномерно для всех х, удовлетворяющих условию x − x0 ≤ h. Каждый член ряда (1.48) есть непрерывная функция верхнего предела. Следовательно,
Y = lim y n
n → ∞
существует и является непрерывной функцией от x.
Докажем, что полученная таким образом функция Y (x) удовлетворяет данному уравнению (очевидно, что начальное условие удовлетворено, так
как Y ( x0 ) = lim yn ( x) = y0 ).
n→ ∞
Возьмем равенство (1.46n):
x
yn (x) = y0 + ∫ f (x, yn−1 )dx
x0
и перейдем к пределу при n → ∞ .
32
Благодаря равномерной непрерывности функции f (x, y ) |
по y в области |
D мы для любого заренее заданного положительного числа ε |
можем найти |
такое δ > 0, что неравенство |
|
| f (x, y′) − f (x, y′′) |< ε |
|
будет выполнено для тех пар точек (x, y′) и (x, y′) области D, для которых выполняется неравенство y′ − y′′ < δ (в силу условия Липшица достаточно
взять δ = ε ). Далее из равномерности стремления последовательности yn
N
к пределу возникает возможность для выбранного δ так подобрать натуральное число n0, чтобы при n – 1 > n0 для всех значений x в интервале (x0 − h, x0 + h) имело место неравенство
|
|
|
|
|
yn −1 ( x) − Y ( x) |
|
< δ . |
|
|||||||||
|
|
|
|
|
|
|
|||||||||||
Сопоставляя оба эти неравенства, мы получаем при n −1 > n0 : |
|
||||||||||||||||
|
|
|
|
f ( x, y n −1 ( x)) − f ( x,Y ( x)) |
|
< ε . |
|
||||||||||
|
|
|
|
|
|||||||||||||
Отсюда следует: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
x |
|
|
x |
|
|
x |
|
|
||||||||
|
|
|
|
|
|
|
|||||||||||
|
∫ f ( x, yn−1)dx − ∫ f (x,Y )dx |
≤ |
∫ |
|
f ( x, yn−1) − f ( x,Y ) |
|
dx |
< ε h . |
|
||||||||
|
|
|
|
||||||||||||||
|
x0 |
|
|
x0 |
|
|
x0 |
|
|
||||||||
Пользуясь произвольностью числа ε, находим |
|
||||||||||||||||
|
|
|
|
lim |
∫x f (x, yn −1 )dx = ∫x f (x,Y )dx . |
|
|||||||||||
|
|
n → ∞ |
x0 |
|
|
|
x0 |
|
|||||||||
Таким образом, переходя к пределу в равенстве (1.46n) при |
n → ∞ , |
||||||||||||||||
получаем тождество |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Y = y0 + ∫x f (x,Y )dx, |
(1.51) |
|||||||||||
|
|
|
|
|
|
x0 |
|
||||||||||
т.е. функция Y (x) |
удовлетворяет интегральному уравнению (1.45). Далее |
функция Y допускает производную по x, так как в правой части тождества (1.51) стоит интеграл от непрерывной функции, допускающий производную по верхнему пределу. Дифференцируя (1.51), получаем
dY = f (x,Y ), dx
т.е. Y(x) удовлетворяет данному уравнению (1.42). Первая часть теоремы доказана.
2. Единственность решения задачи Коши (1.42), (1.43).
Допустим, что кроме решения Y(x) существует еще другое решение Z(x), удовлетворяющее тому же начальному условию: Z (x0 ) = y0 . Без ограничения общности можно предположить, что значения x, для которых
33
Y (x) ≠ Z(x) , находятся вправо от x0 в любой близости от точки x0 (иначе за точку x0 мы взяли бы ту точку, в любой близости которой значения Y (x) и Z(x) перестают быть равными, или заменили бы x на – x ).
Покажем, что наше допущение приводит к противоречию. Возьмем
некоторое малое постоянное число |
ε >0 ; по предположению, в замкнутом |
||||
интервале x0 ≤ x ≤ |
x0 + ε не всюду |
Y = Z ; следовательно, положительная |
|||
функция |
|
Y ( x) − Z ( x) |
|
достигает в этом интервале в некоторой точке ξ своего |
|
|
|
наибольшего значения θ >0, причем не может быть ξ = x0 , так как при x = x0
обе функции Y |
и Z |
равны между собой. Мы имеем: |
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
Y ( x) = y0 + ∫x f (x,Y )dx, |
Z ( x) = y0 + ∫x f (x, Z )dx. |
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
||||
|
|
Вычтем |
оба |
тождества |
одно |
из |
другого, придавая |
x значение ξ |
||||||||||||||||||||||||
и используя при оценке разности условие Липшица: |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
ξ |
|
|
|
|
|
|
ξ |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
Y (ξ ) − Z (ξ ) |
|
= θ ≤ ∫ |
|
f ( x,Y ) − f ( x, Z ) |
|
dx ≤ N ∫ |
|
Y ( x) − Z ( x) |
|
dx. |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
||||
|
|
Мы только усилим последнее неравенство, если под интегралом заменим |
||||||||||||||||||||||||||||||
разность |
|
|
Y ( x) − Z ( x) |
|
|
|
ее |
наибольшим |
значением |
θ |
и |
если |
интеграл |
|||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||||||||
распространим |
|
на |
|
промежуток (x0 , x0 + ε ) |
вместо промежутка (x0 ,ξ ). Получим |
|||||||||||||||||||||||||||
|
|
|||||||||||||||||||||||||||||||
θ ≤ εθ |
откуда |
по |
|
предположению |
θ > 0 |
, |
находим |
что |
1 ≤ |
Nε , |
что |
|||||||||||||||||||||
|
N , |
|
|
|
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
|||||
противоречиво, |
поскольку ε можно выбрать сколь угодно малым, |
например |
||||||||||||||||||||||||||||||
ε < |
1 |
|
Следовательно, допущение существования второго решения z, |
|||||||||||||||||||||||||||||
N . |
||||||||||||||||||||||||||||||||
|
||||||||||||||||||||||||||||||||
удовлетворяющего |
тому |
же |
начальному |
условию, |
что |
|
|
|
Y, |
приводит |
||||||||||||||||||||||
кпротиворечию. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, при выполнении условий теоремы существует единственное решение задачи Коши (1.42), (1.43).
Продолжение решения задачи Коши (1.42), (1.43)
Существование решения задачи Коши (1.42), (1.43) доказано только |
|||||||||||
для отрезка I = [x0 − h, x0 + h]. |
Покажем, |
что найденное решение может быть |
|||||||||
продолжено на всю область Q, |
где функция |
f (x, y ) |
определена, непрерывна |
||||||||
и удовлетворяет условию Липшица по y. |
|
|
|
||||||||
Действительно, пусть |
|
x |
(1) |
= x |
0 |
|
+ h |
|
|
|
|
|
|
0 |
|
|
, |
|
|
|
|||
y(x0(1) )= y0(1) , причем (x0(1) , y0(1) ) |
Q . |
||||||||||
Найдем прямоугольник D : |
|
x − x(1) |
|
≤ |
|
a , |
|
|
|
||
|
|
|
|
|
|
||||||
1 |
|
|
|
0 |
|
|
1 |
|
|
|
|
|
|
|
y − y(1) |
|
≤ b : D |
Q. |
|
||||
|
|
|
|
|
|||||||
|
|
|
|
0 |
|
|
1 |
1 |
|
|
34
Пусть M1 = max f (x, y) , тогда по теореме Коши – Пикара существует
|
|
D1 |
y(x |
(1) )= y |
|
|
|
|
|
|
|
|
|
|
(1) |
|
|
|
(1.52) |
||
|
|
|
|
|
0 |
0 , |
I1 |
= x0 |
− h1; x0 |
|
являющееся решением задачи Коши на отрезке |
+ h1 , где |
|||||||||
|
|
|
|
|
|
|
|
(1) |
(1) |
|
|
|
b1 |
|
|
|
|
|
|
|
|
h1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
= min a1, |
M1 |
. |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
||
|
Середина отрезка I1 совпадает с концом I в этой точке? оба решения |
совпадают и равны y0(1) . В силу единственности оба решения будут также
совпадать в общей части отрезков I и I1. Таким образом, найденное решение является продолжением ранее полученного на отрезке I решения на вторую половину отрезка I1.
Пусть теперь x0(2) = x0(1) + h1, y(x0(2) )= y0(2) и точка (x0(2) , y0(2) ) Q. Мы можем по начальному условию определить решение на отрезке I2, имеющего общую половину с I1, т.е. продолжить решение на вторую половину отрезка I2 и т.д.
Аналогичное рассуждение может быть проведено и для убывающих значений x.
Таким образом, можно продолжить решение до границы области Q. Пример. С помощью метода последовательного приближения найти
решение задачи Коши
y′ = y, y(0) =1.
x
Эта задача Коши эквивалентна уравнению y(x) =1+ ∫ y(τ )dτ .
0
На его основе строим последовательные приближения: y0 =1,
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y1 =1 + ∫ dτ = 1+ x, |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
x |
2 |
|
|
|
|
|
|
|
y2 =1 + ∫(1 + τ )dτ = 1+ x+ |
|
, |
|
|
|
|
|
|
|||||||||
2 |
|
|
|
|
|
|
|||||||||||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
||
y |
=1 |
+ |
1 + τ + |
τ |
|
dτ = |
1+ |
|
x+ |
x |
+ |
x |
|
, |
|||
x |
2 |
|
3 |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
||||
3 |
|
|
∫0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
3! |
|
|||||
|
|
|
|
2 |
|
|
|
|
|
|
|
...
yn =1 + x + x2 + ... + xn , 2! n!
....
Очевидно, что lim yn (x) = ex = y(x) .
n→ ∞
35
Теорема существования решения задачи Коши (1.42), (1.43) (Теорема Пеано)
|
|
|
Пусть |
|
функция |
f (x , y ) |
непрерывна |
|
в |
прямоугольнике |
||||||||||
D : |
|
|
x − x0 |
|
≤ a, |
|
y − y0 |
|
≤ b . Пусть, далее, |
M = max |
|
f (x, y) |
|
. Тогда задача Коши |
||||||
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D |
|
|
|
|
|
|
(1.42), (1.43) |
|
имеет по |
крайней мере |
одно |
решение, |
определенное для |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
x |
[x0 − h, x0 + h], где h = min a, |
|
. |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
M |
|
|
|
|
|
|
|
|
|
|
|
|
Для доказательства |
этой теоремы |
мы |
применим |
теорему Арцела и |
принцип неподвижной точки Шаудера. Введем предварительно следующие определения.
Определение. Последовательность функций
|
|
|
|
|
f1(x), f2 (x), ..., fn (x), ... |
|
||||
называется равномерно ограниченной на [a,b], если существует число K > 0 |
||||||||||
такое, что |
для всех |
натуральных |
n |
и любого |
x [a,b] справедливо |
|||||
неравенство |
|
|
|
|
|
|
|
≤ |
|
|
|
|
|
|
|
|
f n ( x) |
|
K. |
|
|
|
|
|
|
|
|
|
|
|||
Определение. Представленная выше последовательность непрерывных |
||||||||||
на [a,b] функций называется равностепенно непрерывной, если для любого |
||||||||||
ε >0 существует δ (ε ) > 0 |
такое, что для любых x1, x2 |
[a,b], удовлетворяющих |
||||||||
неравенству |
|
x1 − x2 |
|
< δ , и любого натурального n выполняется неравенство |
||||||
|
|
f n ( x1 ) − f n ( x2 ) < ε .
Теорема Арцела. Из всякой последовательности равномерно ограниченных и равностепенно непрерывных на [a,b] функций можно выбрать равномерно сходящуюся подпоследовательность.
Эта теорема является критерием компактности в пространстве C [a, b] непрерывных на [a, b] функций.
Ее доказательство приведено, например, в [11].
Теорема Шаудера. Если вполне непрерывный оператор F отображает ограниченное замкнутое выпуклое множество B банахова пространства Y на свою часть, то существует неподвижная точка этого отображения, т.е. такая точка y B , что Fy = y.
Доказательство этой теоремы, а также соответствующие определения также приведены в [11].
Перейдем теперь к доказательству теоремы Пеано, воспользовавшись теоремами Арцела и Шаудера.
Уравнение (1.52) вместе с начальным условием эквивалентно интегральному уравнению
36
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y ( x) = y0 + ∫x f (τ , y (τ ))dτ |
. |
|
(1.53) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x0 |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рассмотрим оператор F, определенный равенством |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Fy = y0 + ∫x f (τ , y(τ ))dτ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
на шаре B : |
|
|
|
y − y0 |
|
|
|
≤ b пространства C[x0 − h, x0 + h]. Покажем, что оператор F |
||||||||||
|
|
|
|
|||||||||||||||
вполне непрерывен на этом шаре. |
{yn (x)}, |
|
|
|||||||||||||||
|
|
Прежде |
всего, если последовательность |
принадлежащая шару |
||||||||||||||
|
|
y − y 0 |
|
|
|
≤ b, равномерно сходится к функции |
y(x) , очевидно, принадлежащей |
|||||||||||
|
|
|
|
|||||||||||||||
|
|
тому же |
|
шару, то в силу непрерывности функции f (x, y ) |
имеем |
f (x, y (x)) → |
||||||||||||
→ f (x, y(x)) |
равномерно на [x0 − h, x0 + h] . Отсюда вследствие |
n |
||||||||||||||||
возможности |
предельного перехода под знаком интеграла при равномерной сходимости
получаем, что Fyn → |
Fy , т.е. оператор F непрерывеннашаре |
|
|
|
y − y 0 |
|
|
|
≤ b . |
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||
Далее для любого элемента y(x) шара |
|
|
|
|
y − y 0 |
|
|
|
≤ b имеем |
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
Fy( x) |
|
≤ |
|
y0 |
|
+ |
|
|
|
∫x f (τ , y(τ ))dτ |
|
|
≤ |
|
|
y0 |
|
+ Mh. |
(1.54) |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Если x1 и x2 |
– две любые точки отрезка [x0 − h, x0 + h] , то |
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
Fy( x2 ) − Fy( x1 ) |
|
≤ |
|
|
|
|
∫2 f (τ , y(τ ))dτ |
≤ |
M |
|
x2 − x1 |
|
. |
(1.55) |
|||||||||||||||||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Неравенства (1.54) и (1.55) в силу теоремы Арцела показывают, что |
||||||||||||||||||||||||||||||||||||||||||||||||
оператор F преобразует |
шар |
|
|
|
y − y0 |
|
|
|
≤ b |
|
|
|
в |
относительно |
|
компактное |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
множество. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Покажем, наконец, что оператор F преобразует этот шар в себя. Действительно,
Fy ( x) − y |
|
= max |
x |
f (τ , y(τ ))dτ |
≤ |
Mh ≤ M |
b |
= b |
|
0 |
∫ |
||||||||
|
|||||||||
|
x |
|
|
|
M |
. |
|||
|
|
x0 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
Таким образом, оператор F удовлетворяет всем условиям теоремы Шаудера. Поэтому существует неподвижная точка этого оператора, т.е. такая функция y(x) , что
y( x) ≡ y0 + x∫ f (τ , y(τ ))dτ .
x0
Это равенство равносильно двум равенствам:
dy(x) = f (x, y(x)), y(x0 ) = y0, dx
чем и завершается доказательство теоремы Пеано.
37
Глава 2. Уравнения первого порядка, не разрешенные относительно производной
2.1. Уравнения первого порядка n-й степени |
|
Рассмотрим уравнение первого порядка |
|
F (x, y, y′) = 0 . |
(2.1) |
Будем предполагать, что левая часть уравнения (2.1) является |
|
многочленом относительно y ′ степени n, причем коэффициенты |
этого |
многочлена в некоторой области D есть непрерывные функции от x и y, допускающие по y непрерывные частные производные.
F (x, y, y′) ≡ |
A (x, y )y′n + A |
(x, y )y′n −1 + ... + A (x, y) = 0 . |
(2.2) |
|
|
n |
n −1 |
0 |
|
Уравнение вида (2.2) называется уравнением первого порядка n-й |
||||
степени относительно |
y ′. Затем допустим, что в области D коэффициент |
|||
при старшей степени |
y ′, т.е. |
An (x, y ) , нигде не обращается в нуль. Тогда по |
основной теореме алгебры уравнение (2.2) для всякой пары значений x, y в рассматриваемой области имеет n решений y ′ (действительных или мнимых). По теореме о неявных функциях каждое из действительных решений является непрерывной функцией от x и y и имеет конечную частную
производную |
|
∂ y′ |
, если для рассматриваемых значений x, y, y′ |
имеем |
∂ F |
≠ |
0 |
|||||
|
|
|
||||||||||
|
|
|
|
|
|
|
∂ y |
|
∂ y′ |
|
||
(т.е. если для данных значений x и y уравнение (2.2) относительно |
y ′ |
не |
||||||||||
имеет |
кратного корня). Таким образом, во всякой области, в которой |
|||||||||||
|
|
∂ F |
|
≥ |
a > 0 , |
разрешенное относительно y ′ уравнение (2.2) |
удовлетворяет |
|||||
|
|
|
||||||||||
|
|
∂ y′ |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
||
условию Липшица. |
|
|
|
|
||||||||
|
|
|
|
|
Итак, пусть в области D уравнение (2.1) имеет k действительных |
|||||||
решений ( k ≤ |
|
n ): |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
y′ = f1 (x, y ), y′ = f 2 (x, y ), ..., y′ = f k (x, y ). |
|
(2.3) |
Каждое из уравнений (2.3) есть уравнение такого вида, для которого доказано существование единственного решения (теорема Коши – Пикара), проходящего через точку (x0, y0); следовательно, уравнение (2.1) допускает в точности k интегральных кривых, проходящих через данную точку (x0, y0) области D.
Пример 1.
y′2 + yy′ − x2 − xy = 0.
38
Разрешая относительно y ′ (или разлагая левую часть на множители), получаем два уравнения: 1) y′ = x , 2) y′ = − y − x ; обе правые части однозначны и непрерывны во всей плоскости xy; соответственно, имеем два решения
1) y = x2 + C ; 2) y = Ce−x − x +1.
2
Пример 2.
y′2 + y2 −1 = 0 .
|
|
Разрешая относительно |
y′, |
|
|
: |
|
y′ |
= ± |
1 |
− |
y2 |
; |
в полосе |
− |
1 |
< |
y |
< |
1 |
||||||||||||
|
|
|
имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
плоскости xOy оба рассматриваемых значения |
y ′ действительны и различны, |
|||||||||||||||||||||||||||||||
при |
|
y |
|
> 1 значения |
y ′ |
мнимы, |
и при |
y = ±1 |
|
два корня уравнения для |
y ′ |
|||||||||||||||||||||
|
|
|
||||||||||||||||||||||||||||||
совпадают. Рассмотрим |
полосу |
−1 < y < +1; |
через |
|
каждую точку (x0 , |
y0 ) |
||||||||||||||||||||||||||
проходят два решения |
y′ |
: |
для одного из них значение производной в этой точке |
|||||||||||||||||||||||||||||
0 |
|
|||||||||||||||||||||||||||||||
y′ = + 1 − y2 , для другого |
y′ = − |
1 − y2 |
. В |
|
данном |
|
случае |
в уравнении, |
||||||||||||||||||||||||
0 |
0 |
|
|
|
|
0 |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
разрешенном относительно |
y′, |
переменные разделяются и мы имеем |
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
|
: |
|
|
|
|
|
|
||
|
|
|
|
|
1) |
|
|
dy |
|
= dx |
; 2) |
|
|
dy |
|
= −dx . |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
1 − y2 |
|
|
1− y2 |
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
Получаем решения: 1) arcsin y = x + C , 2) |
arcsin y = −x + C . Откуда |
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
π |
|
|
|
|
|
|
π |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y = sin(x + C), |
− |
|
|
< x + C < |
|
|
; |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
2 |
2 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
y = sin(− x + C), |
|
|
|
π |
|
|
|
|
|
|
|
π |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
− |
|
|
< −x + C < |
|
|
. |
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
2 |
2 |
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Возьмем |
начальные условия |
|
|
|
x = x0 , y = y0 ; для |
|
определения |
||||||
произвольного постоянного C получаем уравнения: |
|
|
|
|
|||||||||
y0 = sin (± x0 + C ), |
± x0 + C = arcsin y0 , |
C = ± x0 + arcsin y0 . |
|||||||||||
Решения, удовлетворяющие нашему начальному условию: |
|||||||||||||
y = sin[(x − x0 ) + arcsin y0 ], |
− |
π |
|
< x − x0 |
+ arcsin y0 < |
π |
; |
|
|
||||
|
2 |
||||||||||||
|
|
|
|
2 |
|
|
|
|
|
|
|
||
y = sin[− (x − x0 ) + arcsin y0 |
], |
− |
π |
< −x + x0 + arcsin y0 < |
π |
. |
|||||||
|
|
||||||||||||
|
|
|
|
2 |
|
|
2 |
|
|||||
Первое |
семейство |
решений |
представляет |
семейство |
дуг синусоид, |
взятых от точки минимума до точки максимума (исключая концы); все кривые получаются из одной, например
y = sin x, − |
π |
< x < |
π |
, |
|
|
|||
2 |
2 |
|
39
переносом параллельно оси Ox; аналогично, второе семейство есть семейство дуг синусоид, взятых от максимума до минимума; они получаются переносом параллельно оси Ox из кривой
y = sin(− x), − |
π |
< x < |
π |
. |
|
|
|||
2 |
2 |
|
Через каждуюточку полосы −1 < y < 1 проходитпооднойкривойкаждого семейства – возрастающая криваяпервогосемейства иубывающая– второго.
Расширим теперь область, в которой рассматриваются кривые, добавив обе граничные прямые y = ±1 , на которых условия Липшица не выполняются
(у нас |
f (x, y) = ± 1− y2 ; производная ∂ f = |
|
y |
2 |
становится бесконечной |
||||||
|
∂ y |
1 |
− y |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|||
при y → |
±1 ), но правые части уравнений остаются непрерывными. Теперь нам |
||||||||||
|
|
|
|
|
|
|
π |
π |
|
||
нет необходимости ограничивать изменение аргумента пределами |
− |
|
, |
|
. |
||||||
2 |
2 |
||||||||||
|
|
|
|
|
|
|
|
|
Замечая, далее, что sin (− x + C ) = sin (x + π + C ), мы можем представить общее
решение одной формулой:
y = sin (x + C ), − ∞ < x < +∞ .
Это семейство полных синусоид, получаемых из одной параллельным переносом вдоль оси x; чтобы получить каждую один раз, надо ограничить область изменения произвольного постоянного С, например 0 ≤ C≤ π2 . В точках
максимума |
одной полной синусоиды кривая первого семейства |
переходит |
вкривую |
второго, в точке минимума имеет место обратный |
переход |
(см. пунктирные кривые на рис. 2.1). Но эти переходы совершаются на прямых y = ±1 , где условия Липшица не выполнены; поэтому, применяя теорему Коши– Пикара, мы не можем связать убывающую и возрастающую дугу синусоиды в одну интегральную кривую; теорема существования дает две отдельные ветви, проходящиечерезкаждуюточкуполосы(сплошныелиниинарис. 2.1).
y
I
0 |
x |
|
II |
Рис. 2.1
40