高通滤波算法

songzy41

一般二阶的差分方程可表示为：
          b(0)+b(1)Z^(-1)+b(2)Z^(-2)
H(Z)=-------------------------------------------
          1+a(1)Z^(-1)+a(2)Z^(-2)
在进行差分方程滤波时，可如下方式进行：
w(2)=0; w(3)=0;
v(2)=0; v(3)=0;
for k=1 : N
w(1)=x(k);
v(1)=b0w(1)+b(1)w(2)+b(2)w(3)-a(1)v(2)-a(2)v(3);
w(3)=w(2); w(2)=w(1);
v(3)=v(2); v(2)=v(1);
y(k)=v(1);
end

frankneu

谢谢所有人的指点，特别感谢songzy41老师的帮助，现在此问题已基本搞定。
主要原因是我所选取的差分方程的参数精度不够，现在改用LongG类型的参数，问题就解决了，谢谢各位了！！！！！！！！

[ 本帖最后由 frankneu 于 2007-1-12 11:13 编辑 ]