菲涅尔衍射程序，能够获得所需图像的菲涅尔衍射图样

z0=5;%范围2——10 x0=0;%范围-8——8 ROWS=256; COLS=256; colormap(gray(255)) sigma=1/z0; y0=-x0; y=-12.8; for r=1:COLS, x=-12.8; for c=1:ROWS, fFZP(r,c)=exp(j*sigma*(x-x0)*(x-x0)+j*sigma*(y-y0)*(y-y0)); x=x+.1; end y=y+.1; end max1=max(fFZP); max2=max(max1); scale=1.0/max2; fFZP=fFZP.*scale; R=127*(1+imag(fFZP)); figure(1) image(R); axis square on axis off