Home > atmlab > mie > mie2_abs.m

mie2_abs

PURPOSE ^

Computation of the Absorption Efficiency Qabs of a sphere

SYNOPSIS ^

function result = mie2_abs(eps1, mu1, x)

DESCRIPTION ^

 Computation of the Absorption Efficiency Qabs of a sphere
 for given complex permittivity and permeability ratios,  
 eps1=eps1'+ieps1", mu1=mu1'+imu1", between inside and outside 
 of the sphere for size parameter x=k0*a. 
 Ref. Bohren and Huffman (1983) BEWI:TDD122,
 and my own notes on this topic;
 k0=2*pi./wavelength; % x=k0.*radius;
 Input: eps1, mu1, x
 C. M�zler, July 2002.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

mie2_abs.m

SOURCE CODE ^

0001 function result = mie2_abs(eps1, mu1, x)
0002 
0003 % Computation of the Absorption Efficiency Qabs of a sphere
0004 % for given complex permittivity and permeability ratios,
0005 % eps1=eps1'+ieps1", mu1=mu1'+imu1", between inside and outside
0006 % of the sphere for size parameter x=k0*a.
0007 % Ref. Bohren and Huffman (1983) BEWI:TDD122,
0008 % and my own notes on this topic;
0009 % k0=2*pi./wavelength; % x=k0.*radius;
0010 % Input: eps1, mu1, x
0011 % C. M�zler, July 2002.
0012 
0013 ep2=imag(eps1); mu2=imag(mu1);
0014 nj=100*round(2+x+4*x.^(1/3));
0015 dx=x/nj;
0016 x2=x.*x;
0017 nj1=nj+1;
0018 xj=(0:dx:x);
0019 en=mie2_Esquare(eps1,mu1,x,nj); % E-field
0020 en1=0.5*en(nj1).*x2;     % End-Term correction in integral
0021 enx=en*(xj.*xj)'-en1;    % Trapezoidal radial integration
0022 inte=dx.*enx;
0023 Qabse=4.*ep2.*inte./x2
0024 hn=mie2_Esquare(mu1,eps1,x,nj); % H-field (by duality)
0025 hn1=0.5*hn(nj1).*x2;     % End-Term correction in integral
0026 hnx=hn*(xj.*xj)'-hn1;    % Trapezoidal radial integration
0027 inth=dx.*hnx;
0028 Qabsm=4.*mu2.*inth./x2
0029 Qabs=Qabse+Qabsm;
0030 result=Qabs;

Generated on Mon 15-Sep-2014 13:31:28 by m2html © 2005