www.gusucode.com > optimdemos 工具箱 matlab 源码程序 > optimdemos/filtfun2.m
function [f,g]=filtfun2(xfree,x,xmask,n,h,dcgain,maxbin) %FILTFUN2 Return frequency response norm and roots for DFILDEMO. % Copyright 1990-2006 The MathWorks, Inc. x(xmask) = xfree; h2=abs(freqz(x(1:n), x(n+1:2*n), 128)); f= h2 - h; f = f'*f; % Make sure its stable g=[abs(roots(x(n+1:2*n)))-1];