mie计算代码

mie散射的计算 matlab源代码。计算各种介质颗粒的散射！

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];

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1064  2002-11-01 12:21  Mie-Matlab-Maetzler-v2\Code\besselplot1.m

     文件       1092  2002-11-01 12:21  Mie-Matlab-Maetzler-v2\Code\besselplot2.m

     文件        929  2002-11-01 12:21  Mie-Matlab-Maetzler-v2\Code\besselplot3.m

     文件       1319  2002-11-01 12:18  Mie-Matlab-Maetzler-v2\Code\Besselplot4.m

     文件        754  2002-11-01 12:21  Mie-Matlab-Maetzler-v2\Code\besselplot5.m

     文件        569  2002-11-01 12:21  Mie-Matlab-Maetzler-v2\Code\besselriccati.m

     文件        512  2002-11-01 12:21  Mie-Matlab-Maetzler-v2\Code\epsice.M

     文件        591  2002-11-01 12:21  Mie-Matlab-Maetzler-v2\Code\epswater.m

     文件        430  2002-11-01 12:21  Mie-Matlab-Maetzler-v2\Code\epswater93.m

     文件       1754  2002-11-01 12:18  Mie-Matlab-Maetzler-v2\Code\Mie.m

     文件       1763  2002-11-01 12:18  Mie-Matlab-Maetzler-v2\Code\Mie2.m

     文件       1098  2002-11-01 12:18  Mie-Matlab-Maetzler-v2\Code\Mie2_ab.m

     文件       1037  2002-11-01 12:18  Mie-Matlab-Maetzler-v2\Code\Mie2_abs.m

     文件       1187  2002-11-01 12:18  Mie-Matlab-Maetzler-v2\Code\Mie2_cd.m

     文件       1574  2002-11-01 12:18  Mie-Matlab-Maetzler-v2\Code\Mie2_esquare.m

     文件        659  2002-11-01 12:18  Mie-Matlab-Maetzler-v2\Code\Mie2_S12.m

     文件        911  2002-11-01 12:19  Mie-Matlab-Maetzler-v2\Code\Mie2_tetascan.m

     文件       1121  2002-11-01 12:19  Mie-Matlab-Maetzler-v2\Code\Mie2_xscan.m

     文件        437  2002-11-01 12:20  Mie-Matlab-Maetzler-v2\Code\Mieab_1.m

     文件        537  2002-11-01 12:20  Mie-Matlab-Maetzler-v2\Code\Mieab_2.m

     文件       2159  2002-11-01 12:20  Mie-Matlab-Maetzler-v2\Code\Miecoated.m

     文件       2166  2002-11-01 12:20  Mie-Matlab-Maetzler-v2\Code\Miecoated_ab1.m

     文件       1473  2002-11-01 12:22  Mie-Matlab-Maetzler-v2\Code\Miecoated_ab2.m

     文件       1474  2002-11-01 12:21  Mie-Matlab-Maetzler-v2\Code\Miecoated_ab3.m

     文件       1265  2002-11-01 12:21  Mie-Matlab-Maetzler-v2\Code\Miecoated_prscan.m

     文件       2240  2002-11-01 12:22  Mie-Matlab-Maetzler-v2\Code\Miecoated_rain3.m

     文件       2009  2002-11-01 12:22  Mie-Matlab-Maetzler-v2\Code\Miecoated_rain4.m

     文件       2023  2002-11-01 12:22  Mie-Matlab-Maetzler-v2\Code\Miecoated_rain5.m

     文件       2179  2002-11-01 12:23  Mie-Matlab-Maetzler-v2\Code\Miecoated_rain6.m

     文件       2182  2002-11-01 12:23  Mie-Matlab-Maetzler-v2\Code\Miecoated_rain7.m

............此处省略28个文件信息