Страница 1 из 1

Нужно переделать под БИХ фильтр

Добавлено: 27 дек 2017, 15:38
iliaW
Привет, задача такая: нужно переделать цикл в этом коде под БИХ фильтр, сейчас он сделан под КИХ, и почему то вообще не хочет пропускать синусы через заданное разностное уравнение. Помогите пожалуйста

Fs = 44100;
t = (0:2047)/Fs;
x = sin(2*pi*1000*t)+sin(2*pi*2000*t)+sin(2*pi*3000*t);

b =[10.23 6 7 4 1];
a =[1.04199 2.13426 -1.19019 0.11856 -0.02964];
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) - y(i-j+1) * a(j);
end
n = n + 1;
end
figure; plot(y);