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