LTPDA_FQ2RI2 Convert frequency/Q pole/zero representation into real and imaginary conjugate pairs. Returns [ri conj(ri)] ri= ltpda_fq2ri2(f0, Q) M Hewitson 26-01-07 $Id: ltpda_fq2ri2.html,v 1.1 2007/06/08 14:15:10 hewitson Exp $
0001 function ri= ltpda_fq2ri2(f0, Q) 0002 0003 % LTPDA_FQ2RI2 Convert frequency/Q pole/zero representation into real 0004 % and imaginary conjugate pairs. Returns [ri conj(ri)] 0005 % 0006 % ri= ltpda_fq2ri2(f0, Q) 0007 % 0008 % M Hewitson 26-01-07 0009 % 0010 % 0011 % $Id: ltpda_fq2ri2.html,v 1.1 2007/06/08 14:15:10 hewitson Exp $ 0012 % 0013 0014 if(nargin==0) 0015 disp('usage: ri = ltpda_fq2ri2(f0, Q)') 0016 return 0017 end 0018 0019 % a = Q - 0.5 ; 0020 % re = 2*pi*f0 / sqrt(1+tan(a)^2); 0021 % im = re * tan(a); 0022 0023 w0 = 2*pi*f0; 0024 re = -w0/(2*Q); 0025 im = w0*sqrt(4*Q*Q-1)/(2*Q); 0026 G = 1/(w0*w0); 0027 tri = complex(re,im); 0028 ri = [tri conj(tri)]'; 0029 0030 % END