home *** CD-ROM | disk | FTP | other *** search
- function [f,g]=filtfun2(xfree,x,xmask,n,h,dcgain,maxbin);
- %FILTFUN2 Return frequency response norm and roots for DFILDEMO.
- 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];
-