用于获取单粒子散射特性、散射系数的Mie理论计算程序

可以完整的仿真计算出单粒子或单气泡的mie散射结果。包括散射系数、散射相函数、效率因子等等，是一个非常好用的程序代码。matlab代码，很友好，初学者可以拿来直接使用。

m=1/1.33;           %气泡相对周围介质的折射率 
% m=10^（-6）;
% x=100*pi;               %气泡半径a=20um时的气泡尺度参数
x=500*pi;                  %气泡半径a=100um时的气泡尺度参数
nsteps=361;                 %计算步长如3601表示360度/（3601-1）=0.1度

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（mxu）;

        SL（j）= real（a（1j）‘*a（1j））;      %左侧部分散射（0-pi）

        SR（j）= real（a（2j）‘*a（2j））;         %右侧部分散射（pi-2pi）

    end;

% y=[teta;SL+SR]‘;                    %计算散射角在0-pi上的散射
% y=[teta+pi;SR（nsteps:-1:1）]‘;     %计算散射角pi-2pi上的散射
y=[teta teta+pi;SL+SR SL（nsteps:-1:1）+SR（nsteps:-1:1）]‘;     % 计算散射角在0-pipi-2pi上的散射


% plot（y（:1）log（y（:2）））;                 %画图直角坐标系0-pi上的散射
polar（y（:1）log（y（:2）））            %画图极坐标系0-pipi-2pi上的散射
% % text（0.520‘\phi=0^{o}‘）;
title（sprintf（‘Mie angular scattering: m=%g+%gi x=%g‘m1m2x））;
% % 
xlabel（‘Scattering Angel/（rad）‘）;
% % ylabel（‘Scattering Intensity lg（I）‘）;
% % xlim（[0 pi]）;ylim（[-5 30]）;

% result=（log（y（:2）））‘; 

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

     文件       1130  2017-06-28 16:22  Mie_ab.asv

     文件       1126  2017-06-23 11:00  Mie_ab.m

     文件       1357  2019-04-28 09:57  Mie_I10.asv

     文件       1357  2019-04-28 09:58  Mie_I10.m

     文件        435  2009-12-12 21:16  Mie_pt.m

     文件       1050  2017-06-23 11:16  Mie_S12.asv

     文件        272  2015-05-30 16:02  para_out.m

     文件        702  2015-05-30 10:40  Mie_tetascan.m

     文件        997  2017-06-23 11:19  Mie_S12.m

     文件       1196  2019-04-28 11:03  main.m

----------- ---------  ---------- -----  ----

                 9622                    10