瑞利索末菲衍射计算伯松亮斑

% 本程序计算圆形高斯光束入射到位于 束腰处的 圆盘产生的衍射及伯松亮斑
close all;clear all;
w = 0.25;                    % 入射光束腰，束腰半径为0.25mm
lamda = 0.632e-3;           % 入射激光波长
k = 2*pi/lamda;             % 波数
q = 1i*pi*w^2/lamda;        % 束腰处的q参数
N = 500;                    % 光强采样精度
x = linspace（-0.50.5N）;   % 观察光斑的坐标范围
z = 10;                     % 衍射传播的距离
E0 = exp（-1i*k*x.^2/2/q）;   % 束腰处的光场复振幅分布

[X Y] = meshgrid（xx）;      % 坐标网格化
[Ex Ey] =meshgrid（E0E0）;   % 光场二维化
E1 = Ex.*Ey;                % 束腰处的二维光场
figure;imagesc（abs（E1）.^2）; % 圆形高斯光束
colormap hot;axis image;

Disk = one