mie_tetascan.m

来自「MieMatlabMaetzler.zip 非常全功能的Mie散射程序 ma」· M 代码 · 共 24 行

function result = Mie_tetascan(m, x, nsteps)

% Computation and plot of Mie Power Scattering function for given 
% complex refractive-index ratio m=m'+im", size parameters x=k0*a, 
% according to Bohren and Huffman (1983) BEWI:TDD122
% C. M鋞zler, May 2002.

nsteps=nsteps;
m1=real(m); m2=imag(m);
nx=(1:nsteps); dteta=pi/(nsteps-1);
teta=(nx-1).*dteta;
    for j = 1:nsteps, 
        u=cos(teta(j));
        a(:,j)=Mie_S12(m,x,u);
        SL(j)= real(a(1,j)'*a(1,j));
        SR(j)= real(a(2,j)'*a(2,j));
    end;
y=[teta teta+pi;SL SR(nsteps:-1:1)]'; 

polar(y(:,1),y(:,2))
title(sprintf('Mie angular scattering: m=%g+%gi, x=%g',m1,m2,x));
xlabel('Scattering Angle')
result=y;