![](/user_photo/_userpic.png)
книги из ГПНТБ / Чепиль А.М. Конспект лекций по элементам теории случайных процессов
.pdfИз теории рядов Фурье и курса ТЭРЦ известно, что одни детерминированные функции имеют дискретный (линейчатый) спектр, а другие — непрерывный (сплошной).
Дискретным спектром обладают периодические и почтипериодические функции *, непрерывным — непериодические функции. Случайные процессы также бывают с дискретными и непрерывными спектрами. Рассмотрим классы таких стацио нарных процессов.
§ 1. СТАЦИОНАРНЫЕ ПРОЦЕССЫ С ДИСКРЕТНЫМ СПЕКТРОМ
Отметим сразу, что стационарный случайный процесс бу дет иметь дискретный спектр только в том случае, когда в его природе заложена «скрытая периодичность» (например, си нусоидальные токи и напряжения со случайными амплитуда ми и фазами). Простейшим примером стационарного случай ного процесса с дискретным спектром является случайная гар моника
g (t) = A cos (a t -Ь ©) ( а — const) .
Спектр ее будет состоять из одной линии.
Возникает вопрос: какой же длины построить эту линию, при условии, что амплитуда А является случайной величиной? Ранее мы вычислили характеристики случайной гармоники
пц — 0 |
и |
Кг (ъ) =■ 52 cos ore, где <з2 = DA . |
Если £ (t) |
— |
случайный ток, то, как известно из физики, |
средняя мощность, развиваемая гармоникой на единичном со противлении в течение периода, пропорциональна квадрату амплитуды гармоники. Покажем, что средняя мощность слу чайной гармоники пропорциональна ее дисперсии. Действи
тельно, |
пусть |
£ (t) — случайный ток |
и |
сопротивление |
|
R = 1 Ом, а |
|
О |
|
|
|
£ (t) — m% = \ (t) — флуктуационная часть этого |
|||||
тока. Тогда полная энергия флуктуационной |
части процесса |
||||
определяется формулой |
|
|
|||
|
|
|
N = b (t ) R = b ( t ) , |
|
|
* Функция f(t) |
называется почти-периодической, |
если она представима |
|||
|
|
оо |
|
|
|
в виде /(<) |
= |
$] |
Ап cos (<o„i+<p„), тде “ л — не кратные некоторому са. |
||
|
|
я=о |
|
|
|
Такой функцией является, например, колебание, периодически модулиро ванное по амплитуде.
80
а средняя мощность флуктуационной части l{t) равна мате-
О
матическому ожиданию квадрата £2 (t):
N cp = М Ь (t) = Ds .
Дисперсия случайной гармоники равна о'\ поэтому ее спектр имеет вид (рис. 6).
Итак, для построения амплитудно-частотного спектра ста ционарного случайного процесса достаточно знать его корре ляционную функцию и построить спектр корреляционной функ
ции, |
так как Кг (0) = |
Ог определяет среднюю |
мощность |
|
флуктуационной части процесса. |
|
|||
Займемся выяснением условий, при которых стационарный |
||||
случайный процесс %(() |
имеет дискретный спектр. |
Эти усло- |
||
бия очевидно связа |
|
|
||
ны |
со |
структурой |
|
|
корре'ляцио иной |
( |
|
||
функции |
процесса, |
|
||
ибо |
фактически нам |
|
|
|
надо |
строить спектр |
|
|
корреляцио иной |
|
_ |
|
||
функции. |
Для |
того |
|
|
|
чтобы спектр корре- |
4 |
w |
- и „ |
||
|
|||||
ляционной |
функции |
|
|
|
|
был дискретен, |
до |
|
|
Рис. 6 |
|
статочно, |
чтобы |
она |
|
|
|
была периодической функцией. Тогда ее можно разложить в ряд Фурье, и, следовательно, спектр корреляционной функции будет дискретен. Если корреляционная функция Кг (т) почтипериодическая функция, то и в этом случае ее можно представить в виде конечной или бесконечной суммы гармо ник различных частот, но не целых кратных некоторому чис лу со, а значит спектр корреляционной функции будет дискре тен. Так как корреляционная функция Кг (т) четная и если она периодическая или почти-периодич^ская функция, то всег да можно представить ее в виде
|
Кг W = |
cos ш* х • |
(О |
|
Если |
соизмеримы, |
то правая часть |
равенства (1) бу |
|
дет рядом Фурье для Кг (т) |
и корреляционная функция — |
|||
периодическая. Если |
несоизмеримы, то |
Кг (т) — почти- |
периодическая функция. В разложении (1) коэффициенты при cos а)Ат всегда будут положительны и равны дисперсиям со-
6. Зак. 525. |
81 |
ответствующих гармоник, ибо каждый из этих коэффициентов пропорционален средней мощности соответствующей гармони ки — величине положительной.
Если корреляционная функция представима в виде (1), то ее спектром будет совокупность чисел (о|, ш*) (/г= 1, 2,...).
Итак, мы пришли к следующему заключению: для того что бы стационарный случайный процесс имел дискретный спектр, необходимо и достаточно, чтобы его корреляционная функция была периодической или почти-периодической функцией.
Остается выяснить, каким условиям должен удовлетворять случайный процесс, чтобы его корреляционная функция была периодической или почти-периодической.
С этой целью обратимся к примерам. Пусть
|
£ (/) “ |
У] |
(a* COS |
t + bk sin |
() |
, |
|
|
*=I |
|
|
|
|
где a.k |
и Ъи — некоррелированные случайные величины с ну |
|||||
левыми |
математическими |
ожиданиями |
и |
дисперсиями |
||
D ak = Dbk = о|, |
a |
— неслучайные (не обязательно целые |
||||
кратные некоторому числу ш) |
действительные числа. В этом |
|||||
случае процесс £ (t ) |
является конечной суммой некоррелиро |
ванных случайных гармоник различных частот. Очевидно, что
процесс £ (t) |
стационарен в широком смысле, так как |
||
гЩ (t) |
= М |V |
(a k cos wk t -f bk sin |
l) |
|
U-i |
|
|
П |
|
|
|
= V |
(M ak cos |
t + Mbk sin &k t ) ~ |
0 |
Kz (*„ /2) = Af | (*,) |
£ (Q |
- |
M \v(akcoswktl + bksinwktl)X |
||
|
|
|
U=i |
|
|
X (afe cos wk t2 -f bk sin (ak г?2)| = V |
M [a\ cos co4 tvcos wk t2 + |
||||
+ ak t>k cos wh |
sin |
t2 + |
bk cos |
sin ш* ty + |
|
-j- b\ sin Wft ti sin (0/^2) *= |
П |
<j2 (cos Ык tl cos оЦ (2 4- |
|||
^ |
|||||
|
|
|
ft=J |
|
|
12
|
п |
~f sin Ш* t l sin CD* t2) |
— V al cos со* (t2 — tt) = |
П |
к -1 |
|
|
= Yi |
° * cos “* x * |
к= |
1 |
где x = t2 — /, .
Процесс E (/) стационарен. Его корреляционная функция
является периодической, |
если со* соизмеримы, и почти-перио- |
дической функцией, если |
со* — несоизмеримы. В том и другом |
случае спектр (°1, со*) |
(/г = 1, 2, ... , п) процесса £ (£) дискре |
тен. Он представлен на рис. 7.
Такую же |
картину |
будем наблюдать и в |
том |
случае, |
||
когда процесс |
? (/) представляет собой сумму |
бесконечного |
||||
ряда взаимно некоррелированных случайных гармоник |
||||||
5 (0 |
= ^ |
cos “а * + |
sin “а 0 |
|
(2) |
|
или, более обще, |
|
|
|
|
|
|
I ( 0 |
= |
Щ + |
2 (a* cos a k t + |
bk sin со* t) , |
(3) |
|
|
|
|
k=i |
|
|
|
где m% — постоянная |
составляющая случайного |
процесса. |
83
При сделанных допущениях относительно случайных величин
а к и bk |
корреляционная функция Ке. (т) процесса 5 (0 |
будет |
|||||
суммой |
|
бесконечного |
ряда |
косинусоид с |
частотами |
со* |
|
(6=1,2, |
3,...), т. е. |
|
|
|
|
|
|
|
|
Л'$ (х) |
= V |
о| cos ш* Т . |
|
|
|
|
|
|
ft-1 |
|
|
|
|
/<^(т) |
будет периодической функцией, если все частоты |
од., |
|||||
кратные какой-нибудь одной из них (со*—6ш, |
6=1, 2, |
3,...), |
и почти-периоднческой функцией, если кратности между час
тотами не существует. |
В обоих случаях процесс 5 (t) |
стацио |
||
нарен в широком смысле и имеет дискретный спектр |
|
|||
(о2, «*) |
( 6 = |
1, |
2, з , . . . ) . |
|
Пусть, наконец, процесс |
|
|
|
|
u |
o = |
£ |
5* в'™*' |
(4) |
|
k—— |
|
|
является суммой бесконечного ряда взаимно некоррелирован ных случайных комплексных гармоник, причем все ЖЕ* = О, a /ЭЕ* = о2. Этот процесс стационарен в широком смысле, так
как (t) — 0, а
Л'«(*., /3) = Af ! ( * , ) ! |
(t2) = M |
V |
|
|
Д-J |
к — ~ <хз |
А=—е* |
fr —« |
зависит только от разности аргументов fi—/2= т:. Если все wk соизмеримы, то корреляционная функция представлена рядом
Фурье в |
комплексной форме, ее спектр (^, ш*) |
(6=0, ± 1, |
± 2 , . . . ) |
имеет вид, изображенный на рис. 8. |
|
Рассмотренные примеры показывают, что корреляционная |
||
функция |
всех стационарных случайных процессов, |
которые |
представимы в виде конечной или бесконечной суммы взаимно некоррелированных действительных или комплексных случай ных гармоник, будет периодической или почти-периодической функцией. Такие случайные процессы называют процессами с дискретным спектром.
Можно доказать, что, если корреляционная функция ста ционарного случайного процесса t (0 является периодиче
84
ской или почти-периодической функцией, то такой процесс представим в виде (2) или (3). Другими словами, всякий ста ционарный случайный процесс с дискретным спектром пред ставим в виде конечной или бесконечной суммы взаимно не коррелированных случайных гармоник с нулевыми средними значениями и дисперсиями of.
Заметим, что если корреляционная функция Kz (т) ста ционарного случайного процесса Е; (t) периодическая или поч- ти-периодическая функция, то она не стремится к нулю при
т-> о о , а это значит, что связь между сечениями |
? (£,) и £ (£,) |
не ослабевает при увеличении разности |
Это обстоятель |
ство используют для обнаружения «скрытой периодичности» случайного процесса.
Итак, стационарные случайные процессы, обладающие «скрытой периодичностью», имеют дискретный спектр. После довательность чисел (°l, coft) (£ = 1 ,2 ,3, . .. ) показывает от
носительный вклад каждой отдельной случайной гармоники в образование случайного процесса. Каждое из чисел of про
порционально средней мощности соответствующей гармоники,
а сумма |
всех |
of равна |
дисперсии |
случайного |
процесса |
|||
£ о? = Di |
|
и пропорциональна средней мощности флуктуаци- |
||||||
I: |
|
процесса t |
(/). |
В связи |
с этим последователь |
|||
онной части |
||||||||
ность чисел |
(of, |
шй) |
называют |
энергетическим |
спектром |
|||
данного процесса. |
|
|
|
|
|
|||
С точки зрения корреляционной теории задание среднего |
||||||||
значения |
пц и |
спектра |
(of, t»k) |
полностью определяет ста- |
S5
Ционарный случайный процесс ё (t), ибо по спектру легко вос становить его корреляционную функцию
|
|
K |
t (t) = |
2 |
° 1 COS |
СО* X , |
|
|
|
|
к |
|
|
если |
Ё (0 — вещественный случайный процесс, и |
|||||
|
|
|
/<, (т) |
= |
2 |
, |
|
|
|
|
|
к |
|
когда |
ё (/) |
— комплексный случайный процесс. |
||||
Полезно |
ввести |
в рассмотрение |
понятие спектральной |
функции стационарного случайного процесса Ё (t). Спектраль ная функция определяется равенством
|
5 Н |
= £ |
. |
(5) |
Спектральная функция пропорциональна средней мощности |
||||
|
О |
|
|
< ш, и явля |
флуктуации Ё (£) на всей полосе частот, где |
||||
ется функцией спектрального |
распределения |
энергии. Эта |
||
функция обладает следующими очевидными свойствами: |
||||
1. 5 |
( — с о ) = 0 . |
|
|
|
2 . 5 |
( 4 - о о ) = D ? . |
|
|
|
3. Спектральная функция — неубывающая функция.
4. Спектральная функция S (ш) непрерывна слева
^ (ша — 0) = S (wft) .
§2. СТАЦИОНАРНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
СНЕПРЕРЫВНЫМ СПЕКТРОМ
Пусть стационарный в широком смысле случайный процесс Ё (t) не принадлежит к классу случайных процессов с дискрет ным спектром. Тогда его корреляционная функция Ki (х) не будет периодической или почти-периодической, и, следователь но, она не представима в виде суммы ряда косинусоид на всей оси Ох. Предположим, что процесс Ё (t) СК непрерывен. Тог да, как известно, его корреляционная функция Ац (т) будет непрерывной в обычном смысле при всех г. Если корреляци онная функция абсолютно интегрируема на всей оси Ох и в каждом конечном промежутке (а < х <Ь) удовлетворяет усло
86
виям Дирихле (достаточным условиям разложения функции й ряд Фурье на данном интервале), то она может быть пред ставлена интегралом Фурье. В силу четности корреляционной функции, она представима косинус-интегралом Фурье
К , (*) = А I cos сот с/ш i К(. (и) cos ши du . |
(1) |
7Г |
|
ОО
Введем обозначение
s5 (ш) = — /<£ (и) cos mi du . |
(2) |
о
Теперь равенство (1) |
можно записать так: |
|
|
со |
|
АГе (т) = |
J se (ш) cos шт d u . |
(3) |
|
о |
|
Равенства (2) и (3) представляют собой два взаимно об ратных косинус-преобразований Фурье для функций К е (х) и Se (со). Функцию
, Л |
2 |
(т) cos ШТ dx |
( 4 ) |
SE (ш) = |
--- |
естественно назвать действительной спектральной плотностью амплитуд.
Определение. Стационарные случайные процессы, корре ляционные функции которых представимы интегралом Фурье, называются случайными процессами с непрерывным спектром.
Хотя мы определили стационарный случайный процесс с непрерывным спектром как процесс, у которого корреляцион ная функция представима интегралом Фурье, можно доказать, «то функция (т) будет корреляционной функцией стацио нарного в широком смысле случайного процесса £ (/) тогда и только тогда, когда она' представима интегралом Фурье. От сюда следует, что все СК непрерывные случайные стационар ные процессы — суть процессы с непрерывным спектром и са ми допускают представление, аналогичное интегралу Фурье. Приведем без доказательства следующую важную теорему.
87
Теорема. Всякий СК непрерывный случайный и стационар- > ный в широком смысле процесс представим в виде
со
Ъ(t) = т* + I" {Zj (ю) cos cot + Z2 (со) sin ш/} d<n , (5) 6
где случайные процессы Z,(u>) и Z2(co) взаимно некоррелированы, имеют нулевые математические ожидания и удовлетво ряют условию
М {Z; (со + Дш) — Zl (to)}2 = |
М {Z2 (со + Аш) — Z2 (со)}2 . |
||
В этом случае формулу (5) |
называют спектральным разло |
||
жением процесса £ (t). |
|
и |
Z2 (со), входящие в формулу |
Случайные процессы Z, (ш) |
|||
(5), могут быть определены с помощью равенств |
|||
Zj (ш) = Нш |
|
|
cos Ы dt |
Т->* |
|
|
|
и |
|
т |
|
|
|
|
|
Z2 (ш) = lim |
2тс |
j* |
5 (t) sin wt dt . |
т-+- |
-T |
|
|
|
|
|
Вместо полного доказательства теоремы приведем лишь нестрогие наводящие соображения.
Так как корреляционная функция h\ (т) удовлетворяет условиям Дирихле в любом конечном промежутке, то ее мож но разложить в ряд Фурье в интервале (—Т, Т). В силу чет ности корреляционной функции, ее ряд Фурье будет иметь вид:
Ки (О - |
£ Dk cos Uik X , |
(6) |
|
А=О |
|
где |
|
|
|
К * |
|
“ а |
т |
|
и |
|
|
Dk |
COS (i)j х rfx , |
|
88.
Можно доказать, что корреляционная функция |
/<= (т) — |
положительно определенная функция * и поэтому все |
0. |
Равенство (6) справедливо только в интервале (—Т, |
Т), так |
как Ki (т) — непериодическая функция. На основании рассуждений предыдущего параграфа можно предполагать, что
равенство (6) определяет корреляционную функцию |
неко |
|||
торого стационарного |
случайного |
процесса |
с дискретным |
|
спектром. Обозначим этот процесс |
а его корреляцион |
|||
ную функцию — Л^(т). |
Тогда процесс Цг (£) |
представим в |
||
виде |
|
|
|
|
£г (t) = ПЧ + |
Ч cos |
t + bk sin Щ t , |
(7) |
Л= 1
где a k и b k — взаимно некоррелированные случайные величины, причем M ak =* МЬк = 0, D ak = D bk — Dk. Преобразуем ряд Фурье для корреляционной функции, домножив и поделив его
is ■
K l СО= £ |
cos |
(8) |
А=о |
|
|
Так как |
|
|
/гк |
k + |
1 |
ША-Н = ---- j,---- те |
||
ТО |
|
|
ДсОд — u)/, f-i |
шЛ — |
jt |
D Т
и величина — -— представляет собой среднюю плотность энертг
гии, приходящейся на участок Д©4. Обозначим эту среднюю плотность s r (шА) и ряд (6) перепишем так:
/<1 ( т ) |
= |
£ |
Sr ( ш* ) |
c o s |
^ |
Дсо* |
• |
(9 ) |
|
|
|
fc=0 |
|
|
|
|
|
|
|
* Функция f { t ) действительного |
переменного |
t называется |
положи |
||||||
тельно определенной в интервале — |
оо < |
t < оэ, |
если, |
каковы' бы ни былп |
|||||
действительные числа |
tt, |
t2, . : . , t n, |
комплексные |
числа z,t г2, . •.. , г п и |
|||||
целое число п |
|
п |
|
|
_ |
|
|
|
|
п |
|
|
|
|
|
|
|
||
Е |
Е |
/ |
- |
'*) |
2* |
> |
0 . |
|
|
А= 1 ш = 1
89