Дифференциальные уравнения. Учебное пособие-3
.pdf
Теоремаосуществовании иединственностирешениядифференциального уравнения, не разрешенного относительно производной.
Пусть в некотором параллелепипеде с центром в точке (x0 , y0 , p0 )
функция |
F(x, y, p) непрерывна по совокупности переменных вместе |
||||||||||||||
с частными производными |
F , |
|
F , причем выполнены условия: |
||||||||||||
|
|
|
|
y |
|
p |
|
|
|
|
|||||
|
|
F (x , y , p ) 0, |
|
F |
(x , |
y , |
p ) 0. |
(6) |
|||||||
|
|
|
|
|
|
|
|||||||||
|
|
0 |
0 |
0 |
|
|
|
p |
0 |
0 |
0 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
||||
Тогда на некотором сегменте |
|
x x0 |
|
|
Н |
существует |
единственное |
||||||||
|
|
||||||||||||||
решение |
|
уравнения |
(5), |
|
|
удовлетворяющее |
начальному условию |
||||||||
y(x ) y |
0 |
, производная которого |
dy |
(x ) p . |
|
|
|||||||||
|
|
|
|||||||||||||
0 |
|
|
|
|
|
dx |
0 |
0 |
|
|
|||||
Определение. Множествоточекповерхности впространстве(x, y, p), в которых нарушаются условия (6) теоремы, определяют кривые, называющиеся дискриминантными. Ветви этой кривой могут быть решениями уравнения (5). Ветви дискриминантной кривой называются осо-
быми решениями уравнения (5).
ПРИМЕР РЕШЕНИЯ ДУ ПЕРВОГО ПОРЯДКА
1. y (4x y 3)2 . 2. 2xy y ln( y ) . 3. (2x2 y ln y x) y y.
4. y(xy y ) x( y )2 (1 x) .
|
|
|
2 |
|
1 |
|
|
5. y |
y |
4x2 . |
|||||
|
|
||||||
6. Определить область задания уравнения, область существования решения задачи Коши, область существования и единственности реше-
ния задачи Коши, особые решения: y 1 . x2
11
Решение
1. Уравнение y (4x y 3)2 – этоДУ, приводящеесякуравнению
с разделяющимися переменными заменой: |
z 4x y 3, |
dz |
4 |
dy |
. |
||||||
dx |
|
||||||||||
|
y z 4 z2 . Получаем новое уравнение |
|
|
dx |
|||||||
Тогда |
z z2 4 – уравне- |
||||||||||
ние с разделяющимися переменными |
dz |
|
dx. |
Общее решение ДУ |
|||||||
z2 4 |
|||||||||||
имеет вид |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|||
|
arctg(z 2) x C, |
|
|
|
|
|
|
||||
или |
arctg(4x y 1) x C, |
y(x) tg(x C) 4x 1. |
|
|
|
||||||
2. Уравнение 2xy y ln ( y ) – уравнение Клеро. Перепишем ДУ ввиде y ln ( y ) 2xy , продифференцируемпо x исделаем замену пе-
ременных: y (x) p(x) . Тогда получим ДУ: 3 p |
dp |
|
1 |
|
|
|
|
|||||||||||
|
|
|
|
|
|
2x |
. Перей- |
|||||||||||
dx |
|
p |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
||||
дем от функции p(x) |
к обратной функции x( p) [2]: |
|
|
dx |
3p |
2x – |
||||||||||||
|
dp |
p |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
это линейное |
ДУ |
первого порядка. Перепишем |
его |
в |
виде |
|||||||||||||
3p |
dx |
2x |
1 |
|
и решим методом вариации постоянной. Общее реше- |
|||||||||||||
|
p |
|||||||||||||||||
|
dp |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
ние однородного ДУ |
3p |
dx |
2x 0 имеет вид |
x( p) Cp 3 . Ищем ре- |
||||||||||||||
|
||||||||||||||||||
|
|
|
|
|
|
|
dp |
|
|
|
|
|
|
|
|
|
|
|
шение неоднородного ДУ в виде x( p) C( p) p |
2 |
|
|
|
|
|
|
|
|
|
||||||||
3 . Подставим решение |
||||||||||||||||||
в исходное уравнение и получим соотношение для определения неиз-
|
|
1 |
dC( p) |
p 1 . Тогда |
C( p) p |
1 |
|
|
вестной функции |
C( p) : |
3p3 |
3 |
C, |
||||
dp |
||||||||
|
|
|
|
|
|
|
12
|
|
|
|
|
1 |
|
|
2 |
|
а |
x( p) |
|
|
3 |
|
3 |
. Решение уравнения Клеро в параметрическом |
||
C p |
|
|
p |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
2 |
|
|
|
|
|
1 |
|
1 |
|
виде: |
x( p) |
|
|
3 |
|
3 |
, |
y( p) ln p 2 |
|
|
3 |
|
3 |
. |
|||
C p |
|
|
p |
|
|
C p |
|
|
p |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3. Уравнение (2x2 y ln y x) y y – уравнение, нелинейноеотноси-
тельно y(x). Перейдем к функции x( y) : |
2x2 y ln y x y |
dx |
. Перепи- |
||||||||||||
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy |
|
шем ДУ в виде |
yx 2 |
dx |
x 1 |
2y ln y |
– это уравнение Бернулли. |
||||||||||
|
|||||||||||||||
|
|
|
|
dy |
|
|
|
|
|
|
|
|
|||
Замена переменных: z(y) |
1 |
, |
|
dz |
x 2 |
dx |
, тогда ДУ примет вид |
||||||||
|
|
|
|
||||||||||||
|
|
|
|
|
|
x |
|
dy |
|
dy |
|||||
y |
dz |
z 2y ln y |
– неоднородное линейное дифференциальное урав- |
||||||||||||
|
|||||||||||||||
|
dy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
нение первого порядка. Решаем методом вариации постоянной. Реше-
ниеоднородногоДУ y dydz z имеетвид z( y) Cy. Тогдарешениенеод-
нородногоДУ z( y) C( y) y. ПодставляярешениевнеоднородноеДУ, по-
лучаем |
dC( y) |
2 |
ln y |
, или C( y) ln2 |
y C. Общее решение неодно- |
dy |
y |
||||
родного уравнения z( y) y(ln2 y C), |
или x( y) y 1(ln2 y C) 1 . |
||||
4. Уравнение |
y(xy y ) x( y )2 (1 x) – ДУ второго порядка, до- |
||||
пускающее понижение порядка. Проверим выполнение условий однородности [3]: при замене y ky, y ky , y ky уравнение примет
вид ky(kxy ky ) k2 x( y )2 (1 x) , где k – постоянная величина. Уравнение однородно, значит, замена y(x) exp ∫ z(x)dx понижает по-
рядок уравнения на единицу. Вычислим y (x) z (x)exp ∫ z(x)dx ,
13
y (x) z (x) z2 (x) exp ∫ z(x)dx иподставимэтивыражениявисход-
ное уравнение. После приведения подобных получаем ДУ xz z 0 .
Общее решение этого уравнения |
z(x) |
C1 |
, тогда |
|||||||||
x |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
C |
|
|
|
|
|
|
y(x) exp |
∫ |
1 |
dx |
C1x C2 . |
||||||
|
|
x |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
||
5. Уравнение y |
|
y |
2 |
|
1 |
|
|
– это уравнение Риккати. Определим |
||||
4x2 |
||||||||||||
|
|
|||||||||||
частное решение в виде правой части: y1(x) ax . Подстановка в ДУ позволяет определить параметр а: a2 a 14 0 , т. е. a 12 . Тогда частное решение: y1(x) 21x . Общее решение ищем в виде y(x) z(x) y1(x) z 21x . Подставляя решение в уравнение Риккати,
получаем z z2 xz 0 – это уравнение Бернулли. Используя прием из пункта 3, решаем уравнение Бернулли. Общее решение имеет вид
y(x) |
|
1 |
|
1 |
. |
|
|
|
|
|
|
|
|
|
(ln x C)x |
|
|
|
|
|
|
|
|
|
|||||
|
|
2x |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
6. Праваячастьуравнения y |
x2 |
– функция, непрерывнаядлявсех |
||||||||||||
|
||||||||||||||
x ( , |
0) (0, |
) . Общее решение ДУ: |
y(x) |
1 |
C. |
Интегральная |
||||||||
x |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
кривая имеет горизонтальную асимптоту y 0 . Условия теоремы суще-
ствования и единственности выполнены на всей действительной прямой,
кроме |
точки x 0 . |
Обратное уравнение |
dx |
x2 , общее решение |
||
dy |
||||||
|
|
|
|
|
||
x( y) |
1 |
. Прямая |
x 0 – частноерешениеобратногоуравнения. |
|||
y C |
||||||
|
|
|
|
|
||
14
Определить тип дифференциальных уравнений и решить:
1) y2 y x2 sin3 x y3 ctg x ; 2) (1 x2 y)dx x2 ( y x)dy 0 ; 3) ( y )4 4y(xy 2 y)2 ;
4) y (3 ( y )2 y) ( y )4 ;
5) xy 2 |
y cos x 2 y ; |
|
||||||||
6) x y cos |
y |
dx xcos |
y |
dy 0 |
; |
|||||
x |
|
|||||||||
|
|
|
|
|
|
|
x |
|
||
7) yy xy x3 ; |
|
|
|
|
|
|||||
|
2 |
|
|
|
|
y 2 |
|
|||
8) ( y ) |
|
y y |
|
; |
|
|||||
|
|
|
||||||||
|
|
|
|
|
|
x |
|
|||
9) y x3 y 3y ;
10) (2x 3y 1)dx (4x 6y 5)dy 0 ;
11) xy |
|
2y 1 |
|
2 |
; |
|
|
( y ) |
|
||||
12) y ( y )2 |
1; |
|
|
|
||
13) yy y2 ctg x cos x ;
14) (cos x xsin x)ydx (xcos x 2y)dy 0 ; 15) x2 ( y )2 2(xy 2) y y2 0 ;
16) x2 ( y2 y ( y )3 ) 2 y2 y 3xy( y )2 ; 17) y y xy3 ;
18) (sin x y)dy ( y cos x x2 )dx 0 ; 19) y2 x2 ( y )5 xy(( y )2 ( y )3) ; 20) y2 y ( y )3 ;
15
21) y x y (x2 1) xy ;
22) (x3 2xy2 )dx 3x2 ydy xdy ydx ;
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
23) y |
y exp(2x) |
y ; |
||||||||||||||
|
||||||||||||||||
24) ( y )2 2xyy 0 ; |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
4xy |
|
|
|
||
25) y |
8x |
|
y x2 |
1 |
; |
|
||||||||||
|
|
|
||||||||||||||
26) 6x5 ydx ( y4 ln y 3x6 )dy 0 ; |
||||||||||||||||
27) ( yy )3 27x( y2 2x2 ) ; |
||||||||||||||||
28) ( y y 3( y )2 ) y ( y )5 ; |
||||||||||||||||
|
|
|
1 |
|
|
|
|
|
3 |
|
|
|
|
|
||
29) y |
2 |
x |
y ; |
|
|
|||||||||||
|
|
|
|
|
||||||||||||
30) yy x |
1 |
x2 |
y2 2 |
|||||||||||||
|
|
|
|
|
|
|
|
; |
||||||||
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
2 |
|
|
x |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||||
31) ( y )4 |
( y )5 y( y )3 y ; |
|
|
|
|||||
32) (x2 1) y y2 |
2xy 1 0 ; |
|
|||||||
33) 2 |
|
2 |
y 1 x |
4 |
y |
2 |
3 |
dy 0 |
; |
x |
|
|
dx x |
|
|||||
|
|
|
|
|
|
|
|
|
|
34) x2 yy 1 (1 y)xy ; 35) (1 x2 ) y 2xy2 xy ; 36) 2 y x ln y ;
37) yy ( y )2 2xy2 ; 38) 2xy y sin y ;
16
39) y (3 y( y )2 ) ( y )4 ; 40) x yy 2xy ln x 1 0 ;
41) (exp y 2xy)dx (exp y x)xdy 0 ;
42) x(x 1) y y3 |
xy ; |
||||||
43) (2x y 5) y 3x 6 ; |
|||||||
44) y 3 2x y 2 ; |
|
||||||
45) y xy 4 |
y ; |
|
|
||||
46) 5y ( y )2 |
x(x y ) ; |
||||||
47) y sin x y cos x 1; |
|||||||
48) xy 1 exp(x y) ; |
|||||||
|
y |
2 |
1 |
|
( y |
2 |
1)dx xydy ; |
49) x |
|
1 |
|
||||
|
|
|
|
|
|
|
|
50) (2x exp y y4 ) y y exp y .
17
2. СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
dyi |
f |
|
(x , y , y , |
, y ) |
i |
|
. |
(1) |
|
i |
1,n |
||||||||
|
|||||||||
dx |
1 1 2 |
n |
|
|
|
|
|||
|
|
|
|
|
|
|
|||
Определение. Система вида (1) называется нормальной системой дифференциальных уравнений (СНДУ).
Теорема о существовании и единственности решения СНДУ.
Если функции fi определены и непрерывны на открытом множестве Г,
а соответствующие частные производные |
dfi |
, i, j |
|
, тоже непре- |
|||||||||||||||
1, n |
|||||||||||||||||||
dy j |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
рывны на Г, то у системы (1) будет существовать решение |
|
|
|||||||||||||||||
|
|
|
yi |
i (x), |
i |
|
|
, |
|
|
|
|
|
|
(2) |
||||
|
|
|
1,n |
|
|
|
|
|
|
||||||||||
а при наличии начальных условий |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
yi (x0 ) y0 , |
|
i |
|
, |
|
|
|
|
|
|
(3) |
|||||
|
|
|
|
1,n |
|
|
|
|
|
|
|||||||||
это решение будет единственным. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Систему (1) можно представить в векторном виде: |
|
|
|
|
|
||||||||||||||
|
dY |
F(x, y) |
Y ( y , |
y , , y |
n |
)T , |
|
F |
( f , f |
2 |
, , f |
n |
)T . (4) |
||||||
|
|
|
|||||||||||||||||
|
dX |
1 |
2 |
|
|
|
|
|
|
1 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
18
СИСТЕМЫ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
Определение. Система дифференциальных уравнений называется линейной (СЛДУ), если она линейна относительно всех неизвестных функций и их производных:
|
dyi |
n |
|
|
|
|
|
||
|
∑ aij (x) y j |
fi (x), |
i |
1,n |
, |
(5) |
|||
|
dx |
||||||||
|
j 1 |
|
|
|
|
|
|||
Векторная форма записи СЛДУ: |
|
|
|
|
|
||||
|
|
|
dY |
AY F. |
|
|
|
(6) |
|
|
|
|
|
|
|
|
|||
|
|
|
dX |
|
|
|
|
|
|
Если задано начальное условие |
|
|
|
|
|
||||
|
|
|
Y (x0 ) Y0 , |
|
|
|
(7) |
||
торешениебудетединственнымприусловии, чтовектор-функцияF непрерывна на Г и коэффициенты матрицы А: aij (x) – тоже непрерывные
функции.
Введем линейный оператор |
L[y] |
dY |
AY, тогда (6) можно пере- |
||
dx |
|||||
писать в виде |
|
|
|
||
|
|
|
|
||
L y F, |
|
(8) |
|||
если F 0, то операторное уравнение (8) |
называется однородным |
||||
и имеет вид |
|
|
|
|
|
|
L y 0 . |
|
(9) |
||
Так как оператор линейный, то для него выполняются следующие свойства:
1) если Y – решение однородной системы (9), то CY (C const) бу-
дет тоже решением уравнения (9).
2) если Y1, Y2 являютсярешением(9), то Y1 Y2 – тожерешение(9).
19
|
|
|
n |
|
|
|
|
|
Следствие. Линейнаякомбинация ∑ CiYi , |
Ci – const, – решение(9), |
|||||||
если Yi – решения (9). |
|
|
i 1 |
|
|
|
|
|
Yn решений однородной СЛДУ (9) и они ли- |
||||||||
Если даны Y1, Y2 , , |
||||||||
нейно независимы, то линейные комбинации вида |
||||||||
|
n |
|
|
|
|
|
|
|
|
∑ iYi (x) 0 |
|
|
|
(10) |
|||
|
i 1 |
|
|
|
|
|
|
|
только при условии, что все i |
0 . |
|
|
|
|
|
||
Это означает что определитель системы (10): |
||||||||
|
y11 |
y12 |
y1n |
|
|
|
||
|
|
|
||||||
W |
y12 |
y22 |
y2n |
|
|
0 . |
||
|
|
|
|
|
|
|
|
|
|
yn1 |
yn2 |
ynn |
|
|
|
||
Определитель системы (10) называется определителем Вронского
для решений Y1, Y2 , , Yn .
Теорема1. ЕслиопределительВронскогодлялинейнойоднородной системы ДУ (9) с непрерывными на отрезке [a, b] коэффициентами
aij (x) |
равен нулю хотя бы в одной точке x0 [a, |
b] , то решения |
Y1, Y2 , |
, Yn линейно зависимы на этом отрезке и, |
следовательно, |
определитель Вронского равен нулю на всем отрезке. |
|
|
Доказательство. Так как aij (x) непрерывны, то система (9) удовле-
творяет условию теоремы о существовании и единственности решения СЛДУ, следовательно, начальное условие Y (x0 ) 0 определяет един-
ственное решение системы (9). Определитель Вронского в точке x0 равен нулю, тогда существует такая нетривиальная система Сi , для кото-
ройвыполняется C1Y1(x0 ) C2Y2 (x0 ) CnYn (x0 ) 0. Соответствующая линейная комбинация для другой точки x будет иметь вид
n
∑CiYi (x) Y (x) , причем Y (x) – непрерывная вектор-функция, удовле-
i 1
20
