Convert comlpex pole/zero into frequency/Q pole/zero representation. [f0, q]= pri2fq(c) M Hewitson 2003-03-12 $Id: pri2fq.m,v 1.1 2007/04/03 06:50:38 hewitson Exp $
0001 function [f0, q]= pri2fq(c) 0002 0003 % Convert comlpex pole/zero into frequency/Q pole/zero representation. 0004 % 0005 % [f0, q]= pri2fq(c) 0006 % 0007 % M Hewitson 2003-03-12 0008 % 0009 % 0010 % $Id: pri2fq.m,v 1.1 2007/04/03 06:50:38 hewitson Exp $ 0011 % 0012 0013 0014 if(nargin==0) 0015 disp('usage: [f0, q]= ri2fq(c)'); 0016 return 0017 end 0018 0019 a = real(c); 0020 b = imag(c); 0021 0022 f0 = -a*sqrt(b^2/a^2 + 1) / 2 / pi; 0023 q = sqrt(b^2/a^2 +1)/2; 0024 0025 % END