Цифровая фильтрация с нулевой фазой
y = filtfilt(b,a,x)xкак в прямом, так и в обратном направлениях. После фильтрации данных в прямом направлении filtfilt изменяет на противоположную отфильтрованную последовательность и запускает ее обратно через фильтр. Результат имеет следующие признаки:
Искажение нулевой фазы.
Передаточная функция фильтра равна квадрату исходной передаточной функции фильтра.
Порядок фильтра, вдвое превышающий порядок фильтра, указанный в b и a.
filtfilt минимизирует пусковые и завершающие переходные процессы путем согласования исходных условий. Не использовать filtfilt с фильтрами дифференциатора и FIR Гильберта, поскольку работа этих фильтров сильно зависит от их фазовой характеристики.
y = filtfilt(d,x)x, используя цифровой фильтр, d. Использовать designfilt произвести d на основе характеристик частотно-отклика.
Нулевая фазовая фильтрация помогает сохранять признаки в отфильтрованном временном сигнале именно там, где они возникают в нефильтрованном сигнале.
Использовать filtfilt для нулевого фильтра формы волны синтетической электрокардиограммы (ЭКГ). Функция, генерирующая сигнал, находится в конце примера. Комплекс QRS является важной особенностью ЭКГ. Здесь он начинается примерно в момент времени 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')
Повреждение ЭКГ с помощью аддитивного шума. Сбросьте генератор случайных чисел для воспроизводимых результатов. Построить низкочастотный КИХ-эквиппльный фильтр и фильтровать шумовую форму сигнала, используя как нулевую фазу, так и обычную фильтрацию.
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.
Повторите вышеизложенное с помощью фильтра второго порядка Butterworth.
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')
Эта функция генерирует сигнал ЭКГ.
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] Густафссон, Ф. «Определение начальных состояний при фильтрации вперед-назад». Транзакции IEEE ® при обработке сигналов. том 44, апрель 1996 года, стр. 988-992.
[2] Митра, Санджит К. Цифровая обработка сигналов. 2-й ред. Нью-Йорк: Макгроу-Хилл, 2001.
[3] Оппенгейм, Алан В., Рональд В. Шефер и Джон Р. Бак. Дискретно-временная обработка сигналов. 2-я эд. река Верхнее Седло, Нью-Джерси: Прентис Холл, 1999.
designfilt | digitalFilter | fftfilt | filter | filter2