MIE散射理论计算

MIE散射理论计算

function result = besselplot1（n m xmin dx nx）



% Computation and plot of Riccati-Bessel Functions of Order n

% for complex argument z=m*x used in Mie Theory. 

% input: order n refractive index m minimum x value xmin 

% x interval dx number of x values nx.

% C. M鋞zler August 2002



m1=real（m）; m2=imag（m）;

nn=（1:nx）;

x=xmin+dx*nn;

nu=n+0.5;

z=m.*x;

sqx= sqrt（0.5*pi*z）; 

psx = besselj（nu z）.*sqx;

chx = -bessely（nu z）.*sqx;

dpic=psx-i*chx;

r=[real（psx）;-imag（psx）;real（chx）;imag（chx）;real（dpic）;imag（dpic）];

if abs（m2）<0.1

    plot（xr（1:4:））

    legend（‘real（psi_n（mx））‘‘-imag（psi_n（mx））‘‘real（chi_n（mx））‘‘imag（chi_n（mx））‘）

    title（sprintf（‘Riccati-Bessel Functions of Order n=%g for m=%g+%gi‘nm1m2））

    xlabel（‘x‘）

else

    semilogy（xr（1:6:））

    legend（‘real（psi_n（mx））‘‘-imag（psi_n（mx））‘‘real（chi_n（mx））‘‘imag（chi_n（mx））‘‘real（dpic（mx））‘‘imag（dpic（mx））‘）

    title（sprintf（‘Riccati-Bessel Functions of Order n=%g for m=%g+%gi‘nm1m2））

    xlabel（‘x‘）

end;

result=[x;r];

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         512  2002-11-01 12:21  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\epsice.M
     文件         430  2002-11-01 12:21  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\epswater93.m
     文件         591  2002-11-01 12:21  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\epswater.m
     文件        1763  2002-11-01 12:18  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie2.m
     文件        1098  2002-11-01 12:18  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie2_ab.m
     文件        1037  2002-11-01 12:18  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie2_abs.m
     文件        1187  2002-11-01 12:18  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie2_cd.m
     文件        1574  2002-11-01 12:18  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie2_esquare.m
     文件         659  2002-11-01 12:18  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie2_S12.m
     文件         911  2002-11-01 12:19  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie2_tetascan.m
     文件        1121  2002-11-01 12:19  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie2_xscan.m
     文件        1754  2002-11-01 12:18  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie.m
     文件         975  2002-11-01 12:19  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_1.m
     文件        1554  2002-11-01 12:19  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_2.m
     文件        1047  2002-11-01 12:19  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_ab.m
     文件         962  2002-11-01 12:19  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_abcd.m
     文件         742  2002-11-01 12:19  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_abs.m
     文件         953  2002-11-01 12:19  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_cd.m
     文件        1364  2002-11-01 12:19  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_esquare.m
     文件         437  2002-11-01 12:20  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_pt.m
     文件         877  2002-11-01 12:20  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_rain1.m
     文件        1002  2002-11-01 12:20  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_rain1c.m
     文件         958  2002-11-01 12:20  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_rain1d.m
     文件         881  2002-11-01 12:20  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_rain1r.m
     文件        1398  2002-11-01 12:20  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_rain2.m
     文件        2011  2002-11-01 12:20  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_rain3.m
     文件        1628  2002-11-01 12:20  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_rain4.m
     文件        1463  2002-11-01 12:20  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_rain5.m
     文件         603  2002-11-01 12:19  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_S12.m
     文件         726  2002-11-01 12:20  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_tetascan.m
     文件         929  2002-11-01 12:20  PFD-s-etc-Alt2\Scattering-programs-2002\Mie-Matlab-Maetzler-v2\Mie_xscan.m
............此处省略22个文件信息