菲涅尔衍射matlab程序

close all;clear all;warning off all;clc;
%% 基本参数设置

% M = 3072;N = 3072;
% M = 2048;N = 2048;
M = 1024;N = 1024;
% M = 512;N = 512;
%% 空域
Lx = 8*10^-6; Ly = 8*10^-6;

%% 频域
Lxx = -Lx/2:Lx/（M-1）:Lx/2;
Lyy = -Ly/2:Ly/（N-1）:Ly/2;
[mn] = meshgrid（LxxLyy）;

%% 复振幅函数参数
lamda = 365e-9;
k = 2*pi/lamda;

z = 1.0959e-05;


%% 屏幕

    %% 方孔1/2
%     Screen1 = zeros（M/4N）;
%     Screen2 = zeros（M/2N/4）;
%     Screen3 = ones（M/2N/2）;
%     Screen4 = zeros（M/2N/4）;
%     Screen5 = zeros（M/4N）;
%     Screen  = [Screen1;Screen2 Screen3 Screen4;Screen5];
    
    %% 方孔1/4
%     Screen1 = zeros（3*M/8N）;
%     Screen2 = zeros（M/43*N/8）;
%     Screen3 = ones（M/4N/4）;
%     Screen4 = zeros（M/43*N/8）;
%     Screen5 = zeros（3*M/8N）;
%     Screen  = [Screen1;Screen2 Screen3 Screen4;Screen5];

    %% 圆孔衍射
    Screen = zeros（MN）;
    [CxCy] = meshgrid（1:M1:N）;
    C =