Lohman3全息图

PI=3.14159; 

%定义一物体为N*N
K=imread（‘E128.bmp‘）; 
N=128; 
for I=1:1:N 
    for J=1:1:N 
        A（IJ）=1-K（IJ）; 
        B（IJ）=0; 
    end;end; 
figure; 
imshow（A）; 

%为降低全息图的动态范围，乘一随机位相因子% 
for I=1:1:N 
    for J=1:1:N 
        R=rand（11）; 
        B（IJ）=A（IJ）*sin（R*2*PI）; 
        A（IJ）=A（IJ）*cos（R*2*PI）; 
        F（IJ）=A（IJ）+j*B（IJ）; 
    end;end; 

%对物函数做二维FFT变换% 
       F=fft2（F）; 
       A=real（F）; 
       B=imag（F）; 

%计算模与幅角% 
for I=1:1:N 
    for J=1:1:N 
        if A（IJ）~=0      
            A1=atan（B（IJ）/A（IJ））; 
            if A（IJ）>0   A1=A1;end; 
            if A（IJ）<0    
                if B（IJ）~=0   A1=sign（B（IJ））*PI+A1;end; 
                if B（IJ）==0   A1=A1;                end; 
            end;end; 
        if A（I