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

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

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

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

Добавлен: 23.10.2024

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

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

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

50

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

[ГЛ . II

где к — const, и известно, что начальное распределение х (при t = 0) является равномерным в интервале от — л

до л:

p„W = P (x,0) = (

i

И <

“ ’

(

0

при | х | >

л.

ФПК-уравнение для логарифмической плотности вероят­ ности в данном случае имеет вид

д In р

, .

д In р

1

».

к sin х

Q

= к cos х.

at

 

ox

 

Легко проверить, что'это уравнение при указанном началь­ ном условии имеет точное решение вида

Р (xi т)

2я (chт — cos a; sh т)

где т = kt — «безразмерное

время».

 

 

Найдем теперь приближенное решение описанным

методом. Учитывая, что

 

00

 

СО

 

 

 

 

к sin х = к 2

T ^ q r W * * v+1=

2 e(w i)*,wl.

v=0

'

Г

 

v=0

 

 

 

 

/ Л\м

k

1)!

®(2v+x) — аи_... i — (— 1 )

(2v +

и что в данном случае

2V +1

 

 

 

 

 

 

 

А°0 = \ п ± ,

 

=

= 0 (N = 1 ,2 ,.. .) ,

, t) = wn (t' t) = exp (t' t),

по формулам (2.23) последовательно находим

= 1п - ^ + т,

^(2) (f) =

(!

— ехр 2т),

At) (T) =

~

~T exp 2т + T " exp 4 t’

4e) (T) =

U

(TF “ 30 exp 2t + eip

4t

exp 6 t) ‘

На рис. 2.1

для нескольких моментов

времени представ­

лены точное

распределение р (х, т)

и

приближенные


S 2.21

 

Ч А С Т Н Ы Е

СЛУ ЧА И

51

распределения

 

 

 

 

Р(2) (*. Т) =

^

еХР (Т +

А ^

2 ) ’

 

 

Р(4) (*. т) =

2^Г ехР (Т + 4 ~ АЫх2 +

т - 'W * ) ’

 

Pw (х, t) =

- ^

exp (t +

- |- Л(2)а:а +

Л(44 + - L

H(e)z e) .

Видно, что для небольших значений т приближенные рас­ пределения весьма близки к точному, особенно при N = = 4, 6. Легко проверить, что полученные приближенные

Ц20

0,12

Ofl't

0,4

0,2

Рис. 2.1.

распределения совпадают с точным в отношении значении производных до соответствующего порядка (N ) в начале

координат.

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

Итак, данный метод применительно к системам без шу- $ов дает при любом заданном N решение, совпадающее

52

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

[ГЛ. П

с точным в отношении всех производных до N -го порядка

включительно, вычисленных в начале координат. Продолжая рассмотрение частных случаев, заметим, что

значительное упрощение вычислений получается в том случае, когда фундаментальная матрица линейного при­ ближения w — I wil{ I является диагональной, т. е. когда

«каналы» объекта в линейном приближении автономны. В этом случае выражения (2.23) принимают вид

( п

Л, (0 =

Aiwa (0, t) + 2 ^ S

aPPi wu (*'. О dt’,

 

 

 

 

 

 

 

ОР =1

t

п

 

 

 

 

 

 

 

 

 

 

 

A v (t)

=

 

 

(0, t) w jj (0, t) + 2 § 2

а ррц (О +

 

 

 

 

 

 

 

 

 

0

p—1

 

 

 

 

 

n

 

 

 

1

 

 

 

 

 

+

2

a Pii (* ') A P (**)j w a (* '• * )

WJ) (*'> 0 d t \

 

 

 

 

P—1

 

 

 

 

 

 

ЛуА-...s (t) =

 

 

(0, t) Wjj (0, t) . . . wss (0, t) -|-

'~Ax~'

 

 

t

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

■-]- ^ Ш (N

1)

^ a ppv---s

( t ) +

 

 

 

 

о

 

 

P=1

'~N+2

 

 

 

 

 

+ N ^ «

^ ( 0

^

( 0 +

 

(2.24)

 

 

 

 

 

p=i

nTi

 

 

 

 

 

f

2

[a p j^ s { П A pi (O +

• • • +

(<') Лр.(01 +

 

p=i

 

"ТТ’

 

 

 

 

n

 

+

J L

-

2

[Opfc..., ( 0 ^

/

( 0

+

 

 

'

ЛГ-

1

Л

 

A--T

 

 

 

 

 

 

• . . +

Яpi.../ (О Лргз (0 1

~b • •

 

 

 

 

N^1

 

 

 

 

 

 

 

 

. . . +

 

Ы

ъ

[аргЛПАг1::.1(П + ...

 

 

 

N ~ 1 р=1

 

 

W-1

 

 

. . .

+

вру (О Л р ^ .8 ( < ' ) ] ( * ' .

О • • •

( *

' , 0)

N- 1


§ 2.2]

Ч А С Т Н Ы Е СЛУЧА И

53

Вычисления и общий анализ с помощью формул (2.24) настолько упрощаются в сопоставлении с (2.23), что часто выгодно путем преобразования переменных или посред­ ством введения дополнительных фазовых координат (рас­ ширение фазового пространства) приводить объект к виду, при котором wij = 0 при / ф I, т. е. искусственным

путем обеспечивать автономность уравнений линейного приближения.

Если динамическая система линейна (аш = сцкы — = . . . = 0), то выражения для коэффициентов логариф­

мической плотности вероятности свободного движения становятся весьма простыми:

 

П

 

 

•4y/£...s if)

2

- Х ^Vi (0 , t) W\ij (0 , t) ...

w Xs (0 , t).

N

v.l-i.....x= l

_ N~~

 

(2.25)

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

П

 

П

lnpo = Ло h 2

X i Л-- ---

2 A iK x ix k +

i= l

i, k=l

П

+ -g“ 2 A°m XiXkXl + . . .

i, kt

координаты Xi заменить на

n

2 wik(0,t)xk.

Осуществляя такую замену для коэффициентов текущего распределения, получаем выражение (2.25).

Как уже отмечалось, довольно широким и интересным классом систем являются системы без прямых нелинейных связей. Для таких систем


54

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

[Г Л . II

и выражения (2.23) принимают вид

П

Ai (0 = 2 A*u>*i(0,t), v—1

П

 

 

(О= 2

 

(О, о ww (о, о +

 

 

 

 

 

 

v,H-=x

 

 

 

 

 

 

 

 

 

 

n

f

п

 

 

 

 

 

 

 

+

 

2 2

$ [ 2

(О ЛР (*') «\i (*',1) и>м(Г, <)] dt’,

 

 

v, р = 1 О р = 1

 

 

 

 

 

 

 

4 , ^ ( 0

=

 

2

^y-X^vi (О, t) W[4 (О, Q. ■■

 

 

IV

 

V, Ц... Х=1

IV

 

 

 

 

 

 

 

 

 

 

п

(

 

п

 

 

 

• • • i^xs (О» О +

^

$ {-N 2

°ру-^х(<') Ap(t') +

 

 

Л

 

 

v, р,....х=ю

 

p = i

IV4-1

• (2.26)

+

taP!i—х (О Ap4{t') + . . . -jj- йрч...е(У)Арх{У)1+

2

 

Р=1

 

^ V

 

 

 

 

 

IV

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

+

 

 

 

2

1®Ре—X (* ) ^pviJ. (<') +

• • ■

 

 

 

 

 

 

р = 1

IV—1

 

 

 

 

 

 

• •

+

дру...-л (О -^РРХ (01 + •

• •

 

 

 

 

 

 

 

N^1

 

 

 

 

 

 

 

 

 

 

 

 

2!

п

 

 

 

 

 

 

 

■•

+

 

^

ррх (О х!рчц...*|(Г) +

• • •

 

N — \

 

 

 

 

 

 

p=i

 

 

 

 

 

 

• • • +

flpv[i (О Aptp....x(t')]

X

 

 

 

 

 

 

 

 

 

X

u>v'i ( t \ t )

W\4

(Г ,

* ) . . .

w Xs (t',

t)} d t ’,

Выражения (2.26) показывают, что для данного клас­ са систем справедливо указанное выше положение, что неустойчивость невозмущенного состояния объекта спо­ собствует сходимости ряда (2.8) и вызывает уменьшение Aijn ... s (t) с течением времени. Действительно, для

системы, неустойчивой в линейном приближении (по Ля­ пунову), функции wtj (t , t'), вообще говоря, —нара­ стающие функции времени t, а функции w^ (t', t) — убы-


§ 2.2] Ч А С Т Н Ы Е СЛУ ЧА И 55

лающие функции времени. При этом нарастание и убы­ вание может быть не монотонным, однако все wtj (0, t) —►

—►0

при

t —*■ оо.

Отсюда

видно,

что первые

слагаемые

в правых

частях

формул

(2.26)

стремятся к

нулю при

 

о о и притом тем быстрее,

чем выше N. Что касается

последующих

членов в правых

частях

формул

(2.26),

то они могут быть представлены в виде

 

 

 

 

2

 

 

 

 

 

.. .wxs(t',t)dt'.

(2.27)

v,y....х=10

 

 

 

 

 

 

 

 

Из

условия

неустойчивости

следует,

что

u>„i (tr, t)x

Xw^fo' ,t). . .wxa (t',t) — нарастающая (не обязательно мо­ нотонно) функция t'. Она обращается в единицу при V = t и становится сколь угодно малой при t —*■ оо на другом конце интервала интегрирования в (2.27), т. е. при t' = 0. Отсюда следует, что если F^.... х (t') есть убывающая функция £', стремящаяся к нулю при t' -»- оо, то величи­ на (2.27) — также убывающая функция t, стремящаяся к нулю при t -*■ оо. На основе этого и выражений (2.26)

при постоянных или убывающих во времени коэффициен­ тах ant... s последовательно получаем, что A t (t ) , A i} (t),...

. . . , Ань ... s (t) суть убывающие функции времени,

'N

стремящиеся к нулю при t -> оо.

Обратим внимание на следующее свойство свободного движения систем без прямых нелинейных связей. Рассмот­ рим два первых выражения (2.26). Если начальное рас­ пределение центральное симметричное, т. е.симметричное

относительно начала

координат, то А° = 0, A t (t) = 0 и

 

П

 

Ан (t) =

2

(0, t) Wpj (0, t),

V,IA=1

т. e. то же самое, что в линейной системе. Назовем сумму членов

ПП

 

Л о + 2 A ix i

Ч----2~ 2

A ikx ix k

 

i = l

г, k=1

 

ряда (2 .8)

нормальной составляющей логарифмической

плотности

вероятности.

Можно

констатировать, что