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

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

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

Соотношения (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 —

 

tt + 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—vz) 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 — ü ,/ + Т — Vz) 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)

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-го импульса; б(ttn) — смещенная на случайное время 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

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