Коэффициенты рассчитал, вроде понял как по этому алгоритму считать
Коэффициенты для фильтра
частота дискретизации 102400, фильтр режекторный с полосой 4000-5000
такие коэффициенты
%a5 a4 a3 ... a0
a=[0.905475 -8.71601 38.0856 -99.4697 171.921 -205.454 171.921 -99.4697 38.0856 -8.71601 0.905475];
%b5 b4 b3 ... b1
b=[0.819867 -8.04842 35.8696 -95.5494 168.446 -205.333 175.271 -103.451 40.4093 -9.43477];
в matlabe написал функцию для вычисления по этой схеме:
function y=filt(x)
N=1024;
[a,b]=get_k();
yi(1:10)=0;
xi(1:11)=0;
y=[];
for(n=1:N)
xi(1:10)=xi(2:11);
xi(11)=x(n);
yi(11)=sum(a.*xi)+sum(-b.*yi);
y(n)=yi(11);
yi=yi(2:11);
end
end
function [a,b]=get_k
%a5 a4 a3 ... a0
a=[0.905475 -8.71601 38.0856 -99.4697 171.921 -205.454 171.921 -99.4697 38.0856 -8.71601 0.905475];
b=[0.819867 -8.04842 35.8696 -95.5494 168.446 -205.333 175.271 -103.451 40.4093 -9.43477];
end
но что то не то, в этом и вопрос.
Входной сигнал пила с частотой 1500, Объем выборки 1024 вот так выглядит
На выходе такое
1.0e+07 *
Columns 1 through 12
0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
Columns 13 through 24
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 -0.0000 -0.0000
Columns 25 through 36
0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0002 0.0003 0.0005 0.0007 0.0010 0.0015
Columns 37 through 48
0.0023 0.0037 0.0062 0.0108 0.0187 0.0324 0.0554 0.0933 0.1548 0.2537 0.4124 0.6675
Columns 49 through 50
1.0803 1.7534
и т.д. увеличивается в бесконечность.