
книги из ГПНТБ / Красовский, А. А. Фазовое пространство и статистическая теория динамических систем
.pdf60 |
РЕ Ш Е Н И Е Ф ПК -У РА В Н Е Й ИЯ |
[ГЛ . И |
§ 2.3. Статистическая устойчивость невозмущенного состояния.
Равновесные распределения и их устойчивость
Понятие статистической устойчивости невозмущенного состояния, так же как и обычной устойчивости этого со стояния, относится к системам без шумов, т. е. к системам вида
it + ft (ян . . • , х п, 0 = 0 . |
(2.40) |
При рассмотрении статистической устойчивости началь ные условия считаются случайными, характеризуемыми
плотностью распределения р 0 |
(хи . |
. . |
, х п). Невозму |
щенное состояние хх = . . . . = |
х п = |
0 |
обладает полной |
статистической устойчивостью [1.13], [2.5], если при
любом начальном распределении |
р 0 (хи . . . , |
х п) теку |
|||
щее |
распределение |
плотности |
вероятности |
р {хг, . . . |
|
.. ., |
хп, t) стремится с течением времени к центральному |
||||
6-распределению: |
|
|
|
|
|
|
р (хи . . . , х п, t) |
-» |
6 (хи . . . , хп) при t -> |
оо. (2.41) |
|
Здесь |
|
оо, 6 (хи . . . , х п) = 0, |
|||
|
6 (0, . . .., 0) = |
если хотя бы одна из координат отлична от нуля:
|
оо |
^ |
^ б (t'i, . .., яп) dx1 . . . dxn = 1 • |
|
— оо |
Более слабым является требование статистической устойчивости при заданном начальном распределении р 0.
Невозмущенное состояние является статистически устойчивым при заданном начальном распределении р 0 {хи . . . , х п), если при этом распределении текущая плотность вероятности стремится к центральной 6-функции 6 (хг, . . ., х п) с течением времени. Необходимым усло
вием статистической устойчивости невозмущенного со стояния является отрицательная определенность суммы ряда
ПП
2 |
“Ь ~2~ 2 |
AiftXiXfc + . . . |
i = l |
ш i, k—l |
§ 2.3] |
СТА ТИ СТИ ЧЕСКА Я УСТОЙЧИВОСТЬ |
61 |
и стремление с течением времени хотя бы части его коэф фициентов четного порядка к — оо.
Наряду с понятием статистической устойчивости не возмущенного состояния важное значение имеет понятие
устойчивости равновесного распределения вероятности.
Как указывалось выше, стационарным равновесным рас пределением плотности вероятности называется распре деление р (zi, . . . , хп), не зависящее от времени и
удовлетворяющее ФПК-уравнению. Равновесное распре деление может быть устойчивым и неустойчивым, причем здесь можно различать устойчивость в «большом» и «ма лом», как для детерминированных систем. Равновесное распределение устойчиво (асимптотически) в большом, если при любом начальном распределении текущее распределение стремится к равновесному с течением вре мени. Равновесное распределение устойчиво (асимптоти чески) в малом, если при любом начальном распределении, достаточно мало отличающемся от равновесного (норма отличия должна быть задана), текущее распределение с течением времени стремится к равновесному.
Кроме асимптотически устойчивых и неустойчивых равновесных распределений существуют нейтральные или просто устойчивые (не асимптотически) равновесные рас пределения. Отклонения от таких распределений не нара стают, но и не уменьшаются во времени.
Вообще говоря, 6-распределение можно также рас
сматривать как возможный вид равновесного распределе ния, и тогда между статистической устойчивостью состояния и устойчивостью распределения нет принципи альных отличий. Однако равновесными мы будем в даль нейшем называть установившиеся распределения, вы ражаемые обычными, не обобщенными функциями. На хождение равновесных распределений и условий их устойчивости представляет значительный интерес для мно гих динамических систем.
Рассмотрим сначала систему без шумов, т. е. свободное движение динамической системы с аналитическими ха рактеристиками. Коэффициенты текущего распределения для такой системы согласно предыдущему определяются уравнениями (2.19). Общее решение этой бесконечной системы линейных дифференциальных уравнений пред ставлено в виде (2.23).
62 |
РЕ Ш Е Н И Е Ф П К -У РА В Н Е Н И Я |
[ГЛ . I t |
Будем считать исходную динамическую систему ста ционарной, т. е. коэффициенты aih, ai:il, . . . постоянны
ми. Весовые функции линейного приближения в этом случае являются функциями разности двух аргументов: wih {t, О = wik (t — t'). Обратим внимание на первое из
выражений (2.23):
t п п
Л0= Л® 4- ^ 2J a.ppdt' = |
-f- t 2 арр• |
(2.42) |
о Р= 1 |
Р=1 |
|
Для стационарного равновесного распределения все коэффициенты, в том числе А 0, должны быть постоянными.
Отсюда следует, что одним из условий существования ста ционарного равновесного распределения в системе без шумов является равенство нулю следа матрицы а линей
ного приближения:
П
2 аРР = 0. |
(2.43) |
p= i
След матрицы а равен взятой с обратным знаком сумме
корней характеристического уравнения линейного при ближения
| к 1 + а | = 0.
Таким образом, сумма корней этого уравнения должна быть равной нулю:
^1 + ^2 + • • • + к п = 0.
Если хотя бы один корень имеет положительную действи тельную часть, то должен быть один или несколько корней с отрицательной действительной частью. Таким образом, вследствие условия (2.43) корни, не лежащие на мнимой оси комплексной плоскости, располагаются по обе сторо ны этой оси. Отсюда следует, что в случае наличия корней с ненулевой действительной частью все или часть весо вых функций wik (0, t) = wih (— t) неограниченно (экс
поненциально) нарастают с течением времени и согласно (2.23) равновесного распределения быть не может, так как Ai (t), A tj (t), . . . не стремятся к конечным пре
делам.
Приходим к выводу, что равновесное распределение может существовать только при расположении корней
S 2.3] |
СТАТИСТИЧЕСКАЯ УСТОЙЧИВОСТЬ |
63 |
характеристического уравнения линейного |
приближения |
|
на мнимой |
оси. На мнимой оси комплексной плоскости |
могут располагаться мнимые корни и нулевые корни. Мни мые корни порождают незатухающие гармонические со ставляющие весовых функций wih. Вследствие этого,
как видно из (2.23), при наличии мнимых корней измене ние коэффициентов A t (г), A lh(t), . . . носит характер
незатухающих колебаний и стационарное равновесное распределение невозможно. Таким образом, стационарное равновесное распределение в системе без шумов может существовать лишь тогда, когда все корни характеристи
ческого |
уравнения |
линейного приближения нулевые: |
= 0 |
(i = 1, 2,. . . |
, п). В этом случае фундаменталь |
ная матрица линейного приближения есть единичная мат рица и выражения (2.23) с учетом постоянства коэффици
ентов |
аца, cintim, . . . принимают вид |
||
|
|
П |
|
А 0 = |
А-%(t) — Ai -f- 2£ 2 |
арр^ |
|
|
|
р= 1 |
|
|
t |
п |
п |
п
N |
N |
О |
п |
|
|
+ N 2 |
арц ■sAp(t') + |
|
p=i |
'lv+Г' |
|
п |
|
|
п
••Н aPij'...rfAprs (01 “Ь
N - 1
64 Р Е Ш Е Н И Е Ф П К - У Р А В Н Е Н И Я [ГЛ. И
3!
+ (/V -D ( N - г ) |
2 l a PIm |
s^pijk ( О |
"Ь • • • |
" |
|
|
|
|
N- 2 |
|
|
|
• • ' “Ь |
® Р и ' m-^P/rs (0 1 “Ь • • • |
|
|
|
N - 2 |
1 (2.44) |
2! |
laprsApH~ fit') + • |
• • |
|
TV—1 2 |
p = i |
N-l |
|
• • • + flpij^pjci- s(t')]}dt',
'"whT
Из второй группы этих уравнений видно, что для того, чтобы величины A t были постоянными, необходимо и до
статочно, чтобы
П
2 &ppi = 0, i — 1 , 2, . . . , и. p = i
С учетом постоянства И г = И? на основе третьей группы выражений (2.44) приходим к заключению, что для по стоянства A tj необходимо и достаточно выполнение ра
венства
ПП
з2 а ррц ”Ь2 ®pij^p= О-
p=i
Продолжая подобные рассуждения, убеждаемся в справедливости следующего положения. Необходимыми и достаточными условиями существования в стационарной система и-го порядка без шумов равновесного распреде ления вероятностей являются:
—наличие у характеристического уравнения линей ного приближения «-кратного нулевого корня;
—выполнение равенств
П
V—1
§ 2.3] |
|
|
СТАТИСТИЧЕСКАЯ |
УСТОЙЧИВОСТЬ |
65 |
|||||||
6 2 а ррЦ + 2 2 a PijAp — о, |
|
|
|
|
||||||||
|
Р = 1 |
|
|
Р = 1 |
|
|
|
|
|
|
|
|
1 \ |
( N |
-\- |
1) 2 |
а РРО ••• s "4" N |
|
2 |
а ргз ... sA p |
~Ь |
|
|
||
|
|
|
р =1 |
N+2 |
|
|
Р=1 |
ЛГ+1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
2 |
(a P.i/C ...s^p»+ • • • + |
a Pih1 rA p s) + |
|
|
|
||||||
|
Р = 1 |
'~ N ~ ' |
|
|
|
|
'~ N ~ ' |
|
|
|
||
|
|
2 |
П |
( а РЩ ... s-Api} + |
|
+ a Pi] - |
|
|
|
|||
+ |
"/V — '1■ 2 |
• • • |
/"4prs) |
+ |
|
|||||||
|
|
|
P=1 |
N- 1 |
|
|
|
N- 1 |
|
|
(2.45) |
|
|
|
|
3! |
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(/V — 1) (/V ZZT572 |
(a Plm |
... sApijk + • • - |
|
|
|||||||
|
|
|
— 2) |
p=l |
|
|
||||||
|
|
|
N—2 |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
... 4" flpij ... m^prs) 4” |
|
|||
|
|
|
|
|
|
|
|
|
JV-2 |
|
|
|
~b |
ДГ _^ 2 (^Pra-Api] ... f 4- •••4- GpijApIcl ... j) — 0 , |
|
||||||||||
|
|
|
P=1 |
|
|
JV—1 |
|
|
|
A '-l |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
(i, /, &, . «• — |
|
1 T2,. . ., и). |
|
|
|
|||
|
Как |
видно, |
условия |
существования |
стационарного |
равновесного распределения в системе без шумов накла дывают определенные ограничения как на коэффициенты системы, так и на коэффициенты самого стационарного распределения. Любое распределение, удовлетворяющее этим условиям, сохраняется в системе неограниченно долго. В этом смысле такое равновесное распределение можно считать устойчивым (не асимптотически) или ней тральным. Из этих общих условий существования равно весного распределения можно получить целый ряд усло вий для частных случаев.
Условия |
существования симметричного относитель |
но центра |
фазового пространства равновесного распре- |
3 А. А. Красовский
66 |
Р Е Ш Е Н И Е Ф П К - У Р А В Н Е Н И Я |
[ГЛ. II |
деления заключаются в наличии n-кратного нулевого кор ня и выполнении равенств
пп
2 а РРг ~ |
2 |
а РРгз — О? |
Р=1 |
р=1 |
|
N (N + 1) 2 |
аррц —в + |
2 |
(®pjft s-^pi + |
• • |
Р=1 |
N+2 |
Р=1 |
N |
|
|
a Plj ... rA ps) -f- |
|
||
|
|
|
|
(2.46) |
|
3! |
(a Plm ••• s-^pij/c “Ь |
|
|
|
__21 2 |
• |
||
(TV — 1) (TV — 2) |
|
|
|
|
|
P=1 |
|
N - 2 |
|
|
|
|
|
|
|
|
d p ij... mA Pfrs) |
= 0, |
|
|
|
|
N - 2 |
|
- . • — l f 2 , . * ., j i .
Для равномерного в пределах всего фазового пространства распределения вероятностей все коэффициенты A t , A i k ,
A ih[, . . . равны |
нулю. |
В соответствии |
с (2.45) Необходимыми и достаточными |
условиями существования в стационарной системе без шу мов равномерного равновесного распределения вероят ностей являются наличие n-кратного нулевого корня у
характеристического уравнения |
линейного |
приближения |
|
и выполнение |
равенств |
|
|
п |
п |
п |
|
2 a PPi — 0» |
2 а РРг) = 0| • • .) |
2 a PPij---s |
— 0» (2.47) |
Р =1 |
Р =1 |
Р=1 |
|
i, j, s... = 1, 2, . . . , п.
Динамические системы, в которых выполняются условия
П
(2.47) и условие 2 a w — 0, названы выше обобщенно
p = i
консервативными системами. Подклассом таких систем
§ 2.3] |
СТАТИСТИЧЕСКАЯ УСТОЙЧИВОСТЬ |
67 |
являются системы без прямых связей, у которых все коэффициенты с двумя одинаковыми первыми индексами равны нулю (арр = 0, appi = 0, . . . ).
Таким образом, в стационарных обобщенно консерватив ных системах и системах без прямых связей при наличии n-кратного нулевого корня у характеристического уравне ния линейного приближения в отсутствие шумов сущест вует равномерное равновесное распределение вероятно стей. Если в таких системах в начальный момент было равномерное распределение, оно сохраняется неограни ченно долго.
Согласно тем же соотношениям (2.45) необходимыми и достаточными условиями существования в стационарной системе без шумов равновесного нормального распределе ния вероятностей являются наличие n-кратного нулевого корня и выполнение равенств
п п п
2 |
^ppi= |
о» |
® 2 appi3 “i~ 2 2 |
Я'р'и^р = |
о» |
|||
p=i |
|
|
|
p=i |
р=1 |
|
|
|
N |
(N + |
1 ) |
2 |
а ррЦ ...» + |
^ |
2 |
а рИ — » ^ р |
+ |
|
|
п |
35=1 |
P=1 |
'ТК? |
(2.48) |
||
|
|
|
|
|
|
|
|
|
|
+ |
2 |
(а РЙ ... s-4pi + |
• • |
• + |
a Pij ... r-Aps) — О, |
||
|
р=i |
N |
|
|
N |
|
i, j, к,. .. = 1 , 2 , ... , n.
Рассмотрим конкретный пример — вращательное дви жение твердого тела, описываемое уравнениями Эйлера
(1.37)
^ 1 “f~ |
^ 1 2 3 ^ 2 ^ 3 |
~ |
2 ~~Ь |
^ 2 1 3 * ^ 1 ^ 3 = |
* 3 |
|
~~Ь ^ 3 1 2 ^ 1 ^ 2 |
= О , |
|
|
|
|
|
|
|
|
(249) |
где |
хх, х2, |
х 3 — угловые скорости в |
связанных |
осях, |
||||
а 123 = |
т---- » |
Я213 = |
--- J ---“ > |
а 312 |
= |
---- • |
(2.50) |
|
|
J x |
|
J V ‘ |
|
|
J z |
|
3*
68 РЕ Ш Е Н И Е Ф П К -У РА В Н Е Н И Я [ГЛ . 11
Характеристическое уравнение линейного приближения
имеет |
здесь |
трехкратный |
нулевой корень. Условие |
||
з |
|
выполняется. |
Кроме |
того, |
ацЫ = ацыт = |
aPPi = 0 |
|||||
р — Х |
— 0. |
Поэтому равенства |
(2.45) |
принимают вид |
|
= |
з
2 ариАр = °*
р=1
(2.51)
2 {a PrsA pij — f 'I" ■• • "Ь a pijA pkI ... s) — 0 ,
Видно, что равномерное во всем фазовом пространстве рас пределение плотности вероятности является в данном случае равновесным. Кроме того, существует множество других равновесных распределений. Раскроем первые две
группы (N = 2, 3) соотношений |
(2.51) с |
учетом того, |
|
что все коэффициенты |
(г, /, к |
= 1, 2,3), |
кроме (2.50) |
и им симметричных в смысле перестановки двух послед них индексов, равны нулю. Получаем
а123А 1 |
== 0, |
|
|
= |
0, |
Яз12^3 — 0, |
||
2йз12^4з1 |
= 0) |
2a2i3-^2i |
= |
0, |
2a3x24 2 — 0, |
|||
Яхгз-^ц |
& 213А 22 + |
Язхз^ЗЗ |
= |
0, |
2я213-423 — |
0, |
||
|
ai23A i2 "Ь 2а |
1234 з |
= |
0. |
|
|
||
■Если а123 ф |
0, я213 Ф 0, |
а312 Ф 0, |
что имеет |
место при |
различных моментах инерции тела относительно главных
осей, то из этих |
выражений вытекает |
|
|
|
||
А \ = А ° = |
4 = 0 , |
А°п = А \ з = 4 з = 0 , |
(2.52) |
|||
J _ |
j„ — j . |
I |
V |
х |
л» |
|
|
|
J v - |
J x |
л 33 |
0. |
|
4 + ■ |
' "Г |
г |
|
Таким образом, «нормальная составляющая» равновесного распределения угловых скоростей является центральной канонической. В главе I на основе привлечения первых
§ 2.3] |
СТАТИСТИЧЕСКАЯ УСТОЙЧИВОСТЬ |
69 |
интегралов исходной системы уравнений было показано что в данной системе существует множество равновес ных распределений, выражаемое формулой
1пр = |
Т ( / Л а + Jyxt + |
J zx l Л х \ + J\x\ + |
Jlxl), |
(2.53) |
||||
где ¥ — произвольная функция. |
Ясно, |
что |
нормальная |
|||||
составляющая распределения (2.53) имеет вид |
|
|
||||||
С1 ( J x X l Ч' J y X 2 Ч~ J 2 Х л‘ ) |
Ч‘ с 2 (? хх 1 Ч~ J'ilx 2 Ч" |
з) = |
|
|
||||
|
|
_ |
_L А0 г2 |
I- _L /1° -А А |
1 |
Л° г2 |
||
|
|
— |
|
' ~ |
л ггхг i — ~ |
-^зз-чи |
||
где |
2 J x ( Сх + |
C2J X), А°22 = |
2 J у (С, Ч- C2J y), |
|||||
А°п = |
||||||||
|
А% = 2/ г (Сх + С2/ г), |
|
|
|
|
|||
a Cj, С2 — произвольные |
постоянные. Подставляяэти |
|||||||
выражения для Ац, |
Л22, |
в (2.52), получаем |
тождество |
2Сх { J z ~ J y A ~ J x — J z ^ -J y —J х) Ч- 2С2 (J xJ z — J xJ у Ч*
Ч- J XJ у JyjZ 4“ J yjz 4 a;/2) ==
Аналогичную проверку удовлетворения условиям (2.51) можно выполнить для старших членов распределения (2.53), представленного в форме степенного ряда.
Перейдем к рассмотрению равновесных распределений в стационарных системах с шумами. Обозначим коэффи циенты стационарного равновесного распределения, кото рые по определению постоянны, чертой сверху. В соот ветствии с (2 .1 1 ) эти коэффициенты удовлетворяют бес
конечной системе нелинейных алгебраических уравнений
|
2 |
З р ч ( A p q + & p A q ) — — 2 °РР> |
|
|
Р, д=1 |
р=1 |
(2.54) |
|
|
|
|
2 |
р |
Ч 2 ^ Р 1 pqi Ч- A ptA q) — — 2 |
2 a ppii |
Р=1 |
|
р, q— l |
Р=1 |