在
系统的不同平面上,电磁场分量的傅里叶变换是连接空间域和k域的
物理光学建模中的频繁操作。我们介绍一个场所谓的几何区域,在该区域中傅里叶变换可以在不进行积分的情况下得到,总之是以非常有效的数值方式得到。在几何场域中,场由波前相位控制,因此允许我们将稳定相位的概念应用于傅里叶变换积分,我们将所得到的傅里叶变换算法称为几何傅立叶变换,这项技术被证明是快速物理光学的基础支柱。
@8=vFP' wqlcLIJPR 1.光学傅立叶变换
/Fv/oY Z&FkLww 在物理光学中,我们处理电磁场的六个复数场分量(分别为E和H)。在空间域,他们表示为
)RwBg8 2'wr={>W 2yln7[a |e*Gz D 其中
,傅立叶变换到k域定义为
~n[b^b
(2)
'lhP!E_)q 其中,我们使用符号
2yN%~C?$ co5y"yj_ (3)
.3VK;au\\ #(%t*"IY; 方程2中积分的数值评估需要对a和k域中的场进行取样,我们用N表示采样点的数量,所得的离散傅里叶变换构成了N2运算。然而快速傅里叶变换(FFT)算法在N中是线性的,这在
原理上使快速物理光学建模成为可能,但FFT需要
的采样。在光学中,我们通常有强梯度的相位函数,从而导致很大的N值,只有在十分对称的
光学系统中,N才可以很小。因此,尽管FFT在N中是线性的,但是我们很容易在光学上遇到N太大而不能进行快速计算傅里叶变换的问题,这是快速物理光学概念的严重阻碍。
E@w[ LBiowd[ 为了进一步研究,我们用波前相位Ψ将
分解(跳过ω)为
Z-)[1+Hs $]@O/[ (4)
UvkJ?Bu j2|XDOf 对于所有分量都是一样的。 显然,方程 4中的分解是模糊的,其依赖于从源场出发建模中恰当的相位处理方式。由定义
得分解结果
\\jIl3Z [@m[V1D (5)
[6D>f?z J &!B|TS 类似地,我们可以得到
zKLn!b#> Dp*:Q){>E (6)
L@w0N)P<!{ tbq|," 其中波前相位
在k域上。应该提到的是,根据方程 5与 在几何光学上是已知的,然后
,S为光程函数。我们想强调的是,方程 5的分解在物理光学中是更一般和纯粹的数学方法,我们的目标可以表述如下:我们对不通过采样波前相位因素来进行傅里叶变换的技术十分感兴趣,此时Ψ和
是可通过半解析傅里叶变换实现的二次多项式的形式[1]。这里我们想讨论一个概念,适用于一般的波前相位,但在强波前相位近似,它使用稳定相位的概念。
a{h%DpG u(W^Nou/+ 2 几何傅里叶变换理论
KNy`Lj)VPY 稳定相方法的应用在光学中是众所周知的,例如,用于讨论[2]中的
衍射积分。我们将其用于快速计算方程2的傅里叶变换积分。为此,我们假设除临界点附近以外
在通过z的平面内具有比U(ρ,z)高得多的空间频率。 根据稳定相位的概念,直接导致基本方程(跳过z )
d,W/M(S [ S5bj]D (7)
~32Pjk~ P:
n# S % 其中
方程7表示k和p之间的映射,我们假设这个映射是开放、双射和连续的,这意味着它构成了一个同胚,这是波前相位 平滑的数学表达式,并确保k域中的结果场可以在非等距网格上插值。在光学中,当场不在苛性区时,通常满足这种条件,稳定相位的概念也揭示出来
U.: sK* (8)
tNU-2r 由φ(p)的勒让变换
D%-{q>F!gf z^r|3; (9)
lLH$`Wnv EYG"49
c 复函数
vF"c K<FKu $= (10)
E*BSfn&i LZ{YmD&6] 权重因子
取决于φ(p)的二阶导数,该结果通过将空间域中的场值映射到具有附加权重因子的k域来表示傅里叶变换,其仅作为映射本身而依赖于波前相位。因此,傅里叶变换主要执行场分布的几何畸变,我们称之为几何傅里叶变换。
_R)&k%i} N[pZIH5ho= 我们已经开发了一个数值算法来执行几何傅里叶变换。它利用场的混合采样。相比于函数
,波前相位φ(p)本身可以通过少量N(φ)的非等距分布值而
参数化。样条插值的节点是可能的候选项。
=yo{[&Jz RU6KIg{H 而且,我们必须用等距分布的采样点N(U)来处理函数U(p)的采样。一般来说,我们有
,几何傅里叶变换的数值主要基于 中的线性运算,因此速度非常快;N(U)中U值的智能包也可以快速完成,V的采样可以完全避免。总之,当几何傅立叶变换足够精确时,由此产生的数值算法能够实现非常快速的傅里叶变换,对于强波前相位来说就是这种情况。
<"?*zx&