基于MATLAB入射偏振光的全息成像的程序

以下是一个基于MATLAB的入射偏振光全息成像程序的示例：

% 定义全息图参数
N = 512;                % 全息图大小
pixel_size = 1e-6;      % 像素大小
lambda = 633e-9;        % 激光波长
k = 2*pi/lambda;        % 波数

% 定义物体参数
obj_size = 50e-6;       % 物体尺寸
obj = zeros(N);         % 生成一个全0矩阵
x = linspace(-obj_size/2, obj_size/2, N);    % 生成物体坐标
y = x;
[X, Y] = meshgrid(x, y);

% 定义物体形状
r = sqrt(X.^2 + Y.^2);
obj(r < obj_size/4) = 1;    % 在物体中心生成一个圆形

% 定义入射光参数
pol_angle = pi/4;       % 光偏振角度
pol = [cos(pol_angle), sin(pol_angle)];    % 光偏振方向

% 计算物体上每个像素点的复振幅
obj_complex = obj .* exp(1i*k*obj);

% 计算物体傅里叶变换
obj_fft = fftshift(fft2(obj_complex));

% 计算全息图
H = exp(1i*k*pixel_size*(X.^2 + Y.^2)) .* exp(-1i*k*pixel_size*(X*pol(1) + Y*pol(2))) .* obj_fft;

% 计算全息图的逆傅里叶变换
recon = ifft2(ifftshift(H));

% 显示结果
figure;
subplot(2,2,1); imshow(abs(obj), []); title('物体振幅');
subplot(2,2,2); imshow(angle(obj_complex), []); title('物体相位');
subplot(2,2,3); imshow(abs(H), []); title('全息图振幅');
subplot(2,2,4); imshow(angle(H), []); title('全息图相位');

figure;
subplot(1,2,1); imshow(abs(recon), []); title('全息图重构振幅');
subplot(1,2,2); imshow(angle(recon), []); title('全息图重构相位');

这个程序生成了一个直径为50微米的圆形物体，使用了45度偏振的入射光进行全息成像，最后重构了全息图的振幅和相位。您可以根据自己的需求调整程序中的参数和物体形状来生成不同的全息图。