Расчет цифрового фильтра

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

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

Первые 256 отсчетов полученного белого шума

Пропустив заданный белый шум X через фильтр, получили некую его реакцию Y.

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

Практически, используя MathCad, СПМ можно получить следующим образом

где - функция быстрого преобразования Фурье,

- комплексное сопряженное .

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

Полученный спектр плотности мощности после сглаживания окном Гаусса

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

Другие стьтьи в тему

Разработка преобразователя разности фаз в постоянное напряжение
Разработка преобразователя разности фаз двух сигналов в постоянное напряжение со следующими параметрами: · Частота входных сигналов 10кГц - 100 кГц; · Входное напряжение 50мВ - 5В; · Диапазон измерения ∆φ 0 - 180о ...

Разработка измерительного канала температуры на основе бесконтактных методов
Без грамотного построения измерительного канала невозможно построить систему автоматического регулирования и управления технологическим процессом (АСУТП). Данный курсовой проект дает возможность не только в теории, но на практике познать сущность проблемы проектирования измерительных ...