Файл: Красовский, А. А. Фазовое пространство и статистическая теория динамических систем.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) |