Return a,b coefficients for a real pole designed using the bilinear transform. usage: filt = rpole(p, fs) p - pole object fs - the sample rate for the filter M Hewitson 03-04-07 $Id: rp2iir.html,v 1.1 2007/06/08 14:15:07 hewitson Exp $
0001 function [a,b] = rp2iir(p, fs) 0002 0003 % Return a,b coefficients for a real pole 0004 % designed using the bilinear transform. 0005 % 0006 % usage: filt = rpole(p, fs) 0007 % 0008 % p - pole object 0009 % fs - the sample rate for the filter 0010 % 0011 % M Hewitson 03-04-07 0012 % 0013 % $Id: rp2iir.html,v 1.1 2007/06/08 14:15:07 hewitson Exp $ 0014 % 0015 0016 f0 = get(p, 'f'); 0017 w0 = f0*2*pi; 0018 a(1) = w0 / (2*fs + w0); 0019 a(2) = a(1); 0020 b(1) = 1; 0021 b(2) = (w0-2*fs) / (w0+2*fs); 0022 0023 % END