УДК 53.GS7/.GSS
Г. А. Семенов
Выявление порога чувствительности спектрального метода к поиску периодического сигнала в нормальном шуме
G.A. Semenov
The Threshold Response Search of Spectral Analysis for Periodical Signal Detection in Gaussian Noise
Приведены результаты теоретического и экспериментального выявления порога чувствительности спектрального метода анализа к поиску периодического сигнала в нормальном шуме. Описан сравнительный анализ полученных зависимостей.
Ключевые слова.: выделение сигнала на фоне шума, спектральный анализ.
In this paper the results of theoretical and practical threshold response search of spectral analysis for periodical signal detection in Gaussian noise are presented. Comparative analysis of received dependences is described.
Key words: signal detection in noise, spectral analysis.
Задачи обнаружения периодических сигналов в шуме решаются при приеме и обработке сигналов различной природы во всевозможных технических устройствах. Основным инструментом поиска периодических компонент является спектральный анализ.
Настоящая работа посвящена выявлению границ применимости спектрального анализа к детектированию периодических сигналов в нормальном шуме.
Рассмотрим сигнал вида:
ГХр () = х05Ш (о1), 0 < I < аТ
х(^=1 М Т <, < Т (1)
I хш (), аТ < t < Т,
где а << 1 - относительная длительность периодического фрагмента; хШ ^) - случайный сигнал
(шум), имеющий несмещенную статистическую плотность распределения р(х) (нормированную на единицу) со средним квадратическим отклонением
аш .
Оценим значение а (т.е. диапазон значений
а >а р), при котором еще возможно обнаружить
регулярную компоненту в сигнале типа (1), с помощью анализа спектра мощности исходного сигнала. При этом будем считать, что для анализа доступен только дискретный набор значений одной реализации сигнала, т.е. временной ряд х„.
Фурье-образ сигнала (1) можно представить в виде:
А (о)=-^ Г х ^)в-шЖ = л/2п
, аТ , т
= ^= Г хР и)в-шЛ +^= Г хш О)в~шЖ =
' 42ж а
= АР (а) + Аш (о), (2)
где АР (о) и Аш (о) - Фурье-образы регулярной и случайной составляющих соответственно. В дискретном случае (для временного ряда хп)
о t 2пк
1 N
A (ак )=-pr'L xne
\ 2п n=1
T
(З)
Односторонняя спектральная плотность мощности сигнала:
S(а) J2 И(а)|2, а °
(4)
[0, а < 0.
В соответствии с теоремой Парсеваля для регулярного сигнала:
| |хР ()| Ж =| |АР (а) йа = 2/1 Ар (®)| йа,
—да —да 0
(х2р)аГ «(|Ар|2) А®«(|Ар|2)ОТ, (5)
Здесь (хр} и ^|АР- средние квадраты (интенсивности) регулярной составляющей сигнала и ее Фурье-образа соответственно.
В оценке (5) использовалось соотношение неопределенностей Аа аТ > 2п .
Применяя теорему Парсеваля, аналогично для хЩґ) в случае дискретного преобразования Фурье:
(х2ш)(1 — а)Т *(|Аш\2)АО.ш \АШ|2)П, (6)
где величина частотной полосы шума
До 2nN
ДОш =■
T
(7)
Из соотношений (5), (б) получаем уравнение для оценки величины a :
(xp) = (к 01
1 — a, {Ш (Иш I’)
N
(S)
Обозначим относительную интенсивность регулярной составляющей сигнала:
(xp)
V(x-)
= П :
(9)
относительную оценку спектральной плотности регулярной компоненты -
(кО = (sp)
= p2.
(iG)
(lAШ Г) ^ш)
Разрешая (Б) относительно a, > 0 в обозначениях (9)-(iG) получаем:
p
a, =■
1+
в1
в1
JN\ 4г)2 N 2rf N
(ii)
мальное значение
XSD
Max о •2 _ Sp
а
в2 =
(Ш )ln f N )
а
(SШ )
= lnl
(i3)
в
1 +
в 2
в в
n4N\ 4 Nr/2 2 Nn2
i+
lnl N) ь,fN
(i4)
4Nn2 2Nn2
В силу того, что для N >> 1 верно
К1
N
<< 1,
(15)
граничное значение (14), пренебрегая слагаемыми высоких порядков малости, можно переписать в виде:
a, =
П
,l N
N
(1б)
Если вид регулярной компоненты сигнала известен, величину в можно заменить отношением спектральной плотности мощности регулярной составляющей сигнала к спектральной плотности мощности шума на частоте о0, соответствующей главной гармонике регулярной компоненты, так как на практике в алгоритме поиска периодичности используется не усредненное значение (5”^ , а макси-
Sp массива спектральной
плотности исходного сигнала, на которое накладывается условие (12), легшее в основу алгоритма детектирования периодичности дискретного сигнала xn (Шайдук А.М., Останин С.А., Семенов Г.А. Об улучшении чувствительности спектрального анализа // Лазеры. Измерения. Информация. СПб., 2G1G. С. 45):
Sg >а^ш)lnfyj . (12)
Для получения результата в общем виде перепишем выражение для в2, входящего в (11), используя условие (12), с учетом ^ш } = а(Sш) :
Таким образом, мы получили общее выражение для граничного значения относительной длительности периодического фрагмента, при котором еще возможно его обнаружение классическим спектральным методом анализа. Видно, что величина агр зависит только от объема исследуемой выборки
и относительной оценки спектральной плотности мощности регулярной компоненты, т. е. отношения спектральной плотности мощности регулярной компоненты к спектральной плотности мощности шума на частоте, соответствующей главной гармонике регулярного сигнала.
Границы применимости классического спектрального метода анализа для детектирования периодического фрагмента также были установлены посредством численного моделирования с применением алгоритма, основанного на использовании соотношения (12).
Численный эксперимент проводился для временных рядов, соответствующих сигналам вида (1), для выборок объемами N = 104 и N = 105 элементов. Изменяемыми параметрами временных рядов являлись относительная длительность периодического фрагмента а, относительная интенсивность периодического фрагмента ц. Контрольный параметр - двоичный результат проверки условия (12) для массива спектральной плотности мощности исследуемых временных рядов (при этом дополнительно проверялось соответствие найденного значения частоты периодического фрагмента, для которой условие было выполнено, - с истинным моделируемым значением частоты периодического фрагмента). Для каждого конкретного набора параметров ^, а, ц) проверка условия (12) проводилась 104 раз. Количество положительных исходов, отнесенное к общему количеству испытаний, определяло вероятности успешного детектирования периодического фрагмента при заданном наборе контрольных параметров:
2
1
При проведении эксперимента результат детектирования считался успешным, если количество отрицательных исходов не превышало 10%, т.е. значение эмпирически определенной вероятности
Р(N,а,п) > 0,9, (18)
т > 9-103. (19)
На рисунке 1 приведена зависимость Р(а) вероятности обнаружения периодического фрагмента от его относительной длительности для разных значений относительной интенсивности периодического фрагмента п для выборок объемом N = 104 элементов.
Вероятность обнаружения Р(а)
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
Относительная длительность периодического фрагмента а , %
Рис. 1. Зависимость вероятности обнаружения периодического фрагмента от его относительной длительности при разных относительных интенсивностях
Р (Ы,а,п) =----------------------, (17)
v ’ 10000
где т - количество положительных исходов испытания (проверки условия (12) для всех значений массива $(юі)).
Построение зависимости Р(Ы, а) для фиксированных значений п сделало возможным определение связи граничных значений данных параметров а р и пгр , для которых еще возможно проведение
успешного детектирования, и сравнить полученные результаты с аналогичными теоретическими результатами (соотношение (16)).
Линии, соединяющие экспериментальные точки на рисунке 1 (и последующем рисунке 3), проведены для наглядности, смысловой нагрузки не несут; подписи рядом с каждой кривой на рисунке соответствуют значениям параметров, для которых данная кривая была построена, также приведены для большей наглядности.
Анализ рисунка 1 позволил установить, что успешное детектирование (с вероятностью > 0,9) периодического фрагмента с единичной относительной интенсивностью (п = 1) возможно при относительной длительности периодического фрагмента а > 3%, а, например, для периодического фрагмента с половинной относительной интенсивностью (п = 0,5) - при а > 7%.
На рисунке 2 представлены зависимости граничного значения относительной длительности периодического фрагмента, достаточной для его успешного детектирования (в случае выполнения условий (18)-(19)) от относительной интенсивности периодического фрагмента, полученные экспери-
ментально и теоретически в диапазоне 0,3 < п < 1,0 для выборок объемом N = 104 элементов.
Из рисунка 2 видно, что теоретические граничные значения параметров периодического фрагмента, при которых еще возможно его успешное детектирование в шуме, несколько ниже значений, полученных экспериментально (все отличия лежат в пределах 20%), однако сами зависимости при этом имеют аналогичный вид.
Таким образом, с учетом того, что полученные значения параметров, необходимые для детектирования, во-первых, являются граничными, а во-вторых, задаются исходя из требований точности (т. е. в некоторой степени произвольно), можно сделать вывод, что эксперимент адекватно подтвердил результаты, полученные теоретически для данных наборов параметров.
На рисунке 3 приведены результаты численного эксперимента по выявлению границ применимости классического спектрального метода анализа, проведенного для выборок объемом N = 105 элементов.
Относительная длительность периодического фрагмента арр , % 12
— —1 і
\ ,
\
\
Теоретическая зависимость ■ Экспериментальная зависимость
= — (Л = 101)
34 7
0,2 0 Отн 3 0 осител 4 0 ьная т 5 0 ГГЄНСИЕ 6 0 шость 7 0 период 8 0 ИЧЄСКО! 9 'о фраг 1.1 мента г
Рис. 2. Зависимость граничного значения относительной длительности периодического фрагмента от его относительной интенсивности при успешном детектировании
В ер оятно сть обнаружения р(а) 1.0
0.1
0,0
7 = 1, И ь-—1 1——I
п = о, / / ^=0, )
і V
N К >5 II Р 4
7 = о! 7
т/ = 0 ,6/
*=0,
1 к-—
0,9 1,0 1,1 1,2 1,3 1,4 1,5 1,6 1,7 1,8 1.9 2,0 2,1 2,2 2.3
Относительная длительность периодического фрагмента а , %
Рис. 3. Зависимость вероятности обнаружения периодического фрагмента от его относительной длительности при разных относительных интенсивностях
Диапазон изменения относительной интенсивности периодического фрагмента сигнала при этом оставался прежним - 0,3 < п < 1,0; диапазон изменения значений относительной длительности периодического фрагмента, в котором происходили характерные изменения зависимости Р(М, а), был выявлен эмпирически в ходе эксперимента: 1,0% < а < 2,2%.
Из рисунка 3 видно, что успешное детектирование (с вероятностью >0,9) периодического фрагмен-
та с единичной относительной интенсивностью (п = 1) возможно при относительной длительности периодического фрагмента а > 1%, для фрагмента с половинной относительной интенсивностью (П = 0,5) - при а > 2,2%.
Рисунок 4 иллюстрирует зависимости, аналогичные приведенным на рисунке 2, для выборок объемом 105 элементов в диапазоне относительной интенсивности периодического фрагмента сигнала 0,5 < п < 1,0.
Относительная длительность периодического фрагмента , %
'•г
2,1 -
1,7
1 1
Теоретическая зависимость ■ Экспериментальная зависимость
™гр = — (Л = 105)
■ і
0,4 0,5 0,6 0,7 0,8 0,9 1 1,1
Относительная интенсивность периодического фрагмента V
Рис. 4. Зависимость граничного значения относительной длительности периодического фрагмента от его относительной интенсивности при успешном детектировании
Анализ рисунка 4 позволяет сделать вывод, что, как и в случае для выборок объемом N = 104 элементов, характер теоретической и экспериментальной зависимостей полностью аналогичен, причем в данном случае отличие экспериментальных и теоретических значений не превосходит 5,5%.
Таким образом, в ходе проведенного эксперимента было установлено, что в случае выборок объемом N = 105 элементов экспериментальные
результаты также адекватно согласуются с теоретическими результатами, причем с лучшей точностью.
На основе соотнесения теоретических и экспериментальных результатов можно сделать вывод, что полученное выражение (16) позволяет адекватно находить граничные значения чувствительности спектрального метода к детектированию периодических сигналов в нормальном шуме.