离轴全息建模

kkciomp

各位大佬帮忙看看我写的一个离轴全息的程序，不知道错在哪里了
首先是输入一个球面波
lamda = 4.25e-6;
x_min=10e-5;
y_min=10e-5;
D=1;
N=1280;
M=1024;
x=linspace(-(N/2)*x_min,(N/2)*x_min,N);
y=linspace(-(M/2)*y_min,(M/2)*y_min,M);
[X,Y]=meshgrid(x,y);
r=sqrt(D^2+X.^2+Y.^2);
object_wave=exp(1i*2*pi/lamda.*r)./r;




然后加入一个参考波进去





reference_wave=exp(1i*2*pi/lamda.*(x_dot.*cosd(theta)+y_dot.*sind(theta)));
hologram=reference_wave+object_wave;
I_holo=hologram.*conj(hologram);


再重建

function [restruction] = off_restruction(hologram,lamda,D,theta,x_min,y_min,N,M)

x=linspace(-640*x_min,640*x_min,N);
y=linspace(-512*y_min,512*y_min,M);
[X_dot,Y_dot]=meshgrid(x,y);
reference_wave=exp(1i*2*pi/lamda.*(x.*cosd(theta)+D*sind(theta)));
r_hologram=reference_wave.*hologram;
phase_one = exp(1i*2*pi/lamda*D)/1i/lamda/D*exp(1i*2*pi/lamda/2/D.*(X_dot.^2+Y_dot.^2));
phase_two = exp(1i*2*pi/lamda/2/D.*(X_dot.^2+Y_dot.^2));
restruction = fftshift(fft2(r_hologram.*phase_two)).*phase_one;

Alec

是不是显示器分辨率不够，没办法把准确的图像显示出来？看着像是有摩尔纹

dreica

代码很混乱，生成全息图是实数的I_holo为啥重建用复数的r_hologram。标量衍射理论计算可以参考Computational Fourier Optics: A MATLAB Tutorial这本书。

kkciomp

Alec 发表于 2023-3-8 10:06
是不是显示器分辨率不够，没办法把准确的图像显示出来？看着像是有摩尔纹

不太清楚，总感觉复现的有问题，但是一时间又看不出来问题在哪

kkciomp

dreica 发表于 2023-3-8 13:13
代码很混乱，生成全息图是实数的I_holo为啥重建用复数的r_hologram。标量衍射理论计算可以参考Computationa ...

哦对，只能读取到全息图的强度分布，不能用复振幅进行处理！

QEPAS

你看看是不是.这个该用的时候没用，不该用的时候用了，我之前仿真相位恢复算法的时候，这个弄错了，结果差异蛮大的