Файл: Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи.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(a2—x*), |
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, |
?(*,) = О, |
4» (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