用MATLAB软件和液晶光阀实现傅立叶变换计算全息制作及其再现供参习.docVIP

用MATLAB软件液晶光阀实现傅立叶变换计算全息姚雪灿 指导教师 阎晓娜 （上海大学理学院物理系，上海 200444） 摘要：利用MATLAB语言制作了一个迂回相位编码的傅立叶变换全息图，使用电寻址的液晶光阀作为全息图的实时记录介质对得到的傅立叶计算全息图进行光学再现，并对编码过程中加随机相位和不加随机相位后的再现图进行了比较讨论。 关键词：计算全息 傅立叶变换全息 MATLAB 液晶光阀 迂回相位编码 全息制作包括二种方式，光学全息和计算全息。光学全息用光学干涉原理制作，计算全息是用计算机对物波的数学描述编码制作。计算全息可以制作已存在物体的全息图，也可以制作不存在物体的全息图，只要物光波可以用数学。计算全息图要以适合光学再现的尺寸和方式来输出。由于计算全息图上每个抽样单元在微米量级，需要专门的光学缩微照相系统或微光刻系统。在要求较低情况下也可用照相机将显示在计算机屏幕或打印输出的计算全息原图缩拍到高分辨感光胶片上，通过显影、定影等处理得到可用于光学再现的全息图。由于记录介质是照相胶片这就限制了它在实时处理中的应用。近年来，随着高分辨电寻址空间光调制器的发展，像元尺寸在微米量级，像素数超过万的振幅型或相位型空间光调制器已经完全实用化。代替全息干板实现计算全息图的实时输出和再现。本文提出一种利用电寻址液晶光阀作为实时记录介质的计算机制全息图的产生方法，实验结果证明了这种方法的可行性。 Matlab软件实现傅立叶变换计算全息图 傅立叶变换全息记录的复数波面是物光波的傅立叶变换。计算傅立叶变换全息图的制作对物光波抽样、离散傅立叶变换、编码、画图、图像的输出。在制作全息图的过程中,编码是最关键的一步，通过编码把二维光场的复振幅分布变换为全息图的二维透过率分布。迂回相位编码。 设抽样后物光波的复振幅经过离散傅里叶变换后的频谱分布为复数F(m,n), 记为 　F(m,n) = R(m,n)+iI(m,n)F(m,n) = A(m,n)·exp[i(m,n)] (1) 其中, A(m,n)和(m,n)分别代表全息图上各点的幅值和相位A(m,n) =, (m,n) =arctg[I(m,n)/R(m,n)] (2)  由于光学模板的最大透过率为1,所以在编码前还应对A(m,n)的值进行一化,使其最大值为1。假定将物面分为N×N个抽样单元, 抽样间距为δx和δy, 其间距要遵循Nyquist判据。采用罗曼Ⅲ型编码方法,通过改变每个抽样单元内通光孔径的面积来编码振幅,通过改变通光孔径中心与抽样单元中心的位置来编码相位。矩形孔的宽度为Wδx, 其中W为一常数。矩形孔径的高度为Lmnδy,与一化振幅成正比, Pmnδx是孔径中心与单元中心的距离,并与抽样点的位相成正比。孔径参数与复值函数的关系如下=, =/2πK （3） 经过计算,取W =1/2, K =1。 根据以上二元傅里叶变换全息图的实现原理，采用以下的算法思想在MATLAB中进行二元傅里叶变换计算全息图的制作，采用罗曼Ⅲ型编码方式且以字母K为例。其编码如下： b=zeros(128,128); %采样点阵为128X128。 b(8:120,24:40)=1; for M=0:56 b(8+M,(96-M):(112-M))=1; end for N=1:56 b(64+N,(40+N):(56+N))=1; end %K图形用矩阵表示如上。也可选取简单图形F表示。 A=b;a=rand(128,128); aa=exp(i*2*pi.*a); ％给矩阵图一个随机相位，随机相位因子的作用是平滑傅立叶频谱。 AA=double(A).*aa; Afft2=fft2(AA); ％用matlab工具箱对带有随机相位因子的图进行快速傅立叶。 A1=abs(Afft2); B2=angle(Afft2)/(2*pi);%对频谱的相位进行归一化 A1max=max(max(A1)); A1=A1/A1max; %对频谱的幅值进行归一化 s=1; ％设定正方形单元的边长 figure; axis([0 128 0 128]); hold on; for J=1:128 y0=s/2+(J-1)*s; for I=1:128 x0=s/2+(I-1)*s; ％ 计算抽样单元的中心位置 H=A1(J,I)*s; ％矩形高度直接等于归一化的频谱幅值 F1=B2(J,I)*s; ％偏离单元中心的量 W=s/2; ％矩形宽度 if abs (F1)=s/4 x2=x0-W/2;