книги / Методы помехоустойчивого приема ЧМ и ФМ сигналов
..pdfдостаточной точности приближенного уравнения (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), в которой |
влиянием |
шума на время переходных процессов в астатической сисистеме можно пренебречь.
При аналогичном подходе можно показать, что в си стеме ФАПЧ с интегрирующим фильтром время пере ходных процессов с ростом шума незначительно умень шается. Качественное различие влияния шума на эти системы, вероятно, связано с тем, что в астатической системе ФАПЧ на фазовой плоскости имеются широкие инкрементные зоны, в которых накопление энергии из-за шума происходит более интенсивно, чем ее рассеяние в декрементных зонах.
Выводы. Расчет динамических характеристик стоха стических систем ФАПЧ позволил сделать следующие выводы: широкополосный шум приводит к затягиванию среднего времени установления синхронного режима в системах со значительными инкрементными зонами на фазовой плоскости, т. е. в таких, у которых затухание за период биений меняет знак (астатическая система, система ФАПЧ с ПИФ), и практически не изменяет (не