泽尼克多项式(Zernike-Polynomials)-泽尼克系数教程文件

泽尼克多项式（Zernike Polynomials），泽尼克系数

什么是 Zernike Polynomials

通常人们会使用幂级数展开式的形式来描述光学系统的像差。由于泽尼克多项式和光学检测中观测到的像差多项式的形式是一致的，因而它常常被用来描述波前特性（泽尼克，1934）。但这并不意味着泽尼克多项式就是用来拟合检测数据的最佳多项式形式。在某些情况下，用泽尼克多项式来描述波前数据具有很大的局限性。比如说，当需要考虑空气扰动的时候，泽尼克多项式几乎没有什么价值。同样地，我们也无法找到一组合适的泽尼克多项式来描述单点金刚石车削加工（single point diamond turning process）中的制造误差。为了准确地描述圆锥面光学元件（conical optical elements）的对准误差，必须对泽尼克多项式进行修正。盲目地使用泽尼克多项式来表达检测数据只会导致糟糕的结果。

泽尼克多项式是由无穷数量的多项式完全集组成的，它有两个变量，ρ和θ，它在单位圆内部是连续正交的。需要注意的是，泽尼克多项式仅在单位圆的内部连续区域是正交的，通常在单位圆内部的离散的坐标上是不具备正交性质的。

泽尼克多项式具有三个和其他正交多项式集不一样的性质。

⒈泽尼克多项式Z(ρ, θ)可以被化解为径向坐标ρ和角度坐标θ的函数，其形式如下：

Z (ρ, θ) = R ( ρ ) G ( θ )，

这里，关于角度的函数G(θ)是一个以2π弧度为周期的连续函数，并且满足当坐标系

旋转α角度之后，其形式不发生改变，也就是旋转不变性：

G (θ + α ) = G ( θ ) G ( α )

其三角函数集形式如下：

G(θ) = e± i m θ

这里m是任意正整数或0。

⒉泽尼克多项式的第二个性质是径向函数R ( ρ )（Radial Function）必须是ρ的n次多

项式，并且不包含幂次低于m次的ρ方项。

⒊第三个性质是当m为偶数时R(ρ)也为偶函数，m为奇数时，R(ρ)也为奇函数。

径向多项式R ( ρ )可以看作是雅可比多项式（Jacobi polynomials）的特例，记做。它们的正交和归一化性质可由如下式子表示：

上式中的δmn‘是克罗内克符号（Kronecker delta），即当n=n’时，δmn‘=1当n≠n'时，δmn‘=0。并且它具有归一化的性质：

在计算径向多项式时，为了方便起见，我们通常会将其分解成如下形式：

其中的次数为2(n-m)，由下式给出：

《光学原理》下册第9.2.1小节给出了上述径向函数的前几个m，n值的显函数形式。

通常我们会用实数形式的多项式（正弦和余弦函数）来代替复制数多项式，这样的话，波前像差函数W(ρ, θ)的泽尼克展开式就有如下形式：

这里W是平均波前差，An，B nm，C nm是多项式展开系数。由于”0级“项是个常数（或者叫平移项）1并且所有其他的泽尼克项在单位圆区域上的平均值是均为零，∴波前像差函数W的平均值就是这个“0级”项的系数A0，这样，上述公式就等价于：

对于一个回转对称的光学系统来说，物体位于子午面内，因而波前像差相对于yz面是对称的，也就是只有θ的偶函数（余弦项）项是非零项。对于一般情况，波前是不对称的，因而也就是同时包含两种三角函数形式。

泽尼克项

下面给出了48项泽尼克多项式，外加一项常数项。需要注意的是，读者并不需要严格按照下文所示的顺序排列这些泽尼克项，实际上在不同的应用和机构会采用不同的排列顺序。

表中的#0项是个常数或者说是平移项（piston term），这一项的系数也代表了平均光程差；而#1和#2项分别是x和y方向的倾斜项（tilt terms），#3代表了聚焦，因此，#1到#3项代表了波前的高斯或者近轴特性；#4和#5项代表了像散和离焦，#6和#7项代表彗差和倾斜，而#8项代表了3级像差和离焦，也就是说#4到#8项为3级相差项；同样地，#9到#15项代表了5级像差，而#16到#24项代表了7级像差，#25到#35项代表了9级像差，#36到#48项代表了11级像差。

2.1 极坐标形式的泽尼克多项式

2.1 笛卡尔坐标系下的泽尼克多项式

2.1 OSC泽尼克多项式

很多早期的用泽尼克多项式来对干涉图样做计算机分析的工作，是在上个世纪七十年代，由亚利桑那大学光学科学中心（OSC，Optical Sciences Center）的John Loomis进行的。OSC泽尼克多项式采用了n从1到5，以及n=6，m=0的项。n=m=0的常数项（piston term）也用来做干涉图样分析，但是这一项并不包含在泽尼克多项式中。因此，OSC泽尼克多项式包含36项，外加一项平移项（piston term）。这也是在光学设计软件OSLO和Code V中采用的形式。

泽尼克多项式的图形

n=1～6，m=0时的泽尼克多项式的三维图

2像散：2x2-ay2，其中a∈(-4, 4)

3彗差：2ρ2x+ax，其中a∈(-5, 3)；2ρ2x+ay，其中a∈(-4, 4)

4球差和离焦：ρ2(2ρ2+1.3a)，其中a∈(-5, 3)

5前36个泽尼克多项式

泽尼克多项式常用于干涉测试，而光学设计人员用的更多的则是赛德尔像差多项式。

泽尼克多项式和赛德尔像差

波前的初级和3级像差系数可以用泽尼克多项式来表示。我们将波前函数用泽尼克项的Z0～Z8这九项来表示成如下形式：

这些泽尼克项和像差的对应关系见表四。

表四：前9个泽尼克项和像差的对应关系

Z0平移（piston）

Z1

x轴倾斜

Z2y轴倾斜

Z3离焦

Z4像散@ 0° & 离焦

Z5像散@ 45° & 离

焦

Z6

彗差& x轴倾斜

Z7

彗差& y轴倾斜

Z8

球差& 离焦

继续将上述波前函数改写成如下形式：

W(ρ, θ ) = W11 cosθ + W20ρ2

+ W40ρ4 + W31ρ3 cosθ + W22ρ2 cos2θ