"ЖУРНАЛ РАДИОЭЛЕКТРОНИКИ" N 5, 2006 |
УДК 621.396.96
Оценка
вероятности обнаружения сигнала в
многоканальной радиотехнической системе
С. В. Субботин,
Д. Ю. Большаков
ОАО "Научно-исследовательский
электромеханический институт" (ОАО "НИЭМИ")
Получена 20 ноября 2006 г.
Для вычисления вероятности обнаружения сигнала в n каналах радиотехнической системы из m предложена аппроксимация многомерной плотности распределения огибающих аддитивной смеси сигнала и помехи. Эта аппроксимация путем некоторого искажения корреляционных связей между огибающими в смежных каналах системы дает возможность достаточно просто вычислять интеграл от данной плотности распределения. Рассмотрен пример вычисления вероятности обнаружения сигнала в четырехканальной системе размера 2x2.
Рассмотрим приемную радиотехническую систему, предназначенную для измерения временной задержки и доплеровского сдвига частоты сигнала, отраженного от цели. Такая система является многоканальной и имеет в каждом канале дальности набор цифровых фильтров с фиксированными значениями доплеровских частот. Корреляционно-фильтровой матрицей (КФМ) будем называть совокупность приемных каналов для различных дальностей и доплеровских частот. Такая матрица получается, например, путем оптимальной обработки принятого сигнала при использовании алгоритма быстрого преобразования Фурье.
Соседние численные значения элементов КФМ могут отличатся на достаточно малую величину, что приводит к появлению отраженного от цели сигнала в двух и более элементах матрицы. Представляет практический интерес вычисление характеристик обнаружения сигнала в n элементах КФМ из m (m³2) с учетом корреляционных связей между огибающими выходных значений смежных элементов данной матрицы. Знание указанных характеристик позволит определить вероятность того, что сигнал от одной цели занимает n элементов при заданном отношении сигнал/помеха и фиксированном уровне ложных тревог, который обеспечивает идеально функционирующая шумовая АРУ.
Считаем, что на вход приемной системы поступает аддитивная смесь сигнала и флуктуационной помехи. Сигнал является когерентной пачкой дружно флуктуирующих радиоимпульсов с релевским законом распределения огибающих и равномерным распределением начальных фаз.
Мгновенные значения сигнала на выходе линейной части приемной системы имеют нормальное распределение с параметрами (0, ). Флуктуационная помеха в полосе пропускания приемной системы является стационарным нормальным случайным процессом с параметрами (0, ).
Вероятность обнаружения сигнала в n элементах КФМ из m равна [1]
где - совместная плотность распределения огибающих на выходах элементов КФМ,- пороговое значение, которое связано с - вероятностью ложной тревоги в одном элементе КФМ выражением
В дальнейшем отношение сигнал / шум считаем равным
Из выражения (1) следует, что для вычисления вероятности обнаружения сигнала в нескольких элементах КФМ необходимо определить функцию совместной плотности распределения . Для ее получения может быть использована методика, изложенная в [2]. При этом выражение для m – мерной функции плотности распределения огибающих представляют в виде произведения показательных и бесселевых функций. Интегрирование данного произведения при разложении бесселевых функций в ряд по коэффициентам корреляции не позволяет получить соотношения, удобного для вычисления вероятности обнаружения сигнала в n элементах матрицы КФМ из m.
Введем эмпирическую аппроксимацию, которая путем незначительного искажения корреляционных связей между огибающими аддитивной смеси сигнала и помехи в смежных элементах КФМ , дает возможность достаточно просто вычислять интеграл в выражении (1). Представим вероятность обнаружения сигнала в n элементах КФМ из m в следующем виде
В формуле (4) интегрирование по каждой переменной ведется от нуля до бесконечности, а значения порога учитывается в показательных функциях с помощью коэффициентов , которые введены в подынтегральное выражения. Выбор показательных функций основан на предположении, что с ростом значений огибающей и вероятности превышения порога умножение подынтегрального выражения в соотношении (1) на функцию вида
не изменяет характера этой зависимости. Аналогичные рассуждения можно провести для случая не превышения порога в - м элементе КФМ, где соответствующая вероятность с ростом значения огибающей уменьшается и, следовательно, выбирается функция вида
Значение коэффициента в выражениях (5) и (6) может быть найдено из условия совпадения соотношений (1) и (4) для одномерного случая m=1. Согласно предложенной аппроксимации имеем
Для введенных выше статистических характеристик сигнала и помехи распределение огибающей аддитивной смеси сигнала и помехи совпадает с релеевской функцией [1]
при ,Подставим формулу (9) в выражения (7) и (8)
Учитывая, что вероятности и образуют полную группу событий, т.е.
Из соотношений (10)-(12) находим
В формуле (4) заменим плотность распределения огибающих аддитивной смеси сигнала и помехи на совместную плотность распределения квадратурных составляющих этой смеси . Представим совместную плотность распределения в виде 2m-мерной нормальной плотности распределения
где - дисперсия квадратурной составляющей смеси сигнала и помехи,
- детерминант корреляционной матрицы , соответствующей данной плотности распределения,
- алгебраическое дополнение элемента в матрице .
Подставим соотношение (14) в выражение (4)
Введем обозначение
Используя формулу произведения двучленов [3], получаем
где ,
,
,
и в общем виде есть сумма всевозможных произведений чисел , , …. Предположим, что величина , соотношение (17) можно представить в виде
при ,
где число сочетаний из n по i.
Подставим выражение (18) в формулу (15) с учетом соотношения (16), меняя местами в формуле (15) операции интегрирования и суммирования, а также учитывая независимость квадратурных составляющих огибающих аддитивной смеси сигнала и помехи и , в совпадающие моменты времени, получаем
при ,
где .
Из анализа правой части соотношения (19) следует, что выражение вида
соответствует вероятности не превышения порога в (m+i-n) элементах матрицы КФМ. Тогда вероятность обнаружения сигнала в n элементах матрицы КФМ из m может быть определена в виде
при .
Из выражения (20) следует, что для определения вероятности обнаружения сигнала в n элементах КФМ из m необходимо вычислить вероятности непревышения порога для любого m. При n=0 из соотношения (19) получаем
где
В векторной форме выражение (21) имеет вид
где - вектор столбец размерности [1´m],
- положительно определенная матрица размера [m´m].
Интеграл в правой части соотношения (23) равен [3]
Вероятность непревышения порога в m элементах матрицы КФМ с учетом соотношения (24) равна
Подставляя выражение (25) в формулу (20), получаем итоговое соотношение для определения вероятности обнаружения сигнала в n элементах КФМ из m
при ,
где - детерминант корреляционной матрицы , а элементы матрицы определяются по формуле (22).
В качестве примера рассмотрим КФМ, состоящую из четырех элементов – m=4 (см. рис. 1). Для каждого столбца или строки этой матрицы (при m=2) корреляционная матрица, соответствующая совместному распределению квадратурных составляющих огибающих аддитивной смеси сигнала и помехи в совпадающие моменты времени, имеет вид
где - коэффициент корреляции между значениями огибающих аддитивной смеси сигнала и помехи в строке или столбце матрицы.
Рис. 1
Детерминант корреляционной матрицы равен
Вычислим матрицу по формуле (22). Детерминант этой матрицы равен
Подставляя соотношения (28) и (29) в формулу (26), получаем
В выражении (30) при имеем точное равенство
где вероятность непревышения порога в одном элементе КФМ вычисляется по формуле (10).
Аналогично вероятности могут быть определены вероятности непревышения порога для трех и четырех элементов КФМ, соответственно, и . Полученные значения вероятностей непревышения порога - , , необходимо подставить в выражение (20). Например, вероятность превышения порога в двух или трех элементах КФМ их четырех равна
,
.
(32)
Пусть параметры сигнала цели соответствуют положению точки на рис. 1, т.е. отстоят от центра элемента 1 по времени на и частоте на , где и - размеры элемента КФМ. Если опрос элементов КФМ проводится в один и тот же момент времени, то коэффициенты корреляции в этих элементах равны
, , .
На рис. 2 представлены зависимости (сплошные кривые) обнаружения сигнала в трех из четырех и необнаружения сигнала в четырех элементах КФМ от отношения сигнал / помеха (см. формулу (3)) при фиксированной вероятности ложной тревоги в одном элементе. Выбранные примеры характеризуются большим диапазоном изменения вероятностей и удобны для сравнения результатов расчета и цифрового моделирования.
Цифровое моделирование многоканальной системы проведено с помощью метода комплексной огибающей в матричной форме [4]. В основу данного подхода положены метод комплексных огибающих и уравнения метода пространства состояний, которые записываются в матричной форме и позволяют получить основные зависимости в компактном виде, удобном для программирования на ЦВМ. При этом, задав в качестве компонент вектора состояний значения огибающих выходных напряжений элементов КФМ рассматриваемой системы, можно одновременно проанализировать всю совокупность данных элементов.
На рис. 2 представлены результаты цифрового моделирования (пунктирные кривые) при условиях, описанных выше.
Рис. 2.
Как видно из рис. 2 погрешность между результатами расчета и цифрового моделирования не превосходит 10-15%, что свидетельствует об удовлетворительной аппроксимации многомерной плотности распределения огибающих в элементах КФМ на основе соотношения (26).
Список Литературы:
Тихонов В.И. Статистическая радиотехника. – 2-изд., перераб. и доп. – М.: Радио и связь, 1982.
Пугачев В.П. Теория вероятностей и математическая статистика: Учеб. пособие для вузов. – М.: Мир, 1995.
Бронштейн И.Н., Семендяев К.А. Справочник по математике. – 10-изд. стереотипное. – М.: Наука, 1964.
Субботин С.В., Большаков Д.Ю. Цифровое моделирование многоканальных радиотехнических устройств с помощью метода комплексной огибающей в матричной форме. ОАО "Концерн ПВО "Алмаз-Антей" - Вычислительные устройства и ПО РЛС (тематический сборник) 2003, с. 105-109.
Авторы:
Субботин Сергей Валентинович, Большаков Денис Юрьевич zuml@mail.ru
ОАО "Научно-исследовательский электромеханический институт" (ОАО "НИЭМИ")