matlab修正离轴制作全息图与再现

function [Image] = holo2（）
%UNtitleD Summary of this function goes here
%   Detailed explanation goes here
%Modified off-axis reference beam CGH of Burch% 
%读入原始图象 
PI=3.14159; 
Image=imread（‘D:\Matlab\homework\lena2.jpg‘）; 
Image_b=Image（::3）;
%zero=zeros（256256）;
%Image_B=cat（3zerozeroImage_b）;
%Image=Image_B;
%figure
%imshow（Image）;
N=512; 
A=zeros（NN）; 
B=zeros（NN）; 
for I=1:1:256 
    for J=1:1:256 
        ImageNum=double（Image（IJ3））; 
        A（IJ）=ImageNum/255; 
        B（IJ）=0; 
    end;
end;
%为降低全息图的动态范围，乘一随机位相因子% 
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）; 
       Max=max（max（abs（F）））; 
       F=F/Max; 
       A=real（F）; 
       B=imag（F）; 
%对物函数做二维FFT变换% 
%定义载波参数 
alpha=0.5; 
for I=1:1:N 
    for J=1:1:N 
        Xcos=（J-1）/127; 
        A1（IJ）=cos（2*PI*alpha*Xcos）; 
        B1（I