Принимая частоту дискретизации 20 кГц, создайте фильтр Баттерворта четвертого порядка с частотой 3-dB 2,5 кГц. Фильтрация коэффициентов для butter
должен быть константой для генерации кода.
type ButterFilt
function output_data=ButterFilt(input_data) %#codegen [b,a]=butter(4,0.25); output_data=filter(b,a,input_data); end
Используйте фильтр Баттерворта, чтобы lowpass-фильтровать шумную синусоиду.
t = transpose(linspace(0,pi,10000)); x = sin(t) + 0.03*randn(numel(t),1);
Фильтруйте шумную синусоиду с помощью фильтра Баттерворта. Постройте график фильтрованного сигнала.
fx = ButterFilt(x); plot(fx)
Запуск codegen
команда для получения исходного кода C ButterFilt.c
и Файл MEX:
codegen ButterFilt -args {zeros(10000,1)} -o ButterFilt_mex -report
Code generation successful: To view the report, open('codegen/mex/ButterFilt/html/report.mldatx').
Исходный код C включает пять коэффициентов числителя и знаменателя фильтра Баттерворта четвертого порядка в качестве статических констант. Применить фильтр с помощью файла MEX. Постройте график фильтрованного сигнала.
output_data = ButterFilt_mex(x); hold on plot(output_data) hold off