Е.Б. Мундриевская, Омский государственный университет, кафедра математического моделирования
1. Введение
Классической задачей статистической радиотехники является задача обнаружения сигнала на фоне случайных помех. Большинство из известных в настоящее время алгоритмов основано на байесовском подходе. Недостатком этого подхода является сложность получающихся алгоритмов и не всегда обоснованное на практике требование задания априорных распределений. Подобным недостатком избыточности априорной информации обладают и параметрические методы [2, 3, 4].
Особое положение среди алгоритмов обнаружения занимают непараметрические алгоритмы. Они более неприхотливы и универсальны, чем другие алгоритмы. Действительно, параметры обнаруживаемого сигнала могут быть известны неточно. Поэтому применение параметрических методов может быть затруднительно.
Основные виды непараметрических алгоритмов обнаружения рассмотрены в [5], к ним относятся знаковый, ранговый и их модификации.
Далее предложен новый непараметрический метод обнаружения синусоидального сигнала, основанный на изменении корреляционной функции наблюдаемой выборки.
2. Постановка задачи
Предположим, что на выходе некоторой системы наблюдается процесс, про который мы можем сказать следующее:
1. Это широкополосный шум известной верхней частоты
2. Это смесь широкополосного шума с верхней частотой
Требуется определить, какой из этих двух процессов мы наблюдаем (просто шум или шум с сигналом).
3. Предлагаемое решение
Для проcтоты изложения будем полагать, что
Пусть корреляционная функция шума имеет вид:
Корреляционную функцию сигнала запишем как:
Пусть
Величина этого отношения предполагается неизвестной.
Корреляционная функция смеси шума и сигнала:
Известно, что корреляционную функцию некоторого случайного процесса можно представить в виде канонического разложения [6]:
где,
Di - канонические дисперсии.
Очевидно, что
Из условия
Покажем, что для любого
Запишем последнее выражение в развернутом виде:
Из (1) и (4):
Подставим в (7):
Приведем правую и левую часть неравенства к общему знаменателю:
Т. к.
Отсюда:
T. к.
Очевидно, что при
Следовательно, алгоритм обнаружения сигнала в шуме можно строить исходя из вычисления канонических дисперсий наблюдаемой выборки.
4. Алгоритм обнаружения
При моделировании будем пользоваться следующей модификацией алгоритма, предложенного в [1].
Шаг 1. На основе выборки { yk} вычисляем ковариационные коэффициенты R(0)=Eyiy'i , R(1)=Eyi+1y'i , R(2)=Eyi+2y'i, R(3)=Eyi+3y'i .
Шаг 2. Строим матрицы:
Шаг 3. Вычисляем разложение:
Шаг 4. Определяем:
Шаг 5. Вычисляем разложение:
Шаг 6. Вычисляем сумму:
S=e1+e2.
Утверждается, что при любом
Покажем практическую состоятельность этого вывода.
Действительно, зная
Пусть
Из предложенного алгоритма следует: если rsh(k)=0 (k=1,2,3), то
Пусть теперь
yi=yis + yish.
Шум и сигнал независимы друг от друга, поэтому
Поэтому
Rsh+s(k)=Eyi+ky'i=E(yi+ks + yi+ksh)(yis + yish)'=
=Eyi+ks(yis)'+Eyi+ksh(yish)'
Воспользуемся (8):
Rsh+s(k)=Eyi+ks(yis)'=rs(k).
Т.к.
5. Анализ алгоритма обнаружения
Несмотря на то, что мы предполагали