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

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

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

 

Г/2

 

f

j и Л О е - " - ' dt =

 

-Г/2

 

где

Г/2

 

&,1 о

 

j, |*

{t) dt,

 

Г/2

 

Г/2

ип (t) cos к (£>0tdt,

 

I*

 

—Г/2

 

_21_ (вда — І^к)

При

К > О,

1 К * + і Ьпк)

 

(4.6)

при

к < О,

Г/2

bnK = ~Y I* ы„ (t) sin ,<cco0 1 dt.

—Г/2

В случае неограниченной длительности импульсов, когда k-я реализация случайной последовательности импульсов имеет вид

СО

l W (t)= У tnk){ t - t {nk)),

п=я— 0 0

а относительно функций %{k)(t) известно только, что они ограниче­

ны и стремятся к нулю при |-/|-ч-оо, целесообразнее .'представлять в виде ряда не отдельные импульсы &-й реализации процесса, а достаточно длинные отрезки каждой реализации, на интервале разложения (—Г/2, Г/2) совпадающие с соответствующей реализа­

цией рассматриваемого импульсного случайного процесса

(см. па­

раграф 2 .1 ).

 

Тогда, считая, что импульсный случайный процесс эргодический

и характеризуется условием

 

||ß * (t)|d t< оо,

(4.7)

-00

 

усеченную k-ю реализацию этого процесса можно представить ря­ дом Фурье на интервале усечения (—Г/2, Г/2) [127], т. е.

№ { t ) =

 

Ч

E(fe)

1 К(ÖJ t

 

У

£ Г е

 

(4.8)

°Т

 

 

 

 

где

Г/2

 

 

 

 

 

 

 

—Г/2

 

 

 

s

 

 

Г/2

 

Л/

 

і К (00 t W

 

 

 

Ф j' i S (< ) (''“

*

n=г=—ЛN

 

 

 

 

-Г /2

 

ЛГ

 

— i к (Oo

 

 

N

 

n .(*)

 

~

 

 

 

(4.9)

= l i e

 

gx0 (rt Cöo)

n——NS

 

 

 

 

150

и Т= 2 п / ш о — период функции l(Tk) (/), тождественно равной

(t)

на интервале (—Т/2,

Т/2)

гари .веек значениях к.

 

Коэффициенты

в

разложении (4.8) являются случайными

комплексными величинами, так как функции b , T(k ) ( t ) различны при разных к в соответствии со сменной реализацией рассматриваемо­

го случайного процесса.

а для коэффициентов

раз-

Очевидно, что Ш ) = 1іт

Т-+00

 

 

ложения (4.8) справедливы соотношения (см. [127]):

 

1 іш т 1 { і к * ) } = 0 , 1 і т т 1 { [ ^ * )]2}

= 0

(4.10)

Т-*со Т-+со

В соответствии с соотношениями (4.4), (4.5) и (4.8) результат нелинейного преобразования импульсного случайного процесса мо­ жет быть представлен в виде детерминированных функций време­ ни и некоторого числа случайных параметров, определяемых ко­ эффициентами разложения преобразуемого процесса в функцио­ нальный ряд. Для каждой конкретной реализации преобразуемо­ го импульсного случайного процесса набор коэффициентов разло­ жения представляет собой набор детерминированных величин.

Представление импульсного случайного процесса в виде ряда упрощает решения как общей задачи определения функции рас­ пределения вероятностей нелинейных безынерционных преобразо­ ваний, так и частных задач вычисления моментных функций и энер­ гетических характеристик таких преобразований.

Следует заметить, что определение спектрально-корреляционных характеристик нелинейных преобразований случайных процессов задача, как правило, более сложная, чем определение функций раопределения.

В заключение этого параграфа отметим, что в практике проек­ тирования и исследования систем обмена информацией возникают задачи анализа нелинейных преобразований не всего импульсного процесса в целом, а только отдельных параметров импульсов (ам­ плитуды, длительности, временного положения и т. п.). С учетом возможности разделения операторов таких преобразований на не­ линейные безынерционные и инерционные линейные задача ис­ следования статистических характеристик импульсных случайных процессов на выходе таких преобразователей еще более упро­ щается.

Отмеченные особенности исследования нелинейных преобразо­ ваний импульсных случайных процессов подробнее будут рассмот­ рены в последующих разделах этой и других глав.

4.2. ВЕРОЯТНОСТНЫЕ ХАРАКТЕРИСТИКИ ИМПУЛЬСНОГО СЛУЧАЙНОГО ПРОЦЕССА НА ВЫХОДЕ НЕЛИНЕЙНОГО ЭЛЕМЕНТА

Преобразование амплитуды и формы импульсов. Рассмотрим результаты воздействия импульсного случайного процесса, k-я реа­ лизация которого имеет вид

151

|<*>(0 = | ]

lnk) и

t - t p \

< h) J

П — — СО

 

на безынерционный нелинейный элемент с характеристикой у — =\f(x). Вероятностные характеристики |(7) заданы функциями рас­ пределения вероятностей случайных параметров In, хп и tn, а от­ носительно функции f(x) положим, что она ограничена и непрерыв­ на со всеми своими производными для каждого значения х, т. е. при любом значении х = Хо представима рядом Тейлора:

/(*) = / Ы Iх ■<■») Іт (х)\*-*,+ -і- (X— *>) /|!| W l~«. +

или

 

 

 

f(x) = Y^bpix — x0)p,

 

 

(4.11)

p= 0

 

 

 

где

 

 

 

/(P)W|л:=л-0

_ 1

dp/W

(4.12)

P!

(d*)p

Pi

 

Значение аргумента f(x), при котором строится ряд (4.11), в электросвязи называемое рабочей точкой, выбирается, исходя из требуемого режима работы нелинейного элемента (выпрямление, детектирование, усиление, ограничение и т. л.) в соответствии со статистическими характеристиками преобразуемого процесса.

Коэффициенты разложения Ьр в ряд (4.11) при заданных ха­ рактеристике преобразования y=f(x) и положении рабочей точки постоянны. Очевидно, что разным характеристикам f(x) соответ­ ствуют разные наборы коэффициентов {Ьр}, р = 0, 1 , 2 , ... даже при одном и том же значении аргумента хо, и, наоборот, при за­ данной характеристике нелинейного элемента y=f(x) величины f(-^(x)/x0 и, соответственно Ьр, определяются выбором рабочей точ­

ки х0 (точки разложения f(x) в ряід

Тейлора). Пусть, например,

характеристика нелинейного элемента

описывается функцией

X

 

У — f(x) — А j е~ 2*/2 dz, z — (х а)/о.

(4.13)

На рис. 4.5 изображен график функции (4.13) для нескольких зна­ чений параметра о. Как известно [53, 76], функция (4.13) непре­ рывна вместе со всеми своими производными на интервале —о о < < х < о о . Ряд Тейлора для этой функции имеет вид f(x) =

152

—оэ

характеристики

- a \ 2

S'bp(xxo)p,где коэффициенты fr — j L .

£o

P Pl dxp er У 2л

x = x a

и ершима ют разные значения ів разных точках оси абсцисс OM. рис. 4.5):

Очевидно, что положение рабочей точки влияет также на ха­ рактер аппроксимации іf(x), а значит и на характер преобразова­ ний (рис. 4.5).

По отношению к преобразуемому сигналу \ ( t ) рабочая точка обычно выбирается так, чтобы отклонения §(£) от х0 имели как можно меньшую вероятность, т. е.

Р { [ | ( 0 — *о] = е) = min

Этому условию, как известно, удовлетворяет математическое ожидание случайной функции ті{І^(і)}. Практически, рабочая

4) Предположение о непрерывности y = f ( x ) и всех ее производных в любой точке оси X не онижает общности дальнейшего анализа, так как полученные при этом результаты могут быть распространены и на функции, имеющие конеч­ ное число разрывов [112]. Разрыв непрерывности функции f ( x ) или ее произ­

водных при этом могут быть устранены путем введения обобщенных функ­ ций |[:127].

153

точка х0 совмещается с постоянной составляющей преобразуемого процесса:

 

Г/2

____

 

х0 = 1іш іг

f

l {k)(t)dt = l {k){t)

(4.14)

Г-оо Т

J

 

 

 

-Г /2

 

 

Если процесс эргодический, то ^ k)( t ) — rrii{l,(-h'>(t)}, а положение ра­ бочей точки определяемое соотношением (4.14) неизменно во вре­ мени. Для нестационарных или неэргодических процессов вели­ чина х0 также может выбираться равной математическому ожида­ нию нестационарного процесса:

x0(t) = m1{ l ik) (О) = ml%{t)

(4.15)

или, в соответствии с (4.14), равной временному среднему неэргодического процесса. В этом случае положение рабочей точки будет разным в различные моменты времени для нестационарного, но эргодического процесса или разным для различных реализаций неэргодического процесса, и х0, выбранном в соответствии с (4.14).

Соответственно изменениям х0 будут меняться и значения про­ изводных функций tf(x) и коэффициенты Ьр, см. ф-лу (4.12), также будут детерминированными или даже случайными функциями вре­ мени.

Если Хо не зависит от времени, ход процесса х0 на выходе не­ линейного элемента определяется множителями вида хо), где Xo=mi{gk)(t)}. Без потери общности примем Хо = 0. Тогда k-ю реа­ лизацию процесса на выходе нелинейного элемента можно пред­ ставить в виде:

г

ч о =

£>(*)]'

 

Р- 0

 

 

0 Lfi=Еoo

(4.16)

” V

r(‘)

Заметим, что если импульсы преобразуемого процесса не пе­ рекрываются, то

0 при п Ф /,

(4.17)

Поэтому ф-лу (4.16) можно переписать в виде

5(ft) (0 = 5 ]б р

(4.18)

р=0 /

 

154

или, меняя очередность суммирования по индексам р и п, обозна­ чая

 

 

Р

 

 

(4.19)

запишем

 

 

l ik)(t)=

&k)

(4.20)

 

П= — со

■VI Г

В соотношениях

(4.19), (4.20) учтено, что в нелинейных безы­

нерционных системах реакция в любой момент времени определя­ ется только значением воздействующего процесса в тот же момент времени t.

Таким образом, импульсный характер воздействия сохраняет­ ся и інѳперѳкірыіваюіщ'иеся иміпульсы на входе остаются также інѳперекрывающимися на выходе.

Импульсы процесса £(7) на выходе нелинейного элемента опи­ сываются взвешенной суммой соответствующих функций вида [u(t)]p, взятых с весовыми коэффициентами Ьр (4.19). Каждый п-к импульс к-к реализации процесса t,(t) можно представить как ре­

зультат умножения

детерминированной

функции v(t),

которая

тождественно равна нулю вне интервала

на величину І,\к)

сдвига по оси

времени на величину

и изменения

масштаба

времени в т ^ !

раз.

 

 

 

Таким образом,

 

 

 

£ (fe)( 0 =

 

t — t(k)

 

(4.21)

l (nk ) V

%n n C

 

 

 

T(*)

 

 

 

 

C

 

 

Учитывая, что функция u(t), определяющая форму импульса на входе нелинейного элемента, в какой-то момент времени t*

xhk)) принимает максимальное значение, равное единице, из ф-лы (4.19) найдем, что

f(k)

(4.22)

ЪП

р = 0

1> *1пк) — 0> х„к) = 1>

С другой стороны, если принять, что

то

 

оо

 

* ( 0 = Х М “ (0 ]Р.

(4-23)

р = 0

 

Вероятностные характеристики амплитуды импульсов процесса на выходе нелинейного элемента в этом случае легко определяют-

155

ся, если известны вероятностные характеристики амплитуды импульсов процесса на входе. Так, моменты случайной величины £„ на основании (4.22) равны:

ПІг - ' ; Ш. ([ І “ ]'} = /« і р Д Й * '] '

(4.24)

р=0

 

Видно, что для определения момента г-го порядка случайной величины £„ необходима одномерная функция распределения ве­ роятностей (по меньшей мере моментов до рг-го порядка).

Если известна двумерная функция распределения вероятностей случайных амплитуд п-іго и /-го импульсов процесса £(£), то мо­ жно определить двумерные моменты случайных амплитуд п-го и /-но импульсов процесса £(f) на выходе нелинейного элемента:

mrq { Ük) l)

Й‘7

у ь Д б Г У

(4.25)

LP=O

J

L i=o

 

В частности, среднее значение, дисперсия и коэффициент кор­

реляции амплитуд импульсов процесса t,(t) согласно

(4.24) и

(4.25) равны:

 

 

 

оо

 

 

(4.26)

а-. = ті { Ük)} = У,Ьртр{ lhk)},

 

р

 

 

 

{ Ң*1) - <4 _ £ £

Ь,.Ь, {т„ { 6!»} -

 

р= 0 1=0

 

 

(4.27)

 

 

 

(4.28)

Ѵ м 2{ ^>}М 2 {£<*>}

где

(4.29)

р= 0/=0

Таким образом, для определения энергетического спектра ус­ редненной корреляционной функции нелинейного преобразования импульсов случайного процесса необходимо не только среднее зна­ чение, дисперсия и коэффициент корреляции амплитуд импульсов на входе, но и моменты более высоких порядков.

Очевидно, что, зная моменты функции распределения вероят­ ностей случайной амплитуды, можно построить по известным пра­ вилам [53, 81] ряд, сходящийся к самой функции распределения или к характеристической функции.

Рассмотрим теперь способы определения формы импульсов на выходе нелинейного элемента. Как было показано выше, форма

156

импульсов на выходе определяется как взвешенная сумма слагае­ мых вида {u(t)]P. Для нахождения функций [u(t)]p можно предста­ вить u(t) в виде экспоненциального ряда (4.4). Это позволяет представить {u(t)f в виде

Іи ]2 = « 2 (0 =

V £

АКі Ак, А****) 1

 

 

 

к,=0 к2=0

 

 

Тогда более высокие степени функции u(t)

можно записать как

“•<*) =к ,£— 0

к £, — 0

к,=ОІ !• к А ‘ к А

1 ) , ке(* , к ' і С "

 

“ Р (0 = П

X

Л ^ е х р ^ О .

(4-30)

К=1

кх =0

 

 

Для исследований функций up(t), представленных рядом (4.30), можно использовать преобразование Лапласа [25]. В самом деле, обозначая преобразование Лапласа от f(t) как

F(s) = ] n f ) e - Stdt

о

в соответствии с (4.29) и (4 30), найдем

СО

со

 

 

 

(4.31)

 

ЕЕ

А* А*

 

Л « Л *

 

[»-(**.+ 8*,)]

 

к,=0 ,=0

 

Заметим,

что представления (4.30), (4.31) и свойства преоб­

разования Лапласа [25] позволяют анализировать не только ре­ зультаты воздействия импульсных случайных процессов на безы­ нерционной, но и на многие виды инерционных нелинейных систем [104, 127].

В виде ряда Фурье может быть представлена любая функция f(t), удовлетворяющая на промежутке OsS^s^T условиям Дирих­ ле. Учитывая определение импульсного случайного процесса и ус­ ловие ограниченности характеристики нелинейного элемента f(x), функция v(t) аналогично (4.4) может быть представлена рядом

Фурье на том же, что и u(t)

интервале Т, причем Т ^ х п, т. е.

v(t) = у Акеіка°(\

(4.32)

К— — СО

 

157

где

Г/2

A . - T

f » ( O e - " - ' .

 

(4.33)

 

*

J

 

 

 

 

 

—7/2

 

 

 

 

 

Рассмотрим связь

коэффициентов разложения в ряд

Фурье

функции ц(Т) и и(Т).

 

можно записать, что

 

Л

На основании (4.33) и (4.23)

 

оо

 

 

 

оо

 

=

*

Г/2up(Oe-‘““»<Ä=5]6^Ä/,>

(4.34)

 

ЛИ*

J

 

ши

 

 

р = 0

 

- Г /2

 

р= 0

 

где

Г/2

 

 

 

 

 

 

 

 

 

 

Л р

= ^ г

f

ttp ( O e " “

“ , , Ä .

 

(4.35)

 

*

и

 

 

 

 

—Г/2

Для большей наглядности заменим непосредственное определе­ ние величин ЛкР, по ф-ле (4.34) возведением в р-ю степень раз­ ложения (4.5). Тогда

u2(t)

V

скеіка°*

 

ZJ

K

 

/C=—CO

u3(t)

I

i « (Oo*

 

2

oo

со

 

І (*,+*„) (Ooi

 

И

Ъ

Ск'Ck*

 

 

3

/ C j = — oo K 2= — oo

 

 

 

oo

oo

У

cKl сКг сКзе

' («і4«2+Кз) ®о<

 

У

У

 

/С- — — 00 Ко= — 00 / С , = — 00

(4.36)

р

со

оо

up(t) = X

= у

у ... у , KlcK2...S x

Xexp I і £ /c r coo4

г = 1

Соотношения (4.36) позволяют сделать важный для практики вывод, а именно, в результате нелинейного преобразования в соста­ ве сигнала на выходе появляются новые гармонические составляю­ щие (так называемые комбинационные частоты). Пусть, например, в разложении (4.5) только со, си Съ С - с_ 2 отличны от нуля. Тогда, выбирая рабочую точку х0= с0 и считая, что Сі= с_і, с2= с_г, пред­ ставим (4.5) в виде

и (0 — с0 + 2 Сі cos 0 t 2 c2 cos 2 co01,

а возведение {w(Y)—x0] в квадрат дает

и2(t) =

4 cf cos2 (i)01 + 4 cf cos2 (j)01-f 8cx c2 cos a>0t cos 2 oo01 =

=

( cf + Сг) +

4 Cj c2 cos coo ^ + 2 cf cos 2 co01 +

-f 4 q c2 cos 3 co01 +

2cfcos 4 co01.

158

Таким образом, уже квадратичное преобразование импульса u(t) привело к появлению новых гармоник Зсо0 и 4(0<ь а также к уве­ личению уровня постоянной составляющей входного сигнала.

В более общем случае, при бесконечном числе слагаемых в (4.5), после нелинейного преобразования мощность перераспределяется, так что общая ширина спектра выходного сигнала увеличивается по

сравнению с шириной спектра входного сигнала.

Для определения величин А кр

заметим, что экспоненциальный

р

 

множитель ехр{і У Krmt}, определяющий частоту колебаний функ-

г = 2

при изменении кг(г—\, 2 , . . .,

ции uP(t) и множитель ехр{і/соѴ}

р) и к изменяется дискретно. Поэтому при каждом р величина А кр

будет

определяться

как

сумма

произведений

коэффициентов

Скг (Г= Ь 2,

...,

р),

соответствующих

таким

наборам

кг(г= 1 , 2 ,

. . ., р),

для которых

выполняется

 

р

 

кг= к.

Таким об-

равенство V

 

разом,

 

 

 

 

 

 

 

 

г=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о&

 

 

 

 

 

со

 

 

 

 

 

=

У ЬрАйр- Ь 0~\- Ьі с0 +

Ь2 у

с*, С—К1А~Ь3у у

у

cKi cKt cKj+ ...

 

р

 

 

 

 

« , = —«>

 

К , к г

к ,

 

 

 

 

 

 

 

 

 

 

 

 

к,+к23=0'

 

 

 

 

 

П

ч

+•••

 

 

 

 

 

 

 

* 1

Kt

Кр

1

 

 

 

 

 

 

 

 

где

р

 

 

 

 

 

 

 

 

 

 

 

 

V кт = О

 

 

 

 

 

 

 

 

 

 

 

г= 1

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

=

 

 

- Д

Г К ,

 

 

 

(4.37)

 

 

 

 

 

 

 

 

 

р—0

 

 

р—0

К і к2

 

Кр

Г— 1

 

 

 

 

 

где

Р

 

 

 

 

 

 

 

 

 

 

 

 

у кг = к,

 

 

 

 

 

 

 

 

 

 

/■ =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ьр A-кр = £ /ѵ У

1 J-

• •

I J П Ск<-.

 

 

 

 

 

р = 0

 

р= 0 к, к,

 

кр г— 1

 

 

 

 

 

Р

/сл = — /с.

 

 

 

 

 

 

 

 

 

где + у

 

 

 

 

 

 

 

 

 

 

г=1

 

 

 

 

 

 

 

 

 

 

 

Причем

суммирование в

(4.37)

ведется по всем кг от — о о до

 

, удовлетворяющим

равенству

р

кг— к.

 

 

;...

- ) - о о

у

 

 

/•=1

159

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