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

книги / Методы помехоустойчивого приема ЧМ и ФМ сигналов

..pdf
Скачиваний:
16
Добавлен:
12.11.2023
Размер:
13.73 Mб
Скачать

достаточной точности приближенного уравнения (1): г< 1. Соответствующее уравнению (1) стохастическое уравнение системы имеет вид

Л

где £(•/)— дельта-коррелированный процесс единичной интенсивности. При нулевой интенсивности шума D = co, и это уравнение переходит в уравнения для энергии, полученные в [11], а коэффициент вФг(А) =/(Л ) опре­

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

уменьшения энергии А в отсутствие

шума. Значения

Ф \ ( А ) , Ф2(Л),

в, D

будут вычислены

при анализе конкретных систем.

режим

установившимся,

Будем считать синхронный

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

А

Уравнение для среднего времени /(Лгр /Ло) дости­ жения границы Лгр процессом Л (t) получено в [10]. В нашем случае оно имеет вид

Отсюда для верхней границы Лгр= Л в]>Л 0

t (Ло/Л0) = Г

 

Г — —

е*<*> dy

dz,

(3)

J

 

J

еФ2()0

 

 

 

 

До

о

 

 

 

 

 

для нижней Лгр =

Л„ •< Л0

 

 

 

 

t (Л„/Л0) =

Г

Г

Д -

 

е*м dy

dz,

(4)

 

J

J еФ2(у)

 

 

 

Ан U

г

где

о

Решения (3) и (4) удовлетворяют нулевым граничным условиям [10]

t { A BIA0 = А а) = 0; t ( A J A 0 = А„) = 0 .

(5)

Величина, обратная первому интегралу уравнения (2)

л

dAo/dt, характеризует среднюю статистическую скорость изменения энергии А при действии шума, а при отсугствии шума равна скорости изменения А.

Введем среднее время t достижения границ процес­ сом A ( t ) при некоторой плотности вероятности началь­ ной энергии Wo (Ло) :

t J* t

Mo)

W Q (i4o) d A 0, A j ^ Л 0

Л2.

A,

 

 

 

 

 

Если A„ à,

a A„ =

oo, то

среднее время достиже­

ния процессом A ( i )

границы Д <Л0 характеризует дли­

тельность переходного

процесса

в системе

ФАПЧ при

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

Оценка надежности режима синхронизации в стохас­ тических системах ФАПЧ. Выберем положение грани­ цы А. Примем А равным математическому ожиданию

энергии A ( t ) в стационарном режиме:

А = < Л СТ (/)> .

Величина < А СТ(t) > пропорциональна

сумме средних

квадратов координат системы, определяющих точность

синхронизации. Величина

< А ст (/)>

определяется из

решения стационарного

уравнения ФП

[10]

< А СТ ( * ) > =

| A W CT( A ) d A ,

 

(6)

где W CT( A ) = C

ехр

f

4 È 1 (Л) d A

 

D

 

LoJ

ф»(

 

 

2 е Ф 2(Д )

 

 

 

 

 

(А )

 

 

 

 

 

 

 

а константа С — из условия нормировки j

W „ (Л) d A —

 

 

 

 

о

л

= 1. После достижения границы А за время t(A/Ao) следует определить время пребывания системы в сир-

кронном режиме.

Известно,

что превышение

энергией

A (t) величины 2

(выход из

области жесткого

синхро­

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

ната /1 (t) движется по траектории;

начинающейся при

Л=Д

и целиком лежащей в области жесткого синхро­

низма

А <Л„ =2. Среднее значение

длительности этого

интервала 'Назовем средним временем надежной работы

А

системы 1(2/А). Достижение А—2 не всегда приводит к перескоку фазы на 2я или возникновению длительно­

го

переходного процесса. Точка от Л = 2 может вернуть­

ся

в область синхронного режима без превышения

 

А

Лв =2. Время /(Д/2) возврата траекторий суровня Ло=2 на уровень Д является дополнительной характеристикой работы системы в синхронном режиме.

Будем считать, что система находится в синхронном режиме, если ее движение происходит по траектории,

начинающейся в точке Л0 = Д (точка

1,

рис. 1,а) и кон­

чающейся в точке Л = 2 (точка 2),

и

в асинхронном,

если ее движение происходит по траектории, начинаю­ щейся в точке Л = 2 (точка 2, рис. 1, а) и кончающейся в точке Л = Д (точка /).

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

 

Асинхронныйрежим

|

■ —»

â

-------------------- 1

Синхронный

I

I

J

 

________________ |

режим

|_________ j_____

11*

163

ный и асинхронный режимы. Вероятности синхронного

Рс и асинхронного

Р а режимов определим выраже­

ниями

 

 

,с ^

t (2/А)

 

 

?(2 /Д) + * (Л /2 )

 

 

Л

 

Pa =

<(А/2)

(7)

 

< ( Д / 2 ) + Î (2/Д)

 

 

л

л

Таким образом, зная величины t(A/2)

н t(2jA), мож­

но оценить вероятность работы системы в синхронном режиме Рс. Длительность вхождения системы в син­

хронный режим определяется средним временем дости­

жения

марковским

процессом А (t)

при

движении с

уровня

Ао границы

Д= < Л ст( /) > .

Это

л

время t(A/Ao)

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

Оценка влияния шума на динамические характери­ стики астатической системы ФАПЧ. Проведем расчет

среднего времени достижения синхронного режима в астатической системе ФАПЧ и надежности ее работы в этом режиме.

Операторный коэффициент передачи фильтра в цепи обратной связи такой системы имеет вид /С(р) = ( 1 + +pTi)/pT2, где p=d/dt, Ти Т2— постоянные времени

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

метры

фильтра [ 1 1 ]

b%— QçTt

и

a J =

Qc r 2 (2 c —

полоса

синхронизма

бесфильтровой

системы

ФАПЧ),

найдем

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

е H_ 0 I (A),

Ф2 (Л)

в

асимптоти­

ческом

случае [9] ( e = l/] /a i) :

 

 

 

 

 

-(4^/ЗК,) {/<, ( l - T )+(J4-1)£i} + ^

 

 

2;

ф г И ) =

_ (4й,/ЗА'2) VÂJ2{(A - 1 )Е2-

(8 )

 

- ( A - 2 ) A 2} + - L , А > 2;

 

Ф2 (Л) =

(9)

2А (Еа/К2),

А > 2 ;

где Ег = Е ( YÂ]2), К1= К ( / т

Е2 = Е( У 2/Л),

К 2= К { V2/A) -

— полные эллиптические

интегралы 1 -го и Н-го рода

[ 12].

Значение параметра D, характеризующего интенсив­

ность шума, определяется формулой [9] /> = 2 Э Д ;

( 10)

О

где R (т) — нормированная автокорреляционная функция

входного шума | (t) , приведенного к выходу ФД (5?т = = с2 /?(т)); No — спектральная плотность входного шу­

ма на центральной частоте подстраиваемого генератора; хк — время корреляции £ (t) . Оно предполагается мень­ шим всех постоянных времени системы. Параметр D — это отношение мощностей сигнала и шума в полосе 2 С

[4].

л

л

vbküJ

Для

л

расчета величин t(A{Ао), /(Д /2 ), t(2/A) по фор­

мулам

(3) и (4) с учетом

(8 ) и (9)

нужно использо­

вать ЭВМ, так как интегралы, входящие в эти выраже­ ния, не удается выразить в известных функциях. Одна­ ко оказалось, что качественно верная и количественно достаточно точная картина влияния шума на динамику системы получается при аппроксимации (8 ) и (9) их

разложениями в ряды при А < 1

и

2 соответственно.

На

рис.

2 штрихпунктиром

показана

зависимость

Ф](Л)

при

отсутствии

шума.

Ее

аппроксимация пред­

ставлена зависимостью

Ф\(А)

при Dbi= c o .

Сравнение

их показывает, что при больших начальных энергиях А 0

точность определения длительности переходных процес­ сов по такой аппроксимации удовлетворительна. Поэто­ му вместо (8 ), (9) запишем

Ьх А -Ь 1/D, А < 2; - V 2 + 1 / A Л > 2 ;

беличина, обратная средней статистической скорости

уменьшения энергии А

с уровня А 0 в направлении А —

= 2 < Л 0 в соответствии

с (4), будет

определяться вы­

ражением

 

 

 

d t & I A 0) _

— 1

1

( 12)

d A 0

еФДЛо)-

гф ./2 - l/ D )

Как видно из (12), скорость уменьшения энергии в сис­ теме во время переходного процесса по частоте сни­ жается с ростом шума.

 

А

 

 

 

по

(4)

Величина Л(Д/Л0)/йЛ0 при Л0 < 2 находится

с учетом (11):

 

 

 

 

 

d t ( L j A 0) _

1

1 — e -Dü,(2 -A„) ( J

\ D

Y Лп-^ 2 .

d A 0

еЛ0

 

D b , — 2/

 

 

 

 

 

 

 

 

(13)

Зависимости —d A 0l d t ( A 0), рассчитанные по

(12)

и

(13)

для различных D, приведены на рис. 2 сплошными ли­ ниями.

Выражение для средней статистической скорости, оп­

ределяющей увеличение энергии в системе при Л 0

2 и

Л „ < 2 , рассчитанное по формуле (3), имеет вид

 

---- = в М о ( е ° м . _ I)-1, Л0< 2 .

(14)

166

Величина этой скорости определяет время, через кото­ рое в среднем энергия системы достигает величины Лв<!2 . Эта зависимость дана на рис. 2 пунктиром. С рос­

том шума скорость

увеличения

энергии

от

значения

Л0= Д возрастает.

 

 

 

 

 

 

 

 

 

 

 

Расмотрим стационарный

режим. Из

(6 )

 

среднее

значение энергии

< Л СТ(< )> * = 1 / Dbx = А

 

при ма­

лых шумах

D bx >

1 (рис. 2). Нижняя граница

 

А„ = А

определена. Проводя интегрирование в (12), (13)

и (14)

и учитывая

нормировку

Л

 

 

Г —

 

получаем

(t =

(2 С/ 1/ л,) t),

формулы для времени

 

установления

синхронного ре­

жима ЦА/Ао) =it(2fA0) -И(А/2)

 

и величины / ( 2 /Д) :

t (2 /Л.) =

b -

 

(Л„ -

2 );

 

 

 

(15)

 

 

2 С Dbx — 2

 

 

 

 

 

 

m 2 ) =

 

lnD&, — е _206‘

(

 

4£>&t

 

X

с &1

 

D ô i-

2

 

 

Й

L

 

 

 

 

 

 

 

 

 

 

Х ( Ei* (2Dbx) — Ei* (l))

 

 

 

 

 

(16)

< (2/А) =

[El* (2Dbt) -

Ei* (1) -

In 2Dbt],

 

(17)

Sc bt

 

 

 

 

 

 

 

 

 

 

 

где Ei* — интегральная

 

показательная

функция

[12].

Среднее время переходного процесса по частоте //ш при

включении системы со случайной начальной разностью

фаз [13]

определяется величиной

/(2/Ло) и расстройкой

у, если вместо Л<> в формулу (15)

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

значение А 0— (агу2/2) + 1 [13] :

 

 

tf =

а, Ч 2РЬХ \ ( а, у2

 

а , Т2/2 > 1 - (18)

-

1

•'ll!

2 Сbx \ D b 1— 2 j \ 2

 

 

 

 

 

Ha рис. 3, а приведена зависимость от Db{ отношения

средних времен переходных процессов по частоте в аста­ тической системе ФАПЧ при действии на ее входе ши­

рокополосного

шума и без него

{Db{ =

оо). Величина

этого отношения, как следует из

(18), определяется фор­

мулой

 

 

 

 

 

l fml l f =

DbJ(Dbi - 2 ) .

(19)

Затягивание переходных процессов по частоте в та­

кой системе

из-за шума

практически

не проявляется

1 , 1 ) при отношении сигнал/шум в полосе си ­ стемы q> 5 (рис. 3, а). Отношение q определяется соглас­

но [4] как отношение амплитуды сигнала

UQ к среднему

квадрату

напряжения

шума <тш' Я =

{Ус/ аш = V D b ,

При этом

Рс/Р ш = q2l2.

Условие малого влияния шума

на время переходных

процессов в астатической системе

<7>5 соответствует Р

с/Рш> 12,5 в полосе системы. Ве­

личина t f j t f при этом не зависит от начальной рас­ стройки у. Значение дисперсии времени переходного

процесса по частоте о- можно определить по следую-

*/ш

щей формуле, полученной таким же методом [14]:

Как и следовало ожидать, значение

о!

зависит от Hà-

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

и увеличивает­

ся при

ее росте. На рис. 3, б приведены зависимости

/ al

I tfm

от интенсивности шума

для двух

значений

обобщенной

начальной расстройки

а}у212. Из

рисунка

видно, что увеличение интенсивности шума приводит к росту | / v j f j l t f m .

На рис. 4 приведена зависимость среднего времени t

установления синхронного режима в системе, рассчитан­

ная

по

формулам

(15)

ti2ci>//Jaf

 

и (16),

от

начальной

 

расстройки

 

при

 

раз­

 

 

личных

величинах

ин­

 

 

тенсивности

шума

на

 

 

входе

системы.

 

При

 

 

малых

 

расстройках,

 

 

когда

среднее

время

 

 

установления

синхрон­

 

 

ного режима

определя­

 

 

ется, в

основном,

пере­

 

 

ходным

процессом по

 

 

фазе

(16),

t с увеличе­

 

 

нием шума уменьшает­

 

 

ся и оказывается мень­

 

 

ше

среднего

времени

 

20 40 Dôi

переходного

процесса

8 10

при

отсутствии

шума

-2

 

(Db j= o o ) .

Это

объяс­

-4

 

няется

тем,

что

при

-8

 

увеличении

интенсив­

-10

 

ности

шума

граница

-20

 

A = l/D b u

которая

оп­

-40

 

ределяет

 

появление

SO

 

синхронного

режима,

 

Рис. 5

сдвигается

вправо

(см.

 

 

 

рис.

2 ),

и

среднее

вре­

 

 

мя достижения всеми траекториями этой границы стано­ вится меньше.

Используя (16) и (17), рассчитаем вероятность Р с (7), определяющую надежность работы системы в син­ хронном режиме. Для Dbx>3 расчетная формула Рс

имеет вид

р

~ 1 _ Sbbx-\-2 e _ 2 0 i l

c

Dbx- 2

На рис. 5

приведена зависимость ig A = l g (1— A.)

от интенсивности шума и параметров системы. Видим, что при *7 = 5 в полосе системы Р » < 1 0 - , °. При этом среднее время надежной рабочей системы t (2/Л) без пе­

рескоков фазы оценивается по формуле (17).

Отметим, что при Ь\!ах > 1, используя аналогичный

подход в сочетании с квазистационарным методом [9], можно получить следующее выражение для среднего времени переходного процесса по частоте:

Сравнивая его с (18), видим, что влияние шума на ве­ личину tjm во всей области параметров астатической

системы ФАПЧ одинаково. Поэтому зависимость Ьш I */ (см. рис. 3, а) для асимптотических параметров {bitYâi<^\) верна и в квазистациоиарном случае. Одна­

ко

в полосе

системы в

общем

случае

[4] q —

= )/i>& i/(l+m ôi). Поэтому

области

Dbx>20

соответст­

вует

область

(7 > l / 2 0 /(l+ m fri), в которой

влиянием

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

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

Выводы. Расчет динамических характеристик стоха­ стических систем ФАПЧ позволил сделать следующие выводы: широкополосный шум приводит к затягиванию среднего времени установления синхронного режима в системах со значительными инкрементными зонами на фазовой плоскости, т. е. в таких, у которых затухание за период биений меняет знак (астатическая система, система ФАПЧ с ПИФ), и практически не изменяет (не­