УДК 621.391:621.396.96
Д. И. Попов
Рязанский государственный радиотехнический университет
| Инвариантная обработка многочастотных сигналов
Синтезирован алгоритм обработки многочастотного сигнала, инвариантный к значениям допле-ровских сдвигов фаз компонент многочастотного сигнала. Проведен анализ характеристик обнаружения системы, реализующей разработанный алгоритм, позволяющий определить теоретический предел усовершенствования реальных систем данного класса и направления поисков новых систем.
Доплеровская фаза, многочастотные сигналы, обработка сигналов, пассивные помехи
Обработка многочастотных сигналов на фоне коррелированных (пассивных) и некоррелированных помех рассмотрена в [1], [2]. Априорная неопределенность статистических характеристик сигналов и помех приводит к адаптивному построению систем обработки. При этом ввиду значительного превышения сигнала помехой адаптация обычно осуществляется только к параметрам помехи. Оптимальная обработка многочастотного сигнала в этом случае в каждом частотном канале реализуется на основе адаптивного матричного фильтра (АМФ) (к параметрам помехи) и неадаптивного многоканального (по доплеровской фазе сигнала) фильтра, вычисляющего дискретное преобразование Фурье выходных отсчетов АМФ [1]. Принципиальные трудности адаптации к параметрам сигнала по данным исходной выборки преодолеваются в случае использования выходных отсчетов АМФ. Достаточно эффективное подавление помехи в результате матричной обработки открывает возможности для адаптивного накопления сигнала, позволяющего избежать традиционного многоканального по доплеровской фазе сигнала построения системы обработки в каждом частотном канале.
В работе [2] синтезирован квазиоптимальный алгоритм оценивания и соответствующий ему измеритель доплеровской фазы многочастотного сигнала по выходным отсчетам АМФ в каждом частотном канале. Исследование свойств получаемых оценок показало возможность их использования в системах обработки с адаптивным накоплением сигнала, позволяющим сократить число доплеровских каналов в каждом частотном канале или при прежнем числе доплеровских каналов уменьшить расстройку между каналами, исключив межканальные потери.
Требования дальнейшего упрощения систем обработки приводят к задаче синтеза более простых
8
систем, одноканальных по доплеровской частоте (фазе) сигнала. Представляет интерес синтез алгоритма обработки многочастотного сигнала, инвариантного к значениям доплеровских сдвигов фаз компонент многочастотного сигнала. Исследование характеристик системы, использующей такой алгоритм, позволит определить теоретический предел усовершенствования реальных систем данного класса и направления поисков новых систем.
Цель настоящей статьи - синтез и анализ алгоритма обработки многочастотного сигнала, инвариантного к доплеровским сдвигам фаз компонент многочастотного сигнала.
Синтез алгоритма обработки. Рассмотрим обработку Ь компонент многочастотного сигнала, каждая из которых представлена последовательностью N цифровых отсчетов ип1 = хп1 + гу„1,
п = 1, N, I = 1, Ь, комплексной огибающей аддитивной смеси сигнала, пассивной (коррелированной) помехи и собственного шума. Отсчеты следуют через период повторения Т и образуют в одном элементе разрешения по дальности совокупность вектор-столбцов ип = {ип1 }т ( т - символ
транспонирования). Сигнал и помеха являются узкополосными случайными процессами гауссов-ского типа, статистически независимыми в частотных компонентах, что достигается разносом их несущих частот, выбираемым из условия малости длин волн, соответствующих разностным частотам, по сравнению с радиальными размерами цели. Статистические свойства совокупности {и,..., иь т с точностью до параметров корреляционных матриц Щ векторов иI описываются гауссовской совместной плотностью вероятности [1].
Алгоритм оптимальной обработки Ь частотных компонент, состоящих из N цифровых отсче-
© Попов Д. И., 2015
тов Uni (n = 1, N, l = 1, L), определяется на основе отношения правдоподобия, вычисление которого при использовании упомянутого типа плотностей вероятности сигнала и помехи и одной помехи приводит к минимальной достаточной статистике [1]
L
V(фсЬ Фс2, • • •, 9oL) = XV(Фс1 ) = l=1
N , N 2
L
=x
l=1
X e~ik Фс X Wt(n, k)Uni
(1)
к=1 и=1
где фс/ - доплеровский сдвиг фазы сигнала за период повторения Т; (и, к) = м>п (и, к)е'(и-к)фп1 -элементы обратной корреляционной матрицы помехи, причем фп/ - доплеровский сдвиг фазы
помехи за период повторения Т; * - символ комплексного сопряжения.
Алгоритм (1) описывает оптимальную обработку многочастотного сигнала. Внутренняя сумма алгоритма (1) соответствует матричной фильтрации групп из N отсчетов. Преодоление априорной неопределенности параметров помехи основывается на адаптивном байесовском подходе, в соответствии с которым неизвестные величины W (и, к) или м>1 (и, к) и фп/ заменяются их состоятельными оценками. Неопределенность величин фс1 в доплеровском интервале [-л, л]
первоначально предполагает многоканальное когерентное накопление результатов матричной фильтрации. При этом сигнал от движущейся цели из-за различия доплеровских сдвигов фазы частотных компонент будет попадать в различные доплеровские каналы каждого из когерентных накопителей, что исключает объединение выходных величин последних в соответствии с алгоритмом (1). Другим вариантом преодоления априорной неопределенности доплеровских сдвигов фазы сигнала является усреднение алгоритма (1) по данным параметрам. Полагая величины фс1 равномерно распределенными в интервале [-л, л], в результате усреднения (1), исключающего неопределенность этих величин в пределах указанного интервала, находим
V = (2л) Ь х
л л л
х I | ••• | ¥(фсЬ Фс2,•••, ФсЬ)Фс1 ЛФс2 •• •dФсL =
—л —л —л
L я N , ,
= х(2л)-1 j x e~'(k"p)фс1 x
l=1 -л k, p=1
N
x x W*(n,k)UnlWl (r,p)U*rl. n, r=1
Изменяя порядок интегрирования и суммирования, а также учитывая, что
л
(2л)-1 j e-1 (k-p)фсldфсl = sine[(k-p)я] =
sin [(k - p)л] = [1, k = p,
(k - p) л [0, k ф p,
имеем
L N
v= xx
l=1 k =1
N
x W* (n, k)Unl
n=1
(2)
Полученный алгоритм определяет структуру системы обработки, инвариантной к доплеров-ским сдвигам фазы компонент отраженного многочастотного сигнала. При этом межпериодная обработка отсчетов каждой компоненты является комбинированной, т. е. в условиях априорной неопределенности спектрально-корреляционных характеристик помехи распадается на адаптивную когерентную матричную фильтрацию отсчетов и последующее некогерентное суммирование N результатов матричной фильтрации. Завершается обработка суммированием Ь результатов раздельной обработки отсчетов каждой частотной компоненты.
Весовыми коэффициентами адаптивного матричного фильтра являются оценки элементов обратной корреляционной матрицы помехи, вычисление которых в условиях априорной неопределенности в общем случае представляет собой трудоемкую процедуру, усложняющую реализацию алгоритма обработки. При марковских аппроксимациях гауссовской помехи матричный фильтр преобразуется в векторный (одноканаль-ный) адаптивный режекторный фильтр [2].
Адаптивный режекторный фильтр (АРФ) в случае произвольных корреляционных свойств помехи может выполняться с комплексными весовыми коэффициентами [3], определяемыми с помощью адаптивных алгоритмов по максимально правдоподобным оценкам корреляционных параметров помехи [4]. При реализации данных АРФ в цифровом виде предполагается использование комплексных перемножителей, число кото-
2
рых равно порядку фильтра. При этом существенно усложняется структура АРФ, особенно высоких порядков, и повышаются требования к быстродействию арифметических операций для выполнения обработки в реальном масштабе времени. Избежать указанных трудностей можно предварительной компенсацией доплеровского сдвига фазы помехи, обусловленного взаимным перемещением источника мешающих отражений и носителя радиолокатора. В [5] синтезированы алгоритмы оценивания и предложены принципы построения и структурные схемы автокомпенсаторов доплеров-ской фазы пассивных помех с прямой и обратной связями. Режектирование "остановленной" помехи теперь может быть осуществлено фильтром с действительными весовыми коэффициентами, адаптирующимися к корреляционным свойствам помехи на выходе автокомпенсатора [6]-[9].
Анализ системы обработки. Для определения характеристик обнаружения системы, соответствующей синтезированному алгоритму инвариантной обработки, необходимо найти закон распределения статистики V . С этой целью алгоритм (2) представим в виде квадратичной формы
V=!>=Х и и =
/=1 /=1
Ь N
= ЕЕ й (п,к)ип11ик1, (3)
I=1 п,к =1
где - матрица обработки 1-й частотной компоненты, элементы которой имеют вид
N
й, (п, к )=Х Щ (к, п )Щ>, р).
п=1
В соответствии с интерпретацией алгоритма (3) система инвариантной обработки должна осуществлять весовое суммирование всех возможных комбинаций попарных произведений поступающих отсчетов каждой частотной компоненты. Весовые коэффициенты являются элементами матриц обработки , определяемыми по элементам обратных корреляционных матриц помехи, параметры кото -рых при анализе полагаются известными.
Универсальную методику анализа в рассматриваемом случае дает метод характеристических функций [10]. Характеристическая функция величины V определяется следующим образом:
- Ь - Ь ю ю
©V (И) = е'* =П=П ! ... ! Р(и,;
I=1 I =1 —ю —ю
(и, = й^цйи2, ...^N1.
Используя плотности вероятности Р (и,) из работы [1] и величины V, из алгоритма (3) настоящей статьи, найдем:
©V ('') = (2л)— ш П{ае1 Щ: ,=1
| ... | ехр
—2 иГ (( — 2ИЙ1 )и,
(и,
где Щ - матрица, обратная корреляционной матрице Я, вектора и,. Учитывая, что [11]
ю ю / 1 N
| ... | ехр I— --и*тАи] (и = (2л)^ае1 А
—ю —ю
и аег щ = (аег я) окончательно получим
Ь
—1
©V (''') = п {^ щ ^ (( — 2щ) 1т =
Ь
= п^ (/ — 2пяй ) , ,=1
где I = ЩЩ - единичная матрица.
Искомая плотность вероятности определяется при помощи преобразования Фурье
ю
р(V)=^ г ©V(и)е~Ш(И, 2л
(4)
вычисление которого предполагает приведение определителя det (I — 2ИЯ,Й1) в подынтегральном выражении к виду, удобному для интегрирования. С этой целью используем метод собственных значений [10], позволяющий представить характеристическую функцию в виде
Ь Ь N —1
©v (П) = п ©v, 0'') = п п (1 — 2л'ап1 (5)
,=1 ,=1 п =1
где ап, - собственные значения матриц Я,Й,,
п = , , = ттъ.
Интегрированием в (4) с использованием метода вычетов и с учетом (5) находится плотность вероятности р (V). При этом следует учитывать,
что собственные значения ап, = ап, , = 1, Ь, т. е. являются кратными, причем кратность числа ап равна Ь. Выражение для вероятности превыше-
ю
ния порогового уровня обнаружения ^ статистикой V принимает вид
М 1
р ( > V0 )= I р (v) ¿V = Х( Ь —1)! х
vn ]=1
Л
Ь—1
Л а
Ь—1
Ь—1
а и
ехр
а
^ N ( а, ГЬ
п 1 -ак
и У к=1
к Фи
а
и У
(6)
где М - число различных положительных собственных значений матриц Я/0/.
Использование в выражении (6) собственных значений матриц Я/0/ приводит к вычислению вероятности ложной тревоги Е, а собственных значений матриц
Ясп/ (фс1) 01 = [<?Яс1 (фс1) + Яп/ ] 01
(д - отношение "сигнал/помеха") - к вычислению вероятности правильного обнаружения Б.
Результаты расчетов. При расчетах использовались совместные флуктуации сигнала и гаус-совская корреляционная функция помехи
Р1 к) = ехр{-[лрг/ ( - к)] 2Дв}, где Р = Д/Т - нормированная ширина спектра помехи в первом частотном канале; Г1 = ////1 < 1 -
отношение несущих частот 1-го и первого каналов. Характеристики обнаружения систем инва-
-з
риантной обработки при N = 10, Е = 10 и Р = 0.15 приведены на рисунке. Сплошные кривые соответствуют двухчастотной системе (Ь = 2) при г/ = 0.9 и различных значениях доплеровско-го сдвига фазы сигнала в первом частотном кана-
Б
0.75
0.50
0.25
0
_1_
_1_
_1_
-60
-40
-20
0
д, дБ
ле фс1 = фс. Штриховая кривая соответствует одночастотной системе (Ь = 1) при доплеровском сдвиге фазы сигнала для "слепых" скоростей цели, т. е. фс = +2лк , к = 0, 1, 2, • . Из представленных зависимостей следует, что предельная эффективность двухчастотной системы наблюдается при оптимальной скорости цели в одном (первом) частотном канале (фс1 =фс =л). Однако и на бывших "слепых" скоростях выигрыш двухчастотной системы в значении порогового отношения "сигнал/помеха" д по сравнению с одночастотной системой достигает в зависимости от номера "слепой" скорости существенных значений (до нескольких десятков децибел).
Таким образом, в результате статистического синтеза получен алгоритм обработки многочастотного сигнала на фоне пассивных помех, инвариантный к доплеровским сдвигам фаз его компонент. Проведенный анализ соответствующей данному алгоритму системы позволяет установить предельные возможности обнаружения сигнала от цели для этого класса систем и подтверждает эффективность использования систем обработки многочастотных сигналов для борьбы со "слепыми" скоростями цели.
СПИСОК ЛИТЕРАТУРЫ
1. Попов Д. И. Оптимальная обработка многочастотных сигналов // Изв. вузов России. Радиоэлектроника. 2013. Вып. 1. С. 32-39.
2. Попов Д. И. Адаптивная обработка многочастотных сигналов // Изв. вузов России. Радиоэлектроника. 2013. Вып. 6. С. 15-20.
3. Попов Д. И. Синтез и анализ эффективности систем адаптивной междупериодной обработки сигналов на фоне помех с неизвестными корреляционными свойствами // Радиотехника и электроника. 1983. Т. 28, № 12. С. 2373-2380.
4. Попов Д. И. Оценивание параметров пассивных
помех // Изв. вузов. Радиоэлектроника. 2003. Т. 46,
№ 3. С. 71-80.
5. Попов Д. И. Автокомпенсация доплеровской фазы пассивных помех // Цифровая обработка сигналов. 2009. № 2. С. 30-33.
6. А. с. 711849 СССР, МПК6 601Б7/36, 601Б13/52. Устройство для подавления пассивных помех / Д. И. Попов; опубл. 27.11.98. Бюл. № 33.
7. А. с. 875960 СССР, МПК6 601Б7/36, 601Б13/52. Устройство для подавления пассивных помех / Д. И. Попов; опубл. 27.11.98. Бюл. № 33.
8. А. с. 1015757 СССР, МПК6 601Б7/36. Устройство подавления пассивных помех / Д. И. Попов; опубл. 27.11.98. Бюл. № 33.
9. А. с. 1098399 СССР, МПК6 G01S7/36. Устройство адаптивной режекции пассивных помех / Д. И. Попов; опубл. 20.12.98. Бюл. № 35.
D. I. Popov Ryazan State Radio Engineering University
10. Миддлтон Д. Введение в статистическую теорию связи: в 2 т. М.: Сов. радио, 1961-1962. Т. 1. 1961. 782 с.; Т. 2. 1962. 832 с.
Invariant Processing of Multifrequency Signals
The multifrequency signal processing algorithm that is invariant to the values of Doppler phase components of multifrequency signals have been synthesized. The analysis of the detection characteristics of the respective system conducted that allows to determine the theoretical limit of improvement of real systems of this class, and the direction of the search for new systems.
Doppler phase, multifrequency signals, signals processing, clutter Статья поступила в редакцию 30 ноября 2015 г.
УДК 681.514
С. И. Зиатдинов
Санкт-Петербургский государственный университет аэрокосмического приборостроения
Анализ ошибок узкополосного приема частотно-модулированных сигналов
Показано, что ограничение полосы пропускания при приеме частотно-модулированного сигнала приводит к нежелательным изменениям как амплитуды, так и частоты выходного сигнала, значения которых зависят от значения индекса угловой модуляции.
Частотная модуляция, спектр сигнала, узкополосный прием, ошибки
Частотная модуляция (ЧМ) относится к распространенному методу передачи сообщений по информационному каналу. Исследованиям искажений частотно-модулированного колебания при прохождении через колебательную систему посвящено множество работ (см., например, [1]-[3]). Однако в данных работах с целью недопущения динамических искажений, связанных с конечной полосой пропускания линейной системы, рассмотрен лишь случай медленного изменения частоты входного сигнала.
Напротив, в настоящей статье использован спектральный метод анализа прохождения ЧМ-сиг-нала через линейную систему, который можно без ограничений применять как при медленных, так и быстрых изменениях частоты.
ЧМ-сигналы характеризуются столь большим числом спектральных составляющих в используемой полосе частот, что применение спектрального метода сопряжено с большими, иногда непреодолимыми трудностями вычисления. Однако современные вычислительные средства полностью снимают все трудности обработки большого числа спектральных составляющих в выходном сигнале. 12
При ЧМ частота несущего сигнала ю(t) изменяется по закону передаваемого сообщения 5 (t): ю (t) = ю0 + ks (t), где ю0 - частота немо-дулированного несущего сигнала; k - коэффициент пропорциональности.
При этом полная фаза ЧМ-сигнала записывается следующим образом:
t
0 (t ) = J ю (l) dl = o>0t + Д0 (t), (1)
0
t
где Д0 (t) = j ks (l) dl =Дф^) + фо, причем Дф(t) -
0
изменение фазы несущего сигнала вследствие ЧМ; фо - начальная фаза.
В результате для ЧМ несущего сигнала можно записать выражение
u (t) = Um cos [^ot + Дф() + Фо ], (2)
где Um - амплитуда сигнала.
© Зиатдинов С. И., 2015