Home > atmlab > mie > miecoated_tetascan.m

miecoated_tetascan

PURPOSE ^

Computation and plot of Mie Power Scattering function for given

SYNOPSIS ^

function result = miecoated_tetascan(m1, m2, x, y, nsteps)

DESCRIPTION ^

 Computation and plot of Mie Power Scattering function for given 
 complex refractive-index ratios m1,2=m1,2'+im1,2", size parameters 
 x=k0*a, y=k0*b, of coated sphere
 according to Bohren and Huffman (1983) BEWI:TDD122
 C. M�zler, July 2002.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

miecoated_tetascan.m

SOURCE CODE ^

0001 function result = miecoated_tetascan(m1, m2, x, y, nsteps)
0002 
0003 % Computation and plot of Mie Power Scattering function for given
0004 % complex refractive-index ratios m1,2=m1,2'+im1,2", size parameters
0005 % x=k0*a, y=k0*b, of coated sphere
0006 % according to Bohren and Huffman (1983) BEWI:TDD122
0007 % C. M�zler, July 2002.
0008 
0009 nsteps=nsteps;
0010 m1p=real(m1); m1pp=imag(m1);
0011 m2p=real(m2); m2pp=imag(m2);
0012 nx=(1:nsteps); dteta=pi/(nsteps-1);
0013 teta=(nx-1).*dteta;
0014     for j = 1:nsteps, 
0015         u=cos(teta(j));
0016         a(:,j)=miecoated_S12(m1,m2,x,y,u);
0017         SL(j)= real(a(1,j)'*a(1,j))/(pi*y^2);
0018         SR(j)= real(a(2,j)'*a(2,j))/(pi*y^2);
0019     end;
0020 scan=[teta teta+pi;SL SR(nsteps:-1:1)]'; 
0021 
0022 polar(scan(:,1),scan(:,2)),
0023 title(sprintf('Angular scattering of coated sphere: m1=%g+%gi, m2=%g+%gi, x=%g, y=%g',m1p,m1pp,m2p,m2pp,x,y)),
0024 xlabel('Scattering Angle');
0025 result=y;

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