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

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件?????????381??2014-01-15?11:59??toujing.m
?????文件????????5712??2014-01-13?20:51??toujing.prj