菲涅尔转换的数值计算解析.docVIP

菲涅尔转换的数值计算 在这篇论文中，我们通过均匀空间采样所获得的有限数值来处理计算菲涅尔转换积分这一问题。通常简单的数值采样规则能够让使用者对于任何散射距离都能计算出他的分布。这些规则是如何被扩展到利用基于快速傅里叶变换这一算法来提高计算效率的,在这篇论文中也有所解释。这篇论文也把本方法与其他理论方法做了比较。 引言 光的传播和衍射是个经典问题，在当下仍十分重要。麦斯威尔等式为检测光的行为提供了一个严格的框架，他把光描述成一种电磁波。电场和磁场组成部分的行为都应该特别的考虑，例如，当分析数值孔径高的镜头，就如在显微镜中或者是为平面等离子体匹配模式时发现的。然而在许多场合，对于一个给定的物理问题进行模式匹配时，电磁场的全矢量描述是不需要的，一个简单地标量模型可以用在这里。用一个标量模型去表示一个光学场是一个有用的近似表示。在光学制度中，它为光通过孔径（在这里孔径的大小远大于光波长）的衍射提供了一个高度准确的预测，可参见【1】的实验结果。其他相关理论和数值计算在【2,3】中被提出来了。在这篇论文中我们假设一个标量模型对于我们的目的是足够精确的。 在光传输的标量描述中，区分以下两个定义是非常有用的，即傍轴和非傍轴衍射积分。两类衍射积分来自于麦斯威尔等式，参见【4,5】。然而在傍轴处理过程中，在推倒时应做更多的近似处理。因此，傍轴模式更加脱离于真实的衍射基础物理过程，伴随着这个结论还有非傍轴衍射积分被认为是更加的准确。我们参见在【4】的3、4章中写的衍射处理，并且据此把基尔霍夫定律或瑞利-索末菲，和一个有限的输入范围。衍射的真正的物理过程使得这个信号在衍射平面中有一个无限的输出空间范围，我们后边会定义为SE。 （2） 从一个有限的采样值而不是从1中考虑的连续情形来计算衍射分布从根本上改变了这一问题的性质。这一理论计算的信号将重新产生一个来自于能够产生范围为SE的输出平面的点1的分析结果。然而除此之外，重现波的一个无限集产生于衍射平面。这些重现波的特性取决于这一运算是如何运行的。 （3） 对于数值计算衍射积分考虑了两种方法。直接计算衍射积分会产生彼此间隔为的复制平面，在此λ是波长，z是传输距离，是在输入空间域的采样间隔。通过运算在输入域中信号的空间频率没有发生相位移动，并且不是通过尼奎斯特标准采样获得的空间频率也将会传播并会加强在衍射域中的场。基于光谱的运算在衍射平面中产生了一个彼此间隔的重现波的无限集，在此是空间频率域的采样间隔。在光谱运算过程中衍射信号的空间频率成分经常被限制在空间频率域。 图1，衍射过程的图示。一个无线平面波入射进平面1中的孔径中，该平面以空间变量X表示。考虑到如何计算P(x,z)中的复振幅，仅仅是开放孔径的光对P(x,z)中的复振幅有贡献，我们把开放孔径内的光设想成由无穷多个点光源发光组成的，每一个点光源都带有强度值和相位值（图中只画了一部分点光源）。因为这个孔径是被理想平面波照明的，每一个第二级点光源都带有相同的强度值和相位值，尽管当考虑到距离l1和l4时，由于光的路径不同每一点光源将积累不同的相位值。 （4）我们建议的抽样规则遵从以上1-3点：SE比和都小。当这一点确定后，这个衍射信号被认为是被很好抽样的。我们再一次标注，以上所采取的方法不同于文献中所提出来的方法，这是因为我们在思考时没有用到菲涅尔啁啾函数的特性。 在2-5节中，我们得出了以上的点1-4。在第六节中，我们把注意力转移到计算菲涅尔转换的高效计算算法，对于这一点我们用到了快速傅里叶变换。两种算法被明确的命名为：直接算法和光谱算法，其中光谱算法是我们在第4和5节中讨论过的直接和光谱计算的快速实现（基于快速傅里叶变换）。我们的理论结果与文献中的其他分析方法所得的结论做了对比。最后，在第七节中，我们考虑一个重要的案例在这个案例中我们在第4和5节中提出的抽样规则是不合适的：也就是说，在衍射区被抽样时我们希望对这个结果做菲涅尔逆变换，这后一种情况对数字全息照相术有重要的寓意，然后我们总结这篇论文的主要发现。 一个分析解决方案 在随后的介绍中为了表述简单，我们执行一维分析。在图1中，我们简绘了问题的特性。平面1中的场我们以大写字母U（x）表示，并假设它是已知的，我们的任务是推导出衍射场的表达式，即平面2的表示式。通过【4】中的分析，我们现在写出菲涅尔转换的定义式： = (1) 在这里是单色相干光（在时间和空间上的）的波长，z是平面1和平面2的轴距离（参见图1），j=，变量X和x分别代表在平面1和平面2上的空间坐标。为了余下的分析我们将降低开始的相位项，即。 在有些衍射问题中有可能能解出等式1的解析解。由于一些原因这些解决办法非常重要：他们洞察了这个衍射过程，他们也可以作为一种