Home > atmlab > mie > mie_phasefunctions.m

mie_phasefunctions

PURPOSE ^

Computation of Mie phase functions p1 and p2

SYNOPSIS ^

function result = mie_phasefunctions(m, x, u)

DESCRIPTION ^

 Computation of Mie phase functions p1 and p2
 for complex refractive index m=m'+im", 
 size parameter x=k0*a, and u=cos(scattering angle),
 where k0=vacuum wave number, a=sphere radius;
 s. p. 111-114, Bohren and Huffman (1983) BEWI:TDD122
 C. M�zler, May 2002

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

mie_phasefunctions.m

SOURCE CODE ^

0001 function result = mie_phasefunctions(m, x, u)
0002 
0003 % Computation of Mie phase functions p1 and p2
0004 % for complex refractive index m=m'+im",
0005 % size parameter x=k0*a, and u=cos(scattering angle),
0006 % where k0=vacuum wave number, a=sphere radius;
0007 % s. p. 111-114, Bohren and Huffman (1983) BEWI:TDD122
0008 % C. M�zler, May 2002
0009 
0010 nmax=round(2+x+4*x^(1/3));
0011 ab=mie_ab(m,x);
0012 an=ab(1,:);
0013 bn=ab(2,:);
0014 
0015 pt=mie_pt(u,nmax);
0016 pin =pt(1,:);
0017 tin =pt(2,:);
0018 
0019 n=(1:nmax);
0020 n2=(2*n+1)./(n.*(n+1));
0021 pin=n2.*pin;
0022 tin=n2.*tin;
0023 S1=(an*pin'+bn*tin');
0024 S2=(an*tin'+bn*pin');
0025 Q=mie(m,x);
0026 Qext=Q(1);
0027 p1=S1'*S1./(pi*Qext.*x^2);
0028 p2=S2'*S2./(pi*Qext.*x^2);
0029 result=[p1;p2];

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