题目分析利用菲涅尔—基尔霍夫衍射公式对问题进行处理，通过去不同的 z 值观察近场,菲涅尔区，夫琅和费区的衍射条纹。菲涅尔—基尔霍夫衍射公式：由于单色平面波垂直照射，，由于一般考虑光强的相对值,所以分析时取A=1，,并对方孔进行微元，将方孔每个点在光屏上衍射的复振幅进行叠加，即计算，从而得到光屏上的衍射条纹.对于方孔编程时采纳循环的方式，对于光屏则利用矩阵对光屏整体的复振幅,光强等进行分析，加快运算速度。实现代码clear；lambda=6。32e—7;%波长k=2＊pi/lambda;z=550；％光屏位置a=0。01；b=0.01；%方孔尺寸x1=linspace（-0.005，0。005，199）；y1=linspace（-0。005，0.005,199）;％将方孔分为许多小的微元x=linspace(-0.2,0.2，299);y=linspace（-0.2，0.2,299）;%将光屏分为许多小的微元［X，Y]=meshgrid(x，y）;％用二维矩阵描述光屏上点的位置E=zeros（299,299);%用于存放光屏上每个点的复振幅r=zeros(299，299）;%用于存放光屏上没一点对应方孔上点的距离for l=1:199for j=1：199 %计算方孔上的每个点在光屏上的复振幅r=sqrt（（X-x1（l））.^2+(Y—y1（j)）。^2+z*z)； %当前方孔上点到光屏上各个点的距离E=E+（—i/lambda)＊（exp（i＊k*r））.＊(z。/（2＊r.＊r)+1./(2＊r））＊z; %对不同点在光屏上的复振幅进行累加endendE1=conj（E);I=E1。*E; ％求光强% pcolor(y,x,I） %作图% colorbarIx= I（1：299，200）;％光屏中心水平方向光强分布plot（x,Ix)仿真图形夫琅和费区：Z=550m光强分布光屏中心水平方向光强分布菲涅尔区：Z=1m近场区：Z=0。001m
