КОЭФФИЦИЕНТ ЭФФЕКТИВНОСТИ ОКОННОГО СПЕКТРАЛЬНО-СТАТИСТИЧЕСКОГО МЕТОДА ДЕТЕКТИРОВАНИЯ СИГНАЛОВ В ШУМЕ

Научная статья
DOI:
https://doi.org/10.18454/IRJ.2016.50.230
Выпуск: № 8 (50), 2016
Опубликована:
2016/08/18
PDF

Останин С.А.1, Семёнов Г.А.2, Шайдук А.М.3

1Кандидат физико-математических наук, Алтайский государственный университет, ООО "АГФЗ", 2Кандидат физико-математических наук, Санкт-Петербургский государственный аграрный университет, 3Доктор физико-математических наук, Алтайский государственный медицинский университет

КОЭФФИЦИЕНТ ЭФФЕКТИВНОСТИ ОКОННОГО СПЕКТРАЛЬНО-СТАТИСТИЧЕСКОГО МЕТОДА ДЕТЕКТИРОВАНИЯ СИГНАЛОВ В ШУМЕ

Аннотация

Описан способ увеличения чувствительности детектирования коротких периодических сигналов в гауссовском шуме, основанный на разделении исследуемой аддитивной смеси сигнала и шума на фрагменты с использованием временного окна. Получена функциональная зависимость значения универсального спектрального критерия детектирования от количества участков разделения исследуемого временного ряда. Введён коэффициент эффективности оконного спектрально-статистического метода, позволяющий определить оптимальное количество участков разделения аддитивной смеси сигнала и шума, необходимое для успешного детектирования сигнала.

Ключевые слова: оконный спектрально-статистический метод, детектирование сигналов в шуме.

Ostanin S.A.1, Semenov G.A.2, Shayduk A.M.3

1PhD in Physics and Mathematics, Altai state University, 2PhD in Physics and Mathematics, Saint-Petersburg state agrarian University, 3PhD in Physics and Mathematics, Altai state medical university

EFFICIENCY RATIO OF WINDOW SPECTRAL-STATISTICAL METHOD FOR DETECTION SIGNAL IN NOISE

Abstract

The described method of increasing the sensitivity of the detection of short periodic signals in Gaussian noise based on separation of the additive mixture of signal and noise into fragments using a time window. The obtained function of the spectral values a universal criterion for the detection of the number of sections separating the studied time array. Introduced the efficiency factor of the window a spectral-statistical method to determine the optimal number of segments separating an additive mixture of signal and noise necessary for successful signal detection.

Keywords: spectral-statistical window method, detection signal in noise.

С целью повышения чувствительности спектрального детектирования квазипериодических сигналов в гауссовском шуме необходимо исследование возможностей расширения границ его применимости, для чего рассмотрим этапы создания оконного спектрально-статистического метода, основанного на статистическом методе спектрального детектирования [1‑2].

Рассмотрим аддитивную смесь:

12-08-2016 10-21-22 (1)

где среднее квадратическое отклонение шумовой компоненты 12-08-2016 10-23-37Ш(t) равно σШ, относительная длительность квазипериодического сигнала равна α0. Временной ряд, полученный из (1) – числовая последовательность 12-08-2016 10-26-54, представляющая собой сигнал 12-08-2016 10-23-37(t), измеренный в точках 12-08-2016 10-27-36. Число элементов N этой последовательности будем считать достаточно большим, то есть 1. Число членов временного ряда, соответствующего квазипериодической части сигнала равно 12-08-2016 10-30-18. Из анализа критерия обнаружения квазипериодического сигнала статистическим методом спектрального детектирования [2]:

12-08-2016 10-31-00 (2)

и выражения для универсального спектрального критерия [3]: 12-08-2016 10-36-34 (3)

можно увидеть, что допустимая относительная длительность квазипериодической части сигнала  линейно убывает с увеличением его относительной мощности  при неизменном объёме выборки N и пропорциональна величине 12-08-2016 10-37-46 при изменении объёма выборки. Ясно, что если шумовой участок суммарного сигнала статистически однороден на всей временной области его существования, то значение величины η  не зависит от временного интервала, на котором производилось её вычисление.

С учётом того, что величина относительной длительности квазипериодической части зависит от длительности случайной компоненты сигнала, рассмотрим возможные пути увеличения параметра α посредством разбиения исследуемого суммарного сигнала (1) с помощью временного окна на фрагменты (подмассивы в дискретном случае) и последующего построения их спектров и дальнейшего анализа.

В связи с тем, что разбиение сигнала на фрагменты является линейной операцией, рост величины α также будет линейным. Так как при этом уменьшается объём временного ряда, будет увеличиваться граничное значение относительной длительности квазипериодического сигнала  αгр пропорционально множителю 12-08-2016 10-37-46, что частично ухудшит возможности детектирования, однако общая тенденция повышения чувствительности при этом сохранится. Разбиение рассматриваемого временного ряда на подмассивы следует проводить таким образом, чтобы величина шага сдвига временного окна равнялась половине размера окна. Это требуется для того, чтобы в процессе разбиения временного ряда квазипериодическая компонента полностью содержалась бы в одном из подмассивов и при детектировании учитывалась вся её энергия.

Для определения возможностей оконного спектрального подхода далее проведём оценку его эффективности по сравнению со статистическим методом спектрального детектирования.

Разделим массив Χi на (2n-1) подмассивов так, чтобы сдвиг временного окна составлял половину его ширины (объём подмассивов равен размеру временного окна):

12-08-2016 10-43-07 (4)

где Yjkk-тый член j-того подмассива. Значения индексов меняются в диапазонах: 12-08-2016 10-45-27. Величина  12-08-2016 10-46-10 – размер временного окна. При разбиении временного ряда по правилу (4) каждый подмассив представляет собой j-тую строку матрицы:

12-08-2016 10-47-09 (5)

С учётом гипотезы о соотношении длительности квазипериодической компоненты и размера временного окна, принятой ранее, будем считать, что квазипериодический участок сигнала полностью содержится в одном из подмассивов рассматриваемого временного ряда после разделения. То есть теперь для некоторой строки Yj будет проверяться критерий спектрального детектирования статистическим способом, который ранее применялся ко всему временному ряду Χi. Запишем результаты, полученные для статистического метода спектрального детектирования применительно к некоторому единственному подмассиву Yj.

Объём фрагмента сигнала Yj равен, с учётом чего относительная длительность квазипериодического сигнала αY для временного ряда Yj запишется теперь как

12-08-2016 10-50-26 (6)

то есть получаем линейный рост αY пропорционально увеличению количества участков разбиения (параметра n). Пороговое значение длительности квазипериодической компоненты 12-08-2016 10-51-27, при котором теперь реализуемо его успешное обнаружение в нормальном шуме в пределах временного окна, с учётом (1) равно:

12-08-2016 10-52-04 (7)

Значение относительной мощности квазипериодической части сигнала  сохраняется при этом неизменным. Условие успешного детектирования квазипериодической компоненты теперь выражается соотношением:

12-08-2016 10-53-14  (8)

Переходя от (8) к условию для α0 с учётом (6), получаем:

12-08-2016 10-55-27  (9)

то есть пороговое значение относительной длительности квазипериодической части сигнала при его детектировании в нормальном шуме оконным спектральным методом по отношению ко всему сигналу (1) равно:

12-08-2016 10-56-32  (10)

Универсальный спектральный критерий оценки квазипериодического сигнала на предмет возможности его обнаружения в нормальном шуме примет вид:

12-08-2016 10-57-17  (11)

Отношение ξ величин пороговых значений относительной длительности квазипериодической компоненты, требуемых для их успешного обнаружения, соответствующих статистическому методу спектрального детектирования (2) и оконному спектральному методу (11) определим как коэффициент эффективности оконного метода:

12-08-2016 10-58-42 (12)

Коэффициент эффективности ξ определяет число фрагментов разбиения рассматриваемого временного ряда, при котором реализуемо успешное детектирование квазипериодической компоненты в одном из участков. Определим из (12) значение параметра n при фиксированных значениях ξ:

12-08-2016 11-00-13  (13)

Численно решаем уравнение (13) относительно n, находим значение n0, число фрагментов разбиения сигнала при этом определяется как:

12-08-2016 11-01-46   (14)

С целью поиска оптимального значения параметра , соответствующего наиболее эффективному для успешного детектирования квазипериодического сигнала количеству разбиений временного ряда Χi, проведём исследование функции ξ(n) на наличие экстремумов.

12-08-2016 11-03-36 (15)

Из (15) видно, что на рассматриваемом интервале производная ξ(n)  нулей не имеет, всюду положительна, то есть функция ξ(n) возрастает на всей её области определения. При этом ξ(n) имеет область значений:

12-08-2016 11-04-53  (16)

Понятно, что значение параметра 12-08-2016 11-06-05 практически не реализуемо, так как ширина временного окна при этом должна составлять всего два элемента; однако, на всей практически осуществимой области изменения параметра n эффективность оконного метода растёт при увеличении числа фрагментов разбиения сигнала. Верхняя граница параметра n определяется минимальным количеством членов временного ряда, требуемым для корректного задания квазипериодической части и возможностью обнаружения этого сигнала посредством построения его спектра. Также она определяется длительностью квазипериодической компоненты: при завышенном числе фрагментов разбиения сразу несколько подмассивов временного ряда будут составлять квазипериодический участок, что не является необходимым и, более того, потребует выполнения некоторого количества лишних вычислительных операций. Таким образом, применение оконного спектрального метода детектирования при уменьшении относительной длительности квазипериодической компоненты α0 требует увеличения коэффициента эффективности ξ, что в свою очередь требует увеличения числа фрагментов разбиения сигнала.

Главным препятствием здесь является то, что при большом увеличении числа фрагментов разбиения сигнала несколько повышается вероятность ошибочного детектирования, то есть эффективность оконного метода в вероятностном смысле несколько снижается.

Литература

  1. Семёнов Г.А. Поиск сигнала в нормальном шуме методом анализа статистики спектра / Г.А. Семёнов // Известия Алтайского госуниверситета. – 2011. – Ч.1(69) – С. 189 191.
  2. Семёнов Г.А. Выявление порога чувствительности спектрального метода к поиску периодического сигнала в нормальном шуме / Г.А. Семёнов // Известия Алтайского госуниверситета. – 2011. – Ч.1(69) – С.192 196.
  3. Семёнов Г.А. Универсальный спектральный критерий возможности детектирования периодических сигналов в нормальном шуме / С.А. Останин, Г.А. Семёнов, А.М. Шайдук // Труды конференции "Лазеры. Измерения. Информация". – СПб: Изд-во СПб государственного политехнического университета. – 2011. – С.4 5.

References

  1. Semjonov G.A. Poisk signala v normal'nom shume metodom analiza statistiki spektra / G.A. Semjonov // Izvestija Altajskogo gosuniversiteta. – 2011. – Ch.1(69) – S. 189 191.
  2. Semjonov G.A. Vyjavlenie poroga chuvstvitel'nosti spektral'nogo metoda k poisku periodicheskogo signala v normal'nom shume / G.A. Semjonov // Izvestija Altajskogo gosuniversiteta. – 2011. – Ch.1(69) – S.192 196.
  3. Semjonov G.A. Universal'nyj spektral'nyj kriterij vozmozhnosti detektirovanija periodicheskih signalov v normal'nom shume / S.A. Ostanin, G.A. Semjonov, A.M. Shajduk // Trudy konferencii "Lazery. Izmerenija. Informacija". – SPb: Izd-vo SPb gosudarstvennogo politehnicheskogo universiteta. – 2011. – S.4 5.