книги из ГПНТБ / Коновалов Г.В. Импульсные случайные процессы в электросвязи
.pdfСоотношения (3.49) и (3.50) можно преобразовать с учетом им пульсного характера воздействия:
N оо
т с(*) = £ |
2 |
р(Хл = |
*) Yi |
r ( ö |
|
|
(3.49a) |
Re V г(ü#p) |
|
|
|
||||
п=0 х=0 |
г = 1 |
|
|
|
|
||
N |
£оо |
|
|
|
|
|
(3.50a) |
ти (/) = V |
Р (х„ = |
X) V |
Jm Ап_r(top)^ ,» |
|
|||
n=0 v.=0 |
r=1 |
|
|
|
|||
где |
|
|
|
|
|
) |
|
|
|
|
|
— 1 tö„ |
6 * ) |
(3.51) |
|
|
|
|
|
i p |
n-r |
. |
Перейдем теперь к определению моментных функций второго порядка. Согласно (3.14) и (3.36) корреляционная функция слу чайного процесса t,(t) на выходе узкополосной системы
U+т
В, (t, t + x ) = B <.(t, t) = j j B^ (vt, D2) h0 (t — vj) h0 (t 4- X— oa) X
о0
|
XCOS [(0P |
(t — Vx) |
+ |
|
(t |
— üj)] |
|
D |
|
dvt d D |
|
|
|
|
|
||||||||||||||
|
X |
— |
|
D |
|
ф (t |
|
|
|
|
|
|
|
||||||||||||||||
•cos [ü)p (/ - f |
|
|
|
2) + |
ф |
- f T — |
2)] |
|
|
|
2. |
|
|
|
|
|
|||||||||||||
|
С другой стороны, согласно определению корреляционной функ |
||||||||||||||||||||||||||||
ции можно записать для процесса t,(t): |
|
|
|
|
|
|
|||||||||||||||||||||||
В, |
(/, |
t + |
x) = |
|
m1{ 1{к) |
(0 |
l (k) (t + |
x)}. |
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Учитывая |
|
(3.41), |
|
после ряда тригонометрических преобразова |
||||||||||||||||||||||||
ний корреляционную функцию процесса £(7) запишем в виде: |
|||||||||||||||||||||||||||||
в,, (t, |
X) = -L { [£с с ((, т) - f |
ßc s (С т) cos (0Р т — |
|
|
|
|
|||||||||||||||||||||||
|
- [ |
ßCcS^> T) |
|
- |
ßCscV' Т)] Sin CDp Т + |
|
|
|
|
|
|
|
|
||||||||||||||||
|
+ |
[ßc с (t, x) — BCs (t, T)J cos (Op (21+ T) — |
|
|
|
|
|
||||||||||||||||||||||
|
— [Bc cs (t, X) |
i - B , sc (t, T)] sin ü)p (2/ + |
|
T)}. |
|
|
|
|
(3.52) |
||||||||||||||||||||
|
Запишем далее соотношения для автокорреляционных BQc(t, т), |
||||||||||||||||||||||||||||
B , s(t, х) и взаимно корреляционных Bi cs(t,,x), |
ßg sc(t, т) |
функций |
|||||||||||||||||||||||||||
процесса t,(t). |
|
|
|
|
|
|
составляющей £с(0 |
корреляционная |
функция |
||||||||||||||||||||
|
Для |
косинусной |
|
||||||||||||||||||||||||||
Bu (t, Т) |
t |
f+т |
|
|
|
|
|
Da) h0 |
(t |
— |
Dl) |
hü (t + |
T • |
■ V2) X |
|
||||||||||||||
=J |
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
I |
В% |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
о |
|
0 |
|
|
|
|
(Dl, |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
X COS |
|
|
— |
ф |
(t |
— ■D l)] COS [(Op Da — ■ф |
{t |
+ |
X |
— |
D |
dv du2. |
(3.53) |
|||||||||||||||
|
|
|
[(Op D l |
|
|
|
|
|
2)[ i |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
120
После .преобразований можно представить это выражение в виде
t і+ Х —ѵ
B^c(t, |
|
|
-Jа |
|
v , t + T |
— z)Re{h1{v)X |
|
т> |
Оj |
Bl (t — |
|
t—t + xv—v |
|
,----------------- |
i |
(ö„ (T—2) . |
|
J B^(t— V, t + X- |
||
Xh1(v + z)e |
|
' ] d ü d z + - ^ - J |
0—и
—ü — z) Re [Aj (v) Aj (i> -f- г) ei top (2 i 4 - z ) — І (Op ( 2 / + T ) ] cfocfc
или, раскрывая ßg (А ^+т), можно записать
N N оо оо
S" ('’ T, = T S S S S p<z"=,‘' х' " Х) х
n=0 /=0 x=0 X=0
X X
X ^ ^ |
[Re |
j, r, q (®p> |
|
®P^ " Re V /, r 9 (®p> ШР^]- |
(3.55> |
||||||
|
r=l ?=1 |
|
|
|
|
|
|
|
|
|
|
В (3.54) |
введены обозначения: |
|
|
|
|
||||||
A„. |
2., Ч |
- |
“p) = Щ { £<*>r 6}f) T<*V T<*>X |
|
|||||||
|
|
|
____________ |
— i CDp ( |
<<*) |
- |
<<fe) ) I |
(3.56> |
|||
X |
к |
|
V)g,+T (C0p T<*),)e |
|
n’r |
|
,' q ], |
||||
k n. i, |
r, q Ч - |
®P> = |
ml { |
|
б)*», |
|
|
|
|
||
|
|
|
|
|
|
— i Ир ( |
<(*) |
+ |
)) |
(3.57). |
|
X g« 4 |
*<,*>,)£,+, к |
*}%) e |
|
|
|
' ’ / |
|||||
Аналогичные рассуждения приводят к таким выражениям для |
|||||||||||
корреляционной |
функции |
синусной составляющей процесса |
£(/): |
||||||||
|
|
|
t t + X - v |
|
|
|
|
|
|
|
|
В,. s (t, т) = |
|
1 (‘ |
ß g |
(t — |
u, |
t-{- x—v—z) Re (A t (V) Ax (v + z) x |
|
||||
|
|
|
О |
—V |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t t + x —v |
|
|
|
|
X e “p( |
+ Z) ] dv dz |
|
- J |
I |
ß j R |
— |
V, t + x — v — z) X |
|
О—а
|
r , |
, |
І GJn (2U+ 2) |
— І 0)„ (2f+ T ) |
(3.58). |
||
X Re [Ax (o) hx(у X 2) e p |
e |
p |
J dv dz |
||||
ИЛИи |
|
|
|
|
|
|
|
|
|
n=0 j = |
О X—0 >1=0 |
|
|
|
|
X |
Ä. |
|
|
|
|
|
|
|
£ |
[Re A»- ь '• *(ö)p’ ~ Шр) ~ |
Re |
/.'■*(u)p’ |
(3.59> |
||
|
|
||||||
r = ] |
q = z 1 |
|
|
|
|
|
|
121
Взаимно корреляционные функции квадратурных составляющих тоже могут быть найдены по той же методике. Так, для ßg Cs(t, х) получим
t t + t — V
|
|
|
О —V |
|
|
|
|
z) |
X |
|||
|
|
|
І СОт%(т—z) |
|
|
|
|
|||||
X |
|
|
„---------- |
]dvd z— |
|
|
||||||
X Jm [ Іц (ü) ht {V+ |
г) e |
p |
|
|
|
|||||||
|
t <+t—0 |
|
|
|
|
|
|
|
|
|
||
---- ^-j* |
j* |
ByXt — ü ,/ + Т — V— z) Jm [hi (v) hi (v + z)X |
||||||||||
|
о |
—о |
e— i(Opp (21 + |
|
|
|
|
|
|
|
||
X ei (Opp (2 o + z ) |
t) |
]dvdz. |
|
|
|
|
||||||
или с учетом импульсного характера воздействия |
||||||||||||
|
|
|
N |
N |
|
оо |
со |
|
|
|
|
|
ßc с 5(/, т) = |
ТГ 5] |
5] |
(Хя = |
*. X/ = *) X |
|
|||||||
X X |
|
л=0/=0 х=0 \=0 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||
х £ |
5 J |
[Jm v |
/. г,, к - “ р) - |
|
Jm к , і, г,, к > |
—®р)]. |
||||||
Г=1 <7=1 |
|
|
|
|
|
|
|
|
|
|
||
а для ße зс(Т, т) соответственно: |
—v+ t —z)x t |
|
||||||||||
|
|
|
t t + 1 — V |
|
|
|
||||||
|
|
|
0 |
J |
ß |
|
/ + Т — V |
|||||
|
|
|
|
—V |
|^—v’ t |
|
|
|
|
|
||
|
|
|
|
|
|
Шр ( z) ] dv dzz -----J |
J B^{t — v,t |
|||||
X Jm |
[hx (v) hx (V+ |
2) e |
||||||||||
|
|
|
|
|
|
|
|
|
|
0 |
|
— V |
■y + T — z) Jm [hi (V) h± (v -f- z) e |
Mp <2o + z) e |
' Ир (2<+T) ]d ойг |
||||||||||
или |
|
|
N |
N |
|
оо |
оо |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
4..<м |
-Ф 5 ]Е 2 1 И 'ос- . *і = ч х |
|||||||||||
|
|
|
п—О/=0 и=0 Х=0 |
|
|
|
|
|
||||
- S |
S |
[Jm |
/, г |
q(®р> ®р) + |
-*т |
hn, i, r, q(Юр —(Öp)]. |
(3.60)
(3.61)
(3.62)
(3.63)
Г — 1 < 7 = 1
Средняя по множеству мощность процесса £(t) определится как значение корреляционной функции в совпадающие моменты времени, т. е.
122
где Bzc(t, 0), BiS(t, 0), Bs-cs(t, 0) и BQSc(t, 0) определяются ф-ла ми (3.54) — (3.63) при т = 0.
Дисперсия процесса t,(t) на выходе узкополосной линейной си стемы может быть выражена через вторые центральные моменты квадратурных составляющих:
DI (t) = В, ((, 0) — т2к (t) = |
|
|||
= |
Jt |
W 2 с с {t) + |
М2с s (^)] 4——[44г с с (0 — М2; s (0] cos 2 озр t |
|
|
|
£ |
|
|
---- —\М2ССs U) + М2Сsc(0] s'n 2 (Op іf. |
(3.65) |
|||
Здесь величины M2iC(t), M2es(t), |
Cs(t) и M2ssc(t) определяются |
|||
формулами: |
, |
|
M2 Cc{t) = Dc (t) = B. c(t, 0) — m\c(/)
M2c s (/) = Ds (t) = B . s(t, 0) — m2s (t)
|
|
(3.66) |
M2 , cs (t) = B. c s (t, 0) — mlc (t) mu (t) |
||
M 2 t s c W = B ; s c |
° ) — |
(m0 um i c ( 0 |
Соотношения для т2с |
(t), m2s (t) и mlc(t)mls(t) получаются из. |
ф-л (3.47) — (3.51).
Полученные в этом разделе формулы для математического ожи дания, авто- и взаимно корреляционных функций, средней мощ
ности и дисперсии определяют статистические характеристики пе реходного режима узкополосной линейной системы при воздейст вии на нее нестационарного случайного процесса.
Известно [24], что для стационарности импульсного случайного процесса необходимо, чтобы статистические характеристики слу чайных параметров импульсов (групп импульсов), зависили толь ко от взаимного расположения импульсов (групп импульсов), но не зависили бы от самых этих номеров. Точнее говоря, при вы полнении этого условия стационарными будут только апериодиче ские процессы, а процессы с детерминированным тактовым интер валом будут циклостационарными, согласно классификации при нятой в данной работе.
Однако, если длительность импульсов во много раз больше пе риода их повторения, т. е. периода изменения моментных функ ций рассматриваемого нестационарного процесса, то амплитуда этих изменений будет ничтожно мала по сравнению со средним временным значением. Это означает, что импульсный случайный процесс с детерминированным тактовым интервалом и сильно пе
123.
рекрывающимися импульсами практически стационарен (по край ней мере в широком смысле).
При воздействии |
импульсного случайного процесса на узкопо |
лосную систему, для |
которой A(tfn ‘С<йр=2л/7’р и Т р х Т , условия |
стационарности квадратурных составляющих отклика системы вы полняются для моментов времени отдаленных от момента вклю чения системы на время не меньше Тѵ. Тогда соотношения для оценки статистических характеристик процесса на выходе можно
вания импульсов процесса £(7).
При этом нестационарность случайного процесса g(0 оказыва ется обусловленной только переходным режимом самой линейной системы. Это видно, например, из соотношения для корреляцион ной функции процесса £(7). В предположении стационарности входного воздействия %(t)
С+т—V
<+т—о
І Шр ( 2 0 + 2 ) |
(3.67) |
] dvdz. |
о— V
Нетрудно показать [53], что второе слагаемое в выражении (3.67) в моменты времени отдаленные от момента включения не меньше, чем на Твкл^Тр, = л/Дсоп пренебрежимо мало, так как
внутренний интеграл его содержит быстроосциллирующую функ цию. Поэтому для достаточно больших значений времени авто- и
взаимно корреляционные функции процесса £(7) будут определять ся формулами:
со
(Д т) = B ,s(t, т) = I*F%(ю) Re [Г (і со, і)Т (і со, t + т)]
О
оо
cos (со —сор) т d со-----— ГF. (со) Jm [Т (і со, t) Т (і со, t -j- т)]
2п J 6
(3.68)
оо
О
124
sin (ü) — (Dp) X d (О~ 2n |
Kg (©) Jm [T (i со, t) T (i со, / + T)] X |
X cos (tu— (Op) xd (а. |
(3.69) |
Причем в совпадающие моменты времени взаимно корреляци онные функции равны нулю, т. е. значения квадратурных состав ляющих процесса t,(t) не коррелированьи
В установившемся режиме, соответствующем бесконечно отда ленному моменту включения системы, пределы интегрирования по времени (суммирования по N) распространяются до ± оо. В этом случае процесс t,(t) становится стационарным, и корреляционная функция его определяется известной формулой [53]:
В (т) == В- с (х) cos cüp X + |
с s(т) sin Юр т; |
(3.70) |
|
|
со |
|
|
где В' с(т) = —-— (Kg (со) I Т (і со)|2 cos (со — сор)тс(со, |
(3.71) |
||
|
2я J |
|
|
|
о |
|
|
В; с s (т) = ---- |
1Kg (со) IТ (і (o)|2sin (0)— (Op) xd со. |
(3.72) |
|
2л |
J |
|
|
|
о |
|
|
Нетрудно показать, что для дисперсии процесса t,(t) и его квадратурных составляющих в установившемся режиме справед ливо соотношение:
со |
|
£>с = DC= DS = Вс(0) — m2K = ~ J ^ H I M T2(®)d |
(3-73) |
о |
|
где Fjii (со) — непрерывная часть энергетического спектра про цесса £(7). Конкретные примеры воздействия импульсного случай ного процесса на узкополосную линейную систему будут рассмот рены в параграфе 3.5.
3.4, РАСПРЕДЕЛЕНИЕ ВЕРОЯТНОСТЕЙ МГНОВЕННЫХ ЗНАЧЕНИЙ РЕАКЦИИ ЛИНЕЙНОЙ СИСТЕМЫ
Рассмотрим воздействие на линейную систему с постоянными параметрами случайной последовательности б-импульсов вида
Оо |
|
1(0 = £ l „ b ( t - t n), |
(3.74) |
П = — со
оо
где £„■= J u(t) dt — случайная величина, пропорциональная пло-
—00
щади n-го импульса; б(t—tn) — смещенная на случайное время tn дельта-функция единичной площади.
125
Эта задача имеет самостоятельное прикладное значение, а кро ме того, к ней может быть сведено большинство задач о линейных преобразованиях импульсных случайных процессов с произволь ной, в том числе и случайной формой импульсов.
В самом деле, пусть форма импульса задается функцией u(t), а импульсная характеристика системы функцией h(t). Тогда, в соответствии с (3.5) реакция системы на воздействие импульса u(t) будет иметь вид
СО
иі (0 = j и (т) h (t — т) dr
— оо
Функции u(t) и Ui(t), в свою очередь, можно рассматривать как импульсные характеристики некоторых линейных систем, причем функция Ui ( t ) определяет характеристику приведенной линейной системы, образующейся при последовательном включении систем с импульсными характеристиками u(t) и h(t) (рис. 3.4). Переда-
t M(t)
Рис. 3.4. Реакция одиночного колебательного контура на .'воздействие случайной последовательности 6-импульсов:
а) воздействующий .сигнал; б) реакция на 1-й импульс; в) реакция «а 2-й им пульс; г) реакция на 3-й импульс; д) реакция на три импульса
точная функция приведенной системы, как легко показать, опре деляется перемножением передаточной функции исходной линейной системы и преобразованной по Лапласу функции u(t).
Очевидно, что предлагаемая методика легко распространяется и на задачи исследования преобразований импульсных случайных процессов в линейных системах со случайными параметрами.
Вероятностные характеристики процессов %(t) задаются, вооб ще говоря, бесконечномерным распределением совокупности слу
126
чайных величин іп(п= 0, ±1, ±2, ±3, . . ± оо). Однако, что бы избежать бесконечной совокупности случайных величин, вы делим на оси времени интервал (О, Т) настолько большой по сравнению с длительностью h(t), средним периодом следования импульсов Ти и временем корреляции тк, чтобы можно было «е принимать во внимание импульсов, частично попавших в этот ин
тервал.
.Пронумеруем N импульсов, попавших в интервал і(0, Т), числа ми натурального ряда, так что t\<h<tz< . . ,<tn- i < t N. Тогда ре акцию линейной системы на последовательность /Ѵ-импульсов в
момент времени t согласно .(3.3) |
можно записать в виде |
t(f) = Y l nh{t — tn). |
(3.75) |
п= 1 |
|
Исчерпывающей характеристикой процесса \(t) является 2N- мерная функция распределения вероятностей параметров gn и tn
импульсов (п = 1, 2, 3, . . ., N), попавших в интервал |
(О, Т): |
|
t (^Г’ |
XN, У\і У2)•• •» У |
(3.76) |
Ввиду того что t,(t) функционально связан с процессом £(7), пользуясь законами преобразования распределений вероятностей случайных величин при функциональных преобразованиях самих случайных величин, можно определить 2УѴ-мерную функцию рас пределения случайных параметров N импульсов процесса t,(t).
Из соотношения (3.75) следует, что процесс £|(0 ІВ каждый
фиксированный момент времени |
t* |
является случайной величиной, |
|
функционально связанной с |
и |
tn- Эта |
функциональная связь |
выражается суммой N случайных величин |
(п = 1, 2, . . ., N), |
||
в свою очередь, зависящих от |
и tn: |
|
|
£n(t*) = l nh(t* — tn) (п = 1, 2,..., N) |
|
(3.77) |
Для определения функции распределения случайной величины t,(t*) необходимо знать вероятность события В, состоящего в том, что %(t*) находится в интервале (z, z + dz):
Р {В) = Р {z < £ (**) < 2 + dz) =а>, с (а) dz. |
(3.78) |
Это означает, что вероятность события В связана с вероятно стью суммы случайных величин Z,n(t*). Следует заметить, что слу чайные величины t,n(t*) связаны со случайными моментами воз никновения tn через промежуточные функции r\n= h(t*-—tn) (как правило, нелинейные). Поэтому для того, чтобы определись функ
цию распределения wNg(Zi, . . |
zN) совокупности случайных вели |
|||
чин |
необходимо найти 2ІѴ-мернѵю функцию распределения со |
|||
вокупности случайных величин £„ и т^. |
вероятность Р(В), |
|||
Наиболее просто, по-видимому, определить |
||||
а следовательно, и плотность |
вероятностей w t |
(z) через |
характе |
|
ристическую функцию. Согласно [53] характеристическая |
функция |
127
одномерного |
распределения |
вероятностей случайной |
величины |
W * ) |
|
N |
|
|
|
|
|
о, |
00 |
1 ® |
|
ѳк (©)= J... |
Jo»„c (гъ...,Z„)e n=1 d z i . . . d z N , |
(3.79) |
—oo — oo
асама функция распределения w K (z) определяется как обратное преобразование Фурье от Ѳк (со):
СО
wu ^ |
==~ ^ |
J e , с(«в)е~ 1 “ zda» = |
|
|
|
|
|
|
|
—оо |
|
|
|
|
|
|
|
|
|
/ |
N |
\ |
|
|
® |
, |
|
*а ( |
2 |
z* ~z ) |
|
= |
- ^ - J J - .jw tfc (2х> |
zn)e |
^ ”= 1 |
' dzx...dzndu>. |
(3.80) |
||
|
— |
"“ A r |
|
|
|
|
|
или после изменения порядка интегрирования и введения много мерной 6 -функции
— оо
получим окончательное выражение для одномерной функции рас пределения вероятностей линейного преобразования импульсного случайного процесса в виде
|
00 00 |
N |
|
W 1 с |
J-J |
У x n h ( t — уп) — Z X |
|
|
|
П=1 |
|
X w 2N u ( X l ’- " ' X N ’ Л»-* Уы) П dxndyn. |
(3.81) |
||
|
|
/і=і |
|
Из соотношения (3.81) с очевидностью следует, что определение распределения вероятностей линейного преобразования импульсно го случайного процесса связано с громоздкими выкладками, в осо бенности сложна задача явного решения (3.81) при наличии ста тистических связей параметров импульсов на входе системы и про извольном виде импульсной характеристики h(t) линейной сис темы.
Решение поставленной задачи возможно только в некоторых простейших частных случаях или приближенными методами, ос нованными опять-таки на упрощающих предположениях относи тельно вероятностных характеристик входного случайного процес са и вида линейной системы.
128
Отметим, что соотношения (3.70) и (3.81) определяют условные характеристическую 'функцию и плотность вероятности линейного преобразования импульсного случайного процесса для произволь ного, но фиксированного числа импульсов N. Для определения бе зусловной плотности вероятностей необходимо дополнительное ус реднение по N, т. е.
Оо |
|
|
|
wK(z) = 2 |
Р{х = N}wit{z/N), |
(3.82) |
|
N = О |
|
|
|
где Р{%— N) — вероятность того, что случайная величина |
прини |
||
мает значение N. |
характеристическая функция определится |
как |
|
Безусловная |
|||
со |
р (у = |
N) Ѳк (со/ІѴ). |
|
Ѳк (со) = у |
(3.83) |
N = О
Соотношения (3.82) и (3.83) показывают, что определение рас пределения вероятностей линейных преобразований импульсных случайных процессов связано с громоздкими выкладками, а при наличии статистических взаимосвязей параметров импульсов яв ные выражения для w K(x) или Ѳг (со) могут быть найдены только
численными методами (82, 107]. Для некоторых частных случаев, при упрощающих предположениях относительно вероятностных ха рактеристик параметров импульсов и вида h(t) (импульсной ха рактеристики линейной системы), задача решена в [107].
В качестве примера рассмотрим случай статистической неза висимости случайных амплитуды \ п и моментов возникновения tn импульсов процесса (3.74) на входе линейной системы. Плотность вероятностей случайной амплитуды n-го импульса задается функ цией w ^n(x), а плотность вероятностей временного положения «-го
импульса функцией wUn(x), причем
оо |
|
J Wi in (х) dx = 1, |
(3.84) |
— со |
|
т |
|
\w\tn{x)dx= 1. |
(3.85) |
Ь |
|
Интервал времени (О, Т) в (3.83) по-прежнему выбирается та ким, чтобы он значительно превосходил длительность функции h(t) и средний интервал между импульсами (средний период их сле дования) .
Допущение о статической независимости параметров ,и tn означает, что 2/7-мерная функция распределения вероятностей этих параметров является произведением одномерных функций при лю бом п, т. е. (3.76) ів этом случае примет вид
N
w2N(xi,..., XN, уъ ..., уп) = П wn n (xn)wun (yn)-
п= 1
5 - 9 2 |
129 |