numerator.m

来自「本压缩包括最小均方滤波器的源码」· M 代码 · 共 9 行

function b=numerator(h,a,k1)
b=zeros(1,k1);
a1=[0 a];
for i=1:k1
    b(i)=h(i);
    for j=1:i
    b(i)=b(i)+a1(j)*h(i+1-j);
    end
end