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

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

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

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

Добавлен: 23.10.2024

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

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

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

74

 

Р Е Ш Е Н И Е Ф П К -У Р А В Н Е Н И Я

 

[ГЛ. И

имеет действительное решение

 

 

 

Аи =

2ап

А22

2(Z22

•^33

2азз

 

iSn

*^22

~ЛГ’

(2.60)

 

 

 

А12 — Л13 А23 =

0.

 

 

Это решение является единственным и совпадает с реше­ нием для линейной системы. Это вытекает из указанного в § 2.2 общего свойства систем без прямых нелинейных свя­

зей (см. выражения (2.30) — (2.34)), к которым относится данная система. Величины (2.60) обращают первое из уравнений (2.59) в тождество.

Нижняя группа уравнений (2.59) с учетом (2.60) при­ нимает вид

aUk An +

aihiAhk +

a,ik А ц = 0,

„ „..

(*, М

= 1, 2 . ..............

П).

 

Ненулевыми являются лишь коэффициенты aiU с раз­

личными числовыми индексами. Поэтому соотношения (2.61) сводятся к одному (остальные или выполняются тож­ дественно, или совпадают с данным соотношением):

а128-4ц “Ь a312^33 4" а213^-22 =

или в раскрытой форме:

J z

J v ап

I J x

J z

022

J„ — /

X a33

= 0.

V

J .

S n

^

J„

S 22

 

^зз

 

Спектральные плотности S n , S 22, S 33 связаны со спект­ ральными плотностями S x, S v, Sz возмущающих моментов

относительно соответствующих осей соотношениями

S 11

S

J2

 

 

JV

Коэффициенты an , а22, а33 пропорциональны коэффициен­

там моментов демпфирования:

 

йол

-- ■

__

Wc

йи —

 

(DZ

азз - “ JГ

 

 

 

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

условие

существования

нормального

равновесного распределения в рассматриваемой системе


§ 2.3]

СТАТИСТИЧЕСКАЯ УСТОЙЧИВОСТЬ

75

сводится к соотношению

 

 

(■h - J y )•

(Jx - Jz)

V v - J x) - ^

= 0.

Это условие выполняется, в частности, тогда, когда спект­ ральные плотности возмущающих моментов пропорцио­ нальны (коэффициенты пропорциональности одинаковы) соответствующим коэффициентам моментов демпфирова­

ния:

т

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

Для того чтобы равновесное распределение реализо­ валось в действительности, оно должно быть устойчивым или нейтральным (устойчивым не асимптотически). Ис­ следование устойчивости в малом равновесных распре­ делений производится на основе линейных уравнений, получающихся из (2.1 1 ) линеаризацией в окрестности рав­

новесных значений коэффициентов: ft

А Й 0

-------2~ S ‘^ 'р д ^ А р д - Ь Н р Д И д - ^ Н д Д Л р ) = 0 ,

 

Р. 9=1

П

 

 

П

 

А Л

2 ® P iA H p

2 Sp<l (A -^ p g i 4 “

 

Р = 1

P ,9 = l _

_

 

 

+

4~ A q S .A p ^ ) = 0 ,

п

( 2.62)

А Й ; Л — 2 ( а р А А р н 4 “ ® р л А Л р{) —

p = i

пп

2 2 а Р г к ^ А р — S S p q f i A A p q i k 4 " 2 Н p i i t ^ A q 4 "

Р=1

Р, 9=1

 

4- 2ДвДЛр;й -f- K p i & A g]s 4- ^afeAHpi) = О,


76

Р Е Ш Е Н И Е Ф П К - У Р А В Н Е Н И Я

1ГЛ. 11

Нахождение общих условий устойчивости тривиаль­ ного решения этой бесконечной системы линейных урав­ нений затруднено. Обычно приходится ограничиваться укороченной системой, в которой, начиная с некоторого значения N,

Mftl...s —

— •••—1

(2.63)

N

Tv+1

 

Условия устойчивости, найденные таким путем, носят характер приближенных.

Если рассматривать устойчивость нормального цент­ рального равновесного распределения и в (2.63) считать N = 3, то уравнения (2.62) принимают вид

П

*1

АЛ0 -----2~ 2л S p, M p4 ~

р.<?=1

пп

A^i

2

(api

2

&pqA ?j) АА р — 0,

 

р= 1

 

 

4

= 1

 

 

 

п

 

 

 

п

 

 

АЛ it

2

( a Pi

 

2

SpqA qi \ AA p!i

 

г>=i

'

 

4=i

'

(2.64)

 

 

n

 

 

 

n

 

 

 

2

[aPk Ч"

2

qh J АЛр;

 

 

p = l

'

 

 

4=1

'

 

 

 

 

 

 

 

n

2 2 ®piftA^4p — 0) p=i

i, к — 1 , 2, . . . , n.

Устойчивость тривиального решения этой системы урав­ нений полностью определяется устойчивостью тривиаль­ ного решения второй группы уравнений.

Действительно, корни характеристического уравнения третьей группы (А р полагаем при записи этого характери­

стического уравнения равными нулю) равны попарным суммам Xt + Xj корней , . . ., Хп характеристического

уравнения второй группы. Это можно показать различны­ ми способами, в частности, на основе того, что при соот­ ветствующих ограничениях на начальные условия реше­ ние уравнений третьей группы (при ДЛР = 0) равно


§ 2.41

У Р А В Н Е Н И Я МОМЕНТОВ

77

произведению решений второй группы уравнений:

A A lh =

AA iA A h.

Таким образом, условия асимптотической устойчивости

в отношении коэффициентов

A t, A ih (коэффициент А 0

в этом случае обладает не асимптотической устойчивостью) сводятся к размещению корней уравнения

 

 

| X 1 -

ат -

SA | =

0,

(2.65)

где 1 — единичная матрица,

в левой

комплексной полу­

плоскости.

 

 

 

 

 

 

 

 

 

П

 

 

 

Согласно

(2.57)

при

2 аррл = О

 

 

 

 

 

p = i

 

 

 

 

ат А + 1 а + Л SA = 0.

(2.66)

Если агА =

А а,

то из

(2.66) следует SA =

—а — ат

и согласно (2.65) условие устойчивости равновесного нор­ мального распределения сводится просто к условию ус­ тойчивости нулевого решения уравнений линейного при­ ближения.

§ 2.4.

Уравнения моментов

По

определению

моменты N -го порядка равны

 

 

оо

 

 

Мщ..., = ^

^ р (эсъ . . .

,хп, t) х{хк . . . x s dx! .. . dxn =

N

 

—со

 

N

= ^ ■■. ^ехр 1^А0 +

2

А рХР + 4 ~ 2 А****** +• • • ) ><

 

—00

 

р = 1

р, 9=1

X xixk . .. xs йхг . .. dxn. (2.67)

Они могут быть вычислены после определения коэффициен­ тов распределения А р, А Рд, . . . Для нормального рас­

пределения функциональная связь первых и вторых мо­ ментов с коэффициентами А р, А рд выражается в об­

щем виде.

Для распределений произвольного характера (пред­ ставленных степенным рядом) такая связь в общем явном виде не выражена и, вероятно, не может быть выражена.


78

Р Е Ш Е Н И Е Ф П К - У Р А В Н Е Н И Я

[ГЛ. И

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

Уравнения моментов могут быть получены различ­ ными способами. Так, для свободного движения системы

ПП

+ 2 a vsx u +

2 a m x kx i + • • • = 0

( 2 .6 8 )

к=1

к, (=1

 

бесконечная линейная система уравнений моментов полу­ чается простой подстановкой выражений для производ­

ных Xi из (2.68) в тождество

(х гх к ' ' * x s) ~ x ix k •■ •x s “Ь ••* x ix k •••Х8

и применением операции вычисления математического ожидания [1.13], [2.8]. Рассматривая систему с шумами

П

П

Zi + h = xi + 2 aikxk +

2 aikixkxi + •••=£*> (2-69)

k = i

к, 1=1

получим уравнения моментов на основе ФПК-уравнения и понятия характеристической функции.

Как известно [2.6], [2.7], характеристической функцией

называется преобразование Фурье плотности вероятности

ОО

 

g (соь ... ,©„,*) = $•••$/> (хъ . . . , x n,t) X

 

— ОО

 

X ехр I/ (сй!*! + • • • + »„£„)] dx1 . .. dxn. (2.70)

Здесь j = У — 1 — мнимая единица.

(Во избежание

путаницы в данном параграфе буква /

в индексах не ис­

пользуется.) Частные производные характеристической

функции в точке

tOj =

. . . = со„ =

0 пропорциональны

моментам

 

 

 

 

 

 

(s5? C b s

; L

. _ . . = r

^

'

Р '71)

^ •

-

 

” .