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

84 Р Е Ш Е Н И Е Ф П К - У Р А В Н Е Н И Я [ГЛ. II

При условии (2.81), которое, конечно, является только достаточным, но не необходимым, ряд (2.82) равномерно сходится в области G. Равномерно сходятся и ряды, полу­ чаю щ и еся дифференцированием dp/dt, Lp.

Подставляя (2.82) в уравнение (2.78), убеждаемся, что оно удовлетворяется. Выполняется также начальное условие (2.80). Итак, (2.82) действительно является иско­ мым решением.

В случае отсутствия шумов (S ih = 0) в точке равно­ весия (/ = 0) оператор (2.80) вырождается в след яко­

биана

При непрерывных функциях / г в достаточно малой окрест­

ности равновесной точки

П

В данном случае можно принять d = 0, а величину D задать равной наибольшему в области G значению моду­

ля величины (2.84):

D — шах

Vi

дх.

G г=1 г

Из (2.85), (2.82) или (2.84) вытекает приближенное вы­ ражение текущей плотности вероятности в малой окрест­ ности равновесной точки системы без шумов:

71

i=l 1

чии шумов

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

(2.89)

где L — произвольный линейный стационарный оператор по х,

— заданные начальные условия. Если ряд

равномерно сходится в области G с его производными,

входящими в (2.89), то он сходится к решению уравнения (2.89), удовлетворяющему начальным условиям (2.90).

Г Л А В А 111

РЕШЕНИЕ ФПК-УРАВНЕНИЯ МЕТОДОМ РЯДОВ

И СТАТИСТИЧЕСКАЯ ДИНАМИКА КОНКРЕТНЫХ КЛАССОВ СИСТЕМ

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

§3, 1. Примеры статистического исследования конкретных классов систем

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

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

Решение, удовлетворяющее произвольному началь­ ному распределению, приходится строить на основе мето­ да рядов.

Либрационное движение искусственного спутника. Правильнее говорить не об искусственном спутнике, а о твердом теле, размеры которого малы в сравнении с рас­ стоянием до центра тяготения. К таким телам относятся не только искусственные спутники, но и различные ос­ колки космических объектов. Статистическое рассмотрение

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

Учтем рассеяние энергии в либрационном движении в виде квадратичной диссипативной функции (линейного демпфирования). Подобное демпфирование может иметь место за счет вихревых токов, наводимых магнитным по­ лем, за счет разреженной атмосферы, а также за счет жид­ костных демпферов (для искусственного спутника). Ор­ биту будем считать круговой (й = const), однако высоту орбиты полагаем случайной, так что й — const — слу­ чайная величина. Вводя обозначения

ю* = *ы со у = х г , со, = х 3, а = ж4, а' = х ь, а" = ж6,

уравнения свободного либрационного движения запишем в виде(см. главу I)

где * i, Г2, r3 — коэффициенты демпфирования. При бук­

венных обозначениях коэффициентов эти уравнения при­ нимают форму