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

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

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

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

Добавлен: 28.06.2024

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

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

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

2. В том случае, когда

 

КІ(Х,

t) =a(t)

(х—с),

 

 

Kz(x,

t)=b{t){x—c)\

(3.126)

функции ><p (t) и i]) (х, t)

следует выбрать в виде:

 

 

 

t

 

 

 

 

о

 

 

 

 

t

 

 

Ф (X,

t) = In (X - с) + j

[6 (/") - а (?)] dt'.

(3.127)

 

 

о

 

 

Проводя

соответствующие

выкладки, можно

прийти

к соотношениям (3.122), (3.123), так что для новой пере­ менной снова получается нормальный закон распределе­ ния (3.124).

Подбор функций <p(t) и \J)(x, t) необязательно дол­ жен быть подчинен выполнению условия (3.122). Так, в |16] описан метод, который позволяет с помощью фор­

мул

преобразования

свести

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

Ki(x, t),

Kz{x,

t) к виду КІ{ХІ),

Kiixi)

соответственно.

Незави­

симость коэффициентов сноса и диффузии от времени

упрощает задачу

решения

уравнения Фоккера •— План­

ка — Колмогорова

особенно, если

ограничиться опреде­

лением стационарного распределения.

 

Как видно, метод введения новых переменных дает

возможность получить

.решение

уравнения

Фоккера —

Планка — Колмогорова

в

полном

объеме, включая ста­

ционарный и нестационарный режимы. К

сожалению,

число охватываемых этим методом вариантов, огра­ ничено.

Решение уравнения Фоккера — Планка — Колмогоро­ ва в полном объеме возможно также, если применить хорошо известный в математической физике метод 'раз­ деления переменных (метод Фурье) [4, 17]. Этот метод весьма трудоемок, а получаемое с его помощью решение обычно выражается через сложные функциональные ряды. Мы не будем рассматривать сущности этого

метода и укажем лишь на то, что отдельные

примеры

его применения можно найти в работах {7, 9, 12].

Займемся

теперь

вопросом о

получении

с т а ц и о ­

н а р н о г о

решения

уравнения

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

Колмогорова. Необходимым условием для этого являет-

102


ся -независимость коэффициентов сноса и диффузии о.т времени. В стационарном режиме w(x, £ )= ш С т(*) и dwCT(x)/dt = 0. Поэтому в соответствии с (3.90), (3.91) имеем

Обычно 'предполагают, что, на границах,

которые

в данном случае удалены в бесконечность, поток

 

0(—oo) = G( + oo)=0.

(3.129)

Если это так, то интегрирование уравнения (3.128) дает, что в стационарном режиме ноток равен нулю для всех точек области

G (х) = /С, (x) шс т (Х)—^.-*- [К, (х) даот (x)} = 0. (3.130)

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

нию. В переходном режиме подобного равенства

не на­

блюдается.

 

 

 

В задаче о воздействии белого шума на і?С-цепочку

(§ 3.5) в

переходном режиме функция

w(x, t)

имела

крутые спады

(рис. 3.8) и, следовательно, большой гра-

 

dw(xj)

 

 

 

диент

^

• Поэтому в переходном

режиме

диф­

фузионный поток оказывается большим по величине, чем конвекционный и результирующий поток приобретает направление, определяемое диффузионным потоком. Вследствие этого дисперсия процесса x(t) растет, плот­ ность w(x, t) расплывается. Однако свойства системы (RС-цепочки) таковы, что вместе с увеличением коорди-

*' Мы рассматриваем одномерное уравнение, следовательно, и поток должен определяться не через площадку, а через линию — се­ чение одномерного процесса на уровне х.

103


наты х растет скорость КІ{Х,

t) =—ах

детерминирован­

ного 'потока, направленного

п р о т и в

диффузионного

потока. Когда эти потоки сравняются по величине, на­ ступает стационарный режим.

Не следует думать, что в переходном режиме диффу­ зионный поток всегда пересиливает конвекционный. Соотношение между потоками целиком определяется ха­ рактером .начальных условий. Действительно, если на­ чальное распределение Wo(x, to) имеет большую диспер­ сию, чем у стационарного распределения адОт00, то до наступления стационарного режима конвекционный по­ ток больше диффузионного. Очевидно, когда Wo(x, U) =

=WCT (x), переходного режима не существует.

Вернемся к уравнению (3.130). Это линейное диф­ ференциальное уравнение первого порядка для функции

WCT(X):

4т [Kt (x) даст (x)] - 2/С, (x) шс т (x) = 0.

(3.131)

Общее решение этого уравнения легко находится

X

 

w " w = і ш е х р [ 2 j f r ë r d x ' ] '

( 3 - 1 3 2 )

где С — постоянная интегрирования.

Нижний предел х \ в интеграле произволен, поскольку его изменение эквивалентно изменению постоянной инте­ грирования С. Совместный выбор этих величин подчи­ нен условию нормировки. В дальнейшем в большинстве случаев будет показано Хі = 0.

Таким образом, зная коэффициенты сноса и диффу­ зии, можно непосредственно записать формулу (3.132) для стационарного распределения wCT(x). Это свидетель­ ствует о высркой эффективности использования уравне­ ния Фоккера — Планка — Колмогорова.

Приведем простой пример [18] (см. также '[10]). Пусть

процесс x(t)

описывается

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

дифферен­

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

 

 

x =

kx(a--x2)

+ n(t), & > 0 ,

(3.133)

где a=const, n{t) —белый шум. Уравнение (3.133) пред­

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

(3.100), поэто­

му в соответствии с (3.101), )3.102) имеем

 

Ki(x)=kx(a2x*),

K2(x)=Nol2.

(3.134)

104


После подстановки (3.134) в (3.132) получаем

^ с т M = С ехр

-щ- (2а°х~ — х")

(3.135)

где

00

 

 

С-1 =

J ехр ^

(2аV - л?)

 

Интересно отметить, что характер плотности вероятно­ сти (3.135), выражение для которой получено без вся­ ких усилий, оказывается довольно сложным (рис. 3.11).

3.8.Нормальный марковский процесс

 

Снова обратимся к за'даче о воздействии белого шу­

ма

на интегрирующую /?С-цепочку.

В § 3.5

было со­

ставлено уравнение

Фоккера — Планка — Колмогорова

(3.82) для плотности

вероятности

выходного

процесса

х(і),

но решения его не проводилось. Метод А. Н. Кол­

могорова, изложенный в предыдущем параграфе, позво­ ляет восполнить этот пробел. Действительно, соотноше­ ния (3.80), (3.81) являются частными случаями выражений (3.119), поэтому необходимо лишь конкре­ тизировать результат (3.125).

Для рассматриваемого случая из (3.120), (3.121) на­

ходим

 

 

«К*,, t0)

= x0,

?(*,) = О,

(JC.

t) =

xec

 

 

(3.136)

105


Подставляя (3.136) в (3.125) и учитывая

(3.80),

(3.81), получаем

 

( ) = 1 7 S f r e s p { - " ~ У и ' ~ " ' ' ) • < 3 - 1 3 7 '

где

 

а ; ( 0 = < [ 1 - е - 2 " { ' - ' ' ) ] , os = - ^ .

( З Л 3 8 )

Очевидно, (3.137), (3.138) совпадает с (3.83),

(3.73),

(3.74).

Результат (3.137), полученный исходя из общих пред­ посылок, позволяет сделать следующий важный вывод: если коэффициент Кі(х) является линейной функцией х, а коэффициент Кі{х)—постоянная величина, то одно*- мерная плотность распределения марковского процесса

является г а у с с о в о й .

При t—>-оо

процесс

вступает

в стационарный режим, для которого

из (3.137), (3.138)

получаем

 

 

 

тв„ (х) =

1 _ ехр (—-ff

У

(3.139)

Для закона (3.139) характерна независимость пара­ метров от начальных условий, что иллюстрирует эффект

«забывания»

марковским процессом своего

прошлого.

Конечно, плотность вероятности

(3.139)

легко найти и

из соотношения (3.132).

 

 

 

 

Вспомним теперь тот факт, что при дельтообразных

начальных

условиях

решение

уравнения

Фоккера —

Планка — Колмогорова

для одномерной

плотности ве­

роятности w{x, t) совпадает с соответствующим

решени­

ем

уравнения для

плотности

вероятности

перехода

v(x,

t\x0,4о). В связи с этим, не обращаясь к уравнению

(3.56), на основе (3.137), (3.138)

запишем

функцию

плотности вероятности

перехода

 

 

Из (3.140)

следует, что плотность

вероятности перехода

зависит от р а з н о с т и

моментов

времени (t—ta),

следо­

вательно,

процесс x (t)

является

о д н о р о д н ы м

мар-

106