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

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

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

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

Добавлен: 28.06.2024

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

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

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

спектральными плотностями /Ѵо/2, то коэффициенты сно­ са и диффузии определяются соотношениями

КІ{Х, t)=Fi{-x,

/ ) ;

(3.175)

N

 

 

*<i(0 = ^ J ] G f t ( ' ) G i f c ( 0 .

(3.176)

Ограничение, связанное с равенством спектральных плот­

ностей для всех шумов nf(t),

не является существенным.

Из дальнейшего будет ясно,

что наиболее интересный

для практики случай имеет место, когда в системе (3.174) имеется лишь один белый шум n(ty.

Решение многомерного уравнения

Фоккера — План­

ка — Колмогорова представляет собой

более

сложную

задачу,

чем решение

одномерного уравнения.

Однако

в одном

частном, но

важном случае

удается

получить

общее решение многомерного уравнения. Именно: если

коэффициенты сноса — линейные функции

пространствен­

ных координат

 

 

 

 

 

КІ

(t, x) = £ ацхі

(t) + atf

.

(3.177)

а

коэффициенты

диффузии — постоянные

величины

 

 

/ ( „ • ( / ) = const,

(3.178)

то

многомерный

процесс x(t)

является

н о р м а л ь н ы м

процессом [2, 5, 7, 9]. При этом

любой компонент

много­

мерного марковского нормального процесса в стационар­ ном режиме обладает дробно-рациональной спектральной плотностью. Справедливо и обратное утверждение: если нормальный процесс имеет дробно-рациональную плот­ ность, то его можно рассматривать как компонент много­ мерного марковского процесса.

При выполнении условий (3.177), (3.178) система сто­

хастических уравнений

(3.174) примет вид

N

N

 

(3.179)

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

120


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

Как было выяснено выше. (§ 3.8), одномерные мар­ ковские процессы недифференцируемы. В то же время характерной особенностью реальных процессов является их дифференцируемость, вызываемая' инерционностью реальных устройств. В связи с этим возникает вопрос о создании математической модели, которая полнее отра­ жала бы свойства реальных процессов. Такую модель можно создать на основе •многомерных диффузионных процессов. Для этого обратимся к системе (3.179) и ограничимся сначала простейшим случаем N=2. Посто­ янные коэффициенты в (3.179) зададим таким образом,

что система уравнений примет вид:

 

^

= « „ * , + « „ * я + О / г ( 0 ;

(3.180)

 

*щ-=хх.

 

 

(3.181)

Подставляя

(3.181)

в (3.180),

получаем стохастиче­

ское дифференциальное

уравнение

в т о р о г о

порядка

или, обозначив

 

 

 

 

 

 

x2{t) =

x[t),

 

 

? * . - a t l % - a l

t x =

Gn(t).

(3.182)

Процесс x(t)

в уравнении (3.182) описывает выход­

ное колебание динамической системы второго

порядка,

на вход которой

воздействует

белый шум n(t).

Следует

отметить, что x(t), будучи компонентом двумерного диф­ фузионного процесса, сам по себе марковским процессом не является. Действительно, для решения уравнения вто­ рого порядка (3.182) необходимо задавать в начальный момент to не только значение координаты хо, но и ее про­ изводную ^ . Однако задание производной (ско-

роста) равносильно заданию д в у х значений х в близкие

121


моменты времени, например x(t0) = Х о и x(t0—t)

й—Ах,

поскольку

dx/dt =

Um(Ax/x).

 

 

 

t-*0

 

В этом

случае

процесс x(t,) в «будущем»

(t>t0) за­

висит не только от «настоящего» (to), но и от «прошло­

го» (to—г). Следовательно,

марковское

свойство не удов­

летворяется и процесс x(t)

не является

марковским. Из­

ложенные

соображения,

конечно,

относятся

(и даже

в большей

степени) к процессам, описываемым диффе­

ренциальными

уравнениями

более

высокого

порядка.

Из корреляционной теории случайных процессов (см.,

например, (7]) известно, что для определения

спектраль­

ной плотности

стационарной

функции

y(t\),

связанной

с другой функцией z(t) уравнением вида:

* ( 0 = К Sr У (0 + К - g ^ r у {t) + .-.• + ьту (t) =

=Pm{^)y{t),

необходимо разделить спектральную плотность Sz(co) на квадрат модуля полинома, получаемого из Pm(d/dt) за­ меной d/dt на /со, т. е.

5 » < " > = т а т -

( 3 - 1 8 3 >

На основе (3.182), (3.183) легко находим спектральную плотность процесса x(t)

5 * > ) = W ^ T -

( З Л 8 4 )

По аналогии с (3.180), (3.181) "можно составить си­ стему из N' уравнений

dx

 

 

 

 

 

(0;

—= .л+ .2^ + • • •+ ,л +

 

dt

а

а

3

а

G n

 

dt

=

*,;

 

 

 

(3.185)

dx

 

 

 

 

 

 

dt

 

 

 

 

 

которая эквивалентна стохастическому дифференциаль­ ному уравнению N-ro порядка

dNxN

d»-% "N-2~

..-amxN=Gn(t).

(3.186)

dtN

1 1 dtN~l

1 2 dtN~2

 

 

122


Применяя к (3.186)

формулу

(3.183), находим

спектраль­

ную плотность компонента

Хк(і):

 

 

 

 

 

Sx

(«) = •

G*

ff.

 

(3.187)

 

 

 

 

Физически

соотношения

(3.185) — (3.187)

описывают

выходной процесс xN(t)=x(t,)

линейной системы Л/'-го

порядка, на

вход

которой поступает белый

шум

n(t).

В силу линейности системы и нормального

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

белого шума

выходной процесс x(t)

будет

нормальным.

,y(t) H А

x(t)

n(ti

 

yd)

 

x(t)

 

>•

 

 

 

 

о)

 

 

5)

 

 

 

 

 

Рис.

3.13.

 

 

 

 

Но нормальный процесс с

дробно-рациональной

спек­

тральной плотностью является компонентом многомерно­ го марковского процесса, причем общее число компо­ нентов N равно степени полинома, стоящего в знамена­ теле выражения для спектральной плотности (3.187).

Если процесс x(t), сформированный линейной систе­ мой ІѴ-го порядка, подать на вход другой линейной си­

стемы т-го порядка, то выходной

процесс последней

будет компонентой (JN + m)-мерного

нормального мар­

ковского процесса. Это обстоятельство важно в следую­

щем

отношении. Пусть имеется

некоторая

система А

т-го

порядка

(рис. 3.13,а),

на вход которой

поступает

случайное колебание

y(t)

со свойствами, далекими от

свойств белого

шума.

Необходимо

вычислить

основные

характеристики выходного процесса x(t).

Если подобрать такую систему В (рис. 3.13,6) /-го порядка, на выходе которой под воздействием белого шума п(і\) образуется колебание y(t), то x(t) можно рас­ сматривать как компонент (т + 1)-мерного диффузионно­ го процесса и применять к его исследованию теорию мар­ ковских процессов.

Многомерные диффузионные процессы целесообразно использовать и при анализе нелинейных систем [7]. Пусть система А (рис. 3.13) нелинейна и уравнение, связываю-

123


щее ее выход и вход имеет вид:

^

+

с ,

*£*+\..\+ст_х%+f(x)=y{t),

 

 

 

(3.188)

где f(x)

— нелинейная функция.

 

 

 

 

Если выбрать указанным выше способом линейную

систему

В, то процесс y(t) .будет

иметь

дробно-рацио­

нальную плотность и, следовательно, может

рассматри­

ваться

как компонент

/-мерного

диффузионного про­

цесса.

 

 

 

 

 

 

 

 

Обозначая

 

 

 

 

 

 

 

 

 

у _ „

. d x

dli±i.

1l

 

 

 

 

 

5 - =% - =у г + 3 ; -- . '•• •

 

(3.189)

 

 

 

rf/™-1

~~

dt

—Уі+™>

 

 

вместо

(3.188) получаем

уравнение

 

 

 

%*

= - С^'+™ - -

- с

' - f (Уі») + Уі-

(3.190)

которое

в совокупности

с уравнениями

(3.185), (3.189)

составляет

систему

1 + т дифференциальных уравнений

первого 'порядка, в .правую часть одного из которых вхо­

дит белый

шум.

Следовательно,

каждая из функций

Ui(t) ( і = 1

, 2,

1 + т), в том числе, конечно, и x(t) =

=yi+i(t,),

является компонентом

(/ + т)-мерного марков­

ского (но уже не нормального) процесса. Это очень важ­ ное положение. Оно свидетельствует^о том, что задача

определения

закона распределения

функции на выходе

существенно

нелинейных систем

п р и н ц и п и а л ь н о

разрешима для широкого класса

задач, охватывающего

многие практически интересные

случаи.

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

124