книги / Методы помехоустойчивого приема ЧМ и ФМ сигналов
..pdfвыми для начальном подпороговой области, где линей ная модель приема осуществляется с большой вероят ностью.
Следует заметить, что в известных методах анализа спектральная плотность шума на выходе демодулятора характеризуется средним значением числа перескоков на 2п разности фаз колебания опорного 4M гетеродина
и сигнала. Такое определение возможно лишь для пуас соновского потока перескоков.
Отказ от использования линейного приближения и тем более попытка учесть статистическую зависимость между моментами появления скачков фазы колебания 4M гетеродина приводит к существенным математичес ким затруднениям. Кроме того, в настоящее время не существует методов расчета подавления частотной мо дуляции шумом в случае следящего приема. Все это подчеркивает важность экспериментального исследова ния помехозащищенности следящего приема в подпоро говой области, который пока что является единственно возможным методом исследования в случае достаточно слабого сигнала.
В настоящем сообщении приводятся результаты из мерений помехозащищенности синхронно-фазового де тектора и стандартного детектора с близким к идеаль ному ограничителем частотно-модулированных сигналов, используемых для передачи многоканальной телефонии.
Измерения проводились с использованием 4M демо дуляторов, рассчитанных на прием 60-канальной теле фонии с эффективной девиацией 400 кГц на канал.
В эксперименте использовался СФД с пропорцио- нально-интегрирующим фильтром и стандартный 4M де тектор с фильтром П4 с переменной полосой пропуска ния на входе.
Экспериментальные пороговые |
кривые, |
приведенные |
|||
на рис. |
1 , представляют собой |
зависимость отношения |
|||
сигнал/шум на входе в полосе AM |
(рАМ) |
и отношения |
|||
сигнал/шум на |
выходе телефонного канала (Рс/Рт)> |
||||
а также коэффициента подавления |
(Р). |
группового |
|||
При |
наличии |
предыскажений |
спектра |
||
сигнала |
в нижнем телефонном |
канале в подпороговой |
области отношение снгнал/шум будет наихудшим. Ис ходя из этого, измерения проводились для нижнего ка нала со средней частотой 14 кГц. Кривая 1 показывает
отношение сигнал/шум в канале, кривая 2 — подавление
сигнала для СФД.
Существенно отметить, что пороговые кривые 3, 4, 5, 6
для ЧД снимались для значении отношении ширины по лосы пропускания фильтра ПЧ к средиеквадратической девиации многоканального сообщения, равных 7,5; 6 ; 5;
4 соответственно, что, как показано в [2], существенно влияет на помехоустойчивость. При построении порого вых кривых был учтен псофометрический эффект, а так же уменьшение на 4 дБ девиации сигнала в нижнем ка нале при наличии стандартных предыскажений.
Из сопоставлений пороговых кривых видно, что СФД
суменьшением отношения сигнал/шум теряет выигрыш,
апри слабом сигнале проигрывает ЧД по помехоустой чивости. Во время измерений было отмечено, что приме нение перед СФД фильтра ПЧ с переменной полосой пропускания существенно не повышает его помехоустой чивость.
Исключение составляет случай слабого сигнала. В этом случае можно получить такие же результаты как и для ЧД, что в общем очевидно. С другой стороны, при использовании ЧД существует явный оптимум по поло се пропускания фильтра ПЧ [2]. Значения оптимальной
полосы пропускания заметно уменьшаются при сниже нии отношения сигнал/шум, а помехоустойчивость 4M приемника с переменной полосой пропускания в подпо роговой области несколько возрастает.
Этим, в частности, можно объяснить снижение эф фективности применения следящего приема. Однако следует заметить, что и при постоянной полосе пропус кания фильтра ПЧ пороговая кривая СФД при слабом
сигнале |
имеет более |
резкий спад (кривые |
/ и 3 на |
рис. 1 ), |
а подавление |
полезной модуляции |
возрастает |
быстрее, чем следует из известного соотношения для ЧД Р = (1 —е - Р)2> где р — отношение сигнал/шум в полосе 4M.
Изученные выше особенности поведения характери стик СФД в подпороговой области определяются самим принципом следящего приема, где для организации до полнительной фильтрации сигнала от шума оценивается фаза сигнала и формируется копия сигнала с помощью напряжения частотно-модулнроваииого гетеродина. При большом отношении сигнал/шум копня достаточно хоро шо отражает сигнал, а фаза ЧМГ определяется в основ ном фазой, сигнала и относительно слабо зависит от ре ализаций шума на входе демодулятора. Это обеспечива ет достаточную эффективность дополнительной фильт рации. В области слабого сигнала мощность шума на выходе демодулятора становится соизмеримой с мощно стью демодулированного сигнала. В этом случае с уменьшением отношения сигнал/шум заметно возраста ет дисперсия фазы ЧМГ, а также зависимость фазы ЧМГ от реализаций входного шума. В результате начи нает с большей вероятностью отслеживаться шум, при этом эффективность дополнительной фильтрации сни жается. Последнее приводит к более резкому ухудше нию отношения сигнал/шум на выходе следящего демо дулятора.
В заключение отметим, что при достаточно слабом сигнале следящий прием по помехоустойчивости практи чески равноценен приему с использованием стандартно го ЧД, если фильтр ПЧ выбирается оптимальным.
СПИСОК ЛИТЕРАТУРЫ
1. Витерби Э. Д. Принципы когерентной связи. М., «Сов. ра
дио», 1970,
2. Выбор фильтра тракта 4M приемника, работающего в поро
говой области. — В кп.: Вторая научно-техническая |
конференция |
||
по космической радиосвязи. Тезисы докладов. |
М., |
1971, Авт.: |
|
Ю. С. Агапов, В. В. Герасимов, В. М. Дорофеев, |
В. 3. Озерский. |
||
|
УДК |
621.396,621.33 |
|
|
10. |
С. |
АГАПОВ |
ВЛИЯНИЕ НА ИНТЕНСИВНОСТЬ ФЛУКТУАЦИИ ЧАСТОТЫ СТАТИСТИЧЕСКОЙ ЗАВИСИМОСТИ ПЕРЕСКОКОВ ФАЗЫ НА 2л СУММЫ СИГНАЛА И УЗКОПОЛОСНОГО ШУМА
Иллюстрируется влияние статистической зависимости переско ков фазы на флуктуации частоты. Приводится выражение для спектральной плотности флуктуации на нулевой частоте при рас строенном по частоте сигнале. Полученные результаты сравнивают ся с результатами, рассчитанными по методике Райса.
Знание спектральной плотности флуктуаций произ водной фазы на нулевой частоте Уу (0 ) суммы ВЧ сиг
нала и узкополосного шума позволяет во многих практи ческих случаях определить помехоустойчивость демодуля торов 4M сигналов. Для вычисления / у (0) особенно ши роко и результативно используется в настоящее время им пульсная модель (Райса) взаимодействия сигнала с шумом [1]. Однако соотношения, полученные О. С. Рай сом, недостаточно точны при малом отношении сигнал/ шум. В работе [2] результаты для этой области уточне ны путем учета статистической связи между моментами появления перескоков фазы на 2 л. Влияние статистиче
ской зависимости также отражено моделью, предложен ной H. М. Блекмэном [3], в которой рассматриваются нули ВЧ процесса. Обе модели позволяют получить зна чения для / у (0 ). верные при произвольном отношении сиг.нал/шум. Однако авторами работ [ 2 и 3] исследовал
ся немодулированный сигнал, настроенный на централь
ную частоту фильтра |
ПЧ. Ниже |
показано, |
как влияет |
||
учет статистической зависимости |
|
перескоков |
фазы на |
||
значения / у (0) при |
расстройке |
сигнала |
по |
частоте. |
И
Результаты получены путем некоторого преобразований модели нулей, предложенной Блекмэном в [3].
Сделаем некоторое обобщение известных моделей взаимодействия сигнала с шумом. Для этого рассмот рим методы учета статистической зависимости переско ков фазы на 2 я.
Сумму |
сигнала |
и узкополосного шума запишем в |
виде |
|
|
|
•'I ( t) = E ( t) cos к ; / + > ;(/)î,! |
|
где E(t) |
и ф(t) — медленно меняющиеся амплитуда и |
|
фаза; ©о — средняя |
частота энергетического спектра |
|
шума. |
|
|
Спектр шума будем полагать симметричным относи тельно ©о-
Рассмотрим пересечения фазой q>(/) равноотстоящих
уровней |
с шагом |
Д |
радиан |
на интервале времени |
|||||
(—7/2, |
+ TJ2). На рис. 1 |
показана |
|
||||||
векторная |
диаграмма |
процесса |
ц. |
|
|||||
Фаза |
процесса проходит последова |
|
|||||||
тельно уровни |
0 , я/2 , л, —я/ 2 |
или |
|
||||||
О, я/2 , я |
и опять я, но в обратном |
|
|||||||
направлении. |
|
|
|
уровнями |
|
||||
Шаг |
между такими |
|
|||||||
Д =я/4. |
|
|
N r |
как |
разность |
|
|||
Определим |
|
||||||||
между числом пересечении при по |
|
||||||||
ложительном значении |
производной |
|
|||||||
фазы |
ф'(/) |
и |
числом |
пересечений |
|
||||
при отрицательном ее значении (на |
|
||||||||
рис. |
1 |
пересечения |
первого |
рода |
обозначены на |
||||
годографе вектора 1} точками 1. 2, 3 |
и 4, а второго ро |
да —-точкой 5).
С точностью до единицы эта разность определяется приращением фазы ф за время Т, поделенным на вели
чину шага между уровнями, т. е.
7'/2
Вычисляя отношение второго момента от левой и пра вой частей равенства (1) к интервалу времени Т и пе реходя к пределу при Т -»■оо, можно показать, что
спектральная плотность производной фазы ср^ на йулё-
вой частоте |
9д2 |
(2) |
|
/V (0) = lim 1А <(Д /г)2> . |
|
|
г^оо у |
|
Здесь < > |
означает усреднение по множеству. |
|
Аналогичная методика использовалась в работе [3], где рассматривались нули процесса £(f)cos(p(/), что эк
вивалентно |
пересечениям фазой <|> уровней ±я/2. |
Не нарушая общности рассуждений, можно ограни |
|
читься при |
вычислении iNr областью однозначности ф |
при векторном представлении процесса т} (t). Будем так
же полагать |
. |
0 |
,, |
,о\ |
где к — произвольное целое число. |
содер |
|||
Отметим далее, |
что |
если |
производная фазы |
жит постоянную составляющую, то только в этом слу чае среднее значение не равно нулю н на нулевой час
тоте спектра |
производной фазы |
ей соответствует дис |
|
кретная составляющая спектра. |
|
||
Нулевая |
составляющая |
непрерывной части спектра |
|
производной |
фазы F^ (0) |
определяется дисперсией |
|
разности числа пересечений а Ь |
Отсюда, учитывая (2) |
т
и (3), можно записать, что
2*2
(0) = Iim
Т-+ OU Ttf
N T • |
(4) |
Покажем, что рассматриваемая методика определения Fip' (0) охватывает известные модели взаимодействия
сигнала с шумом.
Предположим, что моменты времени пересечений фа* зой каждого отдельного уровня статистически независи мы между собой. При этом дисперсия разности пересе чений каждого n-го уровня otfT будет равна среднему
значению общего числа пересечений Мпт этого уровня
при положительном и отрицательном значениях произ водной фазы ввиду известного свойства пуассоновского распределения:
с1т = < М „ т > , я = 1, 2, 3, .... к.
Кроме того, из (4), которое справедливо при любом значении к, следует, что aNT — к2 а7пТ , т. е. факт пересе
чения фазой одного из уровней с достоверностью, равной
единице, фиксирует пересечение ею всех остальных к— 1 уровней (на рис. 1 годограф вектора YJ изображен
сплошной линией). Действительно, в противном случае фаза должна пересечь данный уровень в противополож ном направлении (на рис. 1 годограф вектора rj изобра
жен штриховой линией). Однако обратное пересечение
уровня, являясь следствием прямого, |
должно |
зависеть |
|||
от него, что исключено принятым допущением. |
|||||
В |
результате мы получили совокупность |
независи |
|||
мых |
скачков фазы |
на 2 я, |
что характеризует |
импульс |
|
ную |
составляющую |
шума |
в модели |
Райса. |
С другой |
стороны, очевидно, что для совокупности независимых
скачков фазы |
на 2 п моменты времени пересечения фа |
зой любого уровня статистически независимы. |
|
В итоге для |
(0) согласно (4) имеем |
/V (0) = 8w3 < М „ > ,
где <Л 4„> — среднее число пересечений н-го уровня в единицу времени.
Очевидно, что условие независимости пересечении может выполняться с большой вероятностью при доста точно хорошем отношении сигиал/шум, если положить Æ=l, а значение уровня принять равным y s + л, где q>s —
фаза сигнала. Таким образом, мы окончательно пришли к модели, предложенной Рансом [1].
Заметим, что в случае расстройки или модуляции по частоте сигнала в модели Рапса уровень, пересечения которого фиксируются, должен изменяться с угловой скоростью, равной расстройке сигнала по отношению к соо, что вносит существенную дополнительную погреш ность при достаточно малом отношении енгнал/шум в расчетах с использованием формул Райса. Действитель но, с уменьшением мощности сигнала до нуля среднее число пересечении этого уровня зависит от расстройки сигнала (остается «память» о сигнале), в то время как очевидно, что влияние модуляции сигнала на флуктуа ции частоты должно исчезать с уменьшением отношения сигнал/шум.
Вернемся к рассмотрению зависимых пересечений. Как показано выше, число уровней на интервале одно значности и их значения могут быть произвольными. Отсюда следует, что/7®' (0 ) и в этом случае определяет
ся скачками, или, точнее говоря, приращениями фазы
на 2 л, но уже с учетом статистической зависимости
между моментами их появления. Заметим, однако, что удачный выбор уровней пересечения в исходной модели может заметно упростить математические выкладки и окончательные выражения для /V (0 ), расчет по кото
рым приводит к одинаковым результатам.
В работе >[2] рассматривается модель, |
где |
фикси |
||
руются |
пересечения фазы <р с |
уровнем, равным |
я при |
|
Д= 2л. |
Момент пересечения |
этого уровня |
и |
знак ф' |
определяется с помощью произведения дельта-функции и ступенчатой функции соответственно от ортогональ ной и противофазной составляющей суммарного про цесса по отношению к сигналу. Это дало возможность учесть статистическую зависимость моментов времени пересечений. Однако окончательные результаты полу чились довольно громоздкими. В работе {3] рассматри вались нули процесса Tj(f), что, как было показано выше, эквивалентно рассмотрению пересечений фазой ф
уровней ± я / 2 при Д =я. |
|
N г выражается через ин |
|||
|
В этом случае значение |
||||
теграл Стильтьееа [3]: |
|
|
|
|
|
|
|
Г/2 |
|
|
|
|
\NT = — |
j* |
sgn Y d sgn X, |
(5) |
|
|
|
- T p t |
|
|
|
где |
Y = E{t) sin <p(t), |
|
X = E (i) cos <p (i). |
||
|
Спектральная плотность производной фазы на нуле |
||||
вой частоте при этом согласно (3] будет равна |
|||||
|
|
|
Г/2 7"/2 |
|
|
|
/V(0) = Hm |
|
j |
J sgnr^gnKaX |
|
|
|
-Г /2 |
-772 |
|
|
|
х - - (да/ |
: , г |
|
|
<б> |
|
ut\ ut2 |
|
|
|
|
где |
Х> = Х(Ь), У г= У (^ ), |
X 2 = X ( t 2), |
Y 2= Y ( t 2). |
||
Некоторая модификация |
последней |
модели позво |
ляет существенно упростить математические преобразо вания и в итоге получить новые результаты.
Запишем |
выражение (6 ) для /V |
(0) в несколько |
ином виде, для чего проинтегрируем |
его по частям по |
|
переменному |
В результате имеем |
|
Здесь |
(0 ) определяется |
пределом |
от смешанного |
|||
второго |
момента разности |
числа |
пересечений |
уровней |
||
0 , л и уровней ± л /2 . |
|
|
|
|
|
|
Рассмотрим для примера случай иемодулированного |
||||||
сигнала |
с расстройкой относительно |
средней |
частоты |
|||
энергетического спектра шума |
U> = IÙS — ш01 |
где |
— |
частота сигнала. Будем полагать, что шум имеет нор мальное распределение с параметрами (0, а). В этом случае
X (t) = |
Ас ( 0 + Ат cos Ш , |
Y (t)=- A s-\-Am sin Ш , |
||||||||||
где |
Ат— амплитуда |
сигнала; |
|
Ас и As — квадратур |
||||||||
ные составляющие узкополосного шума. |
и As(t) |
вза |
||||||||||
Так |
как случайные |
процессы |
Ac(t) |
|||||||||
имно независимы, то |
(7) |
можно записать в виде |
|
|||||||||
|
|
|
|
|
Г/2 |
772 |
|
|
|
|
|
|
|
|
|
|
|
Г/2 -772 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(8) |
Подынтегральное |
выражение |
в |
(8 ) |
представляет |
||||||||
собой |
произведение |
производных |
от корреляционных |
|||||||||
функций |
случайных |
процессов |
после |
нелинейного |
||||||||
sgn-преобразования |
нормальных |
случайных процессов |
||||||||||
Y (t) |
и X(t). Оба сомножителя |
|
зависят |
от |
времени |
|||||||
только через коэффициент корреляции R(h—U) |
и сред |
|||||||||||
ние |
значения & i=/l/;isin Aw/i |
и |
Д2 =Л ,„cos Дто/ 2 |
процес |
||||||||
сов |
Y (J) |
|
и X(t) соответственно. |
Это |
позволяет для |
оп |
ределения /V (0 ) эффективно использовать так назы
ваемый метод производной вычисления корреляционной функции.
К каждому сомножителю в (8 ) можно применить
правило дифференцирования сложной функции двух
Переменных [4], при этом, применяя указанный метод, достаточно просто найти частные производные sg.n-фуик- ции по Д, b! и агДальнейшие выкладки довольно гро
моздки. Окончательный |
результат преобразования (8 ) |
|||||||||||||
можно представить в виде суммы |
|
/ у |
(0 ) = |
/ у |
- (0 ) -f- |
|||||||||
- j - ^ v ( 0 )> |
гДе Дг¥- (0 ) |
— дискретная |
|
составляющая |
||||||||||
спектра ф'. Непрерывная часть ify (0 ) |
будет |
|
|
|||||||||||
|
|
|
Г |
д ) ' 2 |
у |
|
\-\-R9~2R cos àoix |
|
|
|||||
/у..(0)!= 2е~р |
R'2 |
|
|
1-R* |
|
|
|
|||||||
I ------te |
|
|
|
|
|
|
|
|||||||
|
|
|
J |
1 - Л 2 |
|
|
|
|
|
|
|
|||
4 А шe-p |
I | 1 |
|
|
/?'(sin Am |
|
X |
|
|
||||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
—{- / ? 2 — 2R cos Am |
|
|
|||||||||
|
(1 — e |
_ |
1+RЛ 91—2R/? cCOSo s |
A<OTo |
\ |
|
|
|
|
|||||
|
|
P |
|
|
|
» |
|
/<# + |
|
|
||||
-f- 2 Дш2?е- <’ |
Г Г1 !— e -p -------------------------------X |
|
||||||||||||
‘ |
J |
[ |
' |
|
|
|
H - A » 2 - 2 / |
? |
C O S A«>T |
|
||||
|
— OPI |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
/ |
|
1 + ^ - 2 ^ |
cos Amt \П |
|
|
|
|
|
|||||
X ( l - e " P ■-** |
|
|
) ] * , |
|
|
|
0 ) |
|||||||
где т=^2—1\. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Более подробный |
вывод |
соотношения |
(9) |
приведен |
||||||||||
в работе [5]. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Результат для |
/ у |
(0). |
полученный |
в |
[3], |
представ |
||||||||
ляет частный случай |
(9) |
при Д'<а = 0 . В этом случае |
||||||||||||
/V |
(0 ) = 2 е-е |
|
Г —ff- |
|
е_1> •+* <*с. |
|
|
|||||||
|
|
|
|
J |
1 |
- Я |
2 |
|
|
|
|
|
|
|
Был проведен расчет по формуле (9) для гауссовой |
||||||||||||||
характеристики фильтра, |
т. е. при R (т) = e_'Л'5,'’ , |
где |
||||||||||||
г — радиус инерции фильтра. Результаты |
расчета, пред |
|||||||||||||
ставляющие |
собой |
зависимость |
|
от |
р |
нормированной |
||||||||
спектральной |
плотности |
Д(0 )/П |
для |
различных |
значе |
|||||||||
ний нормированной |
расстройки |
у — / ± ы / 2 п г, |
представ- |