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?=