用matlab写出米氏散射求散射光强度的代码

以下是利用Matlab实现米氏散射求散射光强度的代码：

clear all;
clc;

% 输入参数
lambda = input('请输入激光波长（nm）：');
d = input('请输入散射体直径（nm）：');
n = input('请输入散射体折射率：');
c = 2.99792458e8; % 光速，单位m/s
theta_max = 180; % 最大散射角度，单位°

% 计算
theta = linspace(0,theta_max,1000); % 散射角度，单位°
Qext = zeros(length(theta),1); % 散射截面，单位m^2
for i = 1:length(theta)
    x = 2*pi*d*n/lambda*sin(theta(i)*pi/180); % 无量纲散射参数
    Qext(i) = (2/x^2)*abs((cos(x)-1i*sin(x)/x)^2); % Mie散射截面
end

I0 = 1; % 入射光强度
Is = I0*Qext.*sind(theta); % 散射光强度，单位W/m^2/sr
Is = Is*(pi/180)^2; % 转换为W/m^2/steradian
P = pi*(d/2)^2; % 散射体投影面积，单位m^2
Is = Is*P; % 散射光强度，单位W

% 绘图
figure;
plot(theta,Is,'LineWidth',2);
set(gca,'FontSize',14);
xlabel('散射角度（°）');
ylabel('散射光强度（W）');
title(['Mie散射，\lambda=',num2str(lambda),' nm, d=',num2str(d),' nm, n=',num2str(n)]);

以上代码中，输入的参数包括激光波长、散射体直径和折射率等，通过计算得出散射角度和散射光强度，并绘制出散射光强度随散射角度变化的曲线。