Mercurial > hg > camir-aes2014
diff toolboxes/MIRtoolbox1.3.2/AuditoryToolbox/SecondOrderFilter.m @ 0:e9a9cd732c1e tip
first hg version after svn
author | wolffd |
---|---|
date | Tue, 10 Feb 2015 15:05:51 +0000 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/toolboxes/MIRtoolbox1.3.2/AuditoryToolbox/SecondOrderFilter.m Tue Feb 10 15:05:51 2015 +0000 @@ -0,0 +1,11 @@ +% filter = SecondOrderFilter(f, q, fs) +% Design a second order digital filter with a center frequency of +% f, filter quality of q, and digital sampling rate of fs (Hz). + +% (c) 1998 Interval Research Corporation + +function filts = SecondOrderFilter(f,q,fs) +cft = f'/fs; +rho = exp(- pi * cft ./ q'); +theta = 2 * pi * cft .* sqrt(1-1 ./(4*q'.^2)); +filts = [ ones(size(rho)) -2*rho.*cos(theta) rho.^2 ];