Re: Проверка фильтра
Добавлено: 26 дек 2017, 16:27
Нашёл на другом форуме вот такой код:
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));
Можно ли его под мой фильтр переделать?
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));
Можно ли его под мой фильтр переделать?