Файл: Красовский, А. А. Фазовое пространство и статистическая теория динамических систем.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 23.10.2024

Просмотров: 102

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

i 3.4]

 

СТАТИСТИЧЕСКОЕ

ИССЛЕД О ВА Н И Е

СИСТЕМ

143

 

 

 

3!

_21 2

(a Pikl-Apm q ... s

•••

 

 

 

(TV—1) (TV—2)

 

 

TV-2

 

 

 

 

 

 

 

 

 

P=1

 

 

 

 

 

•••+ a p frsA pik ... d)

•••

TV — 1 2

^7>X

 

 

 

 

 

 

 

 

 

 

P=1

 

 

 

x (api/f...T^Prs + • • • + U p f r ^ j A p i k ) ~

 

 

 

(3.112)

 

 

^ V -Г "

 

 

 

JV-1

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

2

 

(apiJc... r^lpe + •••+

aPT/e rA p i)

 

 

 

 

P = 1

'— f T '

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

n

 

 

 

 

 

- TV 2

 

 

 

- 4 - 2 S p [ N ( N - i ) x

 

 

 

p = i

'" Ж Г

 

P=1

 

 

 

 

 

X Appjfc „ s 2iV^4p^4pifc... s "t- (■^pj-^pjct... a ”b ■• •

 

 

 

- i v + r

 

 

JV+l

 

N

 

 

 

 

• • •

+

^ps-4«ift... r)

+

2j

 

 

+

• . •

 

 

~n ~ T

(A p a c A p i^ j

 

 

 

 

 

iv

'

 

31

 

'~ N -l

 

 

 

 

“b -^ p rs-^ p is... f)

 

 

 

 

 

 

 

 

~t~ // у _-i) (/v __ ^ ( А р у д А д г... t -j~ •••

 

 

 

 

 

'-^ T

 

 

 

 

 

JV-2

 

 

 

• • • +

A p!rsA q i , . . ( )

+

• • • +

{Ар-,*... r-Aps +

• • •

 

 

 

 

 

TV-2

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

TV

 

 

 

 

 

 

+

A pk"_sA pi) \ = N (tf +

l > 2 kpV-ppjk...i-

 

 

 

 

'~N~'

J

 

P=1

"л+Г'

 

 

Непосредственной

подстановкой убеждаемся,

что при

S t =

k j a

эти уравнения имеют решение

 

 

 

А о =

0,

A t = 0,

 

=

2aaih,

A tkl =

— 2ааш ,...

 

 

 

. .

Aiki ... s —

2ааш ... g, • .

•»

(3.113)

что соответствует решению (3.111). Для исследования устойчивости этого решения по Ляпунову производим линеаризацию уравнений (3.112) в окрестности это­ го решения, т. е. записываем уравнения первого


144 РЕ Ш Е Н И Е Ф П К -У РА В Н Е Н И Я МЕТОДОМ РЯ Д О В 1ГЛ. I I I

приближения:

П

АЛ0-^Г 2

kpS.App — О,

 

 

 

 

р = 1

 

 

п

 

 

 

 

«

 

 

 

 

 

b * i +

2 kpapiAAp -i-

2 крЬАрр{ =

О,

 

 

Р~ 1

 

 

 

Р —1

 

 

 

П

 

 

 

 

 

 

Д-^ik +

2

 

 

+

apfrA-4j>i) +

 

 

р=1

 

 

 

п

 

 

 

 

п

 

 

 

 

 

 

+

2 2

kpdpi^Ap

— ^ kpAAppili — О,

 

 

p = i

 

 

 

p = i

 

 

A^ifcf...» +

2

(CtPiA'4p)cf...s+

• • •+ apsA-ApiA...r)+

n

 

p = 1

 

 

 

 

V'7 T ^

2!

П

 

 

 

 

 

 

2 ^p(apa-^ylp/m...8+

• • •+

o,prsAAPi]i ^f)A-

+ ~дг __ j

 

 

l

 

' N

- l

 

'

'

 

 

3!

 

П

 

 

 

(3.114)

+

 

 

 

 

 

 

 

(ДГ_1) (ДГ — 2) 2 kp (apikl&Apd...$ + • • •

 

 

 

 

33=1

 

 

 

 

 

 

 

• • ~\~ ^PfrsJ^pik.,.m) +

. - .

 

 

 

 

 

 

 

N—2

 

• • • + ДГ_ 1 2

 

 

 

+ • • •

 

 

 

 

P= 1

N

- l

 

 

 

 

 

 

 

n

 

 

 

. . . -f- ttpZm...s

 

+ 2 kp (ары

rAAps + . .

 

 

JV-l

 

P=I

JV

 

 

 

 

 

 

n

 

 

 

•. • +

 

/1jj+pN 2 kpO*pik "9AAp

 

 

 

lv

 

 

p =i

'jv+i'

 

 

 

 

_ ^

 

+ *) V к AA „

- О

 

 

 

 

,i_j repU31pptit...i

— u >

 

 

 

 

 

 

p=1

'lv+ T

 

Рассмотрим случай малых спектральных плотностей шумов, когда величина а велика. В этом случае члены


§ 3.4]

СТАТИСТИЧЕСКОЕ! И ССЛЕДОВАН ИЕ СИСТЕМ

145

с множителем 3/а в уравнениях (3.114) можно отбросить и записать приближенные уравнения:

AAj Ч- 21 JcpQ’piААр — О,

35=1

А

Ч - 2 к р (а р»А^4рк Ч~ «pfcA^lpi) -1

 

p = i

 

Ч" 2 2 крЩлк^Ар — О,

 

p = i

ААш ...8 Ч" 2

^'p(CtPiA4pfci...8 Ч" • • •

N р = 1

JV

• • -Ч" ®PsAApjj; _{) Ч-

п

^^ ^Р fapik^ Aplm^s +

 

 

р=1

 

JV-1

]■ (3.115)

 

 

 

 

 

 

 

 

• • ЫргААрш...!) Ч-

 

 

 

 

 

JV-l

Ч~

O'V-

3!

_ 2 2 kp{0‘pikAApm...t +

 

1) (N - 2)

p=i

J V - 2

.+ Ор/г«АЧрг,,.1(з) Ч“

• ■•

 

•Ч"

2!

"

^р (apife.../A4pr3 ч~

дг

j 2

 

 

р=1

 

N —1

 

• • ■+

a p lm ...A A p ili) Ч"

2 ^Р (a p ik...r^p s 4 “

 

 

 

 

Р=1

 

 

 

•Ч~ ®рм...вЧр|) Ч- Af ^

^p^pik.,.sA4p — 0 .

 

 

"Ч Г "

p = i

' w+I"

Легко показать, что если при kt > 0 в начале координат располагается экстремум-минимум функции F, то триви­

альное решение системы уравнений (3.115) устойчиво, а значит, устойчиво распределение вероятностей (3.111), (3.113). Аналогично, если при kt < 0 в начале координат


146 РЕ Ш Е Н И Е Ф П К -У РА В Н Е Й ИЯ М ЕТОДОМ РЯ Д О В [ГЛ . III

располагается экстремум-максимум функции F, то три­

виальное решение уравнений (3.115) и распределение ве­ роятностей (3.111), (3.113) устойчивы.

Действительно, пусть функция F имеет в начале коор­ динат экстремум-минимум и k t 0. Квадратичная форма

П

~2~ 2 aikxixk i,k=i

при этом положительно определенна и все корни %i ха­

рактеристического уравнения первой группы уравнений

(3.115)

X -f- kia.li

kw n

. .

knara

 

 

 

 

 

 

 

 

kiOn

X -f- klO.22

. .

V n 2

= 0

(3.116)

 

 

 

 

 

* ialn

fc«4n

.

.

X + knann

 

 

отрицательны:

 

 

 

 

 

(3.117)

h

< 0,

I =

1 ,

2, .

 

Как следствие, все переменные ДА,- (t) первой группы

уравнений (3.115) стремятся к нулю при t —>■оо. Перенесем

п

члены 2 2 kp&pih&Ap второй группы уравнений (3.115) в p = i

правую часть и будем рассматривать как правые части неоднородной системы уравнений. Корни характеристи­ ческого уравнения, соответствующего второй группе урав­ нений (3.115), равны попарным суммам величин (3.117):

^

+ kj,

(3.118)

и отрицательны. Отсюда и из стремления к нулю величин ДАДг) следует, что все AAjft( £ ) 0 при *-> оо. Анало­ гично для N-й группы уравнений (3.115), в которой все члены, кроме содержащих переменные с N индексами

(■Alit...»), перенесены в правую часть, корни характери-

~~~N

стического уравнения равны

-f- Xj -f-. . . -f- Хг (3.119)

IF

и отрицательны. Поэтому все ДA y ... „ (t) стремятся к нулю при t-+- 00.


§ 3.4]

СТАТИСТИЧЕСКОЕ И СС ЛЕД О ВА Н И Е СИСТЕМ

147

Итак, при указанных условиях все величины

N

стремятся к нулю при t оо, что означает устойчивость

(в малом) распределения (3.111), (3.113) при достаточно низких спектральных плотностях шумов.

Для исследования устойчивости распределения в боль­ шом, для определения области притяжения в пространстве коэффициентов распределения, необходимо исследование устойчивости рассматриваемого решения нелинейной сис­ темы уравнений (3.112). Это можно осуществить лишь в конкретных частных случаях и в приближенной форме. Иными словами, такое исследование можно выполнить лишь численным методом.

В качестве примера рассмотрим одномерную систему (п = 1). В этом случае уравнения (3.112) принимают вид

Лг к1а11Л1Sj щ

 

АцАт) = 2&1аш ,

 

]

А п

— 2/гхИцЛц — 2/с1а111И1 —

 

 

 

 

 

 

1?! ( З И П 114 - 2 А111А1 +

Ah) — б / ^ а ц ц ,

 

Аи. ..г — N kp (яхДи.л +

ctm Ап

л

 

 

 

 

"ТГ

 

I T

 

iv^i

 

 

 

(3.120)

• • •

+ “и^дИх)-----j- iSi((iV + 1 )

 

Лх^хх,..!-)-

 

 

 

n+i

 

 

 

 

iv+2

iv+x

 

 

+ ^цЛх1_д + • .. + Ихх^хИх] =

 

 

 

 

 

 

'~N~

 

'n+1<

 

 

 

 

 

 

 

 

 

= N (N -\г 1) &рап...1)

 

 

 

 

 

 

 

 

 

 

JV+2

 

Допустим,

что функция F равна

 

 

 

 

т.

е.

F

= х \ -

0,186 х! +

0,01а;?,

 

 

 

 

 

 

 

 

 

 

 

 

ап = 2,

«хххх =

—0,744,

а ц Ш1

 

= 0,06,

 

все остальные осхх---х =

0.

График этой функции представ-