Применение Lowpass к входному сигналу

Принимая частоту дискретизации 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)

Figure contains an axes. The axes contains an object of type line.

Запуск 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

Figure contains an axes. The axes contains 2 objects of type line.

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