Обнаружение флуктуирующих малоразмерных целей на фоне флуктуирующей пассивной помехи по наблюдениям поляриметрической РСА
А.Н. Детков, А.В. Никоненко, Д.А. Ницак
Аннотация
В перспективных радиолокационных системах, использующих автоматическое обнаружение целей, селекция флуктуирующих малоразмерных целей (ФМЦ) на фоне пассивных помех представляет собой достаточно сложную задачу. Предложена вероятностная модель радиолокационного канала разностной обработки векторных наблюдений поляриметрической радиолокационной станции с синтезированной апертурой антенны (РСА), который позволит гарантированно обнаруживать флуктуирующие по поляризационным признакам малоразмерные цели на фоне флуктуирующей пассивной помехи. Показано, что из всех возможных ФМЦ можно извлечь две специфические (сингулярные) цели - максимально видимую и наиболее скрытую в данной пассивной помехе, что в целом повышает эффективность обнаружения малоразмерных целей.
Ключевые слова: обнаружение целей, флуктуирующая малоразмерная цель, поляриметрическая РСА.
Abstract
In looking radar systems using automatic target detection, selection fluctuating small-sized targets (FST) on the background clutter is a rather complicated task. We propose a probabilistic model of the radar channel difference vector processing observations polarimetric synthetic aperture radar (PolSAR), which will guarantee to detect fluctuating on polarization featured small-sized target to fluctuating background clutter. It is shown that among all possible FST can extract two specific (singular) purpose - maximum visible and most hidden in the clutter, which generally increases the efficiency of detection of small-sized targets.
Keywords: detection of targets, the fluctuating small-sized target, polarimetric SAR.
Введение
Основой теории обнаружения стабильных и флуктуирующих целей на фоне пассивной помехи является теория обнаружения детерминированных и стохастических сигналов на фоне коррелированных помех [1, 2]. В случае векторных наблюдений при обнаружении радиолокационных целей в соответствии с этой теорией используются аналогичные алгоритмы скалярной радиолокации для гауссовской модели отражений [3]. В целом алгоритмы обработки векторных наблюдений являются обобщением более общих математических объектов, для которых, тем не менее, статистическое и физическое содержание итоговых результатов сохраняется. В [4] решена задача оптимальной поляризации антенны идеального радиолокатора с использованием матриц (векторов) рассеяния цели и априорных данных, а также определены граничные функции эффективности обнаружения поляриметрической РЛС некоторых классов флуктуирующих целей на фоне пассивной помехи. Эта же задача решена для обнаружения поляриметрической РСА квазистабильных малоразмерных целей на фоне подстилающей поверхности [5], причём разработке алгоритма обнаружения предшествовал анализ поляризационных портретов искусственных объектов и фона, полученных в результате математического моделирования [6]. Анализ поляризационных характеристик распределенных объектов показывает, что в параметрах поляризации рассеянной объектом волны содержатся вспомогательные информационные признаки, которые необходимо использовать в задачах пространственной (поляризационной) селекции.
Однако во всех этих случаях качество обнаружения флуктуирующей цели зависит от многих интегральных и еще большего числа физических параметров, что существенно затрудняет анализ влияния этих параметров (цели и помехи) на указанное качество. Затруднение определяется и тем, что связь параметров с вероятностями правильного и ложного обнаружения сложна и не однозначна. Так что максимизация какого-то из параметров не ведет или не обязательно ведет к улучшению вероятностных показателей [4].
Цель работы: разработка алгоритма функционирования радиолокационного канала в поляриметрической РСА при многолучевой обработке векторных наблюдений, который бы гарантировано повышал эффективность обнаружения флуктуирующих целей на фоне флуктуирующего фона с учётом априорного анализа поляризационных характеристик целей.
Обобщенная задача обнаружения флуктуирующей цели на фоне пассивной помехи
Данные поляриметрической РСА традиционно принято описывать комплексной матрицей рассеяния цели S [7]
,
согласованные (основные) компоненты которой и в общем коррелированны, в то время как корреляции между основными и кроссовыми компонентами и или и () практически отсутствуют (теоретически строго нулевые) [8, 9]. Анализ изображений РСА отмечает наличие слабой корреляции между основными и кроссовыми компонентами. Причиной такой взаимной корреляции, как правило, являются: не идеальность развязки по кроссовой поляризации между ортогональными каналами, тепловые шумы и неточности математического моделирования. Для взаимного случая вектор-наблюдение принимает вид
,
где - знак транспонирования. Для невзаимного случая, например в бистатических системах, где размерность увеличивается до 4. Если объектом наблюдения РСА является шероховатая поверхность с достаточно большим числом элементарных отражателей, модель вектора имеет многомерное комплексное гауссово распределение [10]
, (1)
где - размерность вектора ; C - комплексная корреляционная матрица , - знак эрмитовой сопряженности; - операция усреднения по ансамблю реализаций. Вещественные и мнимые части координат m - мерного вектора - , , , имеют нулевое математическое ожидание - и удовлетворяют условиям аналитического сигнала:
; ; ; . (2)
Данные РСА в целом удовлетворяют (1), однако основным фактором, затрудняющим дешифрирование радиолокационных изображений, в том числе обнаружение целей, являются спекл-шум, обусловленный интерференционными явлениями при суммировании сигналов от отдельных центров рассеяния («блестящих точек») в пределах каждого элемента разрешения. Наиболее эффективным способом снижения уровня спекл-шума в РСА является некогерентное суммирование изображений, полученных от одних и тех же участков местности с некоррелированными спекл-структурами в пределах интервала синтезирования. При цифровой обработке в РСА технически проще реализовать многоканальную (многолучевую) обработку сигналов по азимуту [11]. Статистические характеристики данных, полученных поляриметрическими РСА с учетом многолучевой обработки существенно отличаются от однолучевых. Выборочная корреляционная матрица, полученная усреднением реализаций N независимых субапертур (числа независимых наблюдений)
, (3)
описывается комплексным распределением Уишарта [12, 10].
При обнаружении флуктуирующей цели вид решающей статистики определяется отношением
флуктуирующий малоразмерная цель радиолокационный
, (4)
где K - матрица когерентности вектора , соответствующего отражению от флуктуирующей цели либо помехи (нижний индекс матрицы - 0 соответствует индексу пассивной помехи, на фоне которой решается задача обнаружения; нижний индекс матрицы - 1 соответствует индексу обнаруживаемой цели). Использование в (4) вместо корреляционных матриц и матриц когерентности объясняется выполнением условия (2). В этом случае решающее устройство приемника должно вычислять и сравнивать с порогом квадратичный функционал
. (5)
Весовая матрица W, определяется формулой [6]
. (6)
Основная статистика алгоритма обработки наблюдений (5) для любого из случаев при имеет существенно более сложное распределение в сравнении с основной статистикой задачи обнаружения стабильных целей [5], причём число обобщенных параметров, определяющих вероятность правильного обнаружения , вероятность ложного обнаружения или значение порогового сигнала, здесь практически всегда более двух.
В [4] установлено, что основная вероятностная характеристика обнаружения флуктуирующих целей на фоне пассивных помех имеет общий вид , где , - некоторые векторы, задающие две системы обобщенных (интегральных) параметров, сложно связанных с исходными (физическими) параметрами отражений от обнаруживаемой и мешающей целей. Эта связь сложна и неоднозначна. Максимизация какого-то из параметров не всегда ведет к улучшению вероятностных показателей.
Вид распределения статистики (5) определяется характеризационными свойствами W и [4]:
1) матрица W является эрмитовой и строго положительно определенной;
2) координаты вектора удовлетворяют (2).
При этих условиях характеристическая функция статистики (5) имеет вид
,
где M - корреляционная матрица вектора ( при наличии в сигнальной компоненты и - при её отсутствии); I - единичная матрица размерности .
Произведение U = MW при эрмитовых М, W не обязательно является эрмитовым, но для положительно определенных М и W это произведение имеет не только вещественные, но и положительные собственные числа (СЧ) [13]. Двум значениям матрицы М соответствуют два значения матрицы U: U0 и U1. Пусть обозначает совокупность {л1, л2,..., лm} собственных чисел одного из вариантов матрицы U. Плотность распределения статистики (5) можно записать в виде [4]
. (7)
Здесь подынтегральная функция задана на комплексной плоскости - области задания переменной - и имеет полюсы строго в нижней полуплоскости. Интеграл (7) вычисляется с помощью вычетов в сингулярных точках. Число полюсов определяется числом L различных собственных чисел. Каждому k-кратному СЧ соответствует вычет вида , где - z-полином степени k -1. В результате f(z) оказывается экспоненциально-степенной функцией вида
. (8)
Для конкретизации (7) необходимо задать распределение СЧ л1, л2,..., лm. Двум значениям матрицы U: U0 и U1 соответствуют две совокупности СЧ и , которые определяют две плотности вероятности распределения статистики (4) и . Первая плотность соответствует фону, вторая - фону и цели. Вероятности и определяются формулами:
; , (9)
где порог подбирается по заданной величине .
Наиболее простыми с точки зрения расчетов (9) с учётом (7) является случай некратных СЧ л1, л2,..., лm (наиболее частый случай, применительно к сложным объектам рассеяния). Тогда все полиномы (8) , являются константами:
;
и интегралы (9) приводятся к простому виду
.
Для случая невзаимной поляриметрической РСА при разных , интеграл (7) принимает вид
, (10)
; ;
; ;
.
При сильной взаимной корреляции основных компонент вектора рассеяния и также при снижении мощности кроссовых составляющих (составляющей) до уровня внутренних шумов аппаратуры возможна трансформация входной статистики. Следствием этого является вырождение матрицы когерентности (появление СЧ с нулевыми элементами) и снижение размерности вектора рассеяния до 2-х, т.е. m = 2. Снижение мощности кроссовых составляющих наиболее вероятно в случае воздействия по основному лепестку диаграммы направленности поляриметрической РСА помехи с некоррелированными пространственными координатами [5].
Для анализа характеристик обнаружения при снижении размерности входной статистики преобразуем формулу (10) и её компоненты к виду:
, (11)
где
; , ;
, , , .
Вид распределения СЧ и определяется свойствами произведения MW, в котором оба множителя задаются матрицами когерентности и . Такого рода связи в [4] предлагается рассматривать на основе свойств так называемой главной матрицы , которая не является эрмитовой, но при положительно определенных и имеет вещественные и положительные СЧ [13]. Матрица G связывает корреляционные матрицы двух флуктуирующих целей (обнаруживаемой и мешающей) и является математическим объектом, достаточным для задания вероятностных распределений в рассматриваемой задаче обнаружения [4].
Разностный алгоритм обнаружения флуктуирующей цели на фоне пассивной помехи
Все многообразие комбинаций цели и помехи и соответствующие вариации качества обнаружения в конечном итоге определяются многообразием возможных распределений СЧ чисел главной матрицы. Перебирая эти распределения можно найти наилучшую комбинацию или совокупность для которой вероятность будет максимальной и другую совокупность, для которой данная вероятность будет минимальной. Чтобы обеспечить прочие равные условия, необходимо, кроме вероятности ложного обнаружения фиксировать след главной матрицы , который, таким образом, становится определяющим показателем при вычислении вероятности правильного обнаружения .
На рис.1 приведены граничные функции, определяющие вероятность обнаружения одного ансамбля флуктуирующих целей на фоне другого. Пространственные характеристики ансамблей заданы трёхмерной входной статистикой () для . При уменьшении размерности входной статистики () наблюдается более пологий характер кривых обнаружения, а также незначительное сужение границ между кривыми и . Эти кривые соответствуют конкретным наборам СЧ, определяющим максимально видимые и наиболее скрытые цели на фоне пассивной помехи.
Принимая во внимание возможное (значительное) различие пространственно-временных интервалов корреляции для многолучевых РСА, а, следовательно, и возможность факторизации (5) запишем
,
где матрица W определяется из (6) путём замены , .
Рис. 1. Граничные функции обнаружения для РСА при N = 1: 1 - , 2 - , = 0.01 (штрих пунктирная), = 0.001 (сплошная), = 0.0001 (штриховая).
При обнаружении флуктуирующей цели задача может осложняться нефакторизуемостью корреляционной матрицы помехи из-за наличия системного шума. Однако системный шум (например, аппаратный шум РСА) в большинстве случаев можно считать слабым, что вполне оправдывается при наличии сильной пассивной помехи.
Когда размеры элемента разрешения соизмеримы или меньше средних размеров обнаруживаемой распределённой цели, происходит практически полное экранирование целью фона. В этом случае весовая матрица (6) преобразуется к виду
. (12)
Весовая матрица (12), по сути, является основой разностной весовой обработки, позволяющей эффективно использовать оба выигрыша, полученных за счёт поляризационной селекции одного ансамбля целей на фоне другого и наоборот. Квадратичный функционал разностной обработки в виде алгебраических сумм с учётом факторизации (12) принимает вид: