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

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

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

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

Добавлен: 23.10.2024

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

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

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

10 П РЕД И С Л О В И Е

В главе IV рассматриваются равновесные распределе­ ния в системах с кусочно-линейными характеристиками, а также ставятся некоторые задачи синтеза.

Показано, что решение ФПК-уравнения для объекта без шумов подобно решению задачи аналитического конст­ руирования оптимальных управлений при определенном виде минимизируемого функционала.

Заключительная глава работы имеет самостоятельное значение и лишь косвенно связана с предыдущим. В ней рассматриваются вопросы предельной точности управле­ ния, обусловленной тепловыми шумами. Здесь обобщены известные результаты для линейных пассивных систем и получены оценки предельной точности стабилизации для систем с обратными связями. Вводится понятие микроуп­ равления для задач управления объектами микроскопи­ ческих (десятки ангстрем — микроны) размеров. Найде­ ны некоторые оценки предельной точности управления та­ кими объектами. Есть основания предполагать, что эти вопросы имеют принципиальное значение и получат широ­ кое развитие в будущем.

Автор выражает благодарность профессору Ю. А. Ко­ четкову за замечания, сделанные при рецензировании ру­ кописи книги.

Г Л А В А I

УРАВНЕНИЕ ФОККЕРА — ПЛАНКА — КОЛМОГОРОВА

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

По-видимому, первой непрерывной динамической систе­ мой, описываемой обыкновенными дифференциальными уравнениями со случайными функциями, рассмотренной с позиций статистической теории, явилась броуновская частица. Физическая теория броуновского движения была создана Эйнштейном и Смолуховским в 1905—1906 гг. [1.1], [1.2]. Еще раньше Башелье получил закон, которому подчиняется положение частицы, совершающей одномер­ ное броуновское движение [1.3]. ФПК-уравнение (или, как его еще называют, диффузионное уравнение, уравнение диффузии вероятностей) на физической основе было выве­ дено в работах Фоккера [1.4], Планка [1.5]. Строгий вывод принадлежит А. Н. Колмогорову [1.7]. Диффузион­ ные процессы в фазовом пространстве, к которым относит­ ся броуновское движение, можно рассматривать как марковские процессы с непрерывными траекториями.

Широкая физическая постановка задачи о статистиче­ ском рассмотрении динамических систем, связанная с ФПК-уравнением, была опубликована А. А. Андроновым, Л. С. Понтрягиным, А. А. Виттом в 1933 г. [1.6]. Матема­ тической теории диффузионных процессов *) посвящены труды крупнейших математиков: А. Н. Колмогорова [1.7], [1.8], П. Леви [1.9] и др. [1.10], [1.11]. Этот раздел совре­ менной математики быстро развивается и содержит глу­ бокие результаты.

В данной главе излагается вывод ФПК-уравнения на основе физических представлений.

*) Как уже отмечалось в предисловии, аспекты строгой мате­ матической теории диффузионных процессов совсем не затраги­ ваются в данной книге.


12

У РА В Н Е Н И Е Ф О К К Е РА — П Л А Н К А — КОЛМ ОГОРОВА [ГЛ . X

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

§ 1.1. Вывод ФПК-уравнения

Мы будем рассматривать динамические системы (объ­ екты), описываемые обыкновенными дифференциальными уравнениями вида

“Ь fi (■X'li

-Tm t) — ii> i = 1 , 2, ..., n.

(1 -1 )

Здесь Xi, ..., xn — координаты, которые называются gбазо­ выми; ..., tn — случайные^ функции времени типа бе­

лых гауссовых шумов с матрицей] спектральных плотно­ стей S = l.S'jfc!, Sik — const, и корреляционной матрицей Д (т) = ||Bik (г) I = ||5{к6(т)||; 6 (т) — дельта-функция.

Относительно функций / ;, которые часто будут имено­ ваться характеристиками исходной системы или объекта,

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

Уравнения в форме (1.1) и их частный вид для свобод­ ного движения (система без шумов)

[Xi + ft (xlt ..., хп, t) = О

(1.2)

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

Начальные значения фазовых координат

(0), ...

..., хп (0) считаются случайными, и их распределение харак­

теризуется плотностью вероятности

р {х1, ..., хя, 0) = р 0 (Xi, ..., Хп),

(1.3)

которая по определению удовлетворяет условию норми­ ровки

i . . . d x n = i.

(1.4)

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


S 1.1

В Ы В О Д Ф П К -У Р А В Н Е Н И Я

13

ние текущей плотности вероятности р — р (хи

..., xn, t)

в фазовом пространстве. Простой и достаточно нагляд­ ный, на «физическом уровне строгости», вывод уравнения получается при использовании представления о «фазовом газе». Подобные построения часто применяются в статис­ тической физике. Будем рассматривать уравнения (1.1) или (1 .2) как уравнения движения частиц некоторого во­

ображаемого фазового газа. Считаем, что частицы фазового газа не могут исчезать или возникать, их общее число во всем фазовом пространстве постоянно.

Изменение числа частиц в любом элементе объема dxy...

... dxn фазового пространства равно разности числа частиц,

поступающих и убывающих из этого элементарного объ­ ема. Понятие фазового газа может быть интерпретирова­ но также следующим образом. Допустим, что имеется очень много (или бесконечно много) совершенно одинако­ вых (с одинаковыми / г, | 4) систем (1 .1 ), отличающихся

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

Выделим в фазовом евклидовом пространстве элемен­ тарный объем в виде n-мерного прямоугольного параллеле­ пипеда с центром в точке (х1, х2, ..., хп) и сторонами dxl , dx2, ..., dxn, параллельными координатным осям.

Движение частиц фазового газа в каждой точке фазово­ го пространства имеет регулярную составляющую, обус­

ловленную функциями fi( x $ .... хп,

t), и случайную со­

ставляющую, вызываемую шумами

(t). Поэтому мгно­

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

так и случайную центрированную (с нулевым математи­ ческим ожиданием) составляющую Др (а^, ..., хя, t).

Математическое ожидание плотности фазового газа рав­ но искомой плотности вероятности в фазовом прост­ ранстве.

Скорость изменения количества частиц в выделенном элементарном параллелепипеде равна сумме потоков час­


14

У РА В Н Е Н И Е Ф О К К Е РА — П Л А Н К А — КОЛМ ОГОРОВА [ГЛ . I

тиц, проникающих через его грани. Полагая стороны па­ раллелепипеда сколь угодно малыми, с точностью до ма­ лых высших порядков записываем это условие неразрыв­ ности, или сохранения:

щ

(* 1 .

• • ; Хп , t) + Ь р (ж х, . .

. ,

х п , t)] d x v

. . d x „

=

 

-----|Тр

*

* *’ ^ i—l*

2

 

^i+1* ‘

 

. . . ,

x n , t) +

 

i =

x j ^

 

 

 

 

 

 

.

 

 

"h

^ 1 )

* • *» ^г—X»

 

"2

^‘^'ii *^i+l»

*

* •)

 

X

£

f*^l>

* •

•» ^ i - 1 ’ *^i

2

 

^ i+ 1’ * *

•» ^тп % \

S i ( O j

X d a ^ . . . da:b l da:i+ 1 . . . d x n

|^P

^ 1

» * * ч

1? Щ H-

*2 "

^Ч+Х» * * * т *^ТИ ^ H“

"T

^ i »

• * •»

 

y d x it ^ i+ i» • * • .

* ] ] X

X £ f%

 

^Ч-Х»

" f 2

^*4>

^i+X) • • *»

^ “Ь S i ( 0 ^

X dXi-.. dXi-xdxi+i- .. dxn| =

n

=% £ - H P + b p ) V i - b ) ] d 1=1 1

Отсюда следует

^

+

Ж

 

= 2 5Г [(p +

- 6i)l =

 

 

•• У■

i* l

г

■ ■

 

 

 

 

 

 

 

 

П

 

 

 

 

 

 

=

2

[ajr (Pf0 ~

 

gj-

+ 9^7 (&Pfд —

(ph)

 

 

i= l

*r.

г

 

, г

t

г

Математические ожидания величин

д&р

dt »

равны нулю (Др, — центрированные слу­

чайные' функции* /|, р — неслучайные функции). В


1 1.11

В Ы В О Д

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

15

соответствии

с этим из

(1.5) вытекает

 

1=1 1 '

^

= S { - йЧЬрЫ + щ (ДМ) - si**) +

[ДМ?}.

 

i = l '

г

t

г

 

г

J

где

М [

) — символ

математического

ожидания.

(1'6)

 

 

Второе из этих соотношений с учетом того, что||1

= О»

можно переписать в виде

 

 

 

 

д&р

 

 

д&р

д&р

 

 

dt

+ Е 9# ( - / , + 1,) = ^ + 2 ^

. - д #

 

 

1=1

дх,

i = l

дх.

 

 

 

 

*

4

 

 

 

 

 

“ Ap E ^ + S k - m 'ia ^ i ^ S i . ^

 

 

 

1=1 г 1=1

i

 

»=1

дх.

ИЛИ

 

 

 

1

 

 

 

 

 

 

 

 

Д р - Д р Ъ щ = Е з|-М [Д М ]

1=1

дх.

(1,7)

 

 

1=1

1=1

 

 

 

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

нений исходной системы (1.1). В каждой заданной точке

2dfi

фазового пространства величина 1=1 dxi является функцией

времени, и дифференциальное уравнение (1.7) в этой’точке можно рассматривать как линейное относительно Ар.

В начальный момент, времени t = 0 плотность фазового

газа считаем равной заданной начальной плотности рас­ пределения и Apt=0 = 0. Решение линейного уравнения

(1.7) при этом условии имеет вид

t

П

П

\ f Ар = ^ д р ( ; , п { 2 Д м [ А р У - s Ь

(i.8)