Daniel@0: % mag = FreqResp(filter, f, fs)
Daniel@0: % Find the frequency response (in dB) of a filter (1x5 vector) at
Daniel@0: % frequency f and sampling rate fs.
Daniel@0: 
Daniel@0: % (c) 1998 Interval Research Corporation  
Daniel@0: 
Daniel@0: function mag=FreqResp(filter,f,fs)
Daniel@0: cf = exp(i*2*pi*f/fs);
Daniel@0: mag = (filter(3) + filter(2)*cf + filter(1)*cf.^2) ./ ...
Daniel@0: 	(filter(5) + filter(4)*cf + cf.^2);
Daniel@0: mag = 20*log10(abs(mag));