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