Конспект МАТАН2
.pdf
|
|
Лекція 25. Лінійні неоднорідні диференціальні рівняння |
|
131 |
|
|
|
|
|
|
|
|
4. Інтегруючи C (x) та C2(x) знаходимо шукані функції C1(x) та C2(x). |
||||
|
|
1 |
|
||
|
Приклад 25.1. |
Знайти загальний розв’язок рівняння y y |
. |
||
|
sin x |
||||
|
|
|
|
|
25.3. Лінійні диференціальні рівняння зі сталими коефіцієнтами і правою частиною спеціального вигляду
Функція вигляду
f (x) e x (Pn(x)cos x Qm(x)sin x)
називають функцією спеціального вигляду (квазімногочленом).
До розв’язання таких рівнянь застосовують метод підбирання частинного розв’язку ЛНДР. Розгляньмо типи рівнянь до яких він застосовний.
1. Рівняння вигляду
y(n) p1y(n 1) ... pny Pm(x),
де
Pm(x) a0xm a1xm 1 ... am.
а) якщо 0 не є коренем характеристичного рівняння, то частинний розв’язок ЛНДР шукають у вигляді
yч.н. Qm(x) A0xm A1xm 1 ... Am,
де Ai,i 0,m, — невизначені коефіцієнти;
б) якщо 0 є коренем характеристичного рівняння кратності r (є «резонанс» порядку r), то частинний розв’язок ЛНДР шукають у вигляді
yч.н. xrQm(x) xr (A0xm A1xm 1 ... Am ),
де Ai,i 0,m, — невизначені коефіцієнти.
2. Рівняння вигляду
y(n) p1y(n 1) ... pny eaxPm(x),
де
Pm(x) a0xm a1xm 1 ... am.
а) якщо a не є коренем характеристичного рівняння, то частинний розв’язок ЛНДР шукають у вигляді
yч.н. eaxQm(x) eax (A0xm A1xm 1 ... Am ),
де Ai,i 0,m, — невизначені коефіцієнти;
Лекція 26. Системи лінійних диференціальних рівнянь |
135 |
2)якби 1-й вид жив окремо, то кількість його особин неперервно б збільшувалась зі швидкістю пропорційною наявній кількості особин: x k1x, k1 0.
3)особини 2-го виду харчуються лише особинами 1-го виду;
4)якби 2-й вид жив окремо, то він би поступово вимирав: y k2y, k2 0.
Розгляньмо тепер випадок, коли види співіснують. Припустімо також, що коефіцієнт k1 зменшиться на величину пропорційну y. Так само, що коефіцієнт
k2 завдяки наявності 1-го виду (харчів) зміниться на величину пропорційну x.
За цих припущень маємо таку систему диференціальних рівнянь:
|
|
k y)x, |
|
x (k |
|||
|
1 |
|
3 |
|
|
|
k4x)y, |
y (k2 |
|||
|
|
|
|
|
|
|
|
де k1,k2,k3,k4 — додатні числа.
26.2. Основні поняття
Для розв’язання багатьох задач математики, фізики і техніки часто потрібно розглядати декілька невідомих функцій. Знаходження цих функцій може привести до кількох ДР, які утворюють систему.
Систему n ДР 1-го порядку, розв’язаних щодо похідної, вигляду
|
|
|
|
|
|
|
|
|
|
|
f (t,x ,...,x |
|
|
), |
||||||
|
|
|
|
|
|
|
|
|
x |
n |
||||||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
1 |
1 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
f2(t,x1,...,xn ), |
|||||||||
|
|
|
|
|
|
|
|
|
x |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
............................. |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f |
(t,x ,...,x |
|
|
) |
|||||
|
|
|
|
|
|
|
|
|
x |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
1 |
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|||||
називають нормальною системою ДР. |
|
|
|
|
|
|
|
|||||||||||||
|
Обмежмося розглядом нормальних систем із двома невідомими функціями: |
|||||||||||||||||||
|
|
|
|
|
|
|
a (t)x |
|
a |
(t)x |
|
|
|
f (t), |
||||||
|
|
|
|
|
|
x |
1 |
2 |
|
|||||||||||
|
|
|
|
|
|
|
1 |
|
11 |
|
|
12 |
|
|
|
1 |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a22(t)x2 f2(t). |
||||||
|
|
|
|
|
|
x2 a21(t)x1 |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Цю систему можна записати у матричному вигляді |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
Ax f , |
|
|
|
|
||||
|
x (t) |
a (t) |
a (t) |
|
|
|
f (t) |
|
|
|
|
|
||||||||
|
1 |
|
|
11 |
|
|
12 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
де x |
|
,A |
|
(t) |
a |
|
|
, f |
|
, або |
||||||||||
|
x |
2 |
(t) |
a |
22 |
(t) |
|
|
|
f (t) |
|
|
|
|
|
|||||
|
|
|
|
21 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
L[x(t)] f (t),
де L dtd A.
136 |
Розділ 3. Диференціальні рівняння |
Якщо стовпець f 0, то систему називають однорідною, інакше — неод-
норідною.
Розв’язком системи ДР називають сукупність із двох функцій x1(t),x2(t)
які справджують кожному з рівнянь цієї системи.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Теорема 26.1. |
Якщо функції aij (t), fi(t),i, j 1,2, |
неперервні |
на |
відрізку |
|||||||||||
|
[a;b], то в досить малому околі кожної точки M0(t0;x10;x20),t (a;b), |
викона- |
||||||||||||||
|
но умови існування та єдиності розв’язку задачі Коші. |
|
|
|
|
|
|
|
||||||||
|
Якщо стовпці x1(t) та x2(t) |
є розв’язками лінійної однорідної системи |
||||||||||||||
|
L[x ] 0, то для довільних сталих |
C |
1 |
та C |
2 |
стовпець C x (t) C x |
2 |
(t) також є |
||||||||
|
розв’язком цієї системи. |
|
|
|
1 |
1 |
2 |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
(про структуру розв’язку однорідної системи). |
Загальним |
||||||||||||
|
Теорема 26.2 |
|||||||||||||||
|
розв’язком |
|
|
|
|
|
||||||||||
|
|
нормальної лінійної однорідної системи |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
dx |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
dt |
Ax |
|
|
|
|
|
|
|
|
||
|
з неперервними на відрізку [a;b] коефіцієнтами aij (t),i, j 1,2 |
є лінійна комбіна- |
||||||||||||||
|
ція двох лінійно незалежних в інтервалі (a;b) розв’язків x1(t),x2(t) системи |
xзаг.одн. C1x1(t) C1x2(t),
де C1,C2 const.
26.3. Метод Ейлера розв’язання однорідної системи диференціальних рівнянь із сталими коефіцієнтами
Шукатимемо розв’язок однорідної системи |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
dx |
|
|
|
|
|
|
|
|
|
у вигляді |
|
|
|
dt |
Ax |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
(t) |
|
|
t |
|
|||
|
|
|
t |
|
|
e |
|
|
||||||
|
|
|
|
1 |
|
|
|
|
1 |
|
|
|||
|
x |
e |
|
, |
|
|
|
|
|
. |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
x |
2 |
(t) |
e |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
Підставляючи ці розв’язок у систему і перетворюючи її, дістанемо |
||||||||||||||
|
|
|
|
|
) a |
|
|
|
|
0, |
|
|
|
|
|
|
(a |
11 |
12 |
2 |
|
|
|
||||||
|
|
|
|
|
1 |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a (a ) 0.
21 1 22 2
Щоб ця система мала ненульовий розв’язок необхідно й достатньо, щоб
Лекція 26. Системи лінійних диференціальних рівнянь |
137 |
|
a11 a12 |
|
0 2 (a |
11 |
a |
22 |
) (a a |
22 |
a a |
21 |
) 0. |
|||||||||||||
|
a21 |
a22 |
|
|
|
|
|
|
|
|
|
11 |
|
12 |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Це рівняння називають характеристичним рівнянням для системи. |
||||||||||||||||||||||||
1. Якщо корені характеристичного рівняння 1 |
та 2 |
дійсні різні, то зна- |
||||||||||||||||||||||
ходять відповідні їм нетривіальні розв’язки 1, 2 |
алгебричної системи і запи- |
|||||||||||||||||||||||
сують загальний розв’язок однорідної системи ДР у вигляді |
|
|
|
|
||||||||||||||||||||
|
|
|
|
C |
|
t |
|
|
|
t |
|
|
|
const. |
|
|
|
|||||||
|
|
|
x |
заг. одн. |
e |
1 C |
2 |
e |
2 ,C ,C |
2 |
|
|
|
|||||||||||
|
|
|
|
|
1 1 |
|
|
|
2 |
|
|
1 |
|
|
|
|
|
|
|
|
2. Якщо корені характеристичного рівняння 1 та 2 комплексно спряжені:
1 i , 2 i ,
то знаходять нетривіальний розв’язок алгебричної системи , який відповідає
1 i . Тоді
|
|
|
|
i |
|
|
|
|
i |
|
||
|
x1(t) |
Re( e |
),x2(t) |
Im( e |
|
|
) |
|||||
|
x |
заг. одн. |
(t) C x (t) |
C x |
(t),C ,C |
2 |
const. |
|||||
|
|
|
1 1 |
|
|
2 2 |
1 |
|
|
|||
3. Якщо корені 1 та 2 |
характеристичного рівняння дійсні і рівні, то зага- |
|||||||||||
льний розв’язок системи шукають у вигляді |
|
|
|
|
|
|||||||
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
A Bt |
|
,A,B,C,D const. |
||||||||
|
|
x(t) |
e |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C Dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x 8y, |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
x |
|
|
|
|
||||
Приклад 26.1. |
|
|
|
|
|
|
|
|
|
|
|
|
Розв’язати систему |
|
x |
y. |
|
|
|
|
|||||
|
|
|
|
y |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
26.4. Метод виключення змінних
Нормальну неоднорідну систему ДР |
|
|
|
|
||
|
|
a (t)x |
|
a (t)x |
|
f (t), |
x |
1 |
2 |
||||
|
1 |
11 |
12 |
1 |
||
|
|
a21(t)x1 |
a22(t)x2 |
f2(t). |
||
x2 |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
можна звести до одного диференціального рівняння 2-го порядку, виражаючи або змінну x2(t) із 1-го рівняння системи, або змінну x1(t) із 2-го рівняння системи.
Приклад 26.2. Розв’язати задачу Коші
|
|
dx |
|
|
x y cost, |
|
|
|
|
dt |
|
dy |
|
|
y 2x cost sint,x(0) 1,y(0) 2. |
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Відповіді |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
139 |
|||
9.2. Збігається. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
10.1. |
a3 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10.2. |
8 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10.4. |
27 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10.5. |
|
3 a2. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
11.1. |
268. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
11.2. S 16 15 . |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
11.3. I 2. |
|
|
|
|
|
|
|
|
|
|
|
12.1. |
|
J(u,v) |
|
|
|
|
1 |
|
. S |
b2 |
a2 |
ln |
|
. |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2v |
|
|
|
2 |
|
||||||||||||||||||||
12.2. S ab. |
|
|
|
|
|
|
|
|
|
|
|
12.3. m 5. |
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
13.1. |
V |
1 . |
|
|
|
|
|
|
|
|
|
|
|
13.2. |
V |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
13.3. |
V |
|
4 abc. |
13.4. m |
(2 |
2) R4. |
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
||||
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
14.1. |
a2 b2 (3a2 4 2b2). |
14.2. |
8a. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
15.1. 8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
15.2. |
0. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
15.3. |
I |
R4 |
. |
|
|
|
|
|
|
|
|
|
|
15.4. I 2 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
16.1. |
u(x,y) ln |
|
x |
|
2 ln |
|
y |
|
x |
C. 16.2. u(x,y,z) x2y 3xz2 |
C. |
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
16.3. A 5. |
|
|
|
|
|
|
|
|
|
|
y |
16.4. |
S |
3a2 . |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
17.1. S 4 a2. |
17.2. S |
2 a2. |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
17.3. |
2 ( |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
2 1). |
18.1. 4 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
15 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
19.1. |
7 . |
19.2. div |
|
3. |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
r |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
19.3. V 72 . |
|
|
|
|
|
|
|
|
|
|
|
19.4. 2 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
19.5. 3V. |
|
|
|
|
|
|
|
|
|
|
|
19.7. 0. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
19.6. 5. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
22.1. x Cearctg y,C const. |
22.2. |
x C ctg x y . |
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
22.3. |
y |
|
|
|
y2 |
1 Cx,y x. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
x |
|
|
|
x2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|