книги из ГПНТБ / Коновалов Г.В. Импульсные случайные процессы в электросвязи
.pdfИтак, соотношения (4.22), (4.23), совместно с (4.33) и (4.37), полностью определяют реакцию нелинейного элемента на воздей ствие случайной последовательности нѳперекірыівающихся им пульсов.
Формулам, (приведенным выше, (можно придать -несколько более общий вид, если увеличить до бесконечности интервал разложения в ряд Фурье. При этом ряд Фурье переходит в интеграл Фурье:
|
|
00 |
|
|
о о |
|
|
u ( f ) = - L - |
j > B(©)e,e<dö, Fu(a)= |
J u(t)e~i a t dt. |
|||||
|
— |
oo |
|
|
— со |
|
|
Формулы (4.5), определяющие ряды Фурье для различных сте |
|||||||
пеней функции u(t), примут вид: |
|
|
|||||
и2(0 = |
J 4 2)(ö)e‘“ ^co. |
|
|
|
|||
uP(t)= j / l P)(©)eit0'd®, |
|
|
(4.38) |
||||
где |
|
|
|
|
|
|
|
|
|
|
oo |
|
|
oo |
|
^u2)(ю) = ~ — [u2(t) é~ 1“ *d a |
= — |
( Vu(©— ©0 X |
|||||
|
2n |
J |
|
2л |
J |
|
|
X Fu(©x) e |
1“ *d ©! d ©. |
|
|
|
|||
|
|
|
|
i 1 |
I 1 * |
s |
|
|
|
|
|
|
|
oo |
oo |
4 rt(co) = |
^ |
|
ijuP(t)e-ia‘dt = ( 2 n r P+1 J... |
]>„(*>- |
|||
|
|
|
-F«(o>pi)П•,© |
|
|
— co |
— oo |
|
P — \ |
p — i |
|
|
|
||
— |
У] |
|
d ( i >r- |
|
(4.39) |
||
|
r=I |
|
r=[ |
|
|
|
Соотношения (4.39) определяют комплексную свертку спек тральной плотности функции u(t). Спектральная плотность выход
ного -импульса v(t) |
определится как преобразование Фурье: |
оо |
|
Л — |
І СО / JJ |
е |
at |
или как взвешенная сумма сверток спектральной плотности:
(4.40)
р=о
160
В заключение этого раздела отметим, что изложенная методика определения случайных амплитуды и формы импульсов нелиней ного преобразования последовательности неперекрывающихся импульсов достаточно легко может быть распространена на слу чайные 'последовательности перекрывающихся импульсов.
Преобразование функции распределения вероятностей длитель ности импульсов. В предыдущем параграфе было показано, что в результате прохождения импульсного случайного процесса через нелинейный элемент изменяется не только амплитуда, но и форма импульсов (а значит, и их спектральная плотность). Это означает, что нелинейное преобразование импульсного случайного процесса изменяет и длительность импульсов как один из параметров, опре деляющих их форму.
Длительность импульсов может быть определена множеством способов {25, 53, 86]. Практически наиболее широко используется определение длительности импульса как времени превышения функцией gW (t), описывающей форму импульса, некоторого за
данного уровня (например, уровня, равного половине амплитуды), т. е. за длительность импульса принимается разность моментов времени между первым пересечением заданного уровня функцией g(ft)(7) (или un(t), что фактически равноценно) с положительной производной и последним пересечением этого же уровня с отри цательной производной. Наряду с приведенным определением до вольно часто длительность импульса определяется как ширина ос нования прямоугольника с площадью, равной площади импульса и высотой, равной его амплитуде, т. е.
00
Возможны, очевидно, и другие способы определения длительности импульса.
Закон преобразования длительности импульсов при воздейст вии импульсного случайного процесса на нелинейный элемент с заданной характеристикой f(x) может, по-видимому, оказаться различным при различных определениях длительности импульсов. В общем случае даже при очень простой функциональной связи процесса на выходе с процессом на входе £,(0—Ш(0] преобразо вание длительности найти оказывается весьма сложно. За исклю чением, конечно, тривиальных случаев, например, £(t)— Kl(t), где К •— постоянная величина.
Можно предположить, что в зависимости от способа определе ния длительности должна изменяться и методика определения за кона распределения ее вероятностей.
Следует отметить, что полученные ранее соотношения, связыва ющие амплитуду и форму импульсов процесса на выходе нелиней ного элемента с амплитудой и формой импульсов входного процес са, позволяют трактовать задачу определения изменений длитель-
6—92 |
161 |
«ости и ее статистических характеристик как задачу о функции рас пределения вероятностей длительности импульсов импульсного случайного процесса с известными статистическими свойствами прочих параметров (например, амплитуды, формы, временного по ложения) импульсов.
Заметим, еще, что при любом определении длительность им пульса можно трактовать как разность моментов времени двух ха рактерных точек функции %n(t) (или ип(і)), описывающей форму импульсов, т. е. фактически длительность импульса есть время пре бывания функции ^n(t) над некоторым уровнем Хо. Таким обра зом, задачу об определении функции распределения вероятностей длительности импульсов можно рассматривать как задачу о рас пределении вероятностей длительности выбросов случайной функ ции in(t) над заданным уровнем х0, причем уровень может быть одинаковым для всех импульсов или разным. В последнем случае
необходимо решать задачу о выбросах случайной |
функции t](t) = |
— £,(t)—Xo(t) над нулевым уровнем, причем x0(t) |
может быть так |
же, как и \(t), случайным процессом, в общем случае статистиче ски или функционально связанными с t).
Исследование статистических характеристик выбросов вообще и их длительности в частности — одна из актуальных задач совре менной теории случайных лроцеосов [53, 86, 104]. Рассмотрим более подробно возможности и особенности применения методики опре деления функции распределения вероятностей длительности выб
росов, предложенной в [74] к определению |
плотности вероятности |
|||
длительности |
импульсов. |
так, что |
\ n ( t — і п ) пересекает ее |
|
Пусть |
функция A ( t ) задана |
|||
снизу вверх в некоторый момент времени t 0, т. е. |
||||
Т) (t) = l n ( t ~ tr) — A(t) = 0 при |
/ = и |
|
||
и |
|
|
|
|
dt— > 0 |
при |
t — to- |
|
(4.42) |
|
|
Здесь 1п(0 — функция, описывающая форму п-го импульса в слу чайной последовательности неперекрывающихся импульсов. Мо мент времени to, очевидно, расположен на оси времени вблизи мо мента t n , так что без потери общности можно принять t o = t n -
Учитывая, что по определению функции %n(t)-*-0 при |f|-»-oo в некоторый момент времени ^о+тп(т„^ 0 ) функция \ n ( t ) снова пе
ресечет уровень A ( t ) , |
но уже сверху вниз, так что |
Л (to + т„) = 0, d dt |
< 0 при t = to + V |
Применим к процессу г\(t) нелинейное преобразование вида
00
(4.43)
Я
162
Функция (4.43) может принимать значения 1,0, —1 в соответвии с изменениями г\(t):
I |
1 |
при |
г)(t) > |
0, |
|
v(t) = | |
0 |
при |
rj (/) = |
0, |
(4.44) |
( — 1 |
при |
т) (0 < |
О, |
|
Проинтегрируем v(t) в пределах от нуля до Т. Интервал интег рирования выбирается несколько большим, чем хп, т. е. таким, что бы в пределах отрезка времени 0, Т функция %n(t) успевала пере сечь уровень A(t) снизу вверх и сверху вниз. Тогда
U+т |
гп+хп |
т |
|
(4.45) |
|
j V(t)dt = |
j' |
<tf— |
J di = 2xn — T. |
|
|
Откуда, |
учитывая |
(4.43), получим |
интегральное |
уравнение |
|
Т_ |
Т оо |
|
|
|
|
|
sinx[g„ (О — 4 (0] d* = |
/(£„, А), |
(4.46) |
||
2 |
2л пО—оо |
|
с функцией \ n(t). |
||
связывающее длительность импульса |
|||||
Используя |
известное [53] правило |
преобразования распределе |
ния вероятностей при функциональном преобразовании случайных
величин, для |
|
рассматриваемого случая, считая обратную функцию |
|||||||
| = ф(т) однозначной, запишем |
|
|
|
|
|||||
Wu (х) = |
дх |
»16 [У = ф (*)], |
|
|
(4.47) |
||||
|
|
|
|
|
|
|
|
||
|
¥ |
|
|
|
|
|
|
|
|
где у = ср(х) — функция, обратная f(ln, А). |
пределах |
||||||||
Бела допустить, |
что функция |
А) определена в |
|||||||
изменения |
| п(0> |
а |
функция |
ХѴ[^ |
непрерывна на |
площади |
|||
|
|
|
|
|
|
X |
|
производ |
|
— о о < х < |
оо, |
|
|
|
и имеет в этой области частную |
||||
ную, то |
|
|
|
|
|
|
|
|
|
дх |
I |
со |
|
sin [у — A] z |
|
|
|
|
|
|
|
д |
|
= 1. |
|
(4.48) |
|||
ду |
О—00 ду |
|
|
|
|||||
|
|
|
|
|
|
||||
и последнее |
|
равенство становится |
очевидным, если учесть, |
что |
|||||
т](Т) = £ті(/0—A(t) |
на интервале (О, |
Т) обязательно проходит |
че |
||||||
рез нуль по условию. С учетом |
(4.48) |
предположение об однознач |
ности функции, обратной |
(4.46), не существенно [74]. Таким обра |
|
зом, |
|
|
wix (х) = |
[у = Ф (*)]. |
(4.49) |
Моменты распределения (4.49) определяются по известной ме
тодике [53].
6* |
163 |
|
Учитывая специфические свойства импульсных случайных про цессов и, в частности, постепенное убывание функции \ n(t), опи сывающей форму п-го импульса, можно несколько иначе опреде лить функцию распределения вероятностей длительности импульса.
Для этого функцию vn(t), представляющую собой прямоуголь ный импульс единичной амплитуды и случайной длительности хп, представим рядом Фурье, «а интервале (0. Г0), причем Т0^ 2 х п-
Тогда
CO
ѵп it) = It |
|
|
— |
V 4 |
sin (2/c -(- 1) (m0 t -f- Qn) |
(4.50) |
J |
X |
IT |
Z j |
(2 к + 1) |
||
— |
00 |
|
|
|
|
|
где coo —2я/7 о.
Случайная фаза Ѳп свидетельствует о случайном характере из менения длительности импульса vn(t) при изменениях параметров
п-го импульса преобразуемого |
по закону |
(4.43) |
процесса |
%(t). |
Пределы изменения случайной |
величины |
Ѳп от |
0 до + я/2 |
при |
г|(7)>0, т. е. при положительных значениях амплитуды преобра
зуемаго импульса I« (ft) и |
при т]'(()<0 при отрицатель |
|
ных значениях амплитуды \ n(t). |
Так как |
случайные события |
О и 5n(t)<0 несовместимы, |
а закон |
преобразования (4.43) |
не меняется при изменении полярности преобразуемого процесса
ц(і), то
|
twin (У) |
dy |
|
|
0 < Ѳ Л< |
It |
|
|
||
|
dz |
|
|
2 |
|
(4.51) |
||||
wie (z) = |
|
|
dzdy |
|
|
|
|
|
||
|
twin (— У) |
|
-f<* < 0. |
фазы Ѳп и |
||||||
Покажем, |
что |
функциональная |
связь случайной |
|||||||
г}(t)—%n(t)—A(t) выражается формулой |
|
|||||||||
Ѳ„ = |
arc sin г) (t)lp, |
|
|
|
|
|
|
(4.52) |
||
где |
р — нормирующий |
множитель, |
выбираемый так, |
чтобы функ |
||||||
ция, обратная |
(4.52), не превышала |
1, т. е. |
|
|||||||
sin Ѳ„ = |т](t)\/p < |
1. |
|
|
|
d ц |
|
|
(4.53) |
||
Определяя |
производную |
модуль которой равен |
||||||||
dy |
|
|
|
|
|
|
7ѳ |
|
|
|
pcosz |
|
|
|
|
|
|
|
|
(4.54) |
|
dz |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
и подставляя |
(4.53), |
(4.54) в |
(4.51) с учетом — psinO =rpsin(—Ѳ) |
|||||||
для закона распределения вероятностей случайной величины |
||||||||||
получим |
|
|
|
|
|
|
|
|
|
|
|
pcosza>i„(psinz),---- —< z < — |
|
||||||||
Wie (2) |
|
|
|
|
|
2 |
|
2 |
(4.55) |
|
- |
|
It |
|
|
Ti |
|
|
|||
|
А |
|
|
|
|
|
|
|||
|
О’ < |
~ |
Т |
’ |
г > Т - |
|
|
|
164
С другой стороны, функциональная зависимость между Ѳи и ис комой случайной величиной хп имеет также довольно простой вид.
В самом деле, учитывая разложение vn(t) в ряд (4.50) и объе диняя (4.50) с (4.45), найдем
|
|
|
|
sin (2/с |
1) (Ü)Qt -f- 0п) dt. |
(4.56) |
|
|
|
|
|
(2 к + 1) |
|
После интегрирования получим |
|
|||||
Т_ _ |
1 |
£ |
4 cos (2к + |
1) (to0 t + Ѳп) |
|
|
2 |
2я |
|
coo (2/с 4-1) |
|
||
I |
|
к=0 |
|
|
T |
|
|
|
|
|
(4.57) |
||
2(о0 [|соот 4- Ѳ„ '|®л|] + |
|
|||||
|
|
|||||
Окончательно длительность импульса |
|
|||||
хп = Т — Ѳл/©0. |
|
|
|
(4.58) |
||
Соответственно функция распределения случайных величин т„ |
||||||
и Ѳп связаны зависимостью |
|
|
||||
Wu(x) = |
ö)e а>іѳ [©о (Т — X)], |
|
(4.59) |
|||
а с учетом ф-лы (4.55) |
можно установить непосредственную связь |
|||||
законов распределения |
случайных величин хп и r\(t): |
|
||||
WIT (X) = |
©О р COS ©о (Т --Х) Wir) [р sin ©о (Т — X)], х > 0, |
(4.60) |
||||
0, |
х < 0. |
|
|
|||
|
|
|
|
Для простоты функцию распределения (4.60) целесообразно представить в виде ряда по ортогональным полиномам. Причем, учитывая, что случайная величина тп может принимать только по ложительные значения, в качестве ортогональных полиномов ис пользуем полиномы Лягерра, интервал ортогональности которых охватывает всю положительную полуось, т. е. (0, сю). Тогда
WIT (X) = £ ске~ * ха й а) (X) = |
ха е- * £ скÜ a) (х). |
(4.61) |
||||
к=0 |
|
|
|
к=0 |
|
|
L™ (х) = і |
еД - а ^ - |
(е |
хX K+“) — обобщенный полином Лягерра. |
|||
«! |
dxK |
|
|
|
|
|
Учитывая условия |
ортогональности полиномов |
Лягерра (53], |
||||
|
|
|
Г (я 4~ к -f- 1) |
при к = /, |
|
|
I е - ххай а) (X) L|a) (X) |
|
|
к! |
(4.62) |
||
|
|
|
||||
о |
|
. |
0 |
при к ф і , |
|
|
где Г(Х) — гамма-функция [112].
165
Для коэффициентов ск разложения (4.61) получим
|
к\ |
со |
|
Ск — |
J■ L*a> (х) Wix (X) dx |
||
Г (а + к + 1) |
(4.63)
РП«
Учитывая формулы для первых полиномов Лягерра, приведен ные в [53] из (4.63) для низших коэффициентов разложения wu (X), получим
Со = ß r( a + l) ’ |
Cl = |
ß r(a + 2) |
j |
+ a ---- jj-'j Wir (V) dv |
||||
|
|
|
|
|
о |
oo |
|
|
ßT(aг + 2 г (1 + a |
ß |
) ’ |
|
|
1) (a + 2) ■ |
|||
„ |
ßГ (aЬгі (a |
|||||||
|
|
|
,llx |
1 |
_ |
|
|
|
2 (« + |
2) V |
v' |
|
|
|
1 |
(a + 1) (a + |
2) — |
2(a + 2) |
^ w |
i x (v)dv = |
|
ßT(a + 3) |
||||
2t |
|
|
|
|
|
(4.64) |
||
ß |
mix + |
J |
|
|
|
|
|
|
|
|
|
|
|
|
где а и р - произвольные постоянные, поэтому их можно выбрать так, чтобы сі = С2—0. Тогда
a = т у о \ — 1
(4.65)
ß = о\ІтХх
где ти и о* — математическое ожидание и дисперсия случайной
длительности тп.
В таком случае функция распределения случайной величины г« (4.63) примет вид:
|
|
|
_ * Щ х |
Шіт (х) = |
Щх |
*Щх у |
„2 |
|
е |
+ |
(4.66)
r = 3
Моменты функции распределения (4.66) можно определить с учетом (4.60) в виде
т к X— I **®в р cos ©о (Т —*х) шіті [р sin ©о (Т — x)^dx |
(4.67) |
166
или, переходя в |
(4.67), |
к новой переменной интегрирования z = |
|
= о)о(Т—х) в виде |
|
|
|
я/2 |
|
|
|
mKX= j (w0T — zj^pcosz^infpsinzj^z, |
(4.67а) |
||
— я / 2 |
|
|
|
где w ІГ)(у) — w ^ |
(у-\-А) |
— связана с совместным распределением |
|
вероятностей случайных |
параметров |
импульсов преобразуемого |
|
процесса %(t) соотношениями, приведенными в параграфе 3.4. |
|||
Заметим, что уровень А, по которому определяется длитель |
ность импульсов, может быть задан не только как постоянная ве личина, но и как произвольная (в том числе случайная) функция времени. Из предыдущего рассмотрения .видно, что ограничения, накладываемые на r j ( 7 ) = g(7)—А состоят в том, что г\(і) должна быть дифференцируемой и иметь нуль на интервале (0, Т).
Если уровень A(t) является случайной функцией, допуская, что известна совместная плотность распределения вероятностей Шг(ё, А) для распределения вероятностей т] имеем
Win(y) = ^ w 2(v,y+ V)dv. |
(4.68) |
А |
|
Таким образом, задача определения закона распределения ве роятностей длительности импульсов сводится к отысканию функ ции, обратной (4.46), или, если воспользоваться промежуточным случайным параметром Ѳп и представлением и.'1х (х) в виде ряда
(4.66), к определению моментов w lx(x) по ф-лам (4.67) или (4.67а).
4.3. СПЕКТРАЛЬНО-КОРРЕЛЯЦИОННЫЕ ХАРАКТЕРИСТИКИ НЕЛИНЕЙНЫХ ПРЕОБРАЗОВАНИЙ ИМПУЛЬСНЫХ СЛУЧАЙНЫХ ПРОЦЕССОВ
Энергетические характеристики процесса на выходе нелинейно го элемента. Общее определение энергетических характеристик им пульсного случайного процесса, данное в гл. 2, позволяет для слу чая воздействия на безынерционный нелинейный элемент случай ной последовательности импульсов записать соотношения: .
|
СО 00 |
|
(4.69) |
|
n«=—о» j——00 |
|
|
|
|
|
|
* (ш) = |
P £>»./<“> |
(4.70) |
|
|
|||
где |
n=-N/=—N |
|
|
/ |
t + x - t * |
|
|
|
|
||
В,, n, j (it, t ) = m x J £< *> |
п( ' Ч т 1~г ) и ( |
T)(ft) i C |
(4.71) |
|
|
/ C |
|
167
к , п, і (со) = тх { £<,*>£<*> ^ к{*]к\ёѵ((>>т£\)Х |
|
|||
_______( Цк)_______ Цк)Л |
(4.72) |
|||
Xgv[(*r}fye~ ' aK " С |
' |
*’ |
|
|
g v (со т<*>.) = |
j V (X) е |
dx. |
|
(4.73) |
|
о |
|
|
|
В (4.69) |
и (4.70) ß g |
(7, |
т), ß* (т), и ß g (со) ■— |
соответственно |
корреляционная, усредненная корреляционная функции и энерге тический спектр процесса t,(t) н-а выходе нелинейного элемента, а Цк\ Д.(*) и /)*)- амплитуда, длительность и временное положение
г-го импульса k-й реализации процесса Z,(t).
Учитывая, что форма (функция v(t)] и статистические характе ристики случайных параметров импульсов процесса tj(t) при -воз действии на нелинейный элемент случайной последовательности неперекрывающихся импульсов могут быть найдены по методике предыдущего параграфа, и полагая для простоты разнородные па
раметры импульсов |
независимыми, выражения |
(4.71) и (4.72) |
||||||
представим в виде |
|
|
|
|
|
|
||
|
>= M |
W |
) |
) |
' ф " Z ; W 2 т |
, ( Z „ , Zj)X |
|
|
|
|
|
|
о |
о |
|
|
|
оо |
00 |
|
|
|
|
|
|
|
ХI JУ(Ѳ" |
(Ѳ/ + |
~ |
) |
W2t |
— |
І — Ѳ/Zy)X |
|
|
— оо — оо |
|
|
|
|
|
|
|
|
X d Ѳ„ d Ѳ/dzn dzj. |
|
|
|
— и , |
n, j) X |
(4.74) |
||
n> j (®) = mt { C^*' Cl*»} Ѳ2, c (со, |
|
|||||||
X j |
Jz„ Zj g v (CO Zn) g |
v (со Zj) w2 X |
c ( Z „ , |
Z /) dzn dZj. |
(4.75 |
|||
0 |
0 |
|
|
|
|
|
|
} |
Причем при ограничениях характеристики y = f(x ) нелинейного элемента, принятых в параграфе 4.2, для первого смешанного мо мента случайных амплитуд п-го и /-го импульсов процесса t,(t) на основании ф-лы (4.25) получим
тг { £<*>} = 2 j \ b Kb i m i { f C T [&<*>]'}. (4.76)
к=0 l—О
Двумерная функция распределения вероятностей длительности ц-го и /-го импульсов процесса £(7) определится в рассматривае
мом |
случае формулой |
|
|
w2Xс (гп, zj) = c\w2xr. (Zjclt Zjlcx), |
(4.77) |
||
где |
Ci — коэффициент, |
зависящий |
от формы импульсов процесса |
£('0 |
(т. е. вида функции u(t)] характеристики (нелинейнаго элемен |
||
та и способов задания |
длительности импульса (О ^щ ^Л ). |
168
Относительно статистических характеристик временного поло жения импульсов процесса t,(t) с учетом безынерционности нели нейного элемента и отсутствия перекрытия импульсов воздейст вующего процесса \(t) можно считать, что
02 t с(®, —©, п, /) — Mi |
02с £ (со —со,«,/). |
(4.78) |
|
Форма |
импульсов, |
зависящая от функции v(t), может |
быть |
найдена из |
соотношения (4.23) или одним из указанных в пара |
графе 4.2 методов.
Таким образом, энергетические характеристики безынерционно го нелинейного преобразования импульсных случайных процессов типа различных формирований неперекрывающихся импульсов могут быть достаточно просто1определены, если для всех случай ных параметров (включая и амплитуду) импульсов известны дву мерные функции распределения вероятностей их возможных зна чений.
В более общем случае, когда длительность импульсов преоб разуемого процесса \(t) неограничена или импульсы, возникаю щие в различные моменты времени, могут накладываться друг на друга для определения энергетических характеристик, целесооб разно воспользоваться представлением усеченных реализаций им пульсного случайного процесса в виде ряда (4.8), а не (4.4) или (4.5), как это делалось в предыдущем параграфе для случая им пульсов ограниченной длительности. При этом удобнее найти сна чала корреляционную функцию ßg (t, т) процесса t,(t) на выходе
нелинейного элемента, а затем, преобразуя по Фурье усредненную корреляционную функцию В* (т), определить и энергетический
спектр F (со) процесса £(7). Можно показать это на примере инерционного нелинейного преобразования.
В качестве примера инерционного нелинейного преобразования импульсного случайного процесса рассмотрим воздействие случайной последовательности им пульсов на типовое радиотехническое звено (рис. 4.4).
Преобразование импульсного случайного процесса в типовом радиотехниче ском звене [53] должно, по-івадимому, привести к некоторому сглаживанию флук туаций мгновенного значения процесса (см. параграфы 3.3 и 3.5), т. е. к уве личению длительности импульсов на выходе преобразующей системы, так что процесс на выходе инерционной нелинейной системы можно рассматривать как импульсный с перекрытием .импульсов.
Для простоты положим, что рассматриваемая нелинейная система полностью разделима и линейные фильтры на нелинейные системы рис. 4.4 можно опреде лять независимо от других.
Пусть на вход линейного фильтра I типового радиотехнического звена воз действует случайная последовательность импульсов ограниченной длительности, k-я реализация которой
(4.79)
169