Home > atmlab > mie > mie2_xscan.m

mie2_xscan

PURPOSE ^

Computation and plot of Mie Efficiencies for given complex

SYNOPSIS ^

function result = mie2_xscan(eps1, mu1, nsteps, dx)

DESCRIPTION ^

 Computation and plot of Mie Efficiencies 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, 
 starting at x=0 with nsteps increments of dx
 a=sphere radius, using complex Mie coefficients an and bn 
 according to Bohren and Huffman (1983) BEWI:TDD122
 result: x, efficiencies for extinction (qext), 
 scattering (qsca), absorption (qabs), backscattering (qb), 
 qratio=qb/qsca and asymmetry parameter (asy=<costeta>).
 C. M�zler, May 2002.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

mie2_xscan.m

SOURCE CODE ^

0001 function result = mie2_xscan(eps1, mu1, nsteps, dx)
0002 
0003 % Computation and plot of Mie Efficiencies for given complex
0004 % 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 % starting at x=0 with nsteps increments of dx
0008 % a=sphere radius, using complex Mie coefficients an and bn
0009 % according to Bohren and Huffman (1983) BEWI:TDD122
0010 % result: x, efficiencies for extinction (qext),
0011 % scattering (qsca), absorption (qabs), backscattering (qb),
0012 % qratio=qb/qsca and asymmetry parameter (asy=<costeta>).
0013 % C. M�zler, May 2002.
0014 
0015 nx=(1:nsteps)';
0016 x=(nx-1)*dx;
0017 for j = 1:nsteps
0018     a(j,:)=mie2(eps1,mu1,x(j));
0019 end;
0020 output_parameters='Real(m), Imag(m), x, Qext, Qsca, Qabs, Qb, <costeta>, Qb/Qsca'
0021 
0022 % plotting the results
0023 epsp=real(eps1);epspp=imag(eps1);
0024 mup=real(mu1);mupp=imag(mu1);
0025 plot(x,a(:,1:6))
0026 legend('Qext','Qsca','Qabs','Qb','<costeta>','Qb/Qsca')
0027 title(sprintf('Mie Efficiencies for eps1=%g+%gi, mu1=%g+%gi',epsp,epspp,mup,mupp))
0028 xlabel('x')
0029 
0030 result=[x a];

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