Нулевая фаза цифровая фильтрация
y = filtfilt(b,a,x)x, и в прямых и в обратных направлениях. После фильтрации данных в прямом направлении, filtfilt инвертирует отфильтрованную последовательность и запускает его назад через фильтр. Результат имеет следующие характеристики:
Нулевое искажение фазы.
Передаточная функция фильтра равняется величине в квадрате исходной передаточной функции фильтра.
Порядок фильтра, который удваивает порядок фильтра, заданного b и a.
filtfilt минимизирует запуск и конечные переходные процессы путем соответствия с начальными условиями. Не используйте filtfilt с дифференциатором и Гильбертовыми КИХ-фильтрами, потому что работа этих фильтров зависит в большой степени от их фазового отклика.
y = filtfilt(d,x)x, с помощью цифрового фильтра, d. Используйте designfilt сгенерировать d на основе спецификаций частотной характеристики.
Фильтрация нулевой фазы помогает сохранить функции в отфильтрованной форме волны времени точно, где они происходят в неотфильтрованном сигнале.
Используйте filtfilt к нулевой фазе фильтруют синтетическую электрокардиограмму (ECG) форма волны. Функция, которая генерирует форму волны, в конце примера. Комплекс QRS является важной функцией в ECG. Здесь это начинается вокруг момента времени 160.
wform = ecg(500); plot(wform) axis([0 500 -1.25 1.25]) text(155,-0.4,'Q') text(180,1.1,'R') text(205,-1,'S')
Повредите ECG с аддитивным шумом. Сбросьте генератор случайных чисел для восстанавливаемых результатов. Создайте КИХ lowpass equiripple, фильтруют и фильтруют шумную форму волны с помощью и нулевой фазы и обычной фильтрации.
rng default x = wform' + 0.25*randn(500,1); d = designfilt('lowpassfir', ... 'PassbandFrequency',0.15,'StopbandFrequency',0.2, ... 'PassbandRipple',1,'StopbandAttenuation',60, ... 'DesignMethod','equiripple'); y = filtfilt(d,x); y1 = filter(d,x); subplot(2,1,1) plot([y y1]) title('Filtered Waveforms') legend('Zero-phase Filtering','Conventional Filtering') subplot(2,1,2) plot(wform) title('Original Waveform')
Фильтрация нулевой фазы уменьшает шум в сигнале и сохраняет комплекс QRS одновременно, это происходит в оригинале. Обычная фильтрация уменьшает шум в сигнале, но задерживает комплекс QRS.
Повторите вышеупомянутое использование фильтра секции второго порядка Баттерворта.
d1 = designfilt('lowpassiir','FilterOrder',12, ... 'HalfPowerFrequency',0.15,'DesignMethod','butter'); y = filtfilt(d1,x); subplot(1,1,1) plot(x) hold on plot(y,'LineWidth',3) legend('Noisy ECG','Zero-Phase Filtering')
Это - функция, которая генерирует форму волны ECG.
function x = ecg(L) %ECG Electrocardiogram (ECG) signal generator. % ECG(L) generates a piecewise linear ECG signal of length L. % % EXAMPLE: % x = ecg(500).'; % y = sgolayfilt(x,0,3); % Typical values are: d=0 and F=3,5,9, etc. % y5 = sgolayfilt(x,0,5); % y15 = sgolayfilt(x,0,15); % plot(1:length(x),[x y y5 y15]); % Copyright 1988-2002 The MathWorks, Inc. a0 = [0,1,40,1,0,-34,118,-99,0,2,21,2,0,0,0]; % Template d0 = [0,27,59,91,131,141,163,185,195,275,307,339,357,390,440]; a = a0 / max(a0); d = round(d0 * L / d0(15)); % Scale them to fit in length L d(15)=L; for i=1:14, m = d(i) : d(i+1) - 1; slope = (a(i+1) - a(i)) / (d(i+1) - d(i)); x(m+1) = a(i) + slope * (m - d(i)); end end
[1] Оппенхейм, Алан V, Рональд В. Шафер и Джон Р. Бак. Обработка сигналов дискретного времени. 2-й Эд. Верхний Сэддл-Ривер, NJ: Prentice Hall, 1999.
[2] Mitra, цифровая обработка сигналов Сэнджита К. 2-й Эд. Нью-Йорк: McGraw-Hill, 2001.
[3] Густафсон, F. “Определяя начальные состояния в прямой обратной фильтрации”. IEEE® Transactions на Обработке сигналов. Издание 44, апрель 1996, стр 988–992.
designfilt | digitalFilter | fftfilt | filter | filter2
