Home > atmlab > mie > mie2_tetascan.m

mie2_tetascan

PURPOSE ^

Computation and plot of Mie Power Scattering function for given

SYNOPSIS ^

function result = mie2_tetascan(eps1, mu1, x, nsteps)

DESCRIPTION ^

 Computation and plot of Mie Power Scattering function for given 
 complex permittivity and permeability ratios eps1=eps1'+ieps1",
 mu1=mu1'+imu1" between inside and outside of the sphere
 and range of size parameters x=k0*a, 
 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 ^

mie2_tetascan.m

SOURCE CODE ^

0001 function result = mie2_tetascan(eps1, mu1, x, nsteps)
0002 
0003 % Computation and plot of Mie Power Scattering function for given
0004 % complex permittivity and permeability ratios eps1=eps1'+ieps1",
0005 % mu1=mu1'+imu1" between inside and outside of the sphere
0006 % and range of size parameters x=k0*a,
0007 % according to Bohren and Huffman (1983) BEWI:TDD122
0008 % C. M�zler, July 2002.
0009 
0010 nsteps=nsteps;
0011 epsp=real(eps1);epspp=imag(eps1);
0012 mup=real(mu1);mupp=imag(mu1);
0013 nx=(1:nsteps); dteta=pi/(nsteps-1);
0014 teta=(nx-1).*dteta;
0015     for j = 1:nsteps, 
0016         u=cos(teta(j));
0017         a(:,j)=mie2_S12(eps1,mu1,x,u);
0018         SL(j)= real(a(1,j)'*a(1,j))/(pi*x^2);
0019         SR(j)= real(a(2,j)'*a(2,j))/(pi*x^2);
0020     end;
0021 y=[teta teta+pi;SL SR(nsteps:-1:1)]'; 
0022 
0023 polar(y(:,1),y(:,2))
0024 title(sprintf('Angular pattern for x=%g, eps1=%g+%gi, mu1=%g+%gi',x,epsp,epspp,mup,mupp));
0025 xlabel('Scattering Angle')
0026 result=y;

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