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

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

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

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

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

 

 

 

ной зависимости

математического ожи­

 

 

 

дания и корреляционной функции) сгла­

 

 

 

живаемого процесса.

 

 

 

 

 

 

 

Это означает, что интервал сглажи­

 

 

 

вания (параметры

системы)

необходимо

 

 

 

выбирать

с учетом

вероятностных

ха­

 

 

 

рактеристик сглаживаемого процесса.

 

 

 

 

Условия

применимости

оператора

 

 

 

скользящего

сглаживания

(3.113)

и

 

 

 

принципы

выбора

параметров

системы,

 

 

 

его реализующей,

детально

рассматри­

 

 

 

ваются в [72, 79]. Как показано в этих

Рис. 3.10. Зависимость коэффици­

работах,

высокая эффективность может

быть получена при сглаживании случай­

ента сглаживания Ко

амплитуды

ных процессов, математическое ожида­

импульсов устройством, реализу­

ние которых является линейной функ­

ющим

оператор текущего средне­

цией времени. В таких случаях интер­

го (а)

и &С-ф,ильтром

(б) от ßT

вал сглаживания

может быть

выбран

 

 

 

достаточно

большим.

Если

оператор

сглаживания воздействует на нестационарный случайный процесс, математичес­ кое ожидание которого изменяется во времени по закону, отличному от линей­ ного, максимальное значение интервала сглаживания ограничивается отрезком времени, на котором возможна линейная аппроксимация закона изменения ма­

тематического ожидания оглаживаемого процесса [76]. При этом

верхняя оценка

эффективности сглаживания выражается соотношением

 

t + T

t + T

2

 

Т ~ 2 j*

I Bg (ult u2) d Uldu2 <

,

(3.125)

l t

t

 

 

2

— максимальное значение дисперсии сглаживаемого

процесса; тк —

где оmax

время корреляции сглаживаемого процесса.

Оператор скользящего сглаживания, как уже выше отмечалось, сравнитель­ но просто реализуется линейными системами с сосредоточенными параметрами. Для достижения более высокой эффективности сглаживания применяют услож­ ненные операторы (см. [79]).

Сглаживание флуктуаций временного положения импульсов. Задача сгла­

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

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

Отметим, что подобный подход к задаче сглаживания флуктуаций времен­ ного положения импульсов основывается на взаимосвязи случайной длительности интервалов между последовательными импульсами + n — t nt n - i и остальными параметрами имп-ульсав —моментами возникновения, длительностью, амплиту­ дой. Это означает, что практически сглаженную последовательность импульсов можно сформировать, определяя среднюю длительность интервалов между по-

140

следовательным« импульсами по выборке из произвольной реализации сгла­ живаемого импульсного случайного процесса £(^). Усредняющие устройства мо­ гут ібыть дискретными (как правило, нелинейными), или аналоговыми (напри­ мер, линейная узкополосная система).

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

Если за истинное значение временного положения импульсов принять мо­ менты прохождения выходного сигнала фильтра через нуль, то в случае совпа­ дения средней частоты следования импульсов соо = 2п/Т с частотой свободных колебаний фильтра сор, получим несмещенную оценку временного положения импульсов относительно момента включения фильтра (т. е. момента возникнове­ ния нулевого импульса), причем эффективность такой оценки определяется ши­ риной полосы пропускания фильтра: чем уже полоса пропускания, гем больше длительность выборки из импульсного случайного процесса, участвующей в фор­ мировании оценки, а следовательно, меньше ее дисперсия [79].

Если частота свободных колебаний фильтра -не равна соо, оценка временного положения импульсов оказывается смещенной на величину расстройки Да>р=

-= (Ор—(ОО-

Применение рассматриваемого способа формирования оценки целесообразно, если интервал следования импульсов не меняется во времени (или меняется незначительно), а флуктуации временного положения импульсов ограничены ин­ тервалом (—7/2, 7/2), т. е. для сглаживания импульсных случайных процессов с детерминированным тактовым интервалом.

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

Учитывая, ято практически амплитуда и фаза процесса на выходе узкопо­ лосной системы — медленно меняющиеся функции времени, и, используя прин­ цип суперпозиции, приемлемый для линейных систем, представим квадратурные составляющие процесса %(t) в виде суммы детерминированной (регуляркой) и случайной компоиент. Тогда огибающая и фаза процесса на выходе фильтра оп­ ределяются формулами:

й к) м

= V

^

w]2+ К* (о + i4Ä) (о]*

(з. 126)

и

 

 

 

 

 

У( k ) = arctg

m ls (t) +

( 4 fe> (О

 

(3.126а)

m i c V ) +

И -'/’ (О

 

 

 

 

 

где

=

m ic(t)

и

(t) = ^ і к'>—m is ft) — случайные флуктуации коси­

нусной и синусной составляющих реакций фильтра соответственно, а математи­ ческие ожидания квадратурных составляющих m lc(t) и m u (t) определяются ф-лам.и (3.49) и (3.50) или (3.49а) и (3.50а).

В дальнейшем будем рассматривать вероятностные характеристики, в част­ ности математическое ожидание и дисперсию случайной фазы процесса £(f). При этом нелинейную зависимость случайной фазы от квадратурных составляющих заменим приближенной линейной. Погрешность, возникающая за счет линеари­ зации, как показано в [38], не -превышает 3—5% и поэтому может не учитывать­ ся при количественной оценке эффективности сглаживания случайных флуктуа­ ций временного положения импульсов узкополосной избирательной системой.

141

Тогда согласно методике статистической линеаризации (76] математическое ожидание случайной фазы

т іу (0 = т 1 { y (k). ß ) }

= arctg mis (t)

 

 

 

(3.127)

 

 

{Лю

(t) ’

 

 

 

 

 

а дисперсия флуктуаций фазы

 

 

 

 

 

 

Dv (0 = Ms b (fe) w} =

 

 

 

 

 

 

 

j

De (t) +[І І 1 ) 2

 

д у \ і д у

I

Mcs (О,

 

Лдх

D s (t) + 2i

 

C

■ дУ І

i x=.

 

д х J \ д у

I х = т 1

 

 

 

y = mis

 

y=m .

 

у—тІС

(3 .1 2 8 )

 

 

 

 

 

 

 

аде Y=Y(Ä)(0; x = m lc( t ) + \ i ck)(ty,

 

(ft).

D s (t) и

M 2Cs(t)

^/= mlsfO+ [x^'г,(t), a D c(t),

вторые центральные

моменты квадратурных составляющих [процесса

t,(t).

 

Частные производные, входящие в

выражение (3.128), равны

 

 

д у

д у

 

 

 

 

 

 

” = - У К * 2 +

У2) - > - ^ - = */(*2+ у * ) .

 

 

 

 

 

дх

 

оу

 

 

 

 

 

 

Подставляя

в

(3.128) значения

частных

производных

при х = т и ( 0

и

y —m u ( t ) , получим выражение для дисперсии в

виде

 

 

 

m 2ls (t) De (t) + т\с (t) D s (t) — 2m lc (t) m ls (t) M 2CS V)

(0 =

(3.129)

[m \c V) +

m\s (<)]*

Дальнейшие преобразования выражения для D y (t) связаны с величинами

т\с (t), mjs (t)j m ic ( t) m l s (t), D e ft),

D s (t), MzcS(t), методика определения кото­

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

т\с (t) =

б е г +

б е г cos 2шР < +

б е г

sin 2шР Г

 

ти (0 =Урег -

б е г cos 2“ Р і

-

 

Урег sin 2©п t ,

 

т іе

(t) mis (<)3 =

б е г sin 2c°p t —

б е г cos 2“ P Г

(3.130)

D c

 

 

 

б л

cos

2wP

(4 +

б л

sin 2c°P Г

 

 

D s (0 = ’б л +

 

 

 

(t)

 

 

■ Jen C0S 2ö)P (0

 

,CJ1sin2o)p

 

 

 

 

 

 

- h

 

 

 

 

^ j e s

( 0

— -■с/слл Sin“ --------2(üpП Ѵ(t)

"J слc

 

COS 2ü)p t ,

 

 

 

В соотношениях (3.130)

приняты

обозначения:

 

 

 

t

t

 

 

 

 

 

 

 

 

 

 

 

б ег =

.

j

mlE (( — °i) ml£

 

~

"a)

(Oi) К Ы

X

 

 

0 0

 

 

 

 

 

 

 

 

 

 

 

 

X c o s [Cöp ( v 2 — Щ ) — Ф

( ü j)

+

ф

(t)2) l d V x d o i ,

 

 

 

 

i

t

 

 

 

 

 

 

 

 

 

 

 

' p e r -

"2~ J J m‘£ ^ — 01) m‘£ ^ ~

 

h° ^

x

 

 

 

0 0

 

 

 

 

 

 

 

 

 

 

 

142

Xcos [Mp (v2 +

Vi) — Ф

( с ц ) — ф

(t>2)] dv,doa,

 

 

 

t t

 

 

 

 

 

X s in

[ö p (a2 +

о г — ф

(Oj) — ф ( о а )] dOidvi,

 

 

 

i

t

 

 

 

 

 

x c o s

[Mp (o2 — ü i) — ф

(i>i) + Ф

(o2) l dv^Vi,

 

 

 

t

t

 

 

 

 

 

 

о b

 

 

 

 

 

XCOS [fflp (o2+ fl) —Ф (fl) — Ф (f2)l * A ,

 

 

 

t

t

 

 

 

 

 

 

о

0

 

 

 

 

 

X sin

[Mp (t>2 +

t)x) — Ф

(t)i) — Ф

(t»2) l dv1du2,

 

(3.131)

Тогда, после подстановки соотношений (3.130)

в (3.129)

и преобразований

для D y (t) получим

 

 

 

 

 

 

 

 

 

 

 

(3.132)

Анализ

соотношения

(3.132)

показывает, что

дисперсия

флуктуаций фаз»

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

Учитывая допущения, которые были приняты при выводе формул для m lc(t), rrtis(t), D c{ t ) D 3(t) и M cs(t) , можно заключить, что в первом приближении дисто­ персия определяется статистикой случайных амплитуд и временного положения импульсов воздействия. Что касается характеристик узкополосного фильтра, на величину дисперсии существенно влияет лишь полоса пропускания Лмп =я/е

(т. е. величина коэффициента затухания а) и частота свободных колебаний фильтра Юр (точнее говоря, величина расстройки Аюр= мр—Мо). Более подроб­ ные характеристики фильтра практически не влияют на величину дисперсии D f .

Заметим

еще, что соотношение .(3.132) так же, как и

более точное соотношение

(3.129),

определяет зависимость D y

от времени,

так

как от времени (от

 

зависят .m ic(t), rriis(t),

D c(t), D s (t),

M cs(t). Причем анализ ф-s

(3.91)—'(3.103) показывает, что в переходном режиме на величину дисперсий влияет и нестационарность воздействующего импульсного .случайного процесса l( t ) и изменение во времени характеристик фильтра (в основном полосы про­ пускания Дюп ), a в установившемся .режиме только нестационарность |(<). Это

означает, что соотношения (3.132) и, конечно, (3.129), позволяют определять и продолжительность .переходного процесса и значение дисперсии D y в установив­

шемся режиме, т. е. при t-*-оо (ІѴ =±оо).

143

Напомним, что так как фильтр узкополосный, его реакция .представляет со­ бой последовательность радиоимпульсов бесконечной длительности. Так что нестационарность процесса £(0> обусловленная только зависимостью от временя случайных величин ѵ п, г, за очет сложения большого числа импульсов в каждый момент времени оказывается в значительной степени сглаженной [25].

Можно показать, что характер изменения дисперсии в переходном режиме определяется функциями fi(N) и fx(N), причем

М*) = Ѵ

£

- а(п+/>7г

л=0 /=0

 

N

N

—а (rt_b/) Тг io)p (n-^-j)T

 

\^і

 

 

(3.133)

n= 0 /=0

Так как установившийся режим соответствует бесконечно отдаленному моменту включения фильтра, то в суммах формул для f i ( N ) и fz(N ) в этом случае следует распространить пределы суммирования до ± оо. При этом указанные суммы можно заменить интегралами

t t

1(N)~ И е~“ {х+у)е_ІСОр(х~у)dxdy'

о о

к 0 — J J е-“

е~‘“Р (Х+У) d x d y .

Замечая, что интеграл для f i ( N ) содержит быстроосциллирующий множитель е—12<арх) можно с пренебрежимо малой погрешностью считать, что в устано­

вившемся режиме k ( N ) = 0 и / рег = ^сл = ^рег= ^сл —0. а дисперсия фазы

0 Ѵ= і сл/24ег-

 

 

 

(3134>

Учитывая, что в установившемся режиме

 

 

 

Вес (к х) = ве$((>т) = J JJ

в ъ (т — г) fto (о) h0 (V +

г) cosсор zdzdv =

 

−−−−−−−−09

 

оо

 

оо

 

оо

 

 

= ~

*j в і (t — г) bh (г) cosсоР zdz =

J

(со)

*j b h (г) coscop г X

 

— оо

 

0

— оо

X costo (т — г) dz,

 

 

 

 

 

00

 

 

 

 

-де b h ( z ) =

J h a( v ) h a ( v - s z ) d v ,

 

 

 

 

 

— oo

oo

 

 

 

 

 

 

 

 

B l (® ) =

(“ ) + F j i (co) = —

j*В j (T z) cosсо ( T z) d со

«0

J/i0(v) cos copt) coseovdv =Re <T (i co)>,

—00

144

со

J (o)cos (OpüXsin шvdv = Jm [ T (i со)],

— OO

получим для

"Р +ЛюПj

 

4 , = \

00(Ѵщ (<о)ТЧсо)4ш =

FBl(u)T*{a)da

УРег = ^ j ^дб(“ ) Г2 (“ ) d“ = - у Os (“ о) Г» ( со;) .

о

Таким образом, в установившемся режиме дисперсия флуктуаций фазы D,f

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

цесса £(£), попадающей в полосу пропускания фильтра, ,к мощности периодиче­ ской составляющей процесса F Rg (сор) , соответствующей периоду следования им-

пульсов:

Юр+До)!

D

1

Тщ (и) d со »

.. Дсоп

(3.135)

у

 

2/д&(“ о) тг (“ о)

 

25 (m0)

 

 

 

 

 

 

шр+ АиП

 

 

 

где Асоп =Г_2і(ио) J

Т г (ы )й (й ■—полоса пропускания фильтра.

 

 

“р—Аап

 

 

 

Следует отметить, что ф-лы і(3.132), (3.134), и особенно (3.129), дают более точное значение дисперсии D y как в установившемся, так и в переходном ре­

жиме при минимуме допущений относительно вероятностных свойств сглажи­ ваемого процесса £(/) и характеристик сглаживающего фильтра. Очевидно, что можно из указанных формул получить и иные, нежели (3.135) приближения (см., например, [27, 70]), Однако приближение (3.435) дает, по нашему мнению, наиболее простые формулы связи статистических характеристик параметров им­ пульсов и полосы пропускания фильтра с величиней дисперсии D y . Так, из

ф-лы (ЗЛ35) следует, что при уменьшении ширины полосы пропускания Ашп пропорционально уменьшается и D y . С другой стороны, при уменьшении Дсо^ увеличивается инерционность фильтра а=я/Амп , а значит и длительность пе­

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

из условия [93]:

е “ОТГ<о,05,

(3.136)

которое, после логарифмирования обеих частей -неравенства можно привести к виду

N Js In 20/аТ г « Дсоп /Гг .

(3.136а)

145

Г л а в а 4

Нелинейные преобразования импульсных случайных процессов

4.1. ОСОБЕННОСТИ АНАЛИЗА. НЕЛИНЕЙНЫХ ПРЕОБРАЗОВАНИЙ ИМПУЛЬСНЫХ СЛУЧАЙНЫХ ПРОЦЕССОВ

Сущность нелинейного преобразования состоит в том, что слу­ чайный процесс l(t), являющийся результатом такого преобразо­

вания, связан с преобразуемым процессом ) в общем случае нелинейным дифференциальным уравнением вида

dm i{t)

t

(4.1)

(dtr

' _

 

где F{x(t)] — некоторый нелинейный функционал своих аргументов, характеризующий оператор преобразования Кнл . Конкретный вид функции F определяется свойствами и параметрами системы, реа­

лизующей оператор Ьнл ■

удается

Во многих практически важных случаях оператор Ьнл

представить в виде

 

Ці) = Ьн л {иі)} = Ьл {Ьн і т } -

(4.2)

Это означает, что система, реализующая оператор (4.2),

состоит

из последовательно включенных нелинейной Іц и линейной L л ча­ стей (рис. 4.1), причем каждый из операторов ін и і л оказывает­ ся более простым, чем исходный оператор Сил . Нелинейная сис­ тема, допускающая представление в виде (4.2), называется при­ водимой [84] или разделимой [86] (рис. 4.16). В противном случае говорят о неприводимой или неразделимой нелинейной системе (рис. 4.1а).

Соотношение (4.1) показывает, что в общем случае значения случайного процесса "Q(t) в каждый момент времени t зависят не только от значений процесса £(7) в тот же момент времени, но и от значений его в другие моменты времени. В таком случае нели­

146

нейная система называется инерционной. Системы, процесс t,(t) на выходе которых полностью определяется значениями процесса %(t) на входе в совпадающие моменты времени, называются бе­ зынерционными.

Нелинейное безынерционное преобразование определяет функ­

циональную связь случайных величин t,(t*)=y и £(/*)=х

в каж ­

дый фиксированный момент времени t*\

 

У = /(*),

(4.3)

Для анализа нелинейных систем существенным является вид функ­

ции f(x).

Как известно [112], эта функция может быть однозначной

или

многозначной,

 

Ф

 

 

гладкой или разрывной

 

 

 

 

т

Нелинейная

№ ы Ишо}

(см. рис. 4.2, 4.3).

 

 

 

 

 

система

 

Следует

отметить,

 

 

 

 

что исследование безы­

^

 

 

 

нерционных

нелиней-

Приводимая нелинейная система

ных систем

значитель- .

Нелинейная інлШѴІ.

Линейная W H J LJ W }}

но проще, поэтому при

Ü,

анализе

стремятся

по

 

часть

 

часть

возможности разделить

 

 

 

 

нелинейную систему на

Рис. 4.1. Нелинейные системы:

части, содержащие

или

а) неразделимая; б) разделимая

только линейные инер­

 

 

 

 

ционные

элементы, или нелинейные безынерционные элементы.

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

Рис. 4.2. Характеристики существенно нелинейных элементов:

а)

с насыщением; б) релейный; в) с насыщением и зоной нечувствительности;

г)

релейный с зоной нечувствительности; д) с зоной нечувствительности и от­

сечкой; е) ограничитель с отсечкой

147

схема подобного оператора преобразования, называемого типовым

радиотехническим звеном [53], приведена на рис. 4.4.

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

Множество различных частных методов основывается на учете особенностей оператора преобразования (безынерционный, одно­ значный, кусочнораз­ рывной и т. д.) илиопе-

Рис. 4.4. Типовое радиотехническое звено:

спектральной

плот-

НОСТЬЮ

И Т ,

ІП .).

И'ОСЛб-

/ — линейный

фильтр

I, 2 — нелинейный элемент,

О'СО'ббНіНОСТИ

3 — линейный

фильтр

//

до б аНіия

б елі-шеиіных

 

 

 

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

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

148

нейных системах, сохраняются в таком случае и при анализе им­ пульсных случайных процессов.

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

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

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

При ограниченной длительности импульсов k-я реализация слу­ чайной последовательности импульсов в соответствии с (1.4) может быть записана в виде

Здесь функции un(t) существуют в несовпадающие интервалы вре­ мени и каждая из них может быть представлена рядом

со

(4.4)

где А пк и s K — действительные или комплексные числа. В частном случае, когда sK=\кыо, где ы0 = 2п/Т, Г-интервал задания функций un(t), ф-ла (4.4) определяет ряд Фурье

СО

(4.5)

П К с

М 0 = V с,

 

где коэффициенты спк могут быть найдены по формулам Эйлера:

149

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