учитывает базовый алгоритм что, в частности, позволяет компенсировать искажения, которые дает базовый алгоритм при повышении разрешающей
СПИСОК ЛИТЕРАТУРЫ
1. Goodman R., Carrara W. Synthetic Aperture Radar Algorithms // Handbook of image and video processing. - London: Academic press, 2000. - P. 749-770.
2. Raney R. K. Radar fundamentals: technical perspective // Principles and Applications of Imaging Radar. Manual of Remote Sensing. V. 2. Eds. F.M. Henderson, A.J. Lewis. - N.Y.: John Wiley & Sons, 1998. - P. 9-130.
3. Хлусов В.А., Лигтхарт Л.П., Шарыгин Г.С. Одновременное измерение всех элементов матрицы рассеяния радиолокацион-
способности [1]. Предложен также математический аппарат для оценки потенциальной разрешающей способности РСА.
ных объектов с использованием сложных сигналов // Радиолокация, навигация, связь: Труды VIII Междунар. научно-техн. конф. - Воронеж: Воронежский госуниверситет, 2002. - Т. 3. -С. 1655-1667.
4. Рытов С.М., Кравцов Ю.А., Татарский В.И. Введение в статистическую радиофизику. Ч. 2. Случайные поля. - М.: Наука, 1978. - 463 с.
5. Фалькович С.Е. Оценка параметров сигнала. - М.: Советское радио, 1970. - 336 с.
УДК 621.396.96
АДАПТИВНЫЙ ФИЛЬТР ПОЛЯРИЗАЦИОННОЙ РЛС СО СЛОЖНЫМИ СИГНАЛАМИ
Г.П. Бабур
Томский государственный университет систем управления и радиоэлектроники E-mail: [email protected]
Представлен адаптивный фильтр поляризационной радиолокационной станции для снижения уровня боковых лепестков сжатого принимаемого сигнала при использовании сложного векторного зондирующего сигнала. Обработка принимаемого сигнала производится в три этапа. Снижение уровня боковых лепестков позволяет дать более точную по сравнению с матричным согласованным фильтром оценку элементов матрицы рассеяния точечных отражателей в составе радиолокационной цели путем изменения в процессе адаптации векторов весовых коэффициентов предлагаемого адаптивного фильтра.
Введение
Радиолокационная цель зачастую представляется набором блестящих точек (точечных отражателей). Необходимо оценить на оси задержек (дальности) их расположение и интенсивность отраженных от них сигналов. В случае поляризационной радиолокационной станции (РЛС) интенсивность отраженных сигналов описываются матрицей рассеяния, четыре комплексных элемента которой характеризуют прямые и перекрестные отражения волн с ортогональными поляризациями [1-3].
Во многих современных РЛС используются сложные зондирующие сигналы, например, фазокодома-нипулированные. Использование сложных сигналов по сравнению с простыми имеет ряд преимуществ: достижение высокого энергопотенциала РЛС при ограниченной мощности излучения, повышенная помехозащищенность и др. Однако, после сжатия сложный сигнал имеет не только пик (информативную часть), но и так называемые боковые лепестки. В случае поляризационной РЛС боковые лепестки обуславливаются не только видом автокорреляционных функций используемых ортогональных сигналов, но и их взаимной корреляцией. Уровнем боковых лепестков (УБЛ) сжатых принимаемых сигналов определяются основные ошибки при оценке матриц рассеяния точечных отражателей в составе радиолокационной цели.
Традиционный способ снижения УБЛ состоит в увеличении длительности используемых сложных сигналов, что бывает не всегда оправдано [3]. Встает задача снижения уровня боковых лепестков при заданной, фиксированной длительности сигналов.
Цель настоящей работы - разработка фильтра, позволяющего снизить УБЛ сигнала с выхода согласованного матричного фильтра поляризационной РЛС при фиксированной длительности (базе) зондирующих сигналов с целью измерения матрицы рассеяния радиолокационных целей.
Описание алгоритма работы адаптивного фильтра
поляризационной РЛС
Элементы матрицы рассеяния в общем случае между собой не зависимы, то мощности прямых и кроссовых компонент сигналов, отраженных от точечных отражателей в составе радиолокационной цели, как правило, между собой не равны, а радиолокационная цель может быть любой. Поэтому невозможно разработать фильтр с фиксированной структурой, который будет оптимален для различных радиолокационных целей, представленных разными совокупностями точечных отражателей. В этом случае имеет смысл разработать адаптивный фильтр с меняющейся структурой и определяющий его работу алгоритм.
Основой предлагаемого алгоритма послужил алгоритм адаптивного сжатия импульсов [4]. Данный алгоритм эффективен при подавлении боковых лепестков, обусловленных автокорреляционной функцией одного зондирующего сложного сигнала, и не работает при использовании векторного зондирующего сигнала в поляризационной РЛС.
Совместная оценка координатных и поляризационных параметров радиолокационного нестационарного во времени объекта возможна только при использовании векторного зондирующего сигнала [3, 5], поэтому зондирующий сигнал и(и)=и(и-Дт), Дт- расстояние между отсчетами, программно описывается в виде двух наборов отсчетов:
(1)
8 (/) =
У (/),
(2)
где 8 (/) =
1 .^11 Л 1 /¡12
11 / 21 Л 1 / 22
8 (/) =
где и1=[и1(1)и1(2)...и1(А0т и и2=[и2(1)и2(2)...и2(Л0т - А(/) = два ортогональных сигнала, состоящие из N отсчетов каждый. Формирование сигналов, описываемых выражением (1), основано на применении сложных ортогональных сигналов.
Адаптивный алгоритм для поляризационной РЛС, так же как и [4] реализуется в три этапа. В качестве первого этапа предлагается использовать согласованную фильтрацию, а именно согласованный матричный фильтр.
Матричный согласованный фильтр (СФ) обеспечивает оптимальный прием векторного сигнала на фоне белого шума, выполняя свертку принимаемого сигнала со своей импульсной характеристикой, соответствующей сопряженному векторному зондирующему сигналу. Описание матричного фильтра для поляризационной РЛС можно найти в работе [5].
Согласованную фильтрацию в матричной форме можно представить как
Принимаемый сигнал у(1) разделяется в антенном тракте РЛС на две ортогональные по поляризации составляющие: у1(/) и у2(/), которые могут быть записаны с помощью выражения (3):
[У 1(1) У 2(1)] = А (/Ж и 2] + №(/) Ч2(/)], (3) где . . . .
ъ(1)=ШШ+1)...ъ(М-1)]т, т(1)=[т(Ш+1)...т(М-1)]т;
- отсчеты аддитивного белого шума в каналах соответственно горизонтальной и вертикальной поляризации. .
Матрица Л(/) в выражении (4) представляет собой совокупность матриц рассеяния Ж-1 элементарных отражателей:
8 (/) 8(1 + 1) ... 8 (/ + N-1)
8 (/-1) 8 (/) '•. ;
; '•. '•. 8 (/ +1)
8(/ - N +1) ... 8(/-1) 8(1)
(4)
каждый элементарный отражатель описывается матрицей рассеяния, четыре комплексных элемента которой характеризуют прямые и перекрестные отражения волн с ортогональными поляризациями
^11 ^12
21 22
Итого отклик согласованного матричного фильтра:
8 (I) =
или
8 (/) =
• А • [и1 и 2 ] +
• [у 1(1) У 2(/)] +
• ШО |?2(/)] (5)
• Щ/) т)2(1)].(6)
оценка матрицы рассеяния
1-ой выборки (индекса ячейки по дальности); единица в верхнем индексе в выражении (2) означает первый этап обработки;
и1=[и1(1)и1(2)...и^Л^ и и2=[и2(1)и2(2)...и2(Л0т - импульсные характеристики фильтров, согласованных с двумя ортогональными сложными сигналами ц и и2, состоящими из N отсчетов каждый;
У(0=[У1(0 У2(0] - содержит два вектора длиной каждый по N отсчетов, выделяемые из непрерывного принимаемого сигнала;
У1(0=[У1(0У1(/+1)...У1(/+^1)]Т,
У2(0=[У2(/)У2(/+1)...У2(/+^1)]Т;
(•)Т - оператор транспонирования;
(•)Н - оператор сопряженного транспонирования.
Таким образом, сигналы на выходе согласованного матричного фильтра, входящего в состав адаптивного фильтра, в виде векторно-матричного умножения могут быть записаны как:
(7)
Единица в верхнем индексе в выражениях (5-7) означает первый этап обработки.
Опишем нахождение векторов весовых коэффициентов адаптивных фильтров во втором и третьем каскадах адаптивного фильтра поляризационной РЛС. Для вычисления векторов весовых коэффициентов адаптивных фильтров сжатия второго каскада используются сигналы с ¡выхода согласованного матричного фильтра ¿¡2, и ¿2, являющиеся
^ 1 5и = и" У 1(/)
Л 1 .12 = и2Н • У 1(/)
1 1 / 21 = и1Н У 2 (/)
Л 1 / 22 = и2Н У 2 (/)
и
оценками элементов матрицы рассеяния радиолокационного объекта. Полагая отсчеты по дальности принимаемого сигнала некоррелированными с шумом, векторы весовых коэффициентов четырех адаптивных фильтров, осуществляющих повторное сжатие входного сигнала, запишем как:
' ^12(/) = р112(С12(/))-'и,
^ ат122(/) = р122(С12(/))-1 и,
' w212 (/) =р212(С22(/))-1и2'
^ 222 (/) = р222(С22(/))-1и2
где ру2(1)=1$1(1)12 - параметр оценки элементов матрицы рассеяния радиолокационного объекта, ¿¿'=1,2.
Матрицы С12(/) и С22(/) определяются как:
N -1
С12(/) = X р112(/ + «Ж,пЦ,Н +
n=-N +1
N -1
+ X р122(/ + „)и2,„и2,",
n=-N +1
N-1
С22(/) = X р212(/ + „)и1,„и1,н„ +
n=-N +1
N -1
+ X р22,(/ + „)и2,„и2,Н,
n=-N+1
где и11Л состоит из элементов сигнала и1, смещенного на п отсчетов, и заполненного нулями остатка. Например, для п=2 имеем и12=[0 0 и1(0)...и1(Ы-3)]т, а для п=-2 и1-1=[и1(2)...и1(Ы-1) 0 0]т. Аналогично и2,„ состоит из элементов сигнала и2, смещенного на п отсчетов, и заполненного нулями остатка; -И<п<К
Сигналы на выходе второго каскада адаптивного фильтра записываются соответственно:
' ~ 2
,511 = w112н (/) • у!(/)
Л 2
^ = w122н (/) • у 1(/)
л 2
5 21 = w 212Н (/) • У 2 (/)
~ 2
522 = w 222Н (/) • У 2 (/)
Для нахождения векторов весовых коэффициентов фильтров сжатия третьего каскада адаптивного фильтра поляризационной РЛС используются уточненные оценки элементов матрицы рассеяния, полученные;, соответственно после второго этапа обработки £21, £22, и
' Wl1з(/) = р11з(С1з(/))-1и1 ^ w12з(/) =р12з(С1з(/))-Ч
' w213(/) =р213(С23(/))-1и2' ^ 223 (/) = р223(С23(/))-1и 2
где ргу3(/)=|£|(/)|2 - параметр оценки элементов матрицы рассеяния радиолокационного объекта, ¿¿=1, 2.
Матрицы С13(/) и С23(/) определяются как:
N -1
СШ = X р113(/ + „)Ul,nUl,яH +
n=-N +1
N -1
+ X р123(/ + п)и2,nU2,H,
п=-N+1
N-1
С23(/) = X Р213(/ + «М,nUl,Н +
n=-N +1
N -1
+ X р22ъ(/ + n)u2,nu2,Н,
n=-N +1
Сигналы на выходе третьего каскада адаптивного фильтра записываются как:
' ~ 3
,511 = w11зн (/) • у 1(/)
Л 3
^ ,512 = Wl2зн (/) • у 1(/)
~ 3
21 = w 213н (/) • У 2(/)
~ 3
5 22 = w 22" (/) • у 2(/)
Описание адаптивного фильтра поляризационной РЛС
Адаптивный фильтр поляризационной РЛС, рис. 1, реализующий адаптивный алгоритм, работает следующим образом. На входы фильтра (Вх. 1 и Вх. 2) поступают сигналы горизонтальной у1(/) и вертикальной у2(/) поляризации. Обработка сигналов производится последовательно в три этапа. На первом этапе происходит сжатие входных сигналов согласованными фильтрами (СФ 1 - СФ 4). В результате существенно увеличивается отношение сигнала к шуму, создаваемому в предыдущих при-емно-усилительных каскадах приемника поляризационной РЛС. Однако сжатые сигналы имеют боковые лепестки, определяющие основные ошибки при оценке элементов матрицы рассеяния радиолокационной цели в поляризационной РЛС.
Сигналы с выходов СФ 1 - СФ 4 используются для адаптации и расчета векторов весовых коэффициентов адаптивных фильтров, реализующих второй этап обработки входных сигналов у1(/) и у2(/). На втором этапе сигналы с выходов СФ 1 - СФ 4 поступают на входы первого и второго вычислительно-управляющих устройств (ВУУ), и используются для вычисления векторов весовых коэффициентов первого - четвертого адаптивных фильтров (АФ). С выхода ВУУ 1 и ВУУ 2 в каждый момент времени сигналы, отображающие вектора весовых коэффициентов фильтров, одновременно поступают на управляющие входы четырех адаптивных фильтров АФ 1 - АФ 4. Кроме максимизации отношения сигнал/шум фильтрами АФ 1 - АФ 4 осуществляется минимизация уровня боковых лепестков сжатого сигнала. На третьем этапе сигналы с выходов адаптивных фильтров АФ 1 - АФ 4 поступают на входы третьего и четвертого вычислительно-управляющих устройств, и используются для вычисления векторов весовых коэффициентов пятого - восьмого адаптивных фильтров. С выхода
1 каск. 2 каск. 3 каск.
Рис. 1. Схема адаптивного фильтра поляризационной РЛС
ВУУ 3 и ВУУ 4 в каждый момент времени сигналы, отображающие эти вектора, одновременно поступают на управляющие входы АФ 5 - АФ 8. Фильтрами АФ 5 - АФ 8 осуществляется дальнейшая минимизация уровня боковых лепестков сжатого сигнала. Сигналы с выходов АФ 5, АФ 6, АФ 7 и АФ 8 являются оценками элементов матрицы рассеяния радиолокационного объекта £ 3Ь £ 32, £ 21 и ^ 22.
Результаты моделирования
Результаты математического моделирования работы адаптивного фильтра поляризационной РЛС со сложными сигналами представлены на рис. 2. Радиолокационный объект описывался пятью отражающими элементами. Векторный зондирующий сигнал описывался двумя ортогональными М-по-следовательностями базой N=31. Матрицы рассеяния точечных отражателей задавались случайным образом в комплексном виде. Сигнал дополнительно содержал аддитивный белый шум с нулевым математическим ожиданием и среднеквадратическим отклонением - 40 дБ по отношению к уровню самого сильного сигнала. После сжатия сигнала среднек-вадратическое отклонение шума составляло - 55 дБ.
Результаты моделирования приводятся в логарифмическом масштабе и представлены в виде нормированных графиков абсолютного значения сигналов на выходах фильтров в составе адаптивного фильтра поляризационной РЛС (рис. 2). На рис. 2, а (сверху вниз) изображены сигналы с выходов СФ1 - СФ 4, на рис. 2, б - сигналы с выходов АФ 1 - АФ 4, на рис. 2, в - сигналы с выходов АФ 5 - АФ 8.
По оси абсцисс отложены отсчеты, соответствующие дальности до радиолокационного объекта (соответствующие задержке принимаемого сигнала относительно излученного) /=г/Дт, где / - текущее время, Дт - интервал дискретизации сигнала. За начало отсчета (/=0) принимается начало излучения зондирующего сигнала. Математическое моделирование работы предлагаемого фильтра осуществлялось в программе Ма1;ЬАВ 7.0.
Сигналы с выходов согласованного матричного фильтра (рис. 2, а), использующегося в качестве первого каскада адаптивного фильтра, имеют высокий уровень боковых лепестков, и слабые эхо-сигналы скрыты в боковых лепестках сильных эхо-сигналов. После второго этапа обработки (рис. 2, б) наблюдается существенное снижение УБЛ, обусловленного сильными эхо-сигналами. При наблюдении сигналов после третьего этапа обработки, с выхода адаптивного фильтра поляризационной РЛС (рис. 2, в), наблюдается снижение УБЛ сжатого сигнала, обусловленного слабыми эхо-сигналами. Для используемой модели снижение УБЛ составило 25 дБ.
Сравнение результатов работы согласованного матричного фильтра (рис. 2, а) и предложенного адаптивного фильтра (рис. 2, в) свидетельствует о существенном подавлении боковых лепестков сжатых сигналов в поляризационной РЛС со сложными сигналами. В результате адаптивной обработки становится возможна оценка слабых эхо-сигналов, которые при согласованной фильтрации скрыты в боковых лепестках сильных эхо-сигналов.
После первого этапа обработки
После второго этапа обработки
После третьего этапа обработки
Ш
-40
-60-
20-40-60-SO 100-I
20 40 60 ео 100 I
a .....
& Л
i
IL
-60-
(0 ...i-, J
.....
I
20 4 D 60 SO 100 f
■60
&
-40
20 40 60 80 100 /
ill Ii А
Iii
íl
III
20 40 во 80 100 /
20 40 60 SO " 100 I
Рис. 2. Результаты работы адаптивного фильтра поляризационной РЛС
Выводы
Предлагаемый адаптивный фильтр позволяет снизить уровень боковых лепестков сжатого сигнала в поляризационной РЛС, обусловленный как ненулевой автокорреляцией сигналов в составе векторного зондирующего сигнала, так и их взаимной корреляцией.
Снижение уровня боковых лепестков сжатого сигнала позволяет дать более точную по сравнению
с согласованным матричным фильтром оценку элементов матриц рассеяния точечных отражателей в составе радиолокационной цели без увеличения длительности зондирующего сигнала путем поэтапной обработки и изменения в процессе адаптации векторов весовых коэффициентов адаптивного фильтра. Для рассмотренной модели снижение уровня боковых лепестков сжатого сигнала составило 25 дБ.
СПИСОК ЛИТЕРАТУРЫ
1. Канарейкин Д.Б., Павлов Н.Ф., Потехин В.А. Поляризация радиолокационных сигналов. - М.: Советское радио, 1966. -440 с.
2. Богородский В.В., Канарейкин Д.Б., Козлов А.И. Поляризация рассеянного и собственного радиоизлучения земных покровов. - Л.: Гидрометеоиздат, 1981. - 279 с.
3. Губин В.А., Коростелев А.А., Мельник Ю.А. Пространственно-временная обработка радиолокационных сигналов. - Л.: Ле-
нинградская инженерная краснознаменная академия им. А.Ф. Можайского, 1970. - 201 с.
4. Blunt S.D., Gerlach K.. Adaptive pulse compression // IEEE 2004 National Radar Conf. - Philadelphia, 2004. - P. 271-276.
5. Хлусов В.А. Совместная оценка координатных и поляризационных параметров радиолокационных объектов // Сибирский поляризационный семинар «СИБПОЛ»: Докл. IV Междунар. научно-практ. семинара. - Сургут, 2004. - С. 67-79.