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