Фильтрация нулевой фазы

Разработайте lowpass Фильтра Баттерворта с частотой 3-dB 1 кГц, чтобы реализовать фильтрацию нулевой фазы на данных, выбранных со скоростью 20 кГц.

type myZerophaseFilt.m
function output  = myZerophaseFilt(input) %#codegen

[B,A] = butter(20,0.314);
output = filtfilt(B,A,input);

end

Используйте codegen для создания файла MEX для myZerophaseFilt.m.

codegen myZerophaseFilt -args {zeros(1,20001)} -o myZerophaseFilt_mex -report
Code generation successful: To view the report, open('codegen/mex/myZerophaseFilt/html/report.mldatx').

Сгенерируйте шумный синусоидальный сигнал как вход в фильтр.

Fs = 20000;
t = 0:1/Fs:1;
comp500Hz = cos(2*pi*500*t);
signal = comp500Hz + sin(2*pi*4000*t) + 0.2*randn(size(t));

Пропустите входные данные с помощью функций MATLAB ® и MEX.

FilteredData = myZerophaseFilt(signal);
MexFilteredData = myZerophaseFilt_mex(signal);

Постройте график компонента 500 Гц и отфильтрованных данных.

tms = t*1000;
plot(tms,comp500Hz)
hold on
plot(tms,MexFilteredData)
plot(tms,FilteredData)
hold off

xlabel('Milliseconds')
ylabel('Amplitude')
axis([0 25 -1.8 1.8])
legend('500 Hz component','MEX','MATLAB')

Figure contains an axes. The axes contains 3 objects of type line. These objects represent 500 Hz component, MEX, MATLAB.

Для просмотра документации необходимо авторизоваться на сайте