clear;clc;close?all;
p=512;h=8;%p是圖像大小，h是孔徑寬度
[XY]=meshgrid（0:1:p-10:1:p-1）;
a=zeros（[pp]）;
for?i=1:p
????for?j=1:p
????????if?（i-p/2）^2+（j-p/2）^2<