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