НАУКА и ОБРАЗОВАНИЕ

УДК.621396.6

НИИ РЭТ МГТУ им. Н.Э. Баумана

rodzivilov@bmstu.ru

golosov@bmstu.ru

ksa@bmstu.ru

Улучшить качество функционирования многоканальных РЛС при воздействии помех точечных источников, вынесенных с цели можно, используя методы компенсации [1, 2]. Эти методы достаточно хорошо описаны в литературе [3, 4, 5] и сводятся, в основном, к весовому суммированию сигналов основного и дополнительных (компенсационных) каналов, причем комплексные весовые коэффициенты должны выбираться таким образом, чтобы обеспечивать максимум отношения сигнал/(шум+помеха).  Если условия функционирования изменяются, то весовые коэффициенты должны быть адаптивными. В настоящее время наиболее широкое распространение получили итеративные алгоритмы адаптации весовых коэффициентов [3], применение которых в бортовых РЛС ограничено (поскольку на интервалах адаптации координаты источников излучения могут изменяться), а также алгоритмы адаптации на основе обращения выборочной ковариационной матрицы [3, 4, 5]. К недостатку последних следует отнести «рассыпание» диаграммы направленности вследствие влияния на выборочную ковариационную матрицу отражений от подстилающей поверхности, принимаемых по основному лепестку диаграммы направленности [3, 4].

В работе рассматривается алгоритм адаптации для импульсно-доплеровской БРЛС на основе непосредственного вычисления весовых коэффициентов [6], свободный от перечисленных выше недостатков. Анализируются характеристики когерентного компенсатора помех на его основе.

Работу когерентного компенсатора помех рассмотрим на следующем примере. Предположим, что бортовая импульсно-доплеровская РЛС имеет два приемных канала: основной канал с узконаправленной диаграммой направленности антенны (с коэффициентом усиления ) и компенсационный канал со слабонаправленной антенной (с коэффициентом усиления). Обнаруживаемая цель  находится на направлении, соответствующем максимуму основного лепестка диаграммы направленности антенны основного канала,  источник помехи (ИП)  смещен на угол  относительно цели (рисунок 1).

Рис. 1. Расположение источника помех и цели относительно БРЛС

Рассматриваемый когерентный компенсатор помех (рис. 2) включает основной (ОК) и компенсационный (КК) каналы, причем каждый из каналов содержит аналоговые модули с коэффициентами усиления , , цифровые приемные устройства (ЦПУ_о, ЦПУ_к), реализующие доплеровскую обработку, с выходов которых на процессор обработки информации (ПОИ) поступают вектора комплексных спектральных отсчетов ,  (размерность их одинакова и равна числу элементов Nразрешения БРЛС по частоте). В ПОИ по спектральным отсчетам ,  непосредственно вычисляется весовой коэффициент (один для двухканальной системы) и осуществляется компенсация помехи в основном канале в частотной области. Результат компенсации в каждом элементе разрешения по частоте отображается вектором  (размерность которого также равна N).

Рис. 2. Упрощенная структурная схема двухканального компенсатора помех

В дальнейшем будем использовать следующую модель сигналов, поступающих с выходов антенных систем основного и компенсационного каналов:

,

.

Здесь: , - соответственно комплексные огибающие полезного сигнала и сигнала помехи, приведенные ко входам приемных устройств и полученные при единичном коэффициенте усиления антенн;  - комплексные огибающие собственных шумов основного и компенсационного каналов соответственно, которые далее будем полагать независимыми, распределенными по нормальному закону, с нулевым математическим ожиданием и уровнями спектральной плотности  и  соответственно;  - комплексная величина, учитывающая различие в коэффициентах передачи сигнала помехи основной и компенсационной антенн:

,

где  и  - коэффициенты усиления антенн компенсационного и основного каналов в направлении на источник помехи.

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

Комплексные отсчеты  () на выходе анализатора спектра приемного устройства основного канала запишем в виде:

,                     (1)

где:  – номер элемента разрешения по частоте (доплеровского фильтра) ();  - число доплеровских фильтров; ,  - амплитуда и фаза сигнала, обусловленного помехой, на выходе -го доплеровского фильтра; - коэффициент, равный единице в том случае, если доплеровский сдвиг частоты полезного сигнала соответствует частоте настройки  - го доплеровского фильтра, и нулю в противном случае; ,  - амплитуда и фаза полезного сигнала на выходе -го доплеровского фильтра; ,  - амплитуда и фаза отсчета на выходе -го доплеровского фильтра, обусловленные собственным шумом (СШ).

В предположении, что коэффициент усиления антенны КК в направлении на источник полезного сигнала много меньше коэффициента усиления антенны ОК, сигнал на выходе анализатора спектра приемного устройства КК можно записать как:

,                                       (2)

где: ,  - амплитуда и фаза составляющей, обусловленной СШ.

Поскольку, как отмечалось выше, анализаторы спектра приемных устройств в основном и компенсационном каналах идентичны, помеху предлагается компенсировать в частотной области, вычитая из выходного отсчета -го доплеровского фильтра основного канала выходной отсчет аналогичного доплеровского фильтра компенсационного канала, умноженный на соответствующий коэффициент . Таким образом, составляющая  (рис.2) на выходе компенсатора будет равна:

.

Помеха будет скомпенсирована, если , для чего необходимо, чтобы . Как видно из (1), (2) оценку требуемого значения  можно получить, сравнивая отсчеты на выходах соответствующих доплеровских фильтров в основном и компенсационном каналах. Для того, чтобы уменьшить влияние собственных шумов и исключить возможное влияние полезного сигнала,  алгоритм компенсации помехи представлен в виде следующей последовательности шагов:

 1. Вычисление первого приближения оценки  ():         .

 2. Расчет коэффициента передачи:                                           .

 3. Выполнение предварительной компенсации помех:        .

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

                                                                          ,

где  - пороговое значение, которое должно обеспечить вероятность ложной тревоги при отсутствии сигнала порядка .

5. Вычисление оценки :                                                          .

6. Расчет результирующего коэффициента передачи:           .

7. Компенсация помех:                                                              .         (3)

Характеристики компенсатора анализировались методом математического моделирования. В частности, на рис. 3 приведен отклик приемника основного канала при воздействии широкополосной помехи (отношение сигнал/(шум+помеха) составляет –5 дБ); на рис. 4 - отклик приемника основного канала после компенсации (отношение сигнал/(шум+помеха) составляет 22.5 дБ, сигнал в 2720-ом доплеровском фильтре). Коэффициент подавление помехи составляет 27.5 дБ. На рисунках по оси абсцисс отложены номера доплеровских фильтров, по оси ординат – относительная амплитуда сигналов в каждом элементе разрешения.

Рис. 3. Отклик приемника основного канала при воздействии широкополосной помехи

Рис. 4. Отклик приемника основного канала после компенсации

На рис. 5 представлен график зависимости коэффициента подавления помехи , определенного как отношение средней мощности помехи на выходе анализатора спектра основного канала до и после компенсации, при разных отношениях мощности помехи к мощности собственного шума (П/СШ).

Рис. 5. График зависимости коэффициента подавления помехи от отношения П/СШ

Рассмотренный выше компенсатор позволяет компенсировать только один источник помех. Рассмотрим возможность применения предлагаемого алгоритма компенсации для систем с достаточно большим числом (M) приемных каналов (рис. 6) для компенсации помех Pисточников, находящихся на разных угловых направлениях.

 Далее, исходя из предложенной схемы,  будем рассматривать потенциальный максимум возможных помеховых воздействий - . На рисунке представлены приемные каналы , , …, , которые включают: антенные системы с диаграммами направленности , , …, ; аналоговые модули с коэффициентами усиления , , …,  и приведенными ко входам приемных устройств напряжениями , , …, ; цифровые приемные устройства (анализаторы спектра) , , …, , формирующие вектора комплексных спектральных отсчетов , , …, ; суммарный канал , в котором поэлементным суммированием всех векторов комплексных спектральных отсчетов вычисляется вектор ; ПОИ, в котором осуществляется компенсация помех в суммарном канале и формирование результирующего вектора спектральных отсчетов  после компенсации.

Рис. 6. Структурная схема М-канального компенсатора помех

Результирующие спектральные отсчеты  выходного вектора  компенсатора помех в этом случае запишем в виде:

, .                          (4)

Здесь  - весовые коэффициенты, которые должны быть выбраны таким образом, чтобы обеспечить максимум отношения сигнал/(шум +помеха). Трудность заключается в том, что мы не имеем возможности наблюдать отдельно сигналы помех каждого из источников.

Представим спектральные отсчеты на выходах приемных устройств в виде:

- для первого компенсационного канала ():

, ;                            (5)

- для второго компенсационного канала ():

, .                           (6)

- для М-го компенсационного канала ():

, .                                   (7)

 - для суммарного канала ():

, ;                           (8)

Здесь , ,  - комплексные амплитуды сигнала на выходах соответствующих приемных каналов для -го элемента разрешения;  - комплексные значения коэффициентов усиления антенн М приемных каналов (первый индекс), обусловленные приходом сигналов P источников помех (второй индекс);  - нормированные (при единичном коэффициенте усиления антенны) комплексные амплитуды сигналов для -го элемента разрешения, обусловленные воздействием -го источника помеховых сигналов; , , …,  - отсчеты, описывающие отклик приемника в соответствующем элементе разрешения на собственные шумы - независимые случайные комплексные величины с нулевым математическим ожиданием и дисперсией  (), где  - полоса пропускания доплеровского фильтра,  - уровень спектральной плотности собственного шума соответствующего канала. В дальнейшем для упрощения будем полагать, что .

Отсчеты , ,…,  () являются наблюдаемыми, а приведенные выше выражения для них представляют собой систему уравнений. Если пренебречь собственными шумами, то из (5) - (8) можно получить:

.

Откуда для -того канала (в предположении, что в этом канале отсутствует полезный сигнал) будем иметь:

.

Следовательно:

.                            (9)

Таким образом, отклик на помеховое воздействие в i-ом элементе разрешения основного канала можно представить как взвешенную сумму откликов на помеховый сигнал в компенсационных каналах и, следовательно, зная весовые коэффициенты можно, в соответствии с (4) осуществить операцию компенсации. В уравнении (9) величины , , …,  известны, однако система не решается. Для определения истинного значения этих коэффициентов дополним это уравнение аналогичными, но сформированными для других () элементов разрешения для получения  матрицы. Получим  уравнений с  неизвестными:

, .

Решение этих уравнений позволяет определить требуемые для компенсации вектора весовых коэффициентов:

     (10)

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

, , .

Характеристики M-канального компенсатора, описанного выше, исследовались методом математического моделирования. В частности, на рис. 7 приведен отклик приемника основного канала трехканального компенсатора при воздействии двух источников широкополосных помех (отношение сигнал/(шум+суммарная помеха) составляет –5 дБ); на рис. 8 - отклик приемника основного канала после компенсации, выполненной в соответствии с выражением (4) (полезный сигнал в 2720-ом доплеровском фильтре). Коэффициент подавления помехи в данном случае составил 24.5 дБ. На рисунках по оси абсцисс отложены номера доплеровских фильтров, по оси ординат – относительная амплитуда сигналов в каждом элементе разрешения.

Рис. 7. Отклик приемника суммарного канала при воздействии широкополосных помех (2 источника)

Рис. 8. Отклик приемника суммарного канала при воздействии широкополосной помехи (2 источника) после компенсации

На рис. 9 представлен график зависимости коэффициента подавления  суммарной помехи нескольких источников многоканальными компенсаторами с числом каналов 2 ÷ 6 при разных отношениях мощности результирующего помехового сигнала к мощности собственного шума. При моделировании учитывалось увеличение коэффициента усиления антенны по суммарному каналу в связи с возрастанием числа каналов.

Рис. 9. Зависимость коэффициента подавления мешающих сигналов нескольких источников для многоканальных компенсаторов

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

Таким образом,  предложенный алгоритм:

- позволяет по данным, полученным на интервале когерентного накопления, обеспечить компенсацию помеховых сигналов нескольких источников;

- дает возможность обеспечить работу компенсатора при наличии отражений от подстилающей поверхности, принимаемых по основному лепестку диаграммы направленности антенны (поскольку последние могут быть отселектированы по частоте);

- может быть распространен на системы с достаточно большим числом приемных каналов (БРЛС с ФАР).

Литература:

1.Перунов Ю.М., Фомичев К.И., Юдин Л.М. Радиоэлектронное подавление информационных каналов систем управления оружием. – М.: Радиотехника, 2003. 415 с.

2.Цветнов В.В., Демин В.П., Куприянов А.И. Радиоэлектронная борьба: радиоразведка и радиопротиводействие. – М.: Радиотехника, 2003.248 с.

3.Монзиго Р.А, Миллер Т.У. Адаптивные антенные решетки.  М.: Радио и связь, 1986. 448 с.

4.Ратынский М.В. Адаптация и сверхразрешение в антенных решетках. 2003г.

5.Обработка сигналов в многоканальных РЛС/ А.П.Лукошкин [и др.]/. М.,: Радио и связь,1983. 328 с.

6.Когерентный автокомпенсатор помех для многоканальных импульсно-доплеровских бортовых радиолокационных станций: Патент № 2374661/ Буров А.С., Голосов П.В., Родзивилов В.А., Загородний В.Г., Колотов С.А. Заявлено 27.11.09.