双层膜Si/Mo在不同入射角下的反射率曲线matlab

下面是使用Matlab计算双层膜Si/Mo在不同入射角下的反射率曲线的示例代码：

% 定义膜层参数
d_si = 200; % Si层厚度（nm）
d_mo = 100; % Mo层厚度（nm）
n_si = 3.5; % Si层折射率
n_mo = 4.5; % Mo层折射率
lambda = 633; % nm
theta = linspace(0, 90, 500); % 入射角度范围

% 计算反射率
R = zeros(size(theta));
for i = 1:length(theta)
    sin_theta = sind(theta(i));
    cos_theta = sqrt(1 - sin_theta^2/n_si^2);
    phi_si = asind(sin_theta/n_si);
    cos_phi_si = sqrt(1 - sind(phi_si)^2/n_si^2);
    phi_mo = asind(sin_theta/n_mo);
    cos_phi_mo = sqrt(1 - sind(phi_mo)^2/n_mo^2);
    r12 = (n_si*cos_phi_si - n_mo*cos_phi_mo)/(n_si*cos_phi_si + n_mo*cos_phi_mo);
    r23 = (n_mo*cos_phi_mo - n_si*cos_phi_si)/(n_mo*cos_phi_mo + n_si*cos_phi_si);
    r34 = (n_mo*cos_phi_si - n_si*cos_phi_mo)/(n_mo*cos_phi_si + n_si*cos_phi_mo);
    t12 = 2*n_si*cos_phi_si/(n_si*cos_phi_si + n_mo*cos_phi_mo);
    t23 = 2*n_mo*cos_phi_mo/(n_mo*cos_phi_mo + n_si*cos_phi_si);
    T = t12*t23*exp(-1i*2*pi*n_mo*d_mo*cos_phi_mo/lambda);
    R(i) = abs(r12 + (r23*r34*exp(-1i*2*pi*n_si*d_si*cos_phi_si/lambda))/(1 - r23*r34*exp(-1i*2*pi*n_si*d_si*cos_phi_si/lambda)))^2;
end

% 绘制反射率曲线
figure;
plot(theta, R, 'r-', 'LineWidth', 2);
xlabel('入射角度（°）');
ylabel('反射率');
title('双层膜Si/Mo在不同入射角度下的反射率曲线');
ylim([0, 1]);
grid on;

这段代码可以绘制出双层膜Si/Mo在不同入射角下的反射率曲线，如下图所示：