RI2FQ Convert comlpex pole/zero into frequency/Q pole/zero representation. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% DESCRIPTION: RI2FQ Convert comlpex pole/zero into frequency/Q pole/zero representation. CALL: [f0, q]= ri2fq(c) VERSION: $Id: ri2fq.m,v 1.2 2008/07/11 17:34:47 ingo Exp $ HISTORY: 03-12-2003 M Hewitson Creation %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0001 % RI2FQ Convert comlpex pole/zero into frequency/Q pole/zero representation. 0002 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0003 % 0004 % DESCRIPTION: RI2FQ Convert comlpex pole/zero into frequency/Q pole/zero 0005 % representation. 0006 % 0007 % CALL: [f0, q]= ri2fq(c) 0008 % 0009 % VERSION: $Id: ri2fq.m,v 1.2 2008/07/11 17:34:47 ingo Exp $ 0010 % 0011 % HISTORY: 03-12-2003 M Hewitson 0012 % Creation 0013 % 0014 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0015 0016 function [f0, q]= ri2fq(c) 0017 0018 if(nargin==0) 0019 disp('usage: [f0, q]= ri2fq(c)'); 0020 return 0021 end 0022 0023 a = real(c(1)); 0024 b = imag(c(1)); 0025 0026 % f0 = -a*sqrt(b^2/a^2 + 1) / 2 / pi; 0027 % q = sqrt(b^2/a^2 +1)/2; 0028 0029 k = -a*sqrt(b^2/a^2 + 1) * 2 * pi; 0030 f0 = 1/k; 0031 q = sqrt(b^2/a^2 +1)/2; 0032 0033 % END