Matlab：实现Fra矩形孔衍射仿真

本文基于光学原理，用Matlab实现夫琅禾费矩形孔衍射现象仿真。

---

Project Code

• 算法核心机理：参看基于光的叠加原理——提出近似分析衍射的一种方法

% 作者：ZQJ
% 日期：2021.1.8 星期五

%***********************模拟矩形孔衍射图样**************************
% 方法：在矩形孔中，均匀取点光源模拟；  该仿真可以理解为：基于夫琅禾费衍射（L的大小不影响图像）
clear,clc;
% 基本参数输入及处理*********************************
% 本例输入：Lambda = 632.8;  a = 0.1;  b = 0.1;  L = 1;
Lambda = input('输入单色光的波长（单位为nm）：');
Lambda = Lambda * 1e-9;       %波长单位转为：m
a = input('输入矩形孔的宽度（单位为mm）：');
a = a * 1e-3;                 %矩形孔宽度单位转为：m
b = input('输入矩形孔的长度（单位为mm）：');
b = b * 1e-3;                 %矩形孔长度单位转为：m
L = input('输入矩形孔到屏的距离（单位为m）：');
r = zeros(51,51);
A = zeros(101,101);             

% 接收屏的参数处理*********************************
Screen_Xmax = 6*Lambda*L/a;              %屏X方向位置(边界)--单位为：m
Screen_Ymax = 6*Lambda*L/b;              %屏Y方向位置(边界)--单位为：m
Screen_x = linspace(-Screen_Xmax,Screen_Xmax,101);      %x取样101份
Screen_y = linspace(-Screen_Ymax,Screen_Ymax,101);      %y取样101份
Screen_xy = cell(101,101);
for x = 1:101
    for y = 1:101
        Screen_xy{x,y} = [Screen_x(x),Screen_y(y)];
    end
end

% 矩形孔点光源的参数处理*********************************
light_x = linspace(-a/2,a/2,51);  %矩形孔宽度上取51份点光源
light_y = linspace(-b/2,b/2,51);  %矩形孔宽度上取51份点光源
light_xy = cell(51,51);
for x = 1:51
    for y = 1:51
        light_xy{x,y} = [light_x(x),light_y(y)];
    end
end

% 各点光源的叠加处理*********************************
for s_x =1:101
    for s_y = 1:101
        for x = 1:51
            for y = 1:51
                r(x,y) = sqrt(sum((light_xy{x,y} - Screen_xy{s_x,s_y}).^2)+L^2);  %孔上的点光源到接收屏某点距离
            end
        end
        phi = 2*pi.*(r-L)./Lambda;                      %各点光源在接收屏某点的相位差
        sumcos = sum(sum(cos(phi)));
        sumsin = sum(sum(sin(phi)));
        A(s_x,s_y) = (sumsin^2 + sumcos^2)/(51*51)^2;   %各点光源在接收屏某点的振幅分布，数值[0,1]区间
    end
end

% 矩形孔衍射光强接收屏分布图*************************
I = A .*2000;                        %振幅增大（便于低的色彩强度值显示）
figure;
image(Screen_x,Screen_y,I);          %设置x和y的像素，显示数值
colormap hot;                        %妆点色彩
xlabel('空间坐标x'),ylabel('空间坐标y');
title('矩形孔衍射光强接收屏分布');
colorbar;

仿真结果图：

---

专栏内容供作者本人或大家学习使用，多多指教 ~