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

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

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

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

Добавлен: 23.10.2024

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

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

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

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

Тогда

 

 

 

п

 

 

 

 

 

{_

 

 

 

 

 

2 Ац\ кк'З-И’З'кк' X

 

 

 

2

 

 

 

г, V

 

 

 

X dxц . .. dxnn — | М | exp

s

Мц',кк'Уц.У

 

 

 

 

 

к, k'=l

 

Таким

образом,

 

 

П

 

 

 

 

 

 

Рх =

С2р ст (*г, .. Хп , t) exp

(jg

2

 

№Ун,укк,)

 

 

 

г,

г', /с,

Я*'=:1

'

 

 

 

 

 

 

(3.85)

 

У ц ’

— 2

-^го, /to, l l ' 3

' i % k i

 

где

 

i, /i=l

 

 

 

 

 

 

 

 

 

 

Рст =

е^Р у

2 Лг/.А'о^л)

 

 

 

'

t, /t=x

 

1

 

— распределение в стационарной (без параметрических шумов) системе. Коэффициенты A i0)k0>ц> определяются

уравнениями (3.83), которые с учетом выражений (3.69) можно представить так:2

го, ко, И'

2 (“ piApo, ко, iv ~\~ йркАро, го, г/')

 

п

р= 1

п

 

 

 

 

 

2 ^ Р Р '.И 'А р р ’, г0 , /СО

2 Sp p ’, <т'А(

 

Р, Р '= 1

 

р, р ', «, <г'=1

 

 

2

^ро, (р (-^QO, кО^РО, io, IV

Аро, ioAqQ, ко, iv)

 

р, 4=1

= Хл'/lji), ко Mki'AiQ' 10,

(3.86)

 

 

где

1 при v = р,

Иу|1 = |

Опри v=^p.

Итак, в рамках рассматриваемой «кубической теории» (аппроксимация совместной логарифмической плотности


§ 3.3] Д И Н А М И К А Л И Н Е Й Н Ы Х СТО Х А СТИ ЧЕСКИ Х СИСТЕМ 129

вероятности полиномом третьей степени) динамика стоха­ стической системы практически полностью определяется путем решения системы уравнений (3.86). Если все п2

параметров являются случайными, то порядок системы линейных дифференциальных уравнений (3.86) с учетом симметрии коэффициентов A i0, h0ilV равен «4/4! и при боль­ ших п весьма велик. Однако обычно число случайных (флук­ туирующих) параметров значительно меньше га2 и боль­ шинство коэффициентов A i0l h0, и' тождественно равно

нулю. Это сразу упрощает систему уравнений (3.86). Для установившегося режима уравнения (3.86) обра­

щаются в линейные алгебраические. Устойчивость реше­ ния уравнений (3.86) определяет устойчивость равновес­ ного распределения в исходной стохастической системе для случая устойчивости этой системы без шумов.

Действительно, пусть система без параметрических шумов устойчива. Тогда при любых начальных условиях величины А * „, ь о стремятся при t —*■оо к определенным

конечным значениям. Если при этом решение системы ли­ нейных уравнений (3.86) устойчиво, то величины И го, ko,w

также стремятся к определенным конечным значениям. Таким образом, при указанных условиях распределение вероятности р х (3.85) стремится с течением времени к

стационарному (не зависящему от времени) распределе­ нию и налицо статистическая устойчивость.

Если же при прочих равных условиях решение урав­ нений (3.86) неустойчиво, то коэффициенты A t0, ho, и' не-

ограниченно нарастают с течением времени. Квадратич­

ная форма

П

2 Мц.'№ уи,У]гк.

i, V, К, fe'=l

квадратичных форм

П

У)У — 2 it), ко,

г, к = 1

положительно определена, и неограниченное нарастание Л го, ft0, вызывает неограниченное нарастание плотности вероятности на периферии фазового пространства. Имеет место статистическая неустойчивость стохастической сис­ темы.

5 А. А. Красовский


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

Следует еще раз отметить, что все это справедиво в

рамках «кубической теории».

 

 

Рассмотрим простой пример.

 

П р и м е р . Стохастическая

система первого

порядка

(гг =

1 ) описывается

уравнениями

 

х\

~f" ^ iA +

 

 

 

Д®ц + a ii,ii^ aii

= £и

или

 

 

 

 

 

 

хм “Ь ^ii-^io

хн хю

^ю,

хп ~1~"^н, н -^ii

^н»

где

З-Ю=

ХП

Х11

~ Д^Ц! £ 10 = El*

 

 

 

Система уравнений (3.86) в данном случае сводится к од­ ному уравнению

■^1о,ю,11

2au.4i0)10lll — а н ,li-^n,Ю)Ю

^u, U-4и, 11

X

 

X-^пчо,1о, — 2510)10Л10, хо^хо, 10, и = 24хо,1о-

(3.87)

Коэффициенты 4 10>10, 4 11)П,

определяемые

согласно из­

ложенному

для линейной стационарной системы, в уста­

новившемся режиме равны

л

 

 

 

 

а

_

1

2аlx

_

1

____2а»1,п

10,10

 

./Июдо —

Л'ю,ю ’

 

11,11

3/и,11

^11,11

Подставляя эти выражения

в (3.87), получаем

 

 

 

 

 

 

 

 

а п

 

 

-<4W, ю , и +

( 2 в ц ■+ И ц ,

н ) ^4 ю , ю , 11

=

4 iio , ю

 

или для установившегося режима

а

 

4 а ц

идо

 

ю, 10, и

 

2ап + ап.и

 

Формула (3.85) в данном случае имеет вид

Рх = сгРстЮ exp

Ми, пУхх), Ун =

^ю, ю, пх1-

Таким образом,

 

 

 

 

 

 

, .

Г1

 

16“il

^

Х*

J

р х - с2рст(жх)exp

|fg 2aii ii (2-n + anii)2

 

 

 

’4

i'll

a2 Z4

c2Pct(xi) exp

alin

. (3.88)

9

s b>,10 (2au +

an,ii)2an,n


§ З.з1 ДИНАМ ИЙА Л И Н Е Й Н Ы Х СТО Х А С ТИ ЧЕСКИ Х СИСТЕМ 131

Эта формула в рамках кубической теории дает исчерпы­ вающее представление о влиянии флуктуации параметра в установившемся режиме.

Случай линейно зависимых флуктуирующих парамет­ ров. В рассмотренной задаче предполагалось, что пара­ метрические шумы хотя и зависимы друг от друга, но их корреляционная матрица неособая. Между тем важным

частным случаем является тот,

когда все флуктуирующие

параметры пропорциональны

одной случайной функции

у] (t):

 

Aaik = пг/Ет) (t),

где n ih — постоянные коэффициенты. Сама случайная

функция

т) (t) определяется как выходная величина не­

которого

стационарного линейного фильтра:

 

Т|-f- (ln . П+1Т1 'р +1-1 1, п+2'^'п+2 Ч- •••Ч~ ®п+1, п+т ^Пгт

^п+1»

#п+2 4“ ®n+2, n+l4 + ^n+2, n+i^n +i Ч- •••+ ®п+2, п+т 3-п+т ~

Еп+2>

•^п+тЧ &п+тг п+ 1^1 “Ь ^пЧт, п+2+г-г2 + •■ *4“^п+т, п+т ^п+т — in+mi

где 1п+1, . . ., |п+т — белые шумы. В данном случае нет

смысла вводить координаты с двойными индексами. Луч­ ше обозначить ц = ж„+1. Тогда уравнения данной стоха­ стической системы запишутся в виде

пп

“Н

^ik^k Ч~

Щк^к^п+1

i — 1 , 2, . . ., Н,

к—1

к—1

 

 

 

 

 

n+m

 

 

 

 

 

•Г«+1 Ч"

®n+i, к^к

~

?n+i> i =

1» 2 , . . .,

Ш,

 

)[=п+1

 

 

 

 

 

ИЛИ

 

 

 

 

 

 

 

n+m

 

n + т

 

 

 

Ч" S

к^к Ч"

S

®г> к, n+ l^k'l'n + l=

( 3 . 8 9 )

 

к=1

 

к=1

 

 

i = 1 , 2, . . ., п f т,

5*


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

где

 

 

a ik

при

п,

i ^ n и

 

ащ — •

при

i > n ,

к ^ > п ,

 

 

0

во всех других

случаях

 

 

к

при

i < ^ n ,

(3.90)

 

 

 

 

 

0

при

 

J

 

 

0

при

i ]> п .

Уравнения (2.14) для данного случая имеют вид

 

n-|-m

 

n-fm

 

n

A-i

2 ®3>i-^i

2

^ P5 (Apqi 4" Apj^Aq) = 2 2 a PPii

 

P=1

 

P, 5=1

 

P=1

 

n-fm

 

 

n-fm

 

A ik

2

(a Pi-^pk 4“ a pk-^pd — 2 2

®Pi)c-^P

 

P=1

 

 

p= lj

 

 

 

n-fm

 

 

 

 

2 *^Р5 (3 A p q ik + 2A pikA g 4" A piA q k ) = 0,

 

 

P> 5=1

 

 

 

 

 

n-\~m

 

 

 

 

Aiki

(ttpiApkl + aPkApU+ dpiApifi)

 

 

p = 1

 

 

 

n+m

2 (a P ik A p i + а р ц А р к + U pklA pi)

P=1 n-fm

2 S p q (QApqiki + Ъ А р А ч т 4 - A piA q ki 4"

P. 5=1

+ АркАдц 4“ ApiAgik) = 0,

(3.91)

При определении логарифмической плотности с точностью до кубических членов следует полагать

A p q ik = Apqiki = •. . = 0 .

Кроме того, при центральном начальном распределении и 2 appt = 0 имеет место решение A t(t) = 0.