Қазақстан республикасы жоғары оқу орындарының Қауымдастығы к. Д. Көлекеев К. Ж. Назарова дифференциалдық теңдеулер алматы, 2012


§25. Орнықтылықты бірінші жуықтау арқылы зерттеу



Pdf көрінісі
бет34/44
Дата18.10.2023
өлшемі1,36 Mb.
#118678
1   ...   30   31   32   33   34   35   36   37   ...   44
Байланысты:
kolekeev-differencialdyk

§25. Орнықтылықты бірінші жуықтау арқылы зерттеу
1.Теңдеулер жүйесінің
i
i
n
dx
f t x
x
i
n
dt
1
( , , ..., ),
1,
=
=
(1)
тыныштық нүктесін 
n
i
x
i
,
1
,
0
=

орнықтылыққа зерттеуді, 
осы нүктенің төңірегінде 
i
f
функцияларының сызықтық құрам-
дарын, мысалы Тейлор формуласымен бөліп
n
i
ij
j
i
n
j
dx
a t x
t x
x
i
n
dt
1
1
( )
( , , ..., ),
1, ,
ψ
=
=
+
=

(2)
сызықтық теңдеулер жүйесінің
n
i
ij
j
j
dx
a t x
i
n
dt
1
( ) ,
1,
=
=
=

(3)
тыныштық нүктесін
n
i
x
i
,
1
,
0
=

орнықтылыққа зерттеумен 
алмастырады.
Мұндағы
(
)
n
i
x
x
t
...,
,
1
1
ψ
функциялары бірінші реттен жо ғары 
шексіз аз, яғни
,
,
1
,
)
(
)
,
(
n
i
x
x
x
t
i
=

γ
ψ
(4)


209
,
0
)
(

x
γ
егер 
2
2
1
...
,
0
n
x
x
x
x
+
+
=

.
Ляпунов теоремасы.
(
)
n
i
x
x
t
,
,
,
1

ψ
функциялары (4) шарт-
ты орындайтын теңдеулер жүйесі (2) қарастырылады.
Егер 
( )
ik
a
матрицасының барлық меншікті мəндерінің нақты 
бөліктері теріс таңбалы болса, онда (2) жүйенің нөлдік шешімі 
асимптотикалық орнықты; ал егер кемінде бір меншікті мəнінің 
нақты бөлігі оң таңбалы болса, онда жүйенің (2) нөлдік шешімі 
орнықсыз.
1-мысал.
Теңдеулер жүйесінің тыныштық нүктесін x=0, y=0 
орнықтылыққа зерттеу керек.
dx
y x x
y
t
dt
dy
x y x
dt
2
2
2
sin ,
.

= − +

⎪⎪


= + −
⎪⎩
Шешуі. Жүйенің сызықты емес мүшелері теореманың шарт-
тарын орындайтындықтан, сызықты жүйенің
dx
x y
dt
dy
x y
dt
,
.

= − +
⎪⎪


= +
⎪⎩
тыныштық нүктесін орнықтылыққа зерттейміз.
Сипаттамалық теңдеуінің
0
1
1
1
1
=



k
k
түбірлері 
2
±
=
k
болғандықтан, жүйелердің тыныштық нүк- 
тесі орнықсыз.
2-мысал.
Жүйенің нөлдік шешімін орнықтылыққа зерттеу 
керек.


210
(
)
x
dx
y e
dt
dy
y
x
dt
3
3
ln 4
,
2
1
1 6 .


=
+
⎪⎪


=
− +

⎪⎩
Шешуі: 
)
1
ln(
)
(
x
x
f
+
=
функциясының Тейлор формуласы 
бойынша, 
(
)
x
x
y e
y e
3
3
ln 4
4
1,


+

+

ал 
,
3
1
3
x
e
x



онда
(
)
x
y e
y
x
3
ln 4
4
3 ;

+


α
)
1
(
)
(
x
x
f
+
=
функциясының Тейлор формуласы бойынша,
x
x
2
1
)
6
1
(
3
1



, онда
.
2
2
6
1
1
2
3
x
y
x
y



+

Бұл жіктеулердің қалдық мүшелері теорема шарттарын орын-
дайды, сонымен бірінші жуықтау жүйесі
dx
y
x
dt
dy
y
x
dt
4
3 ,
2
2 ,

=

⎪⎪


=

⎪⎩
сипаттамалық теңдеуінің 
0
2
2
4
3
=




k
k
түбірлері 
.
2
7
2
1
i
k
±

=
Нақты бөліктері теріс таңбалы, демек, жүйелердің нөлдік 
шешімдері орнықты.
3-мысал.
Жүйенің тыныштық нүктесін 
0
,
0
=
=
y
x
орнық-
тылыққа зерттеу


211
dx
x
y
dt
dy
x y
dt
3
3
9 ,
.

= − +
⎪⎪


= − −
⎪⎩
Шешуі:
Бірінші жуықтау жүйесінің 
dx
y
dt
dy
x
dt
9 ,

=
⎪⎪


= −
⎪⎩
сипаттамалық теңдеуінің
0
1
9
=



k
k
түбірлері жорамал сандар 
i
k
3
±
=
, болғандықтан, бұл əдіспен 
орнықтылыққа зерттелінбейді.
Есеп үшін Ляпунов функциясын 
2
2
9
y
x
+
=
ϑ
құру оңай:
1) 
;
0
)
,
(

y
x
ϑ
2) 
(
)
(
)
(
)
d
x
x
y
y
x y
x
y
dt
3
3
4
4
2
9
18
2
9
0
ϑ
=
− +
+
− −
= −
+

,
0
9
4
4
>
+
y
x
егер 
0
,
0


y
x
болғандықтан, жүйенің тыныш-
тық нүктесі асимптотикалық орнықты.
Бұл мысалдағы бірінші жуықтау жүйесінің координаталар 
бас нүктесі центрі еді. Бастапқы жүйедегі сызықтық емес мү-
шелері, центрді орнықты фокуске айналдырды. 
Жалпы жағдайда осы сияқты күрделірек геометриялық су-
реттер орын алады.
Сызықты емес мүшелер бас нүктенің төңірегінде аз шамалар 
болғанымен, бағыттар өрісін аз да болса өзгертеді. Нүктеден 
(
)
0
0
,
y
x
шығатын траектория бас нүктеден айналғанда сызық-
тық жүйедегіден ығысады. Осындай айналу нəтижесінде барлық 


212
траекториялар бас нүктеге жуықтаса, онда бұл нүкте орнықты 
фокуске айналады; егер траекториялар бас нүктеден қашықтаса, 
онда орнықсыз фокус шығады.
Төңірегіндегі барлық траекториялары спиралдар болатын 
тұйық траекториялар
шектік циклдер
 
деп аталады.
Егер 


t
ұмтылғанда спирал траекториялар шектік циклге 
ұмтылса, онда шектік цикл 
орнықты
; егер шектік циклге жақын 
спиралдар 


t
ұмтылғанда, одан қашықтаса, онда шектік цикл 
орнықсыз
; егер шектік циклдің бір жағындағы траекториялар 
оған жуықтап, екінші жағындағылар одан қашықтаса онда шектік 
цикл 
жартылай орнықты.
11-сурет
орнықты циклді, 
12-сурет
жартылай орнықты циклді көрсетеді.
11-сурет 12-сурет
2. Көпмүшеліктің барлық түбірлерінің нақты 
бөліктерінің теріс болуы
Дифференциалдық теңдеулер жүйелерінің кең жиыны ше-
шімдерін орнықтылыққа зерттегенде сипаттамалық теңдеуінің 
нақты бөліктерінің таңбалары қаралатынын байқадық.
Егер сипаттамалық теңдеу жоғары дəрежелі болса, онда оны 
шешпей-ақ, түбірлерінің барлығының нақты бөліктері теріс таң-
балы болатындығына көз жеткізуге болатын əдістердің орны 
ерекше.
Гурвиц теоремасы
. Коэффициенттері нақты көпмүшеліктің
n
n
n
n
a
z
a
z
a
z
+
+
+
+


1
1
1
...
барлық түбірлерінің нақты бөліктері теріс таңбалы болуы үшін, 


213
Гурвиц матрицасының бас диагоналдық барлық минорларының 
оң таңбалы болуы қажетті жəне жеткілікті.
Гурвиц матрицасының бас диагоналында көпмүшелік коэф-
фиценттері 
n
a
a
a
,...,
,
2
1
орналасқан. Бағандары алма-кезек тек 
тақ немесе жұп коэффициенттерден, 
1
0
=
a
деп, 
k
i
ik
a
b

=
2
фор-
муласымен құрылады. Барлық жетпейтін, яғни индекстері n-нен 
төмен немесе жоғары коэффициенттері 0-мен толтырылады.
Гурвиц матрицасының бас диагоналдық минорлары
,
1
1
a
=
Δ
,
1
2
3
1
2
a
a
a
=
Δ
,
,
0
1
3
4
5
1
2
3
1
3



a
a
a
a
a
a
a
=
Δ
n
n
a
a
a
a
a
a
a
a

0
0
0
0
0
1
3
4
5
1
2
3
1
=
Δ
Гурвиц шарты бойынша 
,
0
...,
,
0
,
0
2
1
>
Δ
>
Δ
>
Δ
n
онда 
n
n
n
a

Δ
=
Δ

1
теңдігінен соңғы шарт 
0
>
n
a
болуы талабымен 
анықталады. 
Гурвиц теоремасының екінші, үшінші, төртінші дəрежелі 
көпмүшеліктерге қолданылуын қарастырайық.
а) 
2
1
2
a
z
a
z
+
+
, Гурвиц матрицасы 
,
0
1
2
1
⎟⎟⎠

⎜⎜⎝

a
a
,
0
1
1
>
=
Δ
a
,
0
2
1
2
>
=
Δ
a
a
.
0
2
>
a
Гурвиц шарттары 
0
,
0
2
1
>
>
a
a
келтіреді. Сипаттамалық 
көп мүшелігі 
2
1
2
a
z
a
z
+
+
болатын теңдеулер жүйесінің ты-
ныш тық нүктесінің асимптотикалық орнықты аумағы 
13-сурет-
тегі
бірінші ширек екендігіне көз жеткіздік.


214
ə)
.
3
2
2
1
3
a
z
a
z
a
z
+
+
+
Гурвиц матрицасы 
,
0
0
0
1
3
1
2
3
1










a
a
a
a
a
онда 
0
,
0
,
0
3
3
2
1
1
>
>

>
a
a
a
a
a
болуы керек. Коэффициенттер 
кеңістігінде бұл теңсіздіктер аумағы 
14-суретте
көрсетілген
13-сурет
14-сурет
б)
4
3
2
3
1
4
2
a
z
a
z
a
z
a
z
+
+
+
+


215
Гурвиц матрицасы 
⎟⎟





⎜⎜





4
2
3
4
1
2
3
1
0
0
0
0
1
0
0
1
a
a
a
a
a
a
a
a
,
Гурвиц шарттары
(
)
.
0
,
0
,
0
,
0
4
4
2
1
3
3
2
1
3
2
1
1
>
>


>

>
a
a
a
a
a
a
a
a
a
a
a
Қарастырылған көпмүшеліктер үшін Гурвиц шарттары ың-
ғайлы, ал жоғарғы ретті көпмүшеліктер үшін қолайсыз. Көп-
мүшелік түбірлерінің нақты бөліктерінің теріс таңбалы болуы-
ның басқа да белгілері қолданылады.
4-мысал.
Жүйенің
dx
y
x
dt
dy
ax by
dt
sin ,

= +
⎪⎪


=
+
⎪⎩
нөлді шешімі
b
a
,
параметрлерінің қандай мəндерінде орнықты?
Шешуі. Жүйенің бірінші жуықтау жүйесі:
dx
y x
dt
,
= +
dy
ax by
dt
,
=
+
ал бұл жүйенің сипаттамалық теңдеуі 
,
0
1
1
=


k
b
a
k
сипаттамалық көпмүшелігі 
a
b
k
b
k

+
+

)
1
(
2
. Онда Гурвиц 
шарттарынан 
b
b
a b
1 0, (
1)(
) 0
− − >
+
− >
берілген теңдеу-
лер жүйесі 
1

<
<
b
a
мəндерінде асимптотикалық орнықты 
дейміз.




Достарыңызбен бөлісу:
1   ...   30   31   32   33   34   35   36   37   ...   44




©emirsaba.org 2024
әкімшілігінің қараңыз

    Басты бет