内容发布更新时间 : 2024/12/24 22:13:45星期一 下面是文章的全部内容请认真阅读。

原图（秦时明月-卫庄）

计算傅里叶变化全息

离散图（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