什么是菲涅尔衍射积分公式？
在波动光学中，菲涅尔衍射（Fresnel Diffraction）描述了光波在传播过程中，当观察屏距离衍射孔径较近时发生的衍射现象。与远场的夫琅禾费衍射不同，菲涅尔衍射区域（近场）包含了更丰富的相位信息，是理解光波传播、透镜成像以及全息术的基础。
菲涅尔衍射积分公式是描述这一现象的核心数学工具。它基于惠更斯-菲涅尔原理，认为波前上的每一点都可以看作是发射球面子波的新波源，空间任意一点的光场分布是所有这些子波在该点相干叠加的结果。
核心概念：近场效应
在菲涅尔衍射区域，波前的曲率不可忽略。这意味着光波在传播过程中，振幅和相位都发生复杂变化，导致衍射图样随距离动态演变，不像夫琅禾费衍射那样保持形状不变。
关键参数：菲涅尔数
菲涅尔数（Fresnel Number, ）是判断衍射区域的关键无量纲参数。当 时，通常属于菲涅尔衍射区；当 时，过渡到夫琅禾费衍射区。
物理意义：二次相位因子
与夫琅禾费衍射相比，菲涅尔衍射积分公式多出了一个二次相位因子 ，这反映了球面波传播的几何特性。
公式推导与数学表达
为了准确计算菲涅尔衍射，我们需要从标量衍射理论出发。假设单色光波波长为 ，波数为 。孔径平面为 ，观察平面为 ，两者距离为 。
1. 基尔霍夫衍射公式
根据基尔霍夫衍射理论，观察点 的复振幅 可以表示为：
U(x,y,z) = frac{1}{jlambda} iint_{Aperture} U_0(x',y') frac{e^{jkr}}{r} K( heta) dx' dy'
其中 是孔径点到观察点的距离， 是倾斜因子。
2. 菲涅尔近似条件
在菲涅尔衍射区域，我们假设 远大于孔径尺寸，满足以下近似条件：
- 分母中的
- 指数项中的 保留二次项：
这些近似被称为菲涅尔近似，它使得积分变得可解析求解或高效数值计算。
3. 菲涅尔衍射积分公式最终形式
经过简化，菲涅尔衍射积分公式可以写为卷积形式：
U(x,y) = frac{e^{jkz}}{jlambda z} e^{jfrac{k}{2z}(x^2+y^2)} iint_{-infty}^{infty} left[ U_0(x',y') e^{jfrac{k}{2z}(x'^2+y'^2)}
ight] e^{-jfrac{2pi}{lambda z}(xx'+yy')} dx' dy'
这个公式表明，菲涅尔衍射场可以通过对孔径场乘以二次相位因子后，进行傅里叶变换，再乘以另一个二次相位因子得到。这一特性为快速算法（如FFT）的应用奠定了基础。
主要应用场景与工程价值
菲涅尔衍射积分公式不仅是理论物理的基石，在现代光学工程中也有着广泛的应用。以下是几个关键领域：
计算机生成全息图（CGH）
在全息术中，记录或重建全息图时，必须精确计算光波从全息图平面到观察者眼睛或相机传感器的传播。菲涅尔衍射模型是CGH计算的核心，因为它能准确模拟近场下的波前传播，生成逼真的三维图像。
例如，在设计用于增强现实（AR）显示的全息光栅时，工程师使用菲涅尔衍射积分来优化光栅周期和深度，以确保特定的衍射效率和视场角。
非相干成像系统分析
在显微镜和望远镜系统中，菲涅尔衍射决定了系统的分辨率和点扩散函数（PSF）。当物镜孔径有限时，光波通过孔径发生衍射，导致理想点源成像为一个艾里斑。理解这一过程有助于设计消色差透镜和优化光学系统性能。
激光束传播与整形
高斯光束在通过光学元件（如透镜、光阑）时的传播行为可以用菲涅尔衍射理论精确描述。在激光切割和焊接中，控制光束在焦深范围内的衍射特性对于保证加工精度至关重要。
无透镜成像与数字全息
数字全息显微镜（DHM）利用传感器直接记录全息图，然后通过数值重建算法（基于菲涅尔衍射积分或角谱法）恢复物体的振幅和相位信息。这种方法无需物镜，适用于生物细胞等大样本的高分辨率三维成像。
数值计算与MATLAB实现示例
由于菲涅尔衍射积分通常没有解析解，数值计算是研究衍射现象的主要手段。下面提供一个基于快速傅里叶变换（FFT）的MATLAB代码示例，用于计算圆形孔径的菲涅尔衍射图样。
算法步骤
- 定义孔径平面网格和孔径函数（如圆形孔径）。
- 计算孔径场的傅里叶变换。
- 乘以二次相位因子和传播相位因子。
- 进行逆傅里叶变换得到观察平面的场分布。
- 计算强度分布并绘图。
MATLAB代码示例
% 参数设置
lambda = 632.8e-9; % 波长 (m)
z = 1; % 传播距离 (m)
k = 2pi/lambda; % 波数
N = 1024; % 网格点数
L = 0.01; % 孔径尺寸 (m)
dx = L/N; % 网格间距
x = (-N/2:N/2-1)dx;
[X,Y] = meshgrid(x,x);
% 圆形孔径
R = sqrt(X.^2 + Y.^2);
aperture = double(R < 0.002); % 半径2mm的圆孔
% 菲涅尔衍射计算 (基于FFT)
U0 = aperture;
% 二次相位因子
phase1 = exp(1jk(X.^2+Y.^2)/(2z));
% 傅里叶变换
U_hat = fftshift(fft2(fftshift(U0)));
% 频域坐标
fx = (-N/2:N/2-1)/(Ndx);
[FY,FX] = meshgrid(fx,fx);
% 二次相位因子 (频域)
phase2 = exp(1jpilambdaz(FX.^2+FY.^2));
% 逆傅里叶变换
U = fftshift(ifft2(fftshift(U_hat.phase2)));
% 最终相位因子
U_final = phase1.U;
% 显示结果
figure;
subplot(1,2,1); imagesc(aperture); title('孔径'); axis image;
subplot(1,2,2); imagesc(abs(U_final).^2); title('菲涅尔衍射强度'); axis image;
colormap(jet);
此代码展示了如何利用菲涅尔衍射积分的卷积形式，通过FFT高效计算衍射图样。用户可以根据需要调整孔径形状、波长和传播距离，模拟不同的衍射场景。
历史沿革与关键里程碑
理解菲涅尔衍射的发展历史，有助于我们把握光学理论的演进脉络。
菲涅尔提出波动理论
奥古斯丁·让·菲涅尔（Augustin-Jean Fresnel）独立发展了光的波动理论，并提出了菲涅尔衍射的积分公式，成功解释了泊松亮斑等现象，有力支持了光的波动说。
基尔霍夫衍射理论
古斯塔夫·基尔霍夫（Gustav Kirchhoff）基于电磁波方程，给出了严格的衍射积分公式，为菲涅尔衍射提供了更坚实的数学基础。
傅里叶光学兴起
帕纳德（P.M. Papas）和利夫（J. Linfoot）等人将傅里叶分析引入衍射理论，揭示了菲涅尔衍射与傅里叶变换之间的深刻联系，为现代光学信息处理奠定了基础。
计算机与激光时代
激光的发明和计算机技术的发展，使得菲涅尔衍射的数值模拟成为可能。计算机生成全息术（CGH）和数字全息术应运而生，极大地拓展了衍射理论的应用范围。
常见问题解答（FAQ）
适用条件主要包括：1) 单色光；2) 标量近似成立（即孔径尺寸远大于波长，且不考虑强偏振效应）；3) 满足菲涅尔近似条件，即传播距离 远大于孔径尺寸，且观察区域不大。具体来说，需满足 。
主要依据菲涅尔数 判断。如果 ，通常使用菲涅尔衍射公式；如果 ，则可以使用更简单的夫琅禾费衍射公式（即傅里叶变换）。在实际仿真中，为了确保精度，即使满足夫琅禾费条件，也可以使用菲涅尔衍射公式，因为它更通用。
二次相位因子 反映了球面波传播的几何相位变化。在菲涅尔近似下，光波从孔径传播到观察屏，不仅经历了平移，还经历了波前的弯曲。这个因子补偿了由于波前曲率引起的相位延迟，使得傅里叶变换能够正确描述衍射场。
最快的方法是基于快速傅里叶变换（FFT）的算法。利用菲涅尔衍射积分的卷积形式，将其转化为傅里叶变换的形式，从而利用FFT的高效性将计算复杂度从 降低到 。此外，可以使用GPU加速（如MATLAB的gpuArray）进一步提���计算速度。