菲涅尔转换的数值计算.docx

菲涅尔转换的数值计算 在这篇论文中， 我们通过均匀空间采样所获得的有限数值来处理计算菲涅尔转换积分这 一问题。 通常简单的数值采样规则能够让使用者对于任何散射距离都能计算出他的分布。 这 些规则是如何被扩展到利用基于快速傅里叶变换这一算法来提高计算效率的 , 在这篇论文中 也有所解释。这篇论文也把本方法与其他理论方法做了比较。 1、 引言 光的传播和衍射是个经典问题， 在当下仍十分重要。 麦斯威尔等式为检测光 的行为提供了一个严格的框架， 他把光描述成一种电磁波。 电场和磁场组成部分 的行为都应该特别的考虑， 例如， 当分析数值孔径高的镜头， 就如在显微镜中或 者是为平面等离子体匹配模式时发现的。 然而在许多场合， 对于一个给定的物理 问题进行模式匹配时， 电磁场的全矢量描述是不需要的， 一个简单地标量模型可 以用在这里。 用一个标量模型去表示一个光学场是一个有用的近似表示。 在光学 制度中，它为光通过孔径 （在这里孔径的大小远大于光波长） 的衍射提供了一个 高度准确的预测，可参见【 1】的实验结果。其他相关理论和数值计算在【 2,3 】 中被提出来了。在这篇论文中我们假设一个标量模型对于我们的目的是足够精确 的。 在光传输的标量描述中， 区分以下两个定义是非常有用的， 即傍轴和非傍轴 衍射积分。两类衍射积分来自于麦斯威尔等式，参见【 4,5 】。然而在傍轴处理过 程中，在推倒时应做更多的近似处理。 因此，傍轴模式更加脱离于真实的衍射基 础物理过程， 伴随着这个结论还有非傍轴衍射积分被认为是更加的准确。 我们参 见在【 4】的 3、4 章中写的衍射处理，并且据此把基尔霍夫定律或瑞利 -索末菲 第一和第二定律归类为非傍轴衍射积分。菲涅尔转换是一个傍轴衍射积分。 在这篇论文中我们用有限的采样数值来研究菲涅尔转换的数值计算问题。 这 是因为数值模拟技术对模拟光传输和为没有解析解的普通衍射问题计算解析解 是相当重要的。菲涅尔转换在数字全息术【 6-12 】、迭代相位检索【 13-15】和估 计相位分布的强度方式的传输中扮演着不可缺少的角色， 估计相位分布的强度方 式的传输在生物图像应用【 16,17 】中起着重要的作用。菲涅尔转换也广泛应用 于针对波束成形【 18,19 】的衍射光学元件的设计中，在计量应用中衍射光斑区 域的去相关也被应用于判定变形， 例如， 在测试时光学粗糙表面的压力、 应变以 及裂缝，【20,22 】，或者是位移和倾斜的检测【 23,24 】。在这所有的应用当中， 菲涅尔转换的数值计算是必要的， 因此完全明白如何执行运算具有显著的理论和 实践重要性。 或许在理论计算过程中遇到的第一个问题会是如何选择一个合适的均匀分 布的抽样值来代替输入平面的复杂标量波场。 在文献中有许多通过检测需要多少 抽样值来处理这个问题的手稿， 因此在检测内核或者是菲涅尔啁啾函数时尼奎斯 特采样率（或者是一个相关的标准）被采用了，参见例子【 25-31 】。然而这个处 理可以等效于准确的理论计算，它往往可以是不必要的数值密集型【 28】。用于 计算输入输出信号的范围大小被发现也影响光学采样率， 然而在文献中也出现了 一些关于理想采样率的意见分歧， 例如【 29-31 】。通常也不是很清楚当违反采样 规则时数值计算是为何变得不准确的。 因为这些因素， 很有必要来重新审视这个 问题，即为衍射积分运算提供一个合适的采样率。 尽管我们的分析是针对于菲涅 尔转换，但是这一结果可以概括到在【29-31】中讨论的傍轴积分中。 在此，我们以一种完全不同的方式来处理采样问题， 即依靠一种对数值衍射 过程的直观的物理描述。这一描述对分析普遍的衍射问题和准确的确定通过采样 所引进的理想误差都是有帮助的。和文献中的其他手稿不同的是我们在决定合适 的采样率时不考虑菲涅尔内核和啁啾函数在计算中的作用。 相反地，我们依靠输 入域的空间和空间频率范围并依靠衍射平面中重现波的位置来定义采样率。 这些 重现波的出现是由于采样操作（平面 1中的复杂信号的离散化，看图1和第3 部分）。这些重现波和他们特性之间的离散化取决于菲涅尔转换是如何运算的。 鉴于我们正在计算的信号的输出范围小于相邻重现波之间的离散这一事实， 我们 认为我们的信号是被合适采样的。 这篇论文的论据基于以下考虑： （1） 在空间和空间频率域有一个有效的有限的支持的信号被定义为使用来 自在【27,28,32-35】中描述的空间-带宽产品的概念。这种信号包含定义的（或 使用者选择的）最大的空间频率， fx，和一个有限的输入范围。衍射的真正的 物理过程使得这个信号在衍射平面中有一个无限的输出空间范围， 我们后边会定 义为S巳 （2） 从一个有限的采样值而不是从1中考虑的连续情形来计算衍射分布从 根本上改变了这一问题的性