【光学】基于matlab夫琅禾费圆孔衍射【含Matlab源码 062期】

一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【光学】基于matlab夫琅禾费圆孔衍射【含Matlab源码 062期】

获取代码方式2：
通过订阅紫极神光博客付费专栏，凭支付凭证，私信博主，可获得此代码。

备注：
订阅紫极神光博客付费专栏，可免费获得1份代码（有效期为订阅日起，三天内有效）；

二、部分源代码

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

% 接收屏的参数处理*********************************
Screen_Xmax = 6*Lambda*L/r;              %屏X方向位置(边界)--单位为：m
Screen_Ymax = 6*Lambda*L/r;              %屏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

% 圆孔点光源的参数处理*********************************
% 该取点在x-y方向不是均匀的，中心选取的点光源数量较多，模拟高斯光
angle = 0:2*pi/90:2*pi*89/90;     %90份旋转角
r_l = linspace(r/25,r,25);        %半径分为25份
light_xy = cell(90,25);
for n = 1:90
    light_x(n,:) = r_l.*cos(angle(n));
    light_y(n,:) = r_l.*sin(angle(n));
end
for x = 1:90
    for y = 1:25
        light_xy{x,y} = [light_x(x,y),light_y(x,y)];
    end
end

% 各点光源的叠加处理*********************************
for s_x =1:101
    for s_y = 1:101
        for x = 1:90
            for y = 1:25
                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)/(90*25)^2;   %各点光源在接收屏某点的振幅分布，数值[0,1]区间
    end
end
 

  

 

  
1
  
2
  
3
  
4
  
5
  
6
  
7
  
8
  
9
  
10
  
11
  
12
  
13
  
14
  
15
  
16
  
17
  
18
  
19
  
20
  
21
  
22
  
23
  
24
  
25
  
26
  
27
  
28
  
29
  
30
  
31
  
32
  
33
  
34
  
35
  
36
  
37
  
38
  
39
  
40
  
41
  
42
  
43
  
44
  
45
  
46
  
47
  
48
  
49
  
50
  
51
  
52
  
53
  
54
  
55
  
56
 

三、运行结果

四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 门云阁.MATLAB物理计算与可视化[M].清华大学出版社，2013.

文章来源: qq912100926.blog.csdn.net，作者：海神之光，版权归原作者所有，如需转载，请联系作者。

原文链接：qq912100926.blog.csdn.net/article/details/112257199