用双相位编码方法，采用matlab编写全息成像程序，给出具体程序代码

双相位编码是一种常见的全息成像技术，在MATLAB中编写全息图像生成程序，可以使用傅里叶变换和复数运算来实现。以下是简化的示例代码，它假设你已经有了原始物体波和参考波的数据：

% 假设obj_wave 和 ref_wave 分别代表物体波和参考波的复杂数组
% 这些通常通过物理实验测量得到

% 双相位编码操作
phase_encoded = obj_wave + ref_wave;

% 将数据转换到频率域（使用快速傅立叶变换）
hologram_fft = fft(phase_encoded);

% 对频域数据取对数，用于增强全息图的对比度
hologram_log = log(abs(hologram_fft));

% 取实部或虚部，因为全息图通常是基于干涉图案的幅度信息
% 这部分取决于你的应用需求
hologram_real = real(hologram_log); % 如果只需要实部

% 求反傅立叶变换恢复到空间域
hologram_recon = ifft(hologram_real);

% 显示全息图
imagesc(-size(hologram_recon, 2)/2:size(hologram_recon, 2)/2, ...
       -size(hologram_recon, 1)/2:size(hologram_recon, 1)/2, hologram_recon);
colorbar; % 添加颜色刻度

%