CC BY
139
В.И. Костылев,
доктор физико-математических наук,
И.П. Гресь,
Воронежский государственный
профессор, Воронежский государственный университет университет
ОБОБЩЕННОЕ ЭНЕРГЕТИЧЕСКОЕ ОБНАРУЖЕНИЕ СЛУЧАЙНОГО ГАУССОВСКОГО СИГНАЛА В ГАУССОВСКОМ
КОРРЕЛИРОВАННОМ ШУМЕ
GENERALIZED ENERGY DETECTION OF RANDOM GAUSSIAN SIGNALS IN CORRELATED GAUSSIAN NOISE
Показано, что применение альтернативного обнаружителя случайного сигнала Чена на фоне коррелированного гауссовского шума позволяет заметно повысить вероятность обнаружения. Тогда как его использование в случае некоррелированного шума, напротив, нецелесообразно.
It was shown, that Chen detector for random signals in correlated Gaussian noise provides better performance. But usage the detector in case of uncorrelated noise is beside the purpose.
Задача оптимального энергетического обнаружения случайного сигнала на фоне шума рассматривалась неоднократно [2, 3]. Было показано, что в случае, если полезный сигнал — гауссовский с некоррелированными отсчетами, и шум — также гауссовский с некоррелированными отсчетами, то оптимальный обнаружитель по критерию Неймана — Пирсона — энергетический [3].
Тем не менее, в статье [1] предлагается альтернативный обнаружитель, который описывается следующим правилом вынесения решения:
где T’ — порог, п — объем выборки, Н1 и Но — гипотезы о наличии и отсутствии полезного сигнала соответственно. Такой обнаружитель отличается от энергетического тем, что модуль входных отсчетов возводится в произвольную положительную степень р , а не в квадрат.
Статистические характеристики случайной величины Ж’ определить сложно, если вообще возможно. В [1] предлагается аппроксимация данной решающей статистики с помощью гамма-распределения, путем подбора параметров последнего. Автор указывает на то, что такая аппроксимация хорошо описывает поведение Ж’ и может быть использована для нахождения порога Т’.
H
(І)
Работа выполнена при поддержке Российского фонда фундаментальных исследований, грант
№ 13-01-00773.
В статье [1] представлены графики (см. рис. 1) зависимостей оптимального значения показателя степенир, соответствующего максимальной вероятности правильного обнаружения Рв, от отношения сигнал /шум (ОСШ) при фиксированных значениях вероятности ложной тревоги Рр и объема выборки п, на этом же рисунке ниже представлены графики зависимостей оптимального значения показателя степени р, соответствующего минимальной вероятности ложной тревоги Рр при фиксированных значениях вероятности правильного обнаружения Рв и объема выборки п.
Рис. 1. Зависимость оптимального значения показателя степенир, соответствующего максимальной вероятности правильного обнаружения, от ОСШ при фиксированных значениях вероятности ложной тревоги и объема выборки (вверху); зависимость оптимального значения показателя степени р, соответствующего минимальной вероятности ложной тревоги от ОСШ при фиксированных значениях вероятности правильного обнаружения и объема выборки (внизу) [1].
Из анализа рисунка можно сделать вывод о том, что различным значениям ОСШ будут соответствовать различные значения показателя степени р. Казалось бы, получен новый обнаружитель, обладающий лучшими характеристиками. Однако несложно показать теоретически, что именно энергетический (т.е. обнаружитель, в котором входные отсчеты возводятся в степень р=2) является оптимальным по критерию Неймана-Пирсона для любого ОСШ.
Пусть в результате наблюдения получен вектор х, представляющий собой значение векторной случайной величины X, принимающей значения из пространства наблюдений Я (п-мерная область, которая объединяет все возможные исходы опыта). Опыт, в результате которого выносится решение, можно интерпретировать следующим образом : у нас имеется точка х из пространства наблюдений Я и необходимо выработать решение уо или у1 в соответствии с правилом вынесения решения. Правило решения сводится к заданию поверхности в пространстве наблюдений Я, которая делит это пространство на две части Я0 и Я1, такая поверхность называется решающей поверхностью. Если х попадает в первую часть пространства наблюдений (Я0), необходимо вынести решение уо, если во вторую (Я1) — у1. Так как х — значения отсчетов случайного
процесса, то каким бы правилом решения мы ни пользовались, мы можем допустить ошибку. Возможны две ошибки:
1. Выносится решение у1, а верна гипотеза Н0 — ошибка первого рода (ложная тревога).
2. Выносится решение у0, а верна гипотеза Н1 — ошибка второго рода (пропуск сигнала).
Условная вероятность ошибки первого рода рассчитывается по формуле
Рр = | Р(х | Н0 ^х. (2)
Я1
Условная вероятность ошибки второго рода — по формуле
Р, = | Р(х|Н,)Л. (3)
Я0
Введем целевую функцию
¥ = Р, + (| Р(х | Н0>А - Рр )1, (4)
Я1
где \ — неопределенный множитель. Мы используем критерий Неймана — Пирсона, по которому вероятность ложной тревоги фиксирована, поэтому выражение в скобках равно нулю. Очевидно, что вероятность ложной тревоги можно переписать в следующем виде
Рр = | Р(х | Н0)й?х = 1 - | Р(х | Н0)ёх. (5)
Я1 Я
Подставим (3) и (5) в (4):
¥ = | Р(х |Н1)ёх + (1 - | Р(х | Н0)ёх-Рр)1. (6)
Я0 Я0
Преобразуем выражение (6)
¥ = \ [Р(х | Нх)-1Р(х | Н0)]ёх + 1(1 - Рр ). (7)
Я0
Выражение в квадратных скобках обозначим как У(х). Необходимо провести решающую поверхность в пространстве наблюдений Я таким образом, чтобы Т было минимально. Поэтому если для некоторой точки х е Я У(х)<0, то включим эту точку в Я1, если же для этой точки У(х)>0, то не будем ее включать в Я1. Получим следующее правило:
У
РхН) >1. (8)
Р(х|Н>) < ^
У
Так как вероятность ложной тревоги фиксирована, параметр \ может быть определен из выражения (8).
В [1 ] рассматривается полезный сигнал, представляющий собой гауссовский дискретный п-мерный случайный процесс с нулевым математическим ожиданием и дисперсией а2 и шум, представляющий собой дискретный гауссовский п- мерный случайный процесс с нулевым математическим ожиданием и дисперсией о2.
В этом случае рассмотренные условные плотности вероятности имеют вид:
1 1 п
Р(х 1Н0) = / 2 ехР{-7^12 Гг!!. (9)
Дх!#,) =
ехр{-
■\]2р(&2 + а2) 2(& + а
(10)
2=1
Подставим (9) и (10) в (8) и прологарифмируем полученное выражение її
П
Е
1 >
2 л, <
Ї0
(11)
где
л =
2&2(а2 + &2) г
а
1п 1-—1п 2
&
\Л
(12)
Таким образом, показано, что оптимальным по критерию Неймана — Пирсона является энергетический обнаружитель, т.е. обнаружитель, у которого показатель степени равен двум независимо от ОСШ, значения вероятности ложной тревоги и объема выборки.
Чтобы еще раз удостовериться в правильности приведенных рассуждений, было проведено статистическое моделирование. Моделирование проводилось в два этапа. Сначала была найдена зависимость порога Т от вероятности ложной тревоги Рр и показателя степени р (при моделировании не использовалась аппроксимация из [1]). Эта зависимость представлена на рис. 2. Для сравнения на этом же рисунке представлена та же зависимость, полученная с помощью предложенной в статье аппроксимации (Г’).
30
25
20
- 115
10
1— т т !
і ! Л
/I // /7
// V
у7
1.5
2.5
3.5
4.5
Рис. 2. Зависимость порога Тот показателя степенир (сплошная линия), соответствующая вероятности ложной тревоги Рр =0.01 и объему выборки п=20, та же зависимость, полученная с помощью аппроксимации распределения Ж’ (пунктирная линия)
1
1
П
2=1
Рис. 3. Разность между указанными на рис. 2 кривыми в зависимости
от показателя степени Уже на данном этапе, из анализа представленных графиков, можно сделать вывод о том, что аппроксимация, предложенная в статье, недостаточно точная. Далее, с использованием найденной зависимости Т(р) и аппроксимации Т’(р) из [1], с помощью моделирования были найдены зависимости вероятности правильного обнаружения Рв от показателя степени р при фиксированной вероятности ложной тревоги Рр=0.01 и объеме выборки п = 20. Результат моделирования представлен соответственно на рис. 4 и 5.
Рис. 4. Зависимость вероятности правильного обнаружения Рв от показателя степенир для различных значений ОСШ. Графики получены без использования аппроксимации из [1]
Рис. 5. Зависимость вероятности правильного обнаружения Рв от показателя степенир для различных значений ОСШ. Графики получены с помощью аппроксимации из [1]
Из указанных графиков можно сделать вывод о том, что приведенные теоретические выкладки верны, оптимальным по критерию Неймана — Пирсона является энергетический обнаружитель (р=2). Также можно сделать вывод о том, что представленные в статье [1] графики не соответствуют действительности. Причиной этому является неточная аппроксимация решающей статистики Ж’ гамма-распределением.
Однако мы считаем, что использование обнаружителя, задаваемого правилом (1), назовем его обобщенным энергетическим обнаружителем, имеет смысл. Теоретически доказано, что при обнаружении дискретного гауссова сигнала с некоррелированными отсчетами с нулевым математическим ожиданием на фоне дискретного гауссова шума с некоррелированными отсчетами и нулевым математическим ожиданием оптимальным по критерию Неймана — Пирсона является энергетический обнаружитель. Но если отсчеты шума коррелированы, энергетический обнаружитель уже не является оптимальным. Кроме того, найти вид оптимального обнаружителя достаточно сложно, необходимо знать информацию о корреляции отсчетов, поэтому на практике может использоваться энергетический, в силу его простоты. Используя обобщенный энергетический обнаружитель, можно добиться повышения вероятности правильного обнаружения при неизменном ОСШ. Заметим, что переход от энергетического обнаружителя к (1) не приводит к существенному усложнению обработки входной смеси сигнала и шума.
Для определения характеристик алгоритма (1) при обнаружении в коррелированном шуме было проведено статистическое моделирование, при котором производилось обнаружение дискретного гауссова сигнала с нулевым математическим ожиданием и некоррелированными отсчетами в дискретном гауссовском шуме с нулевым математическим ожиданием и коррелированными отсчетами. Корреляционная функция — экспоненциальная, параметр а=0.25, 0.5, 0.75. Ее график представлен на рис. 6.
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
л А . 'Д. а = 0.25 а = 0.5 а = 0.75
/ К п
/ ^ \
///
/ ■ \ \
/ / ■''V 1 \
/ / .1 \ ^ \
/" / / \ \
/ / \ \ \
' ^ — - - У ~~ - _
-10
-8
-6
-2
8
10
Рис. 6. Экспоненциальная корреляционная функция при различных значениях параметра а
По результатам моделирования построены графики зависимости вероятности правильного обнаружения Рв от показателя степени р при фиксированных вероятности ложной тревоги Рр=0.01, объеме выборки п=20, ОСШ=4дБ и а=0,25, 0,5, 0,75. Эти графики представлены на рис. 7.
0.9г 0.85 -
0.8
0.75
0.7
0.65
0.6
0.55
0.5
0.45
0.4
• -....
/ /■
.■ /
:• ? : I
71 ; /
Т"7 1
Р0(Р), а = 0.25 . ♦ тах[Р0(р)], а = 0.25 Ро(р),а = 0.5 ▼ тах[Р0(р)]р а = 0.5 ■ Р0(р), а = 0.75 * тах[Р0(р)], а = 0.75
1 /
1 / 1
t Т"
0
12
14
16
18
Рис.7. Зависимость вероятности правильного обнаружения Рв от показателя степенир при ОСШ = 4дБ и различных параметрах а корреляционной функции
Рассмотрим рис. 7 (кривую, соответствующую а=0,25). Максимальная вероятность правильного обнаружения достигается при р=9,1 и равна при этом Рв =0,74. Для сравнения прир=2 (энергетический обнаружитель) — Рв =0,58. Таким образом, удалось заметно повысить вероятность правильного обнаружения выбором показателя степени обнаружителя, существенно не усложнив его конструкцию. При этом выигрыш в ОСШ равен примерно 0,8-0,9 дБ, что составляет 20-22% от ОСШ.
Построим характеристики обнаружителя (рис. 8) при различных значениях р=1, 2, 4, 8 при фиксированных а=0,25, вероятности ложной тревоги Р^=0,01 и объеме выборки п=20. Из рис. 8 видно, что при больших ОСШ значение показателя степени р не играет роли, тогда как при малых ОСШ обнаружитель с р=4 и 8 более эффективен по сравнению с энергетическим.
q
Рис. 8.Зависимость вероятности правильного обнаружения PD от ОСШ, при различных
показателях степени обнаружителя.
Таким образом, показано, что при обнаружении гауссовского некоррелированного сигнала на фоне гауссовского некоррелированного шума оптимальным по критерию Неймана — Пирсона является энергетический обнаружитель, а не обнаружитель предложенный в статье [1]. Кроме того, установлено, что обнаружитель (1) обладает лучшими характеристиками, по сравнению с энергетическим, при обнаружении гауссовского некоррелированного сигнала на фоне гауссовского коррелированного шума.
ЛИТЕРАТУРА
1. Chen Y. Improved Energy Detector for Random Signals in Gaussian Noise. — IEEE Transactions on Wireless Communications, Feb., 2010. — V.9. — P. 558—563.
2. Тихонов В. И., Шахтарин Б.И., Сизых В. В. Случайные процессы. Примеры и задачи. Том 4. Оптимальное обнаружение сигналов. — М.: МГТУ им. Н.Э. Баумана, 2005. — 366 с.
3. Левин Б. Р. Теоретические основы статистической радиотехники. — М.: Радио и связь, 1989. — 656 с.
