НО! при реализации в железе фильтр не получается. Что-бы разобраться, в матлабе делаю фильтрацию, как в железе. В результате- тоже не получился фильтр.
Структуру фильтра взял из хелпа матлаба(на картинке).
Ниже привел свой кусок кода для фильтрации.
Вопрос:
как правильно построить БИХ фильтр?
http://nukeuploads.com/getfile/nVIGg8nNIolEYZUlpTaU/1.gif
Код
OUTs=filter(b1,a1,INs);
ii=3;
OUTs11(1)=0;
OUTs11(2)=0;
OUTs11(3)=0;
Z1=0;
Z2=0;
while ii<=length(INs);
OUTs11(ii)=(INs(ii)*b1(1)+Z1)/a1(1);
Z1=INs(ii)*b1(2)-OUTs11(ii-1)*a1(2)+Z2;
Z2=INs(ii)*b1(3)-OUTs11(ii-1)*a1(3);
ii=ii+1;
end;
ii=3;
OUTs11(1)=0;
OUTs11(2)=0;
OUTs11(3)=0;
Z1=0;
Z2=0;
while ii<=length(INs);
OUTs11(ii)=(INs(ii)*b1(1)+Z1)/a1(1);
Z1=INs(ii)*b1(2)-OUTs11(ii-1)*a1(2)+Z2;
Z2=INs(ii)*b1(3)-OUTs11(ii-1)*a1(3);
ii=ii+1;
end;