Computation and plot of Mie Efficiencies for given complex refractive-index ratio m=m'+im" 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: m', m", x, efficiencies for extinction (qext), scattering (qsca), absorption (qabs), backscattering (qb), qratio=qb/qsca and asymmetry parameter (asy=<costeta>). C. M�zler, May 2002.
0001 function result = mie_xscansmooth(m, nsteps, dx, Dx) 0002 0003 % Computation and plot of Mie Efficiencies for given 0004 % complex refractive-index ratio m=m'+im" 0005 % and range of size parameters x=k0*a, 0006 % starting at x=0 with nsteps increments of dx 0007 % a=sphere radius, using complex Mie coefficients an and bn 0008 % according to Bohren and Huffman (1983) BEWI:TDD122 0009 % result: m', m", x, efficiencies for extinction (qext), 0010 % scattering (qsca), absorption (qabs), backscattering (qb), 0011 % qratio=qb/qsca and asymmetry parameter (asy=<costeta>). 0012 % C. M�zler, May 2002. 0013 nsmooth=Dx/dx; 0014 nx=(1:nsteps)'; 0015 x=(nx-1)*dx; 0016 for j = 1:nsteps 0017 a(j,:)=mie(m,x(j)); 0018 end; 0019 output_parameters='Qext, Qsca, Qabs, Qb, <costeta>, Qb/Qsca' 0020 for j=1:6, 0021 b(:,j)=smooth(a(:,j),nsmooth); 0022 end; 0023 b(:,3)=10000*b(:,3); 0024 % plotting the smoothed results 0025 m1=real(m);m2=imag(m); 0026 plot(x,b(:,1:6)) 0027 legend('Qext','Qsca','10^4*Qabs','Qb','<costeta>','Qb/Qsca') 0028 title(sprintf('Smoothed (Dx=%g) Mie Efficiencies, m=%g+%gi',Dx,m1,m2)) 0029 xlabel('x') 0030 0031 result=b;