Mercurial > hg > camir-aes2014
annotate 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 |
rev | line source |
---|---|
wolffd@0 | 1 % filter = SecondOrderFilter(f, q, fs) |
wolffd@0 | 2 % Design a second order digital filter with a center frequency of |
wolffd@0 | 3 % f, filter quality of q, and digital sampling rate of fs (Hz). |
wolffd@0 | 4 |
wolffd@0 | 5 % (c) 1998 Interval Research Corporation |
wolffd@0 | 6 |
wolffd@0 | 7 function filts = SecondOrderFilter(f,q,fs) |
wolffd@0 | 8 cft = f'/fs; |
wolffd@0 | 9 rho = exp(- pi * cft ./ q'); |
wolffd@0 | 10 theta = 2 * pi * cft .* sqrt(1-1 ./(4*q'.^2)); |
wolffd@0 | 11 filts = [ ones(size(rho)) -2*rho.*cos(theta) rho.^2 ]; |