简介这份数字全息仿真实验资料面向光学信息处理、计算全息方向的初学者与进阶学习者以MATLAB为工具帮助读者理解并复现数字全息从记录到再现的完整流程。压缩包共2个文件包含1个bmp位图与1个m脚本整体约3.42MB其中位图作为原始全息图数据输入脚本则承担读取图像、傅里叶变换、空域滤波与逆变换再现等核心计算。内容围绕全息实验原理、MATLAB编程、零级像消除、波前记录与再现等关键知识点展开读者可借此掌握将光强分布转化为复数场、用高斯滤波去除零级干扰、再经模拟衍射重现物体三维信息的方法。已有2430人学习下载适合希望从代码层面吃透数字全息计算过程、为光学成像与计算全息应用打基础的学习者。1. 数字全息仿真实验从一张干涉图到三维再现MATLAB 能帮你省下多少台光学平台如果你正在做光学课程设计、显微成像算法预研或者单纯想搞明白“全息再现”到底是怎么从一张条纹图里把物体“捞”出来的那数字全息仿真实验几乎是绕不开的一环。它的核心逻辑很直接用 MATLAB 在计算机里模拟物光与参考光的干涉生成全息图再用衍射计算把物光波前重建出来。整个过程不需要真实激光器、隔振平台和光路准直一台笔记本就能跑通从记录到再现的完整链路。适合谁做干涉测量、粒子场分析、微结构三维形貌恢复的从业者以及想用仿真验证算法再上真实光路的人。这一章不展开公式先把“为什么值得做”和“能解决什么”说清楚。2. 全息记录与再现的数学骨架为什么 MATLAB 是首选验证工具数字全息仿真的本质是用离散采样和数值衍射来逼近连续光学过程。MATLAB 之所以在这个方向被大量使用不是因为它快而是因为它的矩阵运算天然贴合波前传播的离散表达而且可视化链路短改一个参数就能立刻看到条纹和再现像的变化。这一章把记录和再现的数学骨架拆开落到可执行的 MATLAB 代码上。2.1 物光、参考光与干涉条纹的生成逻辑全息记录的核心是物光波 (O(x,y)) 与参考光波 (R(x,y)) 的干涉。在仿真里物光通常用一个带相位调制的复振幅表示参考光用平面波或球面波。干涉强度写成[ I(x,y)|O|^2|R|^2O R^O^R ]前两项是直流分量后两项才是携带相位信息的干涉项。很多新手一上来就把物光设成纯振幅物体结果再现出来只有轮廓没有深度这就是忽略了相位调制。常见做法是给物体加一个二次相位或随机相位模拟漫反射或透镜引入的波前弯曲。下面是最小可运行的记录端代码生成一张离轴全息图% 参数设置 N 1024; % 采样点数 lambda 532e-9; % 波长 532nm dx 4e-6; % 像素尺寸 4um k 2*pi/lambda; % 坐标网格 x (-N/2:N/2-1)*dx; [X, Y] meshgrid(x, x); % 模拟物体一个带相位调制的矩形 obj zeros(N); obj(abs(X)0.5e-3 abs(Y)0.5e-3) 1; phase_obj exp(1i * 0.8 * pi * obj); % 相位调制 % 物光加一个二次相位模拟传播 z_obj 0.15; % 物到记录面距离 O phase_obj .* exp(1i*k/(2*z_obj)*(X.^2 Y.^2)); % 参考光离轴平面波 theta 1.5 * pi/180; % 离轴角 R exp(1i*k*(X*sin(theta))); % 干涉强度 I abs(O R).^2; I I / max(I(:)); % 归一化 % 显示 figure; imshow(I, []); title(离轴全息图);这段代码里z_obj控制物光传播距离theta控制离轴角。离轴角不能太小否则再现时零级和正负一级会重叠也不能太大否则条纹频率超过采样极限出现混叠。经验值是让条纹周期至少占 3 到 4 个像素。phase_obj里的系数 0.8 决定相位调制深度改这个值能直接看到再现像对比度的变化。2.2 角谱法再现从全息图到物光波前记录端生成全息图后再现的本质是让全息图作为衍射屏在特定距离上计算衍射场。MATLAB 里最稳的再现算法是角谱法因为它对传播距离没有近轴近似限制近场远场都能用。角谱法的传递函数写成[ H(f_x,f_y)\exp\left(j k z \sqrt{1-(\lambda f_x)^2-(\lambda f_y)^2}\right) ]实现时先对全息图做二维 FFT乘上传递函数再逆变换回来。下面是对上面生成的全息图做再现的代码% 再现距离 z_rec 0.15; % 与记录距离一致 % 频率坐标 fx (-N/2:N/2-1)/(N*dx); [FX, FY] meshgrid(fx, fx); % 角谱传递函数 H exp(1i*k*z_rec*sqrt(1 - (lambda*FX).^2 - (lambda*FY).^2)); H(abs(lambda*FX)1 | abs(lambda*FY)1) 0; % 倏逝波截断 % 再现 U_rec ifft2(fft2(I) .* H); U_rec fftshift(U_rec); % 显示强度 figure; imshow(abs(U_rec), []); title(再现像强度);这里有两个关键点。第一z_rec必须和记录时的z_obj匹配否则再现像会散焦。第二H里的倏逝波截断不能省否则高频噪声会污染再现像。实际调试时如果再现像出现明显振铃先检查频率坐标是否写成了(-N/2:N/2-1)/(N*dx)很多人会误写成(-N/2:N/2-1)/N导致频率尺度错位。2.3 零级项抑制为什么你的再现像总被亮斑淹没离轴全息虽然能把零级和一级在频域分开但零级项能量通常远大于干涉项直接再现时中心会有一个大亮斑把细节压住。常见做法有三种频域滤波、相减法和希尔伯特变换法。频域滤波最直观在 FFT 后把零级对应的低频区域置零再逆变换。下面是一个简单的频域滤波示例% 频域滤波抑制零级 I_fft fftshift(fft2(I)); mask ones(N); center N/2 1; r 30; % 滤波半径需根据条纹频率调整 [XX, YY] meshgrid(1:N, 1:N); mask((XX-center).^2 (YY-center).^2 r^2) 0; I_filtered real(ifft2(ifftshift(I_fft .* mask))); % 用滤波后的全息图再现 U_rec2 ifft2(fft2(I_filtered) .* H); U_rec2 fftshift(U_rec2); figure; imshow(abs(U_rec2), []); title(零级抑制后再现像);滤波半径r是核心参数。太小零级残留太大干涉项被削掉再现像变暗甚至丢失细节。我一般会先看全息图的频谱找到零级峰和一级峰的间距取间距的三分之一到二分之一作为r。如果条纹频率较高r可以适当放大但不要超过一级峰的位置。3. 参数怎么设采样、波长、距离与再现质量的定量关系仿真跑通不难难的是让再现像清晰、无混叠、无伪影。这一章把几个核心参数的定量关系讲清楚并给出可复现的调参流程。3.1 采样条件条纹频率与像素尺寸的硬约束数字全息里最容易翻车的地方就是采样不足。离轴全息的最大条纹频率由离轴角决定[ f_{\max} \frac{2\sin(\theta/2)}{\lambda} ]根据奈奎斯特采样定理像素尺寸 (dx) 必须满足[ dx \frac{1}{2 f_{\max}} ]实际调试时我一般先定波长和像素尺寸再反推最大允许离轴角。比如lambda532nm、dx4um最大条纹频率约1/(2*4e-6)125000 1/m对应最大离轴角约 3.8 度。超过这个角度条纹就会出现混叠再现像里会出现对称的假峰。下面这段代码可以快速检查当前参数是否满足采样条件% 采样条件检查 theta_max asin(lambda/(2*dx)); % 最大允许离轴角 fprintf(最大允许离轴角: %.2f 度\n, theta_max*180/pi); % 如果设定角度超过条纹会混叠 theta_set 1.5 * pi/180; if theta_set theta_max warning(离轴角超过采样极限会出现混叠); end这个检查应该放在参数设置之后、生成全息图之前。很多人跳过这一步结果再现像里出现莫名其妙的对称条纹还以为是算法问题。3.2 再现距离与聚焦如何判断再现像是否在焦再现距离不匹配时再现像会散焦表现为边缘模糊、对比度下降。判断是否在焦最直接的方法是看再现像的梯度或方差。在焦时图像梯度最大方差也最大。可以写一个简单的扫描循环% 再现距离扫描 z_scan linspace(0.10, 0.20, 50); sharpness zeros(size(z_scan)); for i 1:length(z_scan) z z_scan(i); H_scan exp(1i*k*z*sqrt(1 - (lambda*FX).^2 - (lambda*FY).^2)); H_scan(abs(lambda*FX)1 | abs(lambda*FY)1) 0; U_scan fftshift(ifft2(fft2(I_filtered) .* H_scan)); sharpness(i) var(abs(U_scan(:))); end [~, idx] max(sharpness); fprintf(最佳再现距离: %.4f m\n, z_scan(idx)); figure; plot(z_scan, sharpness); xlabel(再现距离 (m)); ylabel(清晰度指标);这个扫描的步长建议先粗后细。先以 1mm 步长扫一遍找到峰值附近再以 0.1mm 步长细化。注意sharpness指标对噪声敏感如果全息图噪声大可以先做一次中值滤波再扫描。3.3 波长与像素尺寸的匹配不同配置下的参数表不同应用场景对波长和像素尺寸的要求不同。下面这张表是我在多个仿真项目中总结的常用配置可以直接作为起点应用场景波长 (nm)像素尺寸 (um)典型再现距离 (m)离轴角 (度)微结构形貌5323.450.10–0.201.0–2.0粒子场分析632.84.650.15–0.301.5–2.5生物细胞成像4052.200.05–0.100.8–1.5大视场测量5326.500.30–0.502.0–3.0这张表里的离轴角都是保守值实际可以适当放宽但不要超过采样极限。如果换用更小的像素尺寸离轴角可以更大但视场会变小需要权衡。4. 避坑与排查数字全息仿真里最容易翻车的 5 个地方这一章记录的是我在实际调试中反复遇到的坑每条按“现象 → 原因 → 解决”写方便对照排查。4.1 再现像出现对称双像现象再现像里物体旁边出现一个对称的虚像两个像关于中心对称。原因这是离轴全息里正负一级没有完全分开的典型表现。当离轴角太小或者滤波半径太大正一级和负一级在频域重叠逆变换后就会同时出现实像和虚像。解决增大离轴角或者减小滤波半径。先看频谱确认正负一级是否分开。如果分开不明显优先调大离轴角但不要超过采样极限。4.2 再现像中心有大亮斑现象再现像中心有一个圆形亮斑细节被完全压住。原因零级项没有抑制干净。离轴全息虽然能把零级和一级在频域分开但如果滤波半径太小零级残留仍然很强。解决增大滤波半径或者改用相减法。相减法的思路是记录两幅全息图一幅有物体一幅没有物体两者相减消掉零级和参考光背景。相减法对振动敏感仿真里可以完美实现实际光路里需要稳定平台。4.3 再现像边缘出现振铃现象再现像边缘出现明暗相间的条纹像水波纹一样。原因频域滤波时用了硬边掩膜相当于在频域加了一个矩形窗逆变换后产生吉布斯振铃。解决把硬边掩膜换成高斯掩膜或汉宁窗。高斯掩膜的标准差取滤波半径的三分之一左右过渡更平滑。下面是一个高斯掩膜的示例% 高斯掩膜替代硬边掩膜 sigma r / 3; mask_gauss 1 - exp(-((XX-center).^2 (YY-center).^2)/(2*sigma^2)); I_filtered_g real(ifft2(ifftshift(I_fft .* mask_gauss)));4.4 再现距离扫描找不到峰值现象扫描再现距离时清晰度曲线平坦找不到明显峰值。原因全息图信噪比太低或者物体本身没有足够的纹理。如果物体是纯相位且相位调制很弱再现像的对比度变化不明显清晰度指标就区分不出来。解决先检查全息图是否有明显条纹。如果条纹对比度低增大相位调制深度或参考光强度。如果物体纹理不足可以在物体上叠加一个随机相位屏增加高频信息。4.5 大尺寸仿真内存溢出现象把采样点数从 1024 提高到 4096 后MATLAB 报内存不足。原因角谱法需要存储多个 N×N 复数矩阵N4096 时单个矩阵就占 256MB加上中间变量很容易超过默认内存限制。解决改用分块处理或者用单精度single替代双精度double。如果只是验证算法先用 1024 跑通再逐步放大。另外fft2默认对复数做双精度变换可以用fft2(single(I))降低内存占用但要注意精度损失。5. 进阶技巧用相位解包裹和三维重建把仿真结果用起来仿真跑通、再现像清晰之后真正有价值的是从相位里恢复三维形貌。数字全息的相位是包裹在 ((-\pi, \pi]) 之间的需要解包裹才能得到连续相位。MATLAB 里可以用unwrap函数做二维解包裹但直接对含噪声的相位用unwrap容易出错。我一般先用质量图引导的解包裹或者用最小二乘法。下面是一个简单的相位解包裹和三维显示流程% 提取再现像相位 phase_rec angle(U_rec2); % 相位解包裹 phase_unwrap unwrap(unwrap(phase_rec, [], 1), [], 2); % 转换为高度 height phase_unwrap * lambda / (4*pi); % 反射式配置 % 三维显示 figure; surf(X, Y, height); shading interp; xlabel(x (m)); ylabel(y (m)); zlabel(高度 (m)); title(三维形貌重建);这里lambda/(4*pi)是反射式配置的转换系数透射式用lambda/(2*pi)。解包裹后的相位如果有整体倾斜说明参考光不是理想平面波可以在频域把一级峰移到中心来消除倾斜。具体做法是在滤波后把一级峰平移到频谱中心再逆变换这样再现像的相位就是平的。另一个进阶方向是提高再现分辨率。角谱法的分辨率受限于数值孔径可以通过多角度照明或合成孔径来提升。仿真里可以生成多幅不同离轴角的全息图分别再现后在频域拼接等效扩大数值孔径。这个思路在显微成像里很常见但拼接时要注意相位对齐否则会出现拼接缝。我自己的习惯是每次改完参数先跑一遍 1024 的快速验证确认再现像和相位都正常再放大到 2048 或 4096 做精细分析。这样能省下大量等待时间也不容易因为内存问题中断思路。希望帮到你。本文还有配套的精品资源点击获取