Спектральный и корреляционный анализ детерминированных сигналов. Конспект лекции: Корреляция, автокорреляция, взаимная корреляция. Свойства автокорреляционной и взаимной корреляционной функции. Ограничения, накладываемые на вид автокорреляционной функции

В теории связи корреляционная теория используется при исследовании случайных процессов, позволяя установить связь между корреляционными и спектральными свойствами случайных сигналов. Часто возникает задача обнаружения одного передаваемого сигнала в другом или в помехах. Для надежного обнаружения сигналов и применяется метод корреляции , основанный на корреляционной теории. На практике оказывается полезным анализ характеристики, дающей представление о скорости изменения во времени, а также длительности сигнала без разложения его на гармонические составляющие.

Пусть копия сигнала u(t - т) смещена относительно своего оригинала u(t) на интервал времени т. Для количественной оценки степени отличия (связи) сигнала u(t) и его смещенной копии u(t - т) используют автокорреляционную функцию (АКФ). АКФ показывает степень сходства между сигналом и его сдвинутой копией - чем больше значение АКФ, тем это сходство сильнее.

Для детерминированного сигнала конечной длительности (финитного сигнала) аналитическая запись АКФ представляет собой интеграл вида

Формула (2.56) показывает, что при отсутствии сдвига копии относительно сигнала (т = 0) АКФ положительна, максимальна и равна энергии сигнала:

Такая энергия [Дж] выделяется на резисторе с сопротивлением в 1 Ом, если к его выводам подключить некоторое напряжение u(t) [В].

Одним из важнейших свойств АКФ является ее четность: В(т) = В(- т). Действительно, если в выражении (2.56) произвести замену переменной х = t - т, то

Поэтому интеграл (2.56) можно представить в другом виде:

Для периодического сигнала с периодом Г, энергия которого бесконечно велика (поскольку сигнал существует бесконечное время), вычисление АКФ по формуле (2.56) неприемлемо. В этом случае определяют АКФ за период:

Пример 2.3

Определим АКФ прямоугольного импульса, который имеет амплитуду Е и длительность т и (рис. 2.24).

Решение

Для импульса вычисления АКФ удобно провести графически. Такое построение показано на рис. 2.24, а - г, где приведены соответственно исходный импульс u(t) = u t сдвинутая на т его копия м т (?) = u(t - т) = м т и их произведение u(f)u(t - т) = uu v Рассмотрим графическое вычисление интеграла (2.56). Произведение u(t)u(t - т) не равно нулю на интервале времени, когда имеется наложение друг на друга любых частей сигнала и его копии. Как следует из рис. 2.24, этот интервал равен х - т м, если временной сдвиг копии меньше длительности импульса. В подобных случаях для импульса АКФ определится как В(т) = Е 2 (т и - |т|) при временном сдвиге копии на текущее время |т| В(0) = = Е 2 т и = Э (см. рис. 2.24, г).

Рис. 2.24.

а - импульс; 6 - копия; в - произведение сигнала и копии; г - АКФ

Часто вводят удобный для анализа и сравнения сигналов числовой параметр - интервал корреляции т к, аналитически и графически равный ширине основания АКФ. Для данного примера интервал корреляции т к = 2т и.

Пример 2.4

Определим АКФ гармонического (косинусоидального) сигнала u(t) = = t/ m cos(co? + а).

Рис. 2.25.

а - гармонический сигнал; б - АКФ гармонического сигнала

Решение

Используя формулу (2.57) и обозначив В п (т) = В(т), находим

Из этой формулы следует, что АКФ гармонического сигнала тоже является гармонической функцией (рис. 2.25, б) и имеет размерность мощности (В 2). Отметим еще один очень важный факт, что вычисленная АКФ не зависит от начальной фазы гармонического сигнала (параметр

Из проведенного анализа следует важный вывод: АКФ практически любого сигнала не зависит от его фазового спектра. Следовательно, сигналы, амплитудные спектры которых полностью совпадают, а фазовые различаются, будут иметь одинаковую АКФ. Еще одно замечание заключается в том, что по АКФ нельзя восстановить исходный сигнал (опять же вследствие утраты информации о фазе).

Связь между АКФ и энергетическим спектром сигнала. Пусть импульсный сигнал u(t) имеет спектральную плотность 5(со). Определим АКФ но формуле (2.56), записав и(С) в виде обратного преобразования Фурье (2.30):

Введя новую переменную х = t - т, из последней формулы получим Здесь интеграл

есть функция, комплексно-сопряженная спектральной плотности сигнала

С учетом соотношения (2.59) формула (2.58) примет вид Функцию

называют энергетическим спектром (спектральной плотностью энергии) сигнала, показывающим распределение энергии по частоте. Размерность энергетического спектра сигнала соответствует величине IP/со) - [(В 2 -с)/Гц].

Учитывая соотношение (2.60), окончательно получим выражение для АКФ:

Итак, АКФ сигнала представляет собой обратное преобразование Фурье от его энергетического спектра. Прямое преобразование Фурье от АКФ

Итак, прямое преобразование Фурье (2.62) АКФ определяет энергетический спектр, а обратное преобразование Фурье энергетического спектра (2.61) - АКФ детерминированного сигнала. Эти результаты важны по двум причинам. Во-первых, исходя из распределения энергии но спектру становится возможным оценить корреляционные свойства сигналов - чем шире энергетический спектр сигнала, тем меньше интервал корреляции. Соответственно, чем больше интервал корреляции сигнала, тем короче его энергетический спектр. Во-вторых, соотношения (2.61) и (2.62) позволяют экспериментально определить одну из функций по значению другой. Часто удобнее вначале получить АКФ, а затем с помощью прямого преобразования Фурье вычислить энергетический спектр. Этот прием широко применяют при анализе свойств сигналов в реальном масштабе времени, т.е. без временной задержки при его обработке.

Взаимокорреляционная функция двух сигналов. Если надо оценить степень связи между сигналами u x (t) и u 2 (t), то используют взаимокорреля- ционную функцию (ВКФ)

При т = О ВКФ равна так называемой взаимной энергии двух сигналов

Значение ВКФ не меняется, если вместо задержки второго сигнала u 2 (t) рассматривать опережение его первым сигналом м,(?), поэтому

АКФ является частным случаем ВКФ, если сигналы одинаковы, т.е. u y (t) = u 2 (t) = u(t). В отличие от АКФ ВКФ двух сигналов В 12 (т) не является четной и необязательно максимальна при т = 0, т.е. при отсутствии временного сдвига сигналов.

Наряду со спектральным подходом к описанию сигналов часто на практике оказывается необходимой характеристика, которая давала бы представление о некоторых свойствах сигнала, в частности о скорости изменения во времени, а также о длительности сигнала без разложения его на гармонические составляющие.

В качестве такой временной характеристики широко используется корреляционная функция сигнала.

Для детерминированного сигнала s (t ) конечной длительности корреляционная функция определяется следующим выражением:

где τ - временной сдвиг сигнала.

В данной главе рассматриваются сигналы, являющиеся вещественными функциями времени, и обозначение комплексного сопряжения можно опустить:

. (1.78)

Из выражения (1.78) видно, что B s (t ) характеризует степень связи (корреляции) сигналаs ( t ) со своей копией, сдвинутой на величину т по оси времени. Ясно, что функцияB s ( t ) достигает максимума при τ = 0, так как любой сигнал полностью коррелирован с самим собой. При этом

, (1.79)

т. е. максимальное значение корреляционной функции равно энергии сигнала.

С увеличением τ функция В 8 (τ) убывает (не обязательно монотонно) и при относительном сдвиге сигналовs (t ) иs (t + τ) на время, превышающее длительность сигнала, обращается в нуль.

Из общего определения корреляционной функции видно, что безразлично, вправо или влево относительно своей копии сдвигать сигнал на величину τ. Поэтому выражение (1.78) можно обобщить следующим образом:

. (1.78)

Это равносильно утверждению, что B s (τ) являетсячетной функцией τ.

Для периодического сигнала, энергия которого бесконечно велика, определение корреляционной функции с помощью выражений (1.129) или (1.129") неприемлемо. В этом случае исходят из следующего определения:

При таком определении корреляционная функция приобретает размерность мощности, причем B Sne р (0) равна средней мощности периодического сигнала. Ввиду периодичности сигналаs( t ) усреднение произведения
или
по бесконечно большому отрезкуТ должно совпадать с усреднением по периодуT 1 . Поэтому выражение (1.79) можно заменить выражением

Входящие в это выражение интегралы суть не что иное, как корреляционная функция сигнала на интервале T 1 . Обозначая ее через B sTl (τ ), приходим к соотношению

Очевидно также, что периодическому сигналу s(t ) соответствует и периодическая корреляционная функцияB s пер (τ). Период функцииB s пер (τ) совпадает с периодомТ 1 исходного сигналаs( t ). Например, для простейшего (гармонического) колебания
корреляционная функция

При τ=0
есть средняя мощность гармонического колебания с амплитудойА 0 . Важно отметить, что корреляционная функция
не зависит от начальной фазы колебания.

Для оценки степени связи между двумя различными сигналами s 1 ( t ) иs 2 ( t ) используется взаимная корреляционная функция, определяемая общим выражением

Для вещественных функций s 1 (t) иs 2 (t)

Рассмотренная выше корреляционная функция В s (τ) является частным случаем функции
, когдаs 1 ( t ) =s 2 ( t ).

В отличие от
взаимная корреляционная функция не обязательно является четной относительно τ. Кроме того, взаимная корреляционная функцияне обязательно достигает максимума приτ = 0.

3 Корреляционный анализ сигналов

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

3.1 Автокорреляционная функция (АКФ)

Автокорреляционная функция сигнала с конечной энергией – это значение интеграла от произведения двух копий этого сигнала, сдвинутых относительно друг друга на время τ, рассматриваемое в функции этого временного сдвига τ:

Если сигнал определен на конечном интервале времени , то его АКФ находится как:

где - интервал перекрытия сдвинутых копий сигнала.

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

Вводимая таким образом мера сходства для сигналов, имеющих форму случайных колебаний вокруг нулевого значения, обладает следующими характерными свойствами.

Если сдвинутые копии сигнала колеблются примерно в такт друг к другу, то это является признаком их схожести и АКФ принимает большие положительные значения (большая положительная корреляция). Если копии колеблются почти в противофазе, АКФ принимает большие отрицательные значения (антисходство копий сигнала, большая отрицательная корреляция).

Максимум АКФ достигается при совпадении копий, то есть при отсутствии сдвига. Нулевые значения АКФ достигаются при сдвигах, при которых не заметно ни сходства, ни антисходства копий сигнала (нулевая корреляция,

отсутствие корреляции).

На рис.3.1 изображен фрагмент реализации некоторого сигнала на интервале времени от 0 до 1 с. Сигнал случайным образом колеблется вокруг нулевого значения. Поскольку интервал существования сигнала конечен, то конечна и его энергия. Его АКФ можно вычислить в соответствии с уравнением:

Автокорреляционная функция сигнала, вычисленная в MathCad в соответствии с этим уравнением, представлена на рис. 3.2. Корреляционная функция показывает не только то, что сигнал похож сам на себя (сдвиг τ=0), но и то, что некоторой схожестью обладают и копии сигнала, сдвинутые друг относительно друга приблизительно на 0.063 с (боковой максимум автокорреляционной функции). В противоположность этому копии сигнала сдвинутые на 0.032 с, должны быть антипохожи дуг на друга, то есть быть в некотором смысле противоположными друг другу.

На рис.33 показаны пары этих двух копий. По рисунку можно проследить, что понимается под похожестью и антипохожестью копий сигнала.

Корреляционная функция обладает следующими свойствами:

1. При τ = 0 автокорреляционная функция принимает наибольшее значение, равное энергии сигнала

2. Автокорреляционная функция является четной функцией временного сдвига .

3. С ростом τ автокорреляционная функция убывает до нуля

4. Если сигнал не содержит разрывов типа δ - функций, то - непрерывная функция.

5. Если сигнал является электрическим напряжением, то корреляционная функция имеет размерность .

Для периодических сигналов в определении автокорреляционной функции тот же самый интеграл делят еще на период повторения сигнала:

Так введенная корреляционная функция отличается следующими свойствами:

Значение корреляционной функции в нуле равно мощности сигнала ,

Размерность корреляционной функции равна квадрату размерности сигнала, например .

Для примера вычислим корреляционную функцию гармонического колебания :

Используя ряд тригонометрических преобразований, получим окончательно:

Таким образом, автокорреляционная функция гармонического колебания является косинусоидой с тем же периодом изменения, что и сам сигнал. При сдвигах, кратных периоду колебания, гармоника преобразуется в себя и АКФ принимает наибольшие значения, равные половине квадрата амплитуды. Сдвиги по времени, кратные половине периода колебания, равносильны смещению фазы на угол , при этом меняется знак колебаний, а АКФ принимает минимальное значение, отрицательное и равное половине квадрата амплитуды. Сдвиги, кратные четверти периода, переводят, например, синусоидальное колебание в косинусоидальное и наоборот. При этом АКФ обращается в нуль. Такие сигналы, находящиеся в квадратуре друг относительно друга, с точки зрения автокорреляционной функции оказываются совершенно не похожими друг на друга.

Важным является то, что в выражение для корреляционной функции сигнала не вошла его начальная фаза. Информация о фазе потерялась. Это означает, что по корреляционной функции сигнала нельзя восстановить сам сигнал. Отображение в противоположность отображению не является взаимно однозначным.

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

Актом проявления сигналом своей свободной воли, независимой от воли создателя (возникновение отдельных реализаций некоторого случайного процесса),

Результатом постороннего насилия над сигналом (введение в сигнал измерительной информации, получаемой при проведении измерений какой либо физической величины).

Аналогичным образом обстоит дело с любым периодическим сигналом. Если периодический сигнал с основным периодом Т имеет амплитудный спектр и фазовый спектр , то корреляционная функция сигнала принимает следующий вид:

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

3.2 Взаимнокорреляционная функция (ВКФ).

В отличие от автокорреляционной функции взаимнокорреляционная функция определяет степень схожести копий двух различных сигналов x(t) и y(t), сдвинутых на время τ друг относительно друга:

Взаимнокорреляционная функция обладает следующими свойствами:

1. При τ = 0 взаимнокорреляционная функция принимает значение, равное взаимной энергии сигналов, то есть энергии их взаимодействия

2. При любом τ имеет место соотношение:

где - энергии сигналов.

3. Изменение знака временного сдвига равносильно взаимной перестановке сигналов:

4. С ростом τ взаимнокорреляционная функция хотя и не монотонно, но убывает до нуля

5. Значение взаимнокорреляционной функции в нуле ничем не выделяется среди других значений.

Для периодических сигналов понятие взаимнокорреляционной функции, как правило, вообще не используется.

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

Статистический анализ электроэнцефалограмм и других результатов регистрации биопотенциалов,

Определение пространственных координат источника сигнала по величине временного сдвига, при котором достигается максимум ВКФ,

Выделение слабого сигнала на фоне сильных статических несвязанных помех,

Обнаружение и локализация каналов утечки информации путем определения корреляции между сигналами радиоэфира в помещении и за его пределами,

Автоматизированное обнаружение в ближней зоне, распознавание и поиск работающих радиоизлучающих подслушивающих устройств, включая мобильные телефоны, используемые как подслушивающие устройства,

Локализация мест утечек в трубопроводах на основании определения ВКФ двух сигналов акустического шума, вызываемого утечкой, в двух точках измерения, в которых расположены датчики на трубе.

3.3 Соотношения между корреляционными и спектральными функциями.

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

Взаимная корреляционная функция, как и любая другая функция времени, может быть подвергнута преобразованию Фурье:

Изменим порядок интегрирования:

Выражение в квадратных скобках можно было бы рассматривать как преобразование Фурье для сигнала y(t), но в показателе экспоненты не стоит знак минус. Это говорит о том, что внутренний интеграл дает нам выражение , комплексно сопряженное со спектральной функцией .

Но выражение не зависит от времени, поэтому его можно вынести за знак внешнего интеграла. Тогда внешний интеграл просто даст нам определение спектральной функции сигнала x(t). Окончательно имеем:

Это означает, что преобразование Фурье для взаимной корреляционной функции двух сигналов равно произведению их спектральных функций, одна из которых подвергнута комплексному сопряжению. Это произведение называется взаимным спектром сигналов:

Из полученного выражения следует важный вывод: если спектры сигналов x(t) и y(t) не перекрывают друг друга, то есть располагаются в различных диапазонах частот, то такие сигналы являются некоррелированными, независимыми друг от друга.

Если положить в приведенных формулах: x(t) = y(t), то получим выражение для преобразования Фурье автокорреляционной функции

Это означает, что автокорреляционная функция сигнала и квадрат модуля его спектральной функции связаны друг с другом посредством преобразования Фурье.

Функция называется энергетическим спектром сигнала . Энергетический спектр показывает, как общая энергия сигнала распределяется по частотам его отдельных гармонических составляющих.

3.4 Энергетические характеристики сигналов с частотной области

Взаимная корреляционная функция двух сигналов связана преобразованием Фурье с взаимным спектром сигналов, поэтому ее можно выразить в виде обратного преобразования Фурье от взаимного спектра:

Теперь подставим в эту цепочку равенств значение временного сдвига . В результате получим соотношение, которое определяет смысл равенства Релея :

то есть интеграл от произведения двух сигналов равен интегралу от произведения спектров этих сигналов, один из которых подвергнут операции комплексного сопряжения.

Это соотношение называется равенством Парсеваля .

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

Это соотношение обладает полной аналогией с равенством Парсеваля.

Функции корреляции сигналов применяются для интегральных количественных оценок формы сигналов и степени их сходства друг с другом.

Автокорреляционные функции (АКФ) сигналов (correlation function, CF). Применительно к детерминированным сигналам с конечной энергией АКФ является количественной интегральной характеристикой формы сигнала, и представляет собой интеграл от произведения двух копий сигнала s(t), сдвинутых относительно друг друга на время t:

B s (t) = s(t) s(t+t) dt. (2.4.1)

Как следует из этого выражения, АКФ является скалярным произведением сигнала и его копии в функциональной зависимости от переменной величины значения сдвига t. Соответственно, АКФ имеет физическую размерность энергии, а при t = 0 значение АКФ непосредственно равно энергии сигнала и является максимально возможным (косинус угла взаимодействия сигнала с самим собой равен 1):

B s (0) = s(t) 2 dt = E s .

Функция АКФ является непрерывной и четной. В последнем нетрудно убедиться заменой переменной t = t-t в выражении (2.4.1):

B s (t) = s(t) s(t-t) dt = s(t-t) s(t) dt = B s (-t).

С учетом четности, графическое представление АКФ обычно производится только для положительных значений t. Знак +t в выражении (2.4.1) означает, что при увеличении значений t от нуля копия сигнала s(t+t) сдвигается влево по оси t. На практике сигналы обычно также задаются на интервале положительных значений аргументов от 0-Т, что дает возможность продления интервала нулевыми значениями, если это необходимо для математических операций. В этих границах вычислений более удобным является сдвиг копии сигнала влево по оси аргументов, т.е. применение в выражении (2.4.1) функции s(t-t):

B s (t) = s(t) s(t-t) dt. (2.4.1")

По мере увеличения значения величины сдвига t для финитных сигналов временное перекрытие сигнала с его копией уменьшается, а, соответственно, косинус угла взаимодействия и скалярное произведение в целом стремятся к нулю:

Пример. На интервале (0,Т) задан прямоугольный импульс с амплитудным значением, равным А. Вычислить автокорреляционную функцию импульса.

При сдвиге копии импульса по оси t вправо, при 0≤t≤T сигналы перекрываются на интервале от t до Т. Скалярное произведение:

B s (t) = A 2 dt = A 2 (T-t).

При сдвиге копии импульса влево, при -T≤t<0 сигналы перекрываются на интервале от 0 до Т-t. Скалярное произведение:

B s (t) = A 2 dt = A 2 (T+t).

При |t| > T сигнал и его копия не имеют точек пересечения и скалярное произведение сигналов равно нулю (сигнал и его сдвинутая копия становятся ортогональными).

Обобщая вычисления, можем записать:

B s (t) = .

В случае периодических сигналов АКФ вычисляется по одному периоду Т, с усреднением скалярного произведения и его сдвинутой копии в пределах этого периода:

B s (t) = (1/Т) s(t) s(t-t) dt.

При t=0 значение АКФ в этом случае равно не энергии, а средней мощности сигналов в пределах интервала Т. АКФ периодических сигналов при этом также является периодической функцией с тем же периодом Т. Так, для сигнала s(t) = A cos(w 0 t+j 0) при T=2p/w 0 имеем:

B s (t) = A cos(w 0 t+j 0) A cos(w 0 (t-t)+j 0) = (A 2 /2) cos(w 0 t).

Отметим, что полученный результат не зависит от начальной фазы гармонического сигнала, что характерно для любых периодических сигналов и является одним из свойств КФ.

Для сигналов, заданных на определенном интервале , вычисление АКФ также производится с нормировкой на длину интервала :

B s (t) = s(t) s(t+t) dt. (2.4.2)

В пределе, для непериодических сигналов с измерением АКФ на интервале Т:

B s (t) = . (2.4.2")

Автокорреляция сигнала может оцениваться и коэффициентом автокорреляции, вычисление которого производится по формуле (по центрированным сигналам):

r s (t) = cos j(t) = ás(t), s(t+t)ñ /||s(t)|| 2 .

Взаимная корреляционная функция (ВКФ) сигналов (cross-correlation function, CCF) показывает степень сходства сдвинутых экземпляров двух разных сигналов и их взаимное расположение по координате (независимой переменной), для чего используется та же формула (2.4.1), что и для АКФ, но под интегралом стоит произведение двух разных сигналов, один из которых сдвинут на время t:

B 12 (t) = s 1 (t) s 2 (t+t) dt. (2.4.3)

При замене переменной t = t-t в формуле (2.4.3), получаем:

B 12 (t) = s 1 (t-t) s 2 (t) dt = s 2 (t) s 1 (t-t) dt = B 21 (-t)

Отсюда следует, что для ВКФ не выполняется условие четности, а значения ВКФ не обязаны иметь максимум при t = 0. Это можно наглядно видеть на рис. 2.4.1, где заданы два одинаковых сигнала с центрами на точках 0.5 и 1.5. Вычисление по формуле (2.4.3) с постепенным увеличением значений t означает последовательные сдвиги сигнала s2(t) влево по оси времени (для каждого значения s1(t) для подынтегрального умножения берутся значения s2(t+t)).

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

Во втором случае, с целью проверки эффективности фильтрации, сначала рассчитывается корреляционная функция полезного эталонного сигнала, представленного числовым рядом, и отфильтрованного сигнала. После чего путем применения прямого дискретного преобразования Фурье к корреляционной функции получают коррелограмму. На полученном графике строят линию критического уровня с учетом ошибки фильтрации с использованием критерия Стьюдента. Эффективность фильтрации определяют визуально: выше критического уровня должны находиться только составляющие спектральной плотности полезного сигнала.

Для большей наглядности и объективности рассчитывается выборочный коэффициент корреляции между числовыми рядами эталонного (исходного полезного) и отфильтрованного сигналов. Коэффициент корреляции может принимать значения в интервале –1…1. Отрицательные значения говорят о том, что эталонный и отфильтрованный сигналы коррелируют в противофазе, т.е. при инверсии отфильтрованного сигнала. В случае если цифровой фильтр обладает хорошей эффективностью фильтрации от помех и шумов, коэффициент корреляции принимает значения, близкие к 1 или –1. Качество разных цифровых фильтров применительно к конкретному сигналу может быть определено путем сравнения рассчитанных коэффициентов корреляции.

Расчет корреляционной функции дискретных сигналов производится следующим образом. Для дискретных сигналов Х(i) и Y(i), i = 1… N выбирается фрагмент массива Y(i), i = 1… N/2 и рассчитывается корреляционная функция

где – величина сдвига в дискретах.

Коррелограмму или спектр корреляционной функции получают путем применения прямого дискретного преобразования Фурье к корреляционной функции:

- действительная часть спектра

;