多层膜光栅的matlab模拟仿真代码

多层膜光栅的matlab模拟仿真代码如下：

% 多层膜光栅的matlab模拟仿真代码
% 定义输入参数
lambda = 632.8e-9; % 波长
theta = 0; % 入射角度
d = 500e-9; % 光栅周期
n1 = 1; % 空气
n2 = 1.5; % 光栅材料
n3 = 1; % 衬底
h1 = lambda/(4*n2); % 线宽
h2 = (d - h1)/2; % 线间距
N = 50; % 层数

% 计算反射率和透射率
R = zeros(1, N); % 反射率
T = zeros(1, N); % 透射率
for i = 1:N
    % 计算每层的折射率
    if mod(i, 2) == 1 % 奇数层
        n = n2;
        h = h1;
    else % 偶数层
        n = n1;
        h = h2;
    end
    % 计算反射率和透射率
    k = 2 * pi * n * h / lambda;
    R(i) = ((n1*cos(theta) - n*cos(asin(sin(theta)/n))) / (n1*cos(theta) + n*cos(asin(sin(theta)/n))))^2 + ((n*cos(theta) - n1*cos(asin(sin(theta)/n))) / (n*cos(theta) + n1*cos(asin(sin(theta)/n))))^2;
    T(i) = 1 - R(i);
end

% 计算透射光强
I = zeros(1, N+1); % 透射光强
I(1) = 1; % 入射光强
for i = 1:N
    I(i+1) = I(i) * T(i);
end

% 画图显示结果
x = linspace(0, d, N+1);
figure
plot(x*1e9, I);
xlabel('位置(nm)');
ylabel('透射光强');

在上述代码中，我们先定义了多层膜光栅的输入参数，包括波长、入射角度、光栅周期、折射率、线宽、线间距和层数等。然后，我们利用计算公式计算每层的反射率和透射率，并利用透射率计算出透射光强。最后，我们将透射光强作为纵坐标，光栅位置作为横坐标，画出透射光强随位置变化的曲线。