Дыферэнцыяльныя і інтэгральныя ўраўненні
Памер: 436с.
Гародня 2000
адвольны вектарслупок, t0 е I.
Падстаўляючы цяпер знойдзеную C(t) у выраз x = X(t)C(t), атрымаем агульнае рашэнне зыходнай сістэмы ў выглядзе
x(t) = X(t )\ Х~' (х )f(x)dx + X(t )С. (31)
'о
Пакажам, што формула (31) змяшчае ўсе рашэнні сістэмы (30) на інтэрвале I. Няхай x = x(t)~ любое рашэнне разглядаемай сістэмы. Выберам t0 е I , і няхай x(t0 ) = х°. Пакажам, што рашэнне, зададзенае формулай (31), пры адпаведным выбары С будзе задавальняць такім жа пачатковым умовам. Падставіўшы дадзеныя пачатковыя ўмовы ў (31), будзем мець х° = X(t0)C. У сілу
177
таго, што X(t) фундаментальная матрыца, адназначна знаходзім С = Х~' (t0)x°. Пры такім выбары С рашэнне (31) будзе супадаць з выбраным рашэннем x = x(t). Адсюль вынікае, што рашэнне сістэмы (30) з пачатковай умовай x(t0) = х° (задачы Кашы) можна запісаць у выглядзе
x(t)= X(t)X'(t0)x° + X(t)\ Х~'(t)f(t)dt ■ (ЗГ) ‘о
Заўвага 1. Аб азначэнні вытворнай і інтэграла ад матрычнай функцыі можна прачытаць у [25], глава V, § 6, с. 116, 118.
Прыклад 1. Рашыць сістэму ўраўненняў dx У , dy х — = x + ~ + t ch t, — = — + y + tsht. dt t dtt
Рашэнне. Складзём пачленна два ўраўненні, атрымаем лінейнае ўраўненне з невядомай функцыяй х + у:
d Г Л
^(x + y) = ^ + ~j(x + y) + te .
Рашэнне адпаведнага аднароднага ўраўнення ёсць x + y = Cjte'. Рашэнне неаднароднага ўраўнення адшукваем у выглядзе х + у = C](t)te'. Функцыю С/^^ызначаем з роўнасці t dCj , te —— = te , гэта значыць C,(t) = t + С,. Такім чынам, dt 7 7
х + у = С^е’+t2 е‘. Адымаючы пачленна з першага ўраўнення дадзенай сістэмы другое, маем
d ( 1}
~Т(ХУ) = \ 7“Т (x~y) + te . dt \ t J
tj
Адкуль знаходзім xy = C2te 1 +t2e '. Улічваючы атрыманыя вынікі, будзем мець
178
х = —(СІе‘ +С2е ' ) + t2ch t, y = —(CIe‘ C2e ‘ ) + t2sh t.
Прыклад 2. Рашыць сістэму ўраўненняў dx ■> dy
t—x + yt, r — = 2x + yt (t^O). at at
x dx
Рашэнне. 3 першага ўраўнення знаходзім у—~+~Адкуль t at
dy х 1 dx d2x . dy
— = — — + y ■ Падстаўляючы гэтыя выразы для у і — у
2 d2x dx dx
другое ўраўненне, атрымаем t‘ —~^ + t — — x = —2х + x + t — або dr dt dt
2 d2x
t —y = 0. У сілу таго, што t ^ 0, з апошняга ўраўнення маем dr
d2x „
—— = 0 і пасля інтэгравання атрымаем х = C] + C2t. Значыць, dt
х dx Cj + C2t C]
y = ~t+~dt=1~ ^2 ~ ^2^^ ■ Агульнае рашэнне
С1
дадзенай сістэмы х = С/ + C2t, у = — + 2С2.
§ 5. Метады інтэгравання лінейных неаднародных сістэм з пастаяннымі каэфіцыентамі
1°. Метад выключэння. Прыменім метад, апісаны ў § 2, да разглядаемай сістэмы (30) з пастаяннай матрыцай A.
Прыклад 1. Рашыць сістэму
1 7 5 5. 1 1 1
Х' =~6Х,+'бХ2~6Хз+1е х^~2Х/~2Х2 + 2Хз''
179
12 1 1 .
Хз ~ зх'~ зХг + зХз + зе '
1 .
Рашэнне. Дыферэнцуем па t першае ўраўненне: xt = ~Х/ + 6
7 5 5 ,
+ ~х2 ~Xj + —е • Пасля замены ў гэтым ураўненні х, і х3 іх 6 6 3
выразамі з другога і трэцяга ўраўненняў сістэмы будзем мець
ЗбХ] = 6х, ЗІХ] х2 +11х3 + 50е'. (32)
Атрыманае ўраўненне прадыферэнцуем па t : 36х\ = бх/ З1х3 +
х2 + 11х3 +50е/,адкуль
216'х/ = Збх, 186х, + 25х/ 41х2 + 19х3 + 322е‘. (33) 3 першага ўраўнення сістэмы і ўраўнення (32) знаходзім х2 і х3 праз t,xI,xl,xl'.
5 1
х2 = —X] + ^Х/ + 25хІ 5е ,
7 1 ,
х, = —х,—хІ + Зх, 5е . *22
(34)
Падставім х2 і х5 ва ўраўненне (33) і атрымаем адносна Xj ураў
ненне трэцяга парадку з пастаяннымі каэфіцыентамі і спецыяльнай правай часткай х; + Х] = 2е'.
Агульнае рашэнне гэтага ўраўнення мае від (гл. § 2.6) х, = = Cj + С2 cos t + С3 sin t + е'. Падставіўшы хІ,хІ,Хі у роўнасці (34), знаходзім
1 1
х2 = 2C] +~(С3 С2)cos t—(C3+ C2)sin t,
1 1
x3 = 3CI ~(C2 + C3)cos / + yfCi _ C3)sin t + e.
2° . Метад варыяііыі адвольных пастаянных. Прыменім метад,
180
апісаны ў папярэднім параграфе, да сістэмы (30), у якой матрыца A(t)^A .
dx dy 1
Прыклад 2. Рашыць сістэму — = у, ~ = х +.
dt dt cost
Рашэнне. Агульнае рашэнне адпаведнай аднароднай сістэмы dx dy
— =у—=—х мае выгляд x = CIcost+C>sint, y = C] sin t+C2 cost dt dt
(гл. § 2). Рашэнне зыходнай сістэмы будзем адшукваць у выглядзе x = Cj(t)cos t + C2(t)sint, y = C!(t)sint + C2(t)cost. Прадыферэнцаваўшы х, у і падставіўшы ўсе выразы ў зыходную сістэму, будзем мець наступныя ўраўненні для вызначэння C](t) і C2(t)'.
C](t)cos t + C2(t)sin t = 0,
1
C,(t)sin t + C2(t)cos t =,
cos t
sin t
адкуль C,=, C2(t) = l. Праінтэграваўшы, атрымаем
cos t
C\(t) = lr\cost\JrCl,C2(t) = t + C2
i, значыць,
x = C, cos t + C2 sin t + cos t ■ ln\cos /| +1 sin t,
У = —C] sin t + C2 cos t sin t ■ ln\cos t\ + t cos t.
3° . Метад нявызначаных каэфіцыентаў. Гэты метад прымяняецца ў тым выпадку, калі функцыі f/t) складаюцца з сумаў і здабыткаў функцый b0 +bIt+...+bmtm, еш, cos fit, sin fit.
Калі fi(t) = ea‘l\Pni(t)cosfiit + Qmi(t)sinfiA^, to частковае рашэнне лінейнай неаднароднай сістэмы адшукваецца па тых жа правілах, што і для аднаго лінейнага ўраўнення з пастаяннымі каэфі
181
цыентамі (гл. § 2.6) са зменамі, якія вызначаюцца наступнай тэарэмай
Тэарэма 7. Частковае рашэнне сістэмы
дзе A =
ajtj = l,n
< і=і^ ;
х = Ах+ Pm(t)e^,
^т(О~ п вектарны паліном ступені т,\\.—
камплексны лік, мае выгляд
^'^Q^«)cv.
прычым Qm+k(t)~ вектарны паліном з кампанентамі ступені, не большай за т + к, а лік к супадае з кратнасцю Ц, калі ц з’яўляецца коранем характарыстычнага ўраўнення (22) і к = 0, калі ц не корань ураўнення (22).
Невядомыя каэфіцыенты паліномаў Q‘m+k(t) (кампанент Qm+kO)) вызначаюцца шляхам падстаноўкі і параўнання каэфіцыентаў падобных членаў, а сапраўдныя рашэнні выдзяляюцца па той жа схеме, што і раней.
Прыклад 3. Рашыць сістэму ўраўненняў dx , dy
— = 2х+у + 2е, — = х + 2у3е4 . dt ' dt
Рашэнне. Знаходзім карані характарыстычнага ўраўнення
адпаведнай аднароднай сістэмы
2\
1
1 2Х
= 0
або X' 4\ +
+ 3 = 0; X] = І,)^ = 3. Агульнае рашэнне аднароднай сістэмы ёсць х = Q е'+ С2 е3<, у =C] е'+ С2 е3'. Паколькі Ц = / з’яўляецца коранем характарыстычнага ўраўнення аднароднай сістэмы, а ]1 = 4 не з’яўляецца коранем гэтага ўраўнення, то частковае рашэнне, згодна з тэарэмай 7, будзем адшукваць у выглядзе
182
x*=(A,+ A2t)e‘ + A3e4', y* = (Bj + B2t)e‘ + B3e4‘. Падстаўляючы гэтыя выразы ў дадзеную сістэму, атрымліваем ураўненні для вызначэння невядомых каэфіцыентаў Ait В,,і = 1,2,3: — A] + А2 — B] — 2, A] + Bj — B2 = 0, A2 + B2 = 0,2 A3 ~ B3 = 0, 2B3 — A3 = 0. 3 гэтых ураўненняў знаходзім A] =0, A2 =1, A3 = —l, B] = 1, B2 = 1, B3 = 2. Частковае рашэнне неаднароднай сістэмы х* = te' e4t ,у* = —(t + l)e' 2е4', а яе агульнае рашэнне (гл. тэарэму 3) запішацца так:
х = С,е' + С2е3‘ + te' е4,.у = С,е' + С2е3‘ ~(t + l)e‘ 2е4'.
Прыклад 4. Устанавіць від частковага рашэння сістэм:
x = 4xy + e3‘(t + sint), х = хy + 8t,
а) j> = x + 2y + te3' cost. б) ў = 5ху + Г +sin2t.
Рашэнне. а) Знаходзім карані характарыстычнага ўраўнення адпаведнай аднароднай сістэмы \1—\2 = 3. У зыходнай сістэме ц = 3, a + ф = 3 + і. 3за таго, што Ц = 5 з’яўляецца двухкратным коранем характарыстычнага ўраўнення (т = 1, к = 2), a a + z'P = 3 + z не з’яўляецца коранем гэтага ўраўнення (т = 1), то частковае рашэнне запішацца ў выглядзе х* =(a/t3 + Ь^2 +c1t + dl)e3‘ +(k1t + l1)e3t sint + fmjt + nje3' cost, y* =(a2t3 + b2t2 +c2t + d2)e3' + (k2t + l2)e3' sin t + (m2t + n2)e3' cos t. б) Карані характарыстычнага ўраўнення тут роўныя XJ 2 = +2і. Для дадзенай сістэмы у = 0 не з’яўляецца коранем характарыстычнага ўраўнення (w = 2), лік a ±ф = ±2і з’яўляецца коранем гэтага ўраўнення (к = 1, т = 0). Таму частковае рашэнне дадзенай сістэмы будзе мець выгляд
х* =a/t2 +b]t + C] + (d, + klt)cos2t + (fl + l!t)sin2t, y* =a2t2 +b2t + c2 +(d2 + k2t)cos 2t + (f2 +I2t)sin2t.
183
4° . Пабудова інтэгравальных камбінацый. Метад Даламбера.
Разгледзім гэты метад у выпадку сістэмы двух ураўненняў dx , dv
— = alx + bIy + fI(t), ~ = а2х + b2y + f2(t). Памножым другое
ўраўненне сістэмы на нейкі лік X і складзём пачленна з першым ураўненнем. Атрымаем
d(xy\y)
—= (al+ka2)x + (bl +Xb2)y + f1(t) + Xf2(t).
Перапішам апошняе ўраўненне ў выглядзе (aj +Х а2 & 0)
d(x + "ky) dt
л f b.+kb2 ]
= (а,+Ха2) х + —~Г^У + /і(0 + ^2(0(35) k ax+ka2 )
■ і b] + kb2
Выберам лік А так, каб = A . Тады (35) прывядзецца да a ] + ka2
ўраўнення, лінейнага адносна х + ку, d(x + ky)
= (а,+)ш2)(х + \у) + f,(t) + Xf2(t), інтэгруючы якое, будзем мець
x + ky = е(а^^С+ \ (f,(t) + \f2(t))e(a^}'dt). (36)
Ь/ + kb)
Калі ўраўненне ;= А мае розныя сапраўдныя карані X, і
а1 + ка2
Х2> то з (36) атрымаем два незалежныя інтэгралы дадзенай сістэмы, а, значыць, інтэграванне сістэмы будзе закончана.
Прыклад 5. Рашыць метадам Даламбера сістэму dx . т , dy ' v — = 2х + 2у + е , — = х5у + е . dt 'dt
Рашэнне. Памнажаем другое ўраўненне сістэмы на X, і складваем пачленна з першым:
— (х + ку)(2к)хУ(25к)уУе' + ке2' dt
184
або
d л ( 2~5^ \
—(x + Xy) = (2h) х + ———yl + e +Xe2. (37)
• at \ — 2 — h J
Выбіраем X так, каб каэфіцыент пры у у квадратнай дужцы быў
роўны X , гэта значыць ——~ = X або I.2 — 31. + 2 = 0, адкуль — 2 — К
Х; =1,1.22Пры Х7 = / з (37) атрымліваем
^ІУ).=.3(і+у)^+г.
dt
Значыць, згодна з формулай (36), будзем мець
х+у = е~3' ((J] + \(е' + е2' )e3,dt\
Пасля інтэгравання знаходзім
x + y = CIe~3t +^е' +^e2t. (38)
Пры Х2 =2 з (37) аналагічна запішам
х + 2у = С2е~4‘ +^е' +^е2'. (39)
Вырашаючы сістэму (38), (39) адносна х і у, атрымаем агульнае рашэнне дадзенай сістэмы
х = 2С1е~3‘ C2e~4t + 0,Зе' + ^е2',
2
у = С^3' + С2е~4'0,05е' + — е2'.
dx
5°. Матрычны метад. У сістэму — = Ax + f (t) увядзём замест х новы слупок невядомых функцый z па формуле
х = eAlz. (40)
Дыферэнцуючы (40) і падстаўляючы ў зыходную сістэму, атрымаем
185
dl
= f(t). Адкуль z(t) = C+^e Ax f (x)dx. Вяртаючыся na
формуле (40) да x, атрымаем
x = e
At
C + \eAxf(x)dT
= eAl
C^\eA('~z) f(x)dx,
<0
дзе C слупок з адвольнымі пастаяннымі элементамі. 3 улікам пачатковай умовы x(t0) — х°, з апошняй формулы атрымаем
С = е Аі° х° , і рашэнне задачы Кашы можна запісаць так: х = еА(1~'о)х° +^еА(,~х>f (x)dx.
*0
Прыклад 6. Рашыць задачу Кашы:
X/ = 5Х] х2 4х3 +е',
• х2 = 12Х] + 5х2 + 12х3, хІ(2) = 0, х2(2) = 1, х3(2) = 0.
х3 = 10х/ Зх2 9х3,
Рашэнне. Уласнымі значэннямі матрыцы А будуць ХІ2 = 1, \3 =1. Уласнаму значэнню Х/2 = 1 адпавядае адзін уласны вектар, таму што ранг матрыцы (А~К12Е) = ( AЕ) будзе роўны 2.