filtfun.m

来自「数字通信第四版原书的例程」· M 代码 · 共 8 行

function [f,g]=filtfun(xfree,x,xmask,n,h,maxbin);
%FILTFUN Returns frequency response and roots for DFILDEMO.
x(xmask)  = xfree;
h2=abs(freqz(x(1:n), x(n+1:2*n), 128));
f= h2 - h;     
% Make sure its stable
g=[abs(roots(x(n+1:2*n)))-1];