Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Коновалов Г.В. Импульсные случайные процессы в электросвязи

.pdf
Скачиваний:
8
Добавлен:
25.10.2023
Размер:
16.24 Mб
Скачать

Итак, соотношения (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

 

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

пО—оо

 

с функцией \ 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

 

 

 

 

 

 

 

 

 

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

 

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

Соседние файлы в папке книги из ГПНТБ