计算傅里叶变换全息

计算傅里叶变化全息原图（秦时明月-卫庄）
离散图（512*512）
计算全息图
再现图
实现程序
%用imread () 函数读入物波函数抽样数据并保存到矩阵A中,这样原始物体就以离散数据形式保存在一个二值矩阵中(矩阵中每个元素非0即1)
Image=imread('zcy6202.jpg');
N=512;
A=zeros(N,N);
B=zeros(N,N);
for I=1:1:256
for J=1:1:256
ImageNum=double(Image(I,J,1));
A(I,J)=ImageNum/255;
B(I,J)=0;
end
end
figure;
imshow(A);
%为降低全息图的动态范围，通过乘随机相位因子的办法来平滑物函数的傅里叶变换谱(起到光学全息中毛玻璃漫射光的作用)
PI=3.14159;
for I=1:1:N
for J=1:1:N
R=rand(1,1);
B(I,J)=A(I,J)*sin(R*2*PI);
A(I,J)=A(I,J)*cos(R*2*PI);
F(I,J)=A(I,J)+j*B(I,J);
end
end
%对物函数做二维快速傅里叶变换（FFT）变换%
F=fft2(F);
Max=max(max(abs(F)));
F=F/Max;
A=real(F);
B=imag(F);
%定义载波参数
alpha=0.5;
for I=1:1:N
for J=1:1:N
Xcos=(J-1)/127;
A1(I,J)=cos(2*PI*alpha*Xcos);
B1(I,J)=sin(2*PI*alpha*Xcos);
end
end
%全息图数据区
for I=1:1:N
for J=1:1:N
Holodata(I,J)=0.5+0.5*(A(I,J)*A1(I,J)+B(I,J)*B1(I,J));
end
end
%绘制全息图，制作灰阶全息图像，灰度等级256 M=512;%定义全息图的大小，
Hologram=zeros(M,M);
S=M/N;%定义每个抽样单元大小为S,S
for I=1:1:N
for J=1:1:N
Xa=(J-1)*S+1;
Xb=J*S;
Ya=(I-1)*S+1;
Yb=I*S;
for Ix=Xa:1:Xb
for Iy=Ya:1:Yb
Hologram(Iy,Ix)=Holodata(I,J);
end
end
end
end
Max=max(max(Hologram));
Hologram=Hologram/Max;
figure;
imshow(Hologram);
%计算傅里叶变换全息的再现
Object=fft2(Hologram);
Object=fftshift(Object);
Object=abs(Object);
Object=1000*Object/max(max(Object)); figure;
imshow(Object);。