Проверка фильтра

Все что касается фильтрации
iliaW
Сообщения: 8
Зарегистрирован: 18 дек 2017, 09:35

Re: Проверка фильтра

Сообщение iliaW »

Нашёл на другом форуме вот такой код:

clear all; close all; clc;

Fs = 44100;
t = (0:2047)/Fs;
x = 5*sin(2*pi*5000*t)+5*sin(2*pi*10000*t)+5*sin(2*pi*15000*t);
sp = abs(fft(x));
plot(sp(1:end/2));

b=[0.8 0.2];
figure; freqz(b);
% b = [ 1 3 -2 4 -3 -2 1 ];

n = 1;
y = zeros(1, length(x) );
% свертка
for i = length(b):length(x)
y(n) = 0;
for j = 1:length(b)
y(n) = y(n) + x(i-j+1) * b(j);
end
n = n + 1;
end

sp = abs(fft(y));
figure; plot(sp(1:end/2));

Можно ли его под мой фильтр переделать?

Ответить