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