Конспект лекций по обыкновенным дифференциальным уравнениям. Ч. 1. Учебное пособие
.pdf
|
|
y(m) (x) y(m 1) |
|
( An)m 1 M |
(x x |
)m |
|
( An)m 1 M |
m |
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
, |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
i |
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m! |
|
|
|
|
|
m! |
|
|
|||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y(m) |
(x) y(m 1) |
|
|
|
M (nA )m |
|
m 1,2,... |
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
(x) |
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
i |
|
|
|
|
|
|
i |
|
|
|
nA |
m! |
|
|
|
|
|
|
|
|
|
|||||||||||
Отсюда следует, что ряд |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
0 |
|
(1) |
|
|
|
0 |
|
(2) |
(1) |
|
|
|
|
|
|
... |
(m) |
(m 1) |
|
|
... |
|||||||||||||||
yi |
yi |
(x) yi |
yi (x) yi |
|
|
(x) |
yi |
(x) yi |
|
(x) |
||||||||||||||||||||||||||
мажорируется сходящимся числовым рядом |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
yi0 |
|
|
M |
|
(nA ) |
|
|
|
|
yi0 |
|
|
M |
|
enA 1 |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
nA m 1 |
m! |
|
|
|
|
|
|
|
|
|
|
|
nA |
|
|
|
|
|
|
|||||||
(9)
(10)
и потому абсолютно и равномерно сходится в интервале x0 , x0 . |
|
Т.к. сумма |
||||||||||||||||||
его |
первых |
(m+1) |
членов |
|
равна |
y(m) (x) , то |
отсюда |
следует, |
что |
для |
||||||||||
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
x x0 , x0 при m m0i будет выполнено неравенство |
|
~ |
(m) |
(x) |
|
, |
где |
|||||||||||||
|
|
|||||||||||||||||||
|
yi (x) yi |
|
||||||||||||||||||
~ |
|
|
|
|
|
|
(m) |
(x) при m . Теперь нужно взять |
||||||||||||
yi (x) |
означает сумму ряда (10), т.е. предел yi |
|
|
|||||||||||||||||
m0 max(m01,...,m0n ) . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Докажем теперь, что система функций (7) удовлетворяет системе |
|||||||||||||||||||
интегральных уравнений (4), т.е. равенства |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
~ |
|
x |
|
~ |
~ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
i 1,...,n |
|
|
(11) |
|
||||||||
|
|
yi (x) yi |
fi t, y1 (t),...,yn (t) dt, |
|
|
|
||||||||||||||
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
являются тождествами. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
~ |
(m) |
|
|
(m) |
(x). |
|
|
|
|
|
|
|
|
|
|
Для доказательства положим yi (x) yi |
|
(x) Ri |
|
|
|
|
|
|
|
|
|||||||||
С использованием соотношения (5) будем иметь |
|
|
|
|
|
|
|
|
|
|
||||||||||
|
~ |
m 1 |
|
|
x |
~ |
|
|
|
|
~ |
|
|
|
(t) dt . |
|
|
|
|
|
|
0 |
|
m |
|
|
|
|
|
m |
|
|
|
|
|
||||||
|
yi (x) Ri |
(x) yi |
fi t, y1 |
(t) R1 |
(t),...,yn (t) Rn |
|
|
|
|
|
||||||||||
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Вычтя из обеих частей этого равенства
31
|
|
|
|
|
|
|
|
|
x |
~ |
~ |
|
~ |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
fi t, y1 |
(t), y2 |
(t),...,yn (t) dt , |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
|
|
|
|
x |
~ |
~ |
|
|
|
|
|
|
||||
получим |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
yi (x) |
yi |
|
fi t, y1 (t),...,yn (t) dt |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(t)] dt. |
|
|
|
|
m 1 |
(x) |
|
|
|
|
|
m |
|
|
|
|
m |
|
|
|
||||||||
|
Ri |
|
fi [t, y1 (t) R1 (t),..., yn (t) Rn (t)] |
fi [t, y1(t),..., yn |
||||||||||||||||||||||
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Однако, в силу условия Липшица и (8) имеем при m m0 |
|
|
||||||||||||||||||||||||
|
|
|
|
|
~ |
|
|
m |
~ |
|
|
|
m |
|
|
|
|
~ |
|
~ |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
fi t, y1 (t) R1 (t),...,yn (t) Rn |
(t) fi t, y1 (t),...,yn (t) |
||||||||||||||||||||||
|
|
|
A |
Rm (t) |
|
... |
|
Rm (t) |
|
nA . |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
1 |
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
(t) dt |
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Поэтому |
|
|
|
(x) |
0 |
|
fi |
|
|
|
|
|
|
|
|
m 1 |
(x) |
|
nA dt (1 nA ) . |
|||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
yi |
yi |
|
t, y1(t),.., yn |
|
Ri |
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x0 |
|
|
Т.к. величина справа может быть сделана сколь угодно малой, а величина слева неотрицательна, то равенства (11) тождественно выполняются. Таким образом мы показали, что решение системы интегральных уравнений (4), а, следовательно, и
решение задачи (1)-(2), при выполнении условий теоремы всегда существует и имеет вид (11).
Осталось доказать, что система функций (7) является единственным решением системы (4). Предположим, что функции
|
|
|
|
|
|
|
|
|
y* (x), y* (x),...,y* (x) |
|
|
|
(12) |
|
|||||
|
|
|
|
|
|
|
1 |
2 |
n |
|
|
|
|
|
|||||
также образуют решение системы (4). Ясно, что функции y* (x) |
непрерывны, |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
т.к. таковы правые |
части уравнений |
системы (4); |
кроме |
того, |
y* (x ) y0 . |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
0 |
i |
Поэтому, |
заменяя (если |
|
это |
необходимо) интервал |
x0 , x0 |
|
на |
меньший |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
интервал x0 , x0 , мы можем считать, что |
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
y* (x) y0 |
|
b, |
i 1,...,n |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
i |
i |
|
|
|
|
|
|
|
|
||
|
* |
* |
|
|
|
* |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(x) |
|
|
D . Поэтому справедливы неравенства (6), т.е. |
|
||||||||||||||
т.е. x, y1 |
(x),.., yn |
D |
|
||||||||||||||||
32
|
|
|
|
|
|
|
f |
(x, y*,..., y* ) |
M , |
i 1,..., n , |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
i |
1 |
|
n |
|
|
|
|
|
|
|
|
а поскольку по предположению |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
yi* (x) yi0 |
fi t, y1* (t),...,yn* (t) dt , |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
то с учетом (6) и (6а) имеем |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
~ |
|
|
|
|
|
x |
~ |
|
~ |
(t) fi |
|
|
(t) dt |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
* |
|
|
|
* |
* |
|
|||||||||||
|
|
yi |
yi |
fi t, y1 |
(t),..., yn |
t, y1 (t),..., yn |
|||||||||||||
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
2M (x x0 ). |
|
|
|
|
|
|
|
|
|
|
||||||||
Далее, т.к. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
~ |
|
|
|
x |
|
|
~ |
|
~ |
|
|
|
|
(t) |
|
|||
|
|
* |
|
|
|
|
|
* |
|
* |
dt |
||||||||
|
yi yi |
|
fi t, y1 (t),...,yn (t) fi t, y1 |
(t),...,yn |
|||||||||||||||
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(6a)
(13)
то, воспользовавшись полученным неравенством и неравенством Липшица,
получим
|
|
|
~ |
* |
|
|
x |
|
|
|
|
|
|
|
|
(x x0 )2 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
yi yi |
2MAn (t x0 )dt 2MAn |
|
|
|
|
|
. |
|
(14) |
|||||||||||
|
|
|
|
|
|
2! |
|
||||||||||||||||
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Снова повторяя ту же процедуру, из (13) и (14) имеем |
|
|
|
|
|||||||||||||||||||
|
|
|
|
~ |
* |
|
|
|
2 (x x0 )3 |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
yi yi |
2M ( An) |
|
|
|
|
|
|
и т.д.; |
|
|
|
|
|||||||
|
|
|
|
|
|
3! |
|
|
|
|
|
|
|||||||||||
на m-м шаге получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
(x x0 ) |
m 1 |
|
|
|
|
|
|
|
m 1 |
|
|||||
|
|
~ |
* |
|
|
|
m |
|
|
|
2M ( An ) |
|
|
|
|
||||||||
|
|
yi yi |
2M ( An) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
(m 1)! |
|
An |
|
|
|
(m 1)! |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
т.к. величина справа |
стремится |
к нулю, а |
слева |
неотрицательна, |
то отсюда |
||||||||||||||||||
* |
|
|
~ |
(x) . Ч. т. д. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
следует, что yi |
(x) yi |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Часто вместо условия Липшица
F(x1) F(x2 ) A x1 x2
вводят более грубое, но обычно легко проверяемое условие существования ограниченной по модулю производной F (x) на отрезке [a,b].
33
Действительно, если при x [a, b] выполняется неравенство F (x) N, то,
пользуясь теоремой о конечных приращениях (формулой Лагранжа), получим
F(x1) F(x2 ) F ( )
x1 x2 ,
где x1 x2 , т.е. [a, b], поэтому |
|
|
|
N , а следовательно, |
||||||||
|
|
|||||||||||
|
F ( ) |
|
||||||||||
|
F(x1) F(x2 ) |
|
N |
|
x1 x2 |
|
. |
|||||
|
|
|
|
|||||||||
Приведенный выше пример |
|
функции |
F(x) | x |, для которой условие |
|||||||||
|
|
|
|
|
|
|
|
при x=0 не существует, говорит о |
||||
Липшица выполняется, но производная F (x) |
||||||||||||
том, что условие F (x) N, является более грубым, чем условие Липшица.
Пример. Пусть имеется д.у. 2-го порядка: y y 0 .
Ранее мы убедились, что решениями этого уравнения являются функции: y1 sin x ,
y2 cos x . |
|
|
|
|
|
|
Положим y y1; y |
|
y2 |
, тогда получим нормальную систему д.у.: |
|||
|
||||||
|
|
|
y |
y |
|
, |
|
|
|
1 |
|
2 |
|
|
|
|
|
y1. |
||
|
|
|
y2 |
|||
Пусть y1(0) 0 y10 ; y2 (0) y20 1;
Эквивалентная система интегральных уравнений имеет вид:
|
|
x |
|
|
|
y1 y10 |
y2 (t)dt, |
|
|
|
0 |
|
|
x |
|
|
|
0 |
y1 (t)dt. |
|
y2 |
y2 |
|
|
|
0 |
Решим эту систему методом последовательных приближений:
x
y1(m 1) y10 y2(m) 0
x y2(m 1) y20 y1(m)
0
x
(t)dt y2(m) (t)dt;
0 |
|
x |
|
(t)dt 1 y1(m) (t)dt; |
m 0,1,2,... |
0 |
|
34
Итак, |
y(0) |
0; y(0) |
1; |
y(1) |
x; |
y(1) |
1; |
y(2) x; |
y(2) |
1 |
x2 |
; |
||||||
|
||||||||||||||||||
|
1 |
2 |
|
|
1 |
|
|
2 |
|
1 |
|
|
|
2 |
2 |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
y(3) |
x |
x3 |
; |
y(3) |
1 |
x2 |
; |
y(4) |
1 |
x2 |
|
x4 |
; |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
1 |
6 |
|
2 |
|
2 |
|
2 |
|
2 |
24 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
При x= /2 получим:
y1(3) ( 2 ) 0.9248; y2(3) ( 2 ) 0.2337.
При x= /4 – соответственно y1(3) ( 4 ) 0.7047; y2(3) ( 4 ) 0.6916; y2(4) 0.7074.
Сравним результаты с точным решением:
y sin x x |
x3 |
|
x5 |
..., |
y |
|
cos x 1 |
x2 |
|
x4 |
... |
|
|
2 |
|
|
|||||||
1 |
3! |
5! |
|
|
2! |
4! |
|
||||
|
|
|
|
|
|||||||
4. Теорема о дифференцируемости решений дифференциальных уравнений, разрешенных относительно производной.
Теорема. Если в окрестности точки (x0 , y0 ) функция f (x, y) имеет
непрерывные производные до m-ого порядка включительно, |
то решение y(x) |
||
уравнения |
|
|
|
|
dy |
f (x, y), |
(1) |
|
dx |
||
|
|
|
|
удовлетворяющее начальному условию y(x0 ) y0 , в некоторой окрестности точки
(x0 , y0 ) имеет непрерывные производные до (m+1) порядка включительно.
Доказательство. |
Подставляя решение y(x) в уравнение |
(1), получим |
||
тождество |
|
|
|
|
|
|
dy |
f (x, y(x)), |
(1a) |
|
|
|
||
|
|
dx |
|
|
и, следовательно, y(x) |
имеет в некоторой окрестности рассматриваемой точки |
|||
непрерывную производную, равную f (x, y(x)). Тогда в силу |
существования |
|||
непрерывных производных функции f, будет существовать непрерывная вторая производная решения
35
d 2 y |
|
f |
|
f dy |
|
f |
|
f |
f (x, y(x)). |
||
|
|
|
|
|
|
||||||
dx2 |
x |
y dx |
x |
y |
|||||||
|
|
|
|
|
|||||||
Если m>1, то, в силу существования непрерывных производных второго порядка,
можно, еще раз дифференцируя тождество (1а), обнаружить существование и непрерывность третьей производной
d |
3 |
y |
|
|
2 |
f2 |
|
|
2 |
f |
|
|
2 |
f2 |
|
f |
|
f |
|
f |
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
2 |
|
f |
|
f 2 |
|
|
f |
|||||||||||
|
|
3 |
x |
x y |
y |
|
x |
y |
|||||||||||||
dx |
|
|
|
|
|
|
y |
|
|
||||||||||||
Повторяя это рассуждение m раз, докажем утверждение теоремы.
4.1 Особые точки дифференциальных уравнений.
Рассмотрим теперь точки (x0 , y0 ) , в окрестности которых решение уравнения (1), удовлетворяющего условию y(x0 ) y0 , не существует, или решение существует, но оно не единственно. Такие точки называются особыми точками уравнения.
Кривая, состоящая из особых точек, называется особой кривой. Очевидно,
что особые точки будут среди тех, в которых нарушаются условия теоремы
существования и единственности решения уравнения (1) (см. разд. 3.2 при n=1).
Первое условие теоремы нарушается, если функция |
f (x, y) имеет разрывы. |
||||||||||||||||||||||||
Однако, если |
переменные x и y равноправны, |
то |
в случае, когда при |
||||||||||||||||||||||
(x, y) (x0 , y0 ) |
функция |
|
f (x, y) |
|
, |
|
уравнение |
(1) |
может |
быть |
заменено |
||||||||||||||
|
|
|
|||||||||||||||||||||||
уравнением |
|
dx |
|
1 |
|
, |
правая |
|
часть |
которого |
|
будет |
непрерывна в |
||||||||||||
|
|
|
|
|
|
||||||||||||||||||||
|
|
dy |
f (x, y) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
точке (x0 , y0 ) , если считать, что |
1 |
|
|
lim |
|
|
1 |
|
0 . |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
f (x0 |
, y0 ) |
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
( x, y) ( x0 |
, y0 ) |
f (x, y) |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Поэтому в тех задачах, в которых переменные x и y равноправны, особые |
|||||||||||||||||||||||||
точки могут быть только среди тех, в которых обе функции f (x, y) |
и |
1 |
|
||||||||||||||||||||||
|
|
||||||||||||||||||||||||
f (x, y) |
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
разрывные.
Пусть уравнение (1) имеет вид:
36
|
|
|
|
|
|
dy |
|
M (x, y) |
, |
|
|
|
|
(2) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
dx N (x, y) |
|
|
|
|
|
|
|
|
|||
где функции M (x, y), N(x, y) – непрерывны. В этом случае функции |
|
M (x, y) |
|
и |
|||||||||||||
|
N (x, y) |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
N (x, y) |
будут одновременно разрывными, лишь в тех точках (x , y ) , |
в которых |
||||||||||||||
|
|
|
|
||||||||||||||
|
M (x, y) |
|
|
|
|
|
|
0 |
0 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
M (x , y ) N(x , y ) 0 и не существует пределов: |
lim |
M (x, y) |
и |
lim |
|
N (x, y) |
. |
|
||||||||
|
|
|
|
|
|||||||||||||
0 |
0 |
0 |
0 |
|
|
|
|
y y0 N (x, y) |
y y0 |
|
M (x, y) |
|
|||||
|
|
|
|
|
|
|
|
|
|
x x |
x x |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
0 |
|
|
0 |
|
|
|
|
Особые точки уравнения (2) можно классифицировать по характеру поведения интегральных кривых в окрестности этих точек. Рассмотрим поведение интегральных кривых в окрестности особой точки указанного выше типа на примере следующего простого уравнения:
dy |
|
ax by |
. |
(3) |
|
|
|||
dx |
|
cx dy |
|
|
По характеру такого поведения можно классифицировать особые точки уравнения
(2).
Пусть |
ad bc 0 |
(т.к. в противном случае |
dy |
const ) (т.е. при этом |
|||||
dx |
|||||||||
|
|
|
|
|
|
|
|
||
|
a |
|
b |
A |
ax by |
A). |
|
|
|
|
c |
d |
cx dy |
|
|
||||
|
|
|
|
|
|
||||
Уравнение (3) эквивалентно следующей системе дифференциальных уравнений:
dx cx dy,dt
dy ax by
dt
с матрицей |
c |
d |
A |
. |
|
|
a |
b |
Найдем особые точки этой системы, для чего рассмотрим различные возможные случаи.
37
1. Корни характеристического уравнения
det( A I ) 2 (b c) bc ad 0 , |
(4) |
||
где I – единичная матрица, различны: 1 |
2 . |
|
|
В этом случае существует линейное преобразование |
|
||
|
|
x |
|
|
B |
|
|
|
y |
|
|
с невырожденной матрицей B , такое, что:
|
|
|
|
|
d |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
B 1 |
|
AB 1 |
|
|
, |
|
|
|||||
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
||||
|
d |
|
|
BAB 1 |
|
|
|
0 |
|
|
|
||||||
Или |
|
|
|
|
1 |
|
|
|
|
, |
|||||||
|
|
|
|
||||||||||||||
|
dt |
|
|
|
|
|
|
0 2 |
|
|
|||||||
|
|
|
|
|
|
|
d |
1 , |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
т.е. |
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
(5) |
||
|
|
|
|
|
|
dt |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
d |
|
2 . |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|||
Таким образом, система (5), а следовательно, и уравнение (3), сводится к уравнению
d |
|
2 |
|
|
|
. |
(6) |
||
d |
||||
|
|
1 |
|
|
Очевидно, первое условие теоремы существования и единственности для этого
уравнения нарушается в точке |
0 0 , 0 0, т.е. |
эта точка является особой. |
||||||||
Уравнение (6) решается разделением переменных: |
d |
|
2 |
d |
, |
|||||
|
|
|
||||||||
|
|
|
|
|
|
|
|
1 |
|
|
откуда |
С |
|
|
|
2 1 . |
|
|
|
|
(7) |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
Из (7) видно, что существование и единственность этого решения зависят от характера корней 1, 2 . Рассмотрим возможные варианты.
1.1. Корни 1, 2 действительные и одного знака.
38
Без ограничения общности можно считать, что 2 1 0 . В этом случае решение существует , , но может быть не единственным. Действительно, все интегральные кривые (7) проходят через особую точку 0 , 0 .
Все они в |
начале |
координат касаются оси , т.к. |
||||||||
|
d |
|
|
|
|
|
|
2 1 |
|
0 . |
|
|
|
|
|
||||||
|
|
|
|
|
||||||
|
|
|
|
|
||||||
|
|
|
2 |
C |
|
1 |
|
|||
|
d |
|
||||||||
|
|
|
|
|
|
|
|
|
||
|
|
|
0 |
1 |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
Кроме семейства (7), куда входит и решение =0, существует также, как видно из
(6), интегральная кривая =0.
Через особую точку 0 , 0 (или в первоначальных координатах x0 0 , y0 0 ) C
проходит бесконечно много интегральных кривых уравнения (3), такая особая точка называется узлом.
Пример 1. Рассмотрим дифференциальное уравнение:
d 2 . d
Его решение C 2 представляет из себя семейство парабол в прямоугольной системе координат и с вершиной в особой точке – узле.
1.2. Корни 1, 2 - действительные разных знаков.
Тогда 2 p 0 , |
и решение (7) имеет вид: C |
|
|
|
p . |
В этом случае через |
||
|
|
|||||||
1 |
|
|
|
|
|
|
|
|
проходят только два решения =0 и =0. |
|
|||||||
особую точку 0 , 0 |
|
|||||||
Остальные интегральные кривые не проходят через начало |
|
|||||||
координат. Такая особая точка называется седлом. |
|
|
|
|||||
Пример 2. Рассмотрим дифференциальное уравнение |
d |
|
. |
|||||
|
||||||||
|
|
|
|
|
|
d |
|
|
Общий интеграл этого уравнения C . Для этого уравнения особая точка 0 , 0
является седлом.
1.3. Корни 1, 2 – комплексно сопряженные:
1 p qi, |
2 p qi. |
39
В этом случае, как видно из (5), мы можем считать, что . Таким образом,
система уравнений (5) сводится к одному уравнению: ddt ( p iq) . Его решение имеет вид:
Cept ei(qt ) , Cept e i(qt ) ,
где C 0 - вещественная постоянная.
Получившиеся новые переменные , принимают комплексные значения при вещественных x и y, поэтому перейдем к вещественным переменным U, V при помощи невырожденного преобразования
UCe pt cos(qt ), 2
VCe pt sin(qt ). 2i
Врезультате будем иметь:
|
U 2 V 2 C2e2 pt , |
|
V |
tg(qt ) , |
(8) |
|||||||
|
|
|
|
|||||||||
|
|
|
|
U |
|
|
|
|
|
|||
|
|
|
|
|
|
|
V |
|
|
|||
|
|
|
|
|
|
|||||||
Или |
U 2 V 2 C1 exp |
p |
arctg |
. |
|
|||||||
|
|
|
||||||||||
|
|
|
|
|
q |
U |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
p |
|
В полярных координатах (U r cos , |
V r sin ) имеем r Ce q . |
|||||||||||
Это семейство логарифмических спиралей в плоскости U, V с асимптотической точкой в начале; все кривые примыкают к началу, но без определенной предельной касательной. Они делают у особой точки (0,0) бесконечное количество оборотов. Такая особая точка называется фокусом.
1.4. Корни 1, 2 чисто мнимые: 1 qi; 2 |
qi . |
|
Из (8) имеем |
U 2 V 2 C 2 . Семейство |
интегральных |
кривых есть семейство замкнутых кривых, окружающих особую точку; через саму особую точку не проходит ни одной интегральной кривой. Такая особая точка называется
40
