Файл: Красовский, А. А. Фазовое пространство и статистическая теория динамических систем.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.