Сызықтық жүйелердің орнықтылығы

Қарапайым жағдайдан бастайық – күй кеңістігіндегі қоздырылмаған ашық жүйеден. Басқару жоқ болған кезде жүйе ашық деп аталады

 

, , (5.5)

 

мұндағы A – n×n-өлшемді, t уақытына тәуелсіз матрица, егер кез келген x(0) үшін кезде болса орнықты деп аталады.

Теорема 5.1. (5.5) жүйесінің орнықтылығы үшін А матрицасының барлық si меншікті мәндері сол жақ жарты жазықтықта жатқаны қажетті және жеткілікті:

 

. (5.6)

 

Бұл жағдайда кез келген үшін сондай тұрақты С = C(A,ν), бар болады, келесі теңсіздікті қанағаттандыратын

 

(5.7)

 

Атап өтейік, (5.6) шартын қанағаттандыратын А матрицасы осыдан әрі қарай гурвицтік немесе орнықты деп аталады. Осыған сәйкес, теріс нақты бөлікті меншікті мәндерді де орнықты деп атайтын боламыз.

Егер А орнықты болса, онда шамасы (матрицаның немесе сәйкес жүйенің) орнықтылық дәрежесі деп аталады, яғни ол орнықты матрицаның меншікті мәнінен жорамал оське дейінгі минималды қашықтық.

Дәлелдеу.Қажеттілігі. (5.6) шарты орындалмайды деп ұйғарайық, яғни кейбір меншікті мәндері, мысалы s1 сондай

 

 

болатын А матрицасы табылады. Егер s1 нақты болса, онда аламыз, мұнда , e1 - нақты вектор және . Онда мұндай бастапқы шартты шешімнің түрі болады және сондықтан болғанда x(t) 0. Егер s1 – комплексті сан болса, онда түйіндес меншікті мән табылады, және оған қоса, егер - меншікті вектор болса, s1 меншікті мәніне жауапты, ал - меншікті векторы үшін, s2 меншікті мәніне жауапты. Бұл жағдайда . аламыз, мұнда , онда x(t) сол екі өлшемдік төменгі кеңістікте қалады: , мұнда келесі дифференциалдық теңдеулермен сипатталады

 

 

белгілеуі арқылы , оңай аламыз, яғни және u > 0 болғанда p(t) 0.

Жеткіліктілігін дәлелдеу (матрицалық экспоненттерді бағалау). А матрицасы диогональдық түрге келтірілетін (мысалы, оның барлық меншікті мәндері әртүрлі нақты мән қабылдайды) болсын деп ұйғарайық, яғни сондай азғындалмаған Т матрицасы бар болады, келесі теңдіктер орындалатын

 

.

 

Бұл жағдайда аламыз

 

.

 

Осы орнықты матрица үшін баға (5.7) дәлелдейді. Бұл жағдайда , сондықтан v = σ деп алуға болады, ал С тұрақтысы айқын түрде жазылады

 

.

 

Матрицалық тұлғада алынған нәтижені матрицалық экспонентаның бағалануы түрінде жазуға болады

 

. (5.8)

 

А матрицасы диогональдық түрге келтірілетін жағдай өте қарапайым, себебі түріндегі айнымалыларды алмастыруды жүргізсек, жүйені келесі түрге келтіруге болады

 

 

немесе басқаша айтқанда, , яғни (5.5) шешімдері

 

 

болатын тәуелсіз n теңдеуге ыдырайды.

Егер және барлық болса, онда барлық si үшін және болады, сондықтан .

Жалпы жағдайда А матрицасы ұқсастық түрлендірумен диогональдық емес, блоктық жордано тұлғасына келтіріледі

 

,

 

және жүйе Ji блоктарына сәйкес төмендегі тәуелсіз m жүйеге жіктеледі. Бір Ji жордано блогы үшін келесі l × l өлшемдік матрицаны аламыз

 

.

 

.

 

Сондықтан, егер болса, онда l – 1 дәрежелі сондай R(t) палиномы табылады, кез келген үшін теңсіздігі орындалатын; мұнда . шамасын болғанда ескере отырып, (5.7) бағалауына келуге болады.

Айта кету керек, (5.6) шарттарын қанағаттандыратын (5.5) жүйесі орнықты болатындығына қарапайым талдаудың да көмегімен көз жеткізуге болады.

Теңдеу (5.5) үшін сипаттамалық теңдеуді келесі түрде жазуымызға болады:

 

. (5.9)

 

Сипаттамалық теңдеу (5.9) ізделінетін s мәнін табуға болады, яғни s1, s2, ... sn – олар сипаттамалық теңдеудің түбірлері болады.

(5.5) теңдеуінің шешімдерін алуымызға болады

 

. (5.10)

 

Басқаша айтқанда, шешімнің уақыттан тәуелділігі экспоненциалды функциямен сипатталады. Сипаттамалық теңдеулерінің түбірлері немесе А матрицасының меншікті мәндері уақытқа тәуелсіз w1, w2, …, wт меншікті векторларды табуға мүмкіндік береді, ал олар өсу немесе кему жылдамдықтарын береді. Сi – алғашқы шарттары арқылы анықталады. Атап өту қажет, жалпы жағдайда s1, s2, ... sn шамалары комплексті мәнді қабылдайды, ал w1, w2, …, wт векторлары x(t) айнымалысының векторлық сипаттамасын ескереді.

Сонымен, орнықтылық проблемасы шешіледі, егер s1, s2, ... sn меншікті мәндері табылатын болса, ал w1, w2, …, wт меншікті векторлардың бұл жағдайда қандай да қасиеттерінің болуына байланыссыз. Шынымен, x уақытқа тәуелді (5.10) қатынасымен анықталатын болғандықтан, жазамыз

 

.

 

Мұнда Resi және Imsi – сәйкесінше si комплексті шаманың нақты және жорамал бөлігі. Осыдан алымыз, Resi < 0 болса, x мәні экспоненциалды уақыт бойынша кемиді. Бұл өшу тербелістерімен модуляцияланады, егер Imsi нольге айналмаса. Осыдан аламыз, қоздырылмаған (стандартты) шешімді (күйді) аламыз, болғанда x = 0 шегіне жететін. Бұл сызықтандырылған жүйенің асимптотика-лық орнықтылығын білдіреді. Екінші жағынан, егер Resi > 0 болса, онда қоздырылған x қозғалысыэкспоненциалды өсетін болады, оған сәйкес x = 0 стандартты шешімі орнықсыз болады.