计算傅里叶变换全息

《计算傅里叶变换全息》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、计算傅里叶变化全息原图（秦时明月-卫庄）离散图（512*512）计算全息图再现图实现程序%用imread()函数读入物波函数抽样数据并保存到矩阵A中,这样原始物体就以离散数据形式保存在一个二值矩阵中(矩阵中每个元素非0即1)Image=imread('zcy6202.jpg');N=512;A=zeros(N,N);B=zeros(N,N);forI=1:1:256forJ=1:1:256ImageNum=double(Image(I,J,1));A(I,J)=ImageNum/255;B(I,J)=0;endendfigure;ims

2、how(A);%为降低全息图的动态范围，通过乘随机相位因子的办法来平滑物函数的傅里叶变换谱(起到光学全息中毛玻璃漫射光的作用)PI=3.14159;forI=1:1:NforJ=1:1:NR=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);endend%对物函数做二维快速傅里叶变换（FFT）变换%F=fft2(F);Max=max(max(abs(F)));F=F/Max;A=real(F);B=imag(F);%定

3、义载波参数alpha=0.5;forI=1:1:NforJ=1:1:NXcos=(J-1)/127;A1(I,J)=cos(2*PI*alpha*Xcos);B1(I,J)=sin(2*PI*alpha*Xcos);endend%全息图数据区forI=1:1:NforJ=1:1:NHolodata(I,J)=0.5+0.5*(A(I,J)*A1(I,J)+B(I,J)*B1(I,J));endend%绘制全息图，制作灰阶全息图像，灰度等级256M=512;%定义全息图的大小，Hologram=zeros(M,M);S=M/N;%定义每个

4、抽样单元大小为S,SforI=1:1:NforJ=1:1:NXa=(J-1)*S+1;Xb=J*S;Ya=(I-1)*S+1;Yb=I*S;forIx=Xa:1:XbforIy=Ya:1:YbHologram(Iy,Ix)=Holodata(I,J);endendendendMax=max(max(Hologram));Hologram=Hologram/Max;figure;imshow(Hologram);%计算傅里叶变换全息的再现Object=fft2(Hologram);Object=fftshift(Object);Objec

5、t=abs(Object);Object=1000*Object/max(max(Object));figure;imshow(Object);