数字全息图再现

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（‘全息再现‘）;
end


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