matlab生成菲涅尔透镜

用matlab生成菲涅尔透镜图，焦距可调

clc
ren=0.532;%波长，单位为微米
f=800000;  %焦距
a=1024;b=1280;  
Z=zeros（ab）;
for i=1:a
    for j=1:b
    Z（ij）=（20）^2*（（i-a/2）.^2+（j-b/2）.^2）/（2*ren*f）;  %20为像素大小，可变
    Z（ij）=mod（Z（ij）1）;
    end
end
Z=uint8（Z.*255）;
figure
imshow（Z）;
[fnamepnameindex]=uiputfile（{‘*.bmp‘}‘图片另存为‘）;
if index
   imwrite（Zfullfile（pnamefname））
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         381  2014-01-15 11:59  toujing.m
     文件        5712  2014-01-13 20:51  toujing.prj