
1503
.pdf
Рис. 1.2. Скорости движения кончика пера по оси x при многократном написании студентом пароля «Петр»
Рис. 1.3. Реализации последовательностей длительностей нажатия клавиш при наборе пароля «Гранитный барьер»
Рисунки позволяют сделать вывод о существовании непрерывных (см. рис. 1.1, 1.2) и дискретных (см. рис. 1.3) СП. Если провести измерения значений непрерывного процесса в равноотстоящих точках через интервал t, получим аналоги дискретного процесса на рис.1.3. Такой процесс перехода от непрерывной к дискретной функции называется дискретизацией, а такого рода дискретная функция получила название случайной последовательности или импульсного

случайного процесса. В более общей формулировке импульсный случайный процесс есть последовательность импульсов, параметры которых являются случайными величинами.
Имея множество возможных реализаций случайного процесса, можно получить исчерпывающую его характеристику –
распределение вероятностей реализаций. При его построении принимается за основу базовое положение: при фиксированном значении аргумента значения функции есть набор значений случайной величины. При дискретизации аргумента СП через интервал t получим систему случайных величин. Каждая из этих величин Хi={xi1, xi2,…xik,…, xim} описывается плотностью распределения вероятностей ω(Хi), i=1,2,…,n (рис. 1.4). Если Х1, Х2,…,Xn независимы, дифференциальный закон распределения вероятностей СП определяется соотношением
n |
, |
|
|
X Xi |
n , |
t 0. |
|
i 1 |
|
|
|
Когда «условие независимости» не соблюдается, выражения для распределения вероятностей СП X становятся громоздкими. Исключением служит соотношение, получаемое при выполнении условия «нормальности» составляющих вектора X :
X |
|
1 |
|
|
|
|
1 |
n |
A |
(xik mxi) (xjk |
mxj) |
|
|||||
|
|
|
|
exp |
|
, (1.1) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
n |
|
|
|
|
|
|
|
ij |
|
|
|
|
|
|
||
|
xi |
|
(2 ) |
n |
A |
|
2Ai, j 1 |
|
|
xi |
xj |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
i 1
где mxi – математическое ожидание случайной величины Xi; xi – среднеквадратичное отклонение этой величины; А – определитель n-
го порядка, составленный из коэффициентов корреляции случайных
величин в i м и j –м сечениях rij. По определению:
mxi M xik ; xi M(xik mxi )2;
rij M (xik mxi) (xjk mxj)/ xi xj ,

Рис. 1.4. Графики, поясняющие процедуру построения распределения вероятностей СП
Aij= |
1 r12.......r1n |
|
|
|
r21 1........r2n |
, |
rii rjj |
1, |
|
|
rn1 rn2......1 |
|
|
|
где знак М – математическое ожидание; Aij – алгебраическое дополнение элемента rij, определяется через А: Аij = (-1)i+j A/ , A/ – определитель, получаемый из А путем вычеркивания i-строки и j-го столбца.
В выражение (1.1) входят математические ожидания случайных величин Xi и коэффициенты корреляции между ними. Поэтому, если из физических соображений можно сделать заключение о нормальности каждой из указанных величин, то для описания СП достаточно определить эти параметры.
Полагая в формуле (1.1) n=1 и n=2, получим частные соотношения для одной и двух случайных величин:
Х1 х112 exp х1к2 тх21х1 2 ,

|
|
|
|
1 |
|
|
|
1 |
|
|
х |
т |
х1 |
2 |
||||
Х |
1 |
,Х |
2 |
|
|
|
|
exp |
|
|
|
|
|
1к |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|||||||||||||||
|
|
|
|
2 |
|
|
|
21 r2 |
|
|
2 |
|
|
|
||||
|
|
|
|
|
2 х1 х2 |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
1 r12 |
|
12 |
|
|
|
|
х1 |
|
|
|
2r |
х |
т |
х |
т |
|
х |
к |
т |
2 |
|
||
1к |
х1 |
2к |
х2 |
|
|
2 |
х2 |
|
. |
(1.2) |
||
|
|
|
|
|
|
|
|
|
||||
12 |
|
х1 х2 |
|
|
|
|
|
х2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
При ∆х→0 число градаций Х из ряда (1,2,…,к,…,т) стремится к бесконечности и приведенные выражения переходят в непрерывные распределения вероятностей для одной или двух случайных величин. При выводе частных формул (1.2) использовались значения алгебраических дополнений:
A ( 1)2 |
1 1; |
A ( 1)3 |
r ; |
A ( 1)3 |
r ; A ( 1)4 |
1. |
|
11 |
|
12 |
21 |
21 |
12 |
22 |
|
Наряду с нормальными СП следует охарактеризовать так называемые марковские процессы, являющиеся частным видом СП, которые широко используются в разных научно-прикладных направлениях (радиотехнике, автоматике, теории надежности и массового обслуживания, физике, биологии, медицине и др.).
Рассмотрим процесс х t , для которого в последовательные моменты времени t1 t2 ... tn-1 tn его значения определены: x1 xt1 , x2 xt2 ,...,xn x tn . Процесс называется марковским, если условные вероятности
p xn / xn-1,...,x1 |
x1,x2,...,xn |
|||||
x ,x |
2 |
,...,x |
n-1 |
|
|
|
1 |
|
|
|
зависят лишь от последнего значения xn-1, т.е. если справедливо равенство
p xn / xn-1,...,x1 p xn / xn-1 .
Тогда для процессов этого вида справедливы выражения:

x1,x2,...,xn p xn /xn-1 x1,x2,...,xn-1 ,x1,x2,...,xn-1 p xn-1- /xn-2 x1,x2,...,xn-2 ,
- - |
- |
- |
- |
- |
- |
- |
- |
- |
- |
х1,x2 p x2 / x1 x1 .
Подставляя каждое последующее выражение в предыдущее, получим
x1,x2,...,xn p xn / xn-1 p xn-1,/ xn-2 ... p x2 / x1 x1 .
Таким образом, марковский процесс определяется одномерной начальной плотностью распределения вероятности x1 и
вероятностями перехода его из i-го в ј-е состояние.
В литературе можно найти и такое утверждение: любой процесс можно рассматривать как марковский, если все параметры из «прошлого», от которых зависит «будущее», включить в настоящее. Так, если за настоящее состояние системы принять «исправна», то процесс не марковский, потому что вероятность ее отказа в предстоящее время τ зависит от продолжительности ее работы и даты ее последнего ремонта. Но оба параметра (время работы и дата последнего ремонта), включенные в «настоящее», позволяют оценить вероятность перехода в неисправное состояние в будущем и расценивать последовательность переходов как марковский процесс. Подобный подход позволяет решить ряд практических задач.
Задача
Распределение вероятностей случайного процесса Х (t) описывается выражением
|
|
1 |
|
|
|
x |
2 |
|
|
X,t |
|
|
exp |
|
|
e2 t t |
. |
||
|
|
|
|
|
|||||
|
|
2 |
|
|
2 |
|
|
||
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
Определить математическое ожидание m t |
и дисперсию x2 t |
||||||||
|
|
|
|
|
|
|
|
x |
|
процесса. Изобразить на графике область процесса Х (t) в границах ±3 x t , а также построить график X,t для нескольких моментов времени: t1 0, t2 1 , t3 2
.

Решение
Приведем представленное выражение для плотности распределения вероятностей к виду
|
|
1 |
|
|
2 |
|
e |
2 t |
e |
t |
|
||
X,t |
|
|
exp |
x |
|
e2 t |
|
exp t |
|
|
|||
|
|
|
|
|
|
|
e t |
||||||
|
|
2 |
|
2 |
|
e |
2 t |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
x2 |
|
1 |
|
|
|
x2 |
|
||||||||
|
|
|
|
|
exp |
|
|
|
|
|
|
|
|
|
exp |
|
|
|
. |
|
|
|
|
|
|
|
|
|
2 |
|
|||||||||
|
|
|
|
|
|||||||||||||||
|
|
2 exp t |
|
|
2e |
2 t |
|
|
2 x t |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
2 x |
t |
Имеем выражение для нормального закона распределения вероятностей. В каждом сечении рассматриваемого случайного процесса математическое ожидание случайной величины равно
m t =0, а дисперсия |
x2 t e 2 t изменяется во времени. Связи |
x |
|
между случайными величинами отсутствуют (коэффициенты корреляции равны нулю, определитель А имеет значимые элементы только по диагонали rij 0 при i≠j и rij 1 при i=j).
В обозначенных точках аргумента имеем:
t1 0, x (0) 1; t2 1/ , x (1/2) e 1 0,37; t3 (2/ ) e 2 0,137.
Этих данных достаточно, чтобы построить требуемые графики, представленные на рисунке.

а)
Иллюстрация поведения (X,t) и x (t) во времени (начало)
б)
Иллюстрация поведения (X,t) и x (t) во времени (окончание)

1.2. Числовые характеристики случайных процессов
Существует большое число задач, для решения которых оказывается достаточным использование числовых характеристик плотности распределения вероятностей. В отличие от числовых характеристик случайных величин (математического ожидания, дисперсии и др.), представляющих собой числа, характеристики СП являются функциями.
На рис. 1.4 видно, что если соединить точки mхi (значение математических ожиданий в сечениях случайной функции), полученная кривая при ∆t→0 будет одной из искомых функций mx t :
mх t Х t Х dx.
-
Она и есть математическое ожидание СП – неслучайная функция, около которой группируются ее конкретные реализации.
Аналогичным образом определяется дисперсия СП– неслучайная функция x2 t , описываемая при ∆t → ∞ выражением
2
x2 t Х t mx t dx
ихарактеризующая разброс его реализаций относительно среднего. Математическое ожидание и дисперсия – важные, но
недостаточные характеристики для описания основных свойств СП. В этом можно убедиться, обратившись к рис. 1.5, на котором изображены реализации двух случайных функций Х(t), Υ(t).
Рис. 1.5. Сопоставление двух СП с одинаковыми дисперсиями и математическими ожиданиями
Очевидно, внутренняя структура обоих СП совершенно различна, но это не отражает ни математическое ожидание, ни дисперсия. Для описания динамики изменения СП вводится специальная характеристика – корреляционная (другой встречающийся термин автокорреляционная) функция. Она характеризует степень сходства между сечениями процесса, взятого в момент tj и удаленного от него на расстояние ij i j t, ∆t – интервал между сечениями.
Из физических соображений следует, что при малом шаге ij
сходство соседних сечений будет больше для процесса на рис. 1.5, б, нежели для процесса на рис. 1.5, а. В первом случае за малый интервал времени реализации не успевают заметно измениться – сходство с сечением процесса в tj велико. Обратное заключение можно сделать, анализируя СП на рис. 1.5, а. Поэтому процессы Х(t) и Y(t) должны иметь различные корреляционные функции, определяемые как
K |
|
t |
|
|
|
M |
|
o |
|
o |
|
|
|
K |
t |
|
|
|
|
o |
|
o |
|
|
|
|
|
|
j |
ij |
|
X t |
j |
X t |
j |
; |
j |
ij |
M Y t |
j |
Y t |
j |
ij |
, (1.3) |
|||||||||||
|
x |
|
|
|
|
|
|
ij |
|
y |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
o |
|
tj X |
tj mx tj , |
o |
tj |
|
|
|
X tj |
ij mx tj |
ij – |
||||||||||||||
где |
|
X |
X |
ij |
центрированные реализации СП, М – знак математического
o o
ожидания. Аналогично определяется Y tj Y tj ij .
Таким образом, корреляционной функцией СП называется неслучайная функция двух аргументов, выраженная через математическое ожидание произведения центрированных СП при тех же аргументах.
Функция, полученная аналогично, но для разных СП: Х(t), Υ(t):
K |
|
t ,t |
|
o |
|
o |
|
взаимокорреляционной |
В |
j |
M Y t |
X t |
, i,j=1,2,…, называется |
||||
|
i |
|
i |
|
j |
|
||
|
|
|
|
|
|
|
|
|
функцией. Она характеризует степень сходства i-го сечения одного СП с различными сечениями другого СП.
При ij 0, |
K ti,tj M |
o |
2 |
|
x2 ti и необходимость в |
X ti |
|
||||
|
|
|
|
|
|
дисперсии, как отдельной характеристике СП, отпадает: достаточно знать математическое ожидание и корреляционную функцию.
Так как функция K tj,tj ij не зависит от последовательности, в
o |
o |
при ее |
которой используются сомножители X(tj) и |
X t ij |
вычислении, то она симметрична относительно своих аргументов.
На практике часто пользуются нормированной корреляционной функцией
r tj,tj ij K tj,tj ij / x tj x tj ij ,
которая представляет собой последовательность коэффициентов корреляции случайных величин X tj и X tj ij . При ij 0
r ti,ti / x2 ti 1.
Информацию о СП, которую дает корреляционная функция, можно получить и через так называемую спектральную плотность.
В теории сигналов широко используется преобразование Фурье функции времени х(t):
|
j t |
|
|
1 |
|
- jѓЦt |
|
x t F e |
|
d ; |
F |
|
x t e |
dt, |
(1.4) |
|
2 |
||||||
- |
|
|
|
- |
|
|
|
где F( ) называют спектром амплитуд, |
которые приписываются |
||||||
элементарным функциям |
ej t ; |
ω – круговая частота 2 /T ; |
Т – |
период синусоидальной функции.
Аналогичную запись можно привести и для корреляционной
функции. При фиксированном |
tj функция K tj, ij будет |
зависеть |
||
только от интервала ij . По аналогии с (1.4) при ∆t→0 запишем: |
||||
|
|
|
||
K tj,tj Sj e jω d ; |
(1.5) |
|||
|
- |
|
||
Sj |
|
1 |
K tj,tj e- jω d . |
|
|
|
|||
|
|
2 - |
|
|
|
(1.6) |
|