function?main（）
????clear;
????clc;
?
????p1?=?double（imread（‘G:\quanxi1.bmp‘））;
????[NM]?=?size（p1）;
????z=1;
????w=0.6328e-6;
????k=2*pi/w;
????v=3.2e-6;

????[f]?=?fresnel（p1?M?Nv??z?w）;???????????????????
????imshow（abs（f）[]）;
????title（‘全息再現(xiàn)‘）;
end


function?[f1]=?fresnel（f0?M?Nv1?z?lambda）
????dx0=v1;
????dy0=v1;
????k=2*pi/lambd