Файл: Красовский, А. А. Фазовое пространство и статистическая теория динамических систем.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)а:а + |
Л(4)г4 + - 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! |
п |
|
|
|
|
|
|
|
|
• |
■• |
+ |
|
^ |
1®ррх (О х!рчц...*|(Г) + |
• • • |
|
||||||
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) |
нормальной составляющей логарифмической |
||
плотности |
вероятности. |
Можно |
констатировать, что |