从‘相机人’到模糊图像：用Matlab角谱法理解光的衍射与成像

从‘相机人’到模糊图像用Matlab角谱法理解光的衍射与成像当一束光穿过显微镜下的透明细胞样本时我们看到的往往是一片模糊——因为细胞本身几乎不吸收光线。但奇妙的是随着观察距离的变化这些透明的细胞轮廓竟会逐渐显现。这种现象背后隐藏着傅里叶光学中关于相位与振幅相互转换的深刻原理。本文将带你用Matlab搭建一个光学实验室通过动态仿真揭示光传播过程中那些肉眼看不见的物理机制。1. 复振幅光波的完整身份证在光学世界里光波可以用一个复数U(x,y) A(x,y)·exp(jφ(x,y))完整描述。其中A(x,y)代表振幅分布——光的强弱变化φ(x,y)代表相位分布——光的波前形状这个复数表示就是复振幅它包含了光波的全部信息。有趣的是人眼和普通相机只能检测光强 I |U|² A²相位信息φ在常规观测中隐形了% 创建合成复振幅示例 A im2double(imread(cameraman.tif)); % 振幅用相机人图像 Phi im2double(imread(lena.tif)); % 相位用Lena图像 U0 A.*exp(1j*2*pi*Phi); % 合成复振幅2. 角谱法光传播的频谱视角光在自由空间传播时角谱理论告诉我们传播过程本质上是空间频谱的相位调制。具体数学表达为$$ A(f_X,f_Y;z) A(f_X,f_Y;0) \cdot \exp\left(j2\pi z\sqrt{\frac{1}{\lambda^2}-f_X^2-f_Y^2}\right) $$其中关键点传递函数exp项就是自由空间的传递函数截止频率当f_X² f_Y² ≥ 1/λ²时光无法传播evanescent wave相位调制不同频率分量获得不同相位延迟传播特性数学表达物理意义低频分量√(1/λ²-f²)为实数正常传播高频分量√(1/λ²-f²)为虚数倏逝波截止频率f_c 1/λ分辨率极限% 角谱传播的核心代码 f sqrt(fx.^2 fy.^2); circ_f (f 1/lambda); % 截止频率滤波器 H exp(1j*2*pi*z*sqrt(1/lambda^2 - f.^2)).*circ_f; Uz_fft U0_fft .* H; % 频域相乘3. 相位显形从透明到可见的魔法纯相位物体如生物细胞在常规成像中几乎不可见但随着传播距离增加相位信息会神奇地转化为强度变化。这是因为相位梯度导致光波偏折传播过程中的衍射效应将相位差转换为强度差最终形成纹影图像Schlieren image实验观察初始距离相位物体几乎不可见中等距离出现边缘增强效果远距离形成复杂的衍射图样提示在显微镜成像中这解释了为何需要离焦才能观察到透明样本的轮廓4. 振幅扩散图像模糊的物理本质对于纯振幅物体如相机人图像传播过程表现为经典的模糊现象。其物理机制是高频衰减角谱传递函数抑制高频分量能量扩散点扩散函数导致像素间串扰分辨率下降遵循瑞利判据模糊程度与传播距离z的关系可以用以下经验公式估算$$ \text{模糊半径} \approx \frac{\lambda z}{D} $$其中D是物体的特征尺寸。5. 仿真实验动态观察传播过程让我们用Matlab实现一个交互式仿真% 动态传播演示 figure; for z 0:10:500 % 从0到500mm逐步传播 H exp(1j*2*pi*z*sqrt(1/lambda^2 - f.^2)).*circ_f; Uz ifft2(U0_fft .* H); imshow(abs(Uz), []); title([传播距离 z num2str(z) mm]); pause(0.1); end关键观察点相位信息何时开始影响强度振幅图像的模糊速度与波长关系不同空间频率成分的衰减差异6. 实际应用从仿真到现实这些原理在多个领域有重要应用相位显微术观察透明生物样本全息成像记录和重建光场光学加密利用相位编码信息图像处理模拟光学模糊效果例如在数字全息中我们可以通过记录不同传播距离的衍射图样逆向求解出原始相位分布。