泽尼克多项式ZernikePolynomials,泽尼克系数

1、泽尼克多项式( Zernike Polynomials )，泽尼克系数什么是 Zernike Polynomials通常人们会使用幂级数展开式的形式来描述光学系统的像差。 由于泽尼克多项式和光学检测中观 测到的像差多项式的形式是一致的，因而它常常被用来描述波前特性(泽尼克， 1934)。但这并 不意味着泽尼克多项式就是用来拟合检测数据的最佳多项式形式。 在某些情况下，用泽尼克多项 式来描述波前数据具有很大的局限性。 比如说，当需要考虑空气扰动的时候， 泽尼克多项式几乎 没有什么价值。同样地，我们也无法找到一组合适的泽尼克多项式来描述单点金刚石车削加工(single point diamond

2、turning process )中的制造误差。为了准确地描述圆锥面光学元件(conicaloptical eleme nts)的对准误差，必须对泽尼克多项式进行修正。盲目地使用泽尼克多项式来表达 检测数据只会导致糟糕的结果。泽尼克多项式是由无穷数量的多项式完全集组成的，它有两个变量，p和e,它在单位圆内部是连续正交的。需要注意的是，泽尼克多项式仅在单位圆的内部连续区域是正交的， 通常在单位圆 内部的离散的坐标上是不具备正交性质的。泽尼克多项式具有三个和其他正交多项式集不一样的性质。1.泽尼克多项式z( p , e)可以被化解为径向坐标和角度坐标 e的函数，其形式如下：Z ( p , e )

3、= R (G ( pe ) ，这里，关于角度的函数 g( e是一个以2弧度为周期的连续函数，并且满足当坐标系旋转a角度之后，其形式不发生改变，也就是旋转不变性：e )G(其三角函数集形式如下:g( e)= e i m 9这里m是任意正整数或02. 泽尼克多项式的第二个性质是径向函数R ( p ) Radial Function)必须是 p的n次多项式，并且不包含幕次低于m次的p方项。3. 第三个性质是当 m为偶数时R(曲为偶函数，m为奇数时，R( p也为奇函数。径向多项式 R ( 可以看作是雅可比多项式( Jacobi polynomials)的特例，记做卜送蛀T。它们的正 交和归一化性质可由

4、如下式子表示：j农(P)琦詁上式中的 弗n，是克罗内克符号(Kronecker delta )，即当n=n 时，=1当n工n时,=0。并 且它具有归一化的性质：=1.在计算径向多项式时，为了方便起见，我们通常会将其分解成如下形式:其中的次数为2(n-m),由下式给出:$二0（2打一加 _ $）!s（n-s“一柳一）!光学原理下册第 921小节给岀了上述径向函数的前几个m, n值的显函数形式。通常我们会用实数形式的多项式(正弦和余弦函数)来代替复制数多项式，这样的话，波前像差函数 W(p ,0 )的泽尼克展开式就有如下形式：ir（p,e）=jF+Y40； (M+工 0； (P)PHI( B忻 z

5、 恥+G 汕讷)n=0级“项是个常数(或者叫平移波前像差函数 W的平均值这里W是平均波前差，An , Bnm, Cnm是多项式展开系数。由于”项)1并且所有其他的泽尼克项在单位圆区域上的平均值是均为零，就是这个“ 0级”项的系数A。，这样，上述公式就等价于：0( pH)=出十工出0： (P)十工P) P”迢加审加Y” sin恥)ji-1对于一个回转对称的光学系统来说，物体位于子午面内，因而波前像差相对于yz面是对称的,也就是只有0的偶函数(余弦项)项是非零项。对于一般情况，波前是不对称的，因而也就是同时包含两种三角函数形式。泽尼克项 下面给岀了 48项泽尼克多项式， 外加一项常数项。 需要注意

6、的是，读者并不需要严格按照下文所 示的顺序排列这些泽尼克项，实际上在不同的应用和机构会采用不同的排列顺序。表中的#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项代表了 1

7、1级像差。2.1极坐标形式的泽尼克多项式2.1笛卡尔坐标系下的泽尼克多项式taj2.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中采用的

8、形式。泽尼克多项式的图形n=16，m=0时的泽尼克多项式的三维图编辑版word像散：2x2-ay2,其中 a (-4, 4)彗差：2 px+ax，其中 a (5, 3)；2 px+ay，其中 a (-4, 4)球差和离焦：P 2(2 p1.3a)，其中a ( 3)前36个泽尼克多项式Illi泽尼克多项式常用于干涉测试，而光学设计人员用的更多的则是赛德尔像差多项式。泽尼克多项式和赛德尔像差ZoZ8这波前的初级和3级像差系数可以用泽尼克多项式来表示。我们将波前函数用泽尼克项的 九项来表示成如下形式:W(p7 ff = Zo + Npcos护 + Z2psin01 + Z(2p2 I)+ Z3 co

9、s2fl;十 Zsp2 sin 20*+ Z3p2 - 2)p cos + Z7(3p2 一 2)psin &+ S 亠 +这些泽尼克项和像差的对应关系见表四表四：前9个泽尼克项和像差的对应关系Zo 平移(piston )Zlx轴倾斜Z2y轴倾斜Z3离焦Z4 像散 0& 离焦Z5 像散 45 & 离焦八、Z6 彗差 & x 轴倾斜Z7彗差 & y 轴倾斜Z8球差 & 离焦继续将上述波前函数改写成如下形式：W( p ,0 )i=cWS 0 +Wp+ W40 p 4 + W p cos 0 + Wp coS 0由于这些泽尼克项中与视场无关， 它们并非真正的赛德尔像差。 用干涉测试的方法智能得到单个

10、 视场点的波前数据。 这使得场曲看上去像离焦， 而畸变看上去像倾斜。 因此， 要得到赛德尔像差， 就必须测量一定数量的视场点。我们可以按照初级和 3级像差的形式继续改写上述波前函数，也就是合并同类项，并用波前相差 系数做等价替换，结果如下：w(p , e + Z8移(Piston)+ ( Z i - 2Z6)p cos e +-(Z Z7) p sin 倾斜(Tilt)+ (2Z 3 - 6 Z8 + Z 4COS2 e Z+sin2 e 2)p离焦 + 像散(focus +astigmatism)+3(Z 6cos e+/Zn e3 ) p彗差(coma)+6Z 8 p球差(Spherica

11、l)对上式中做如下恒等变换:a cos aty/az 4- h2cosi Isn编辑版word便可得到如下所示的视场无关的波前相差系数:W(P1 ff) = ZQ-Z3-h za+ Py(ZI-2ZJa + Z2-2Z?)2+p2(2Z.-6Z8+土 2p2s/zJ + Z5 cos3+ 3p./zTTz|+ 6pZgLpistontiltfocusastigmatismcomaspberical表五列岀了上述视场无关像差多项式的度量（Magnitude）,符号和角度（Angle）。注意离焦项的符号选择原则是使得其系数的数值最小，像散符号的选取则相反。表五：用泽尼克多系数表示的3级像差Term

12、DescriptionVaiiiitLdeAikglerihg -2Z屏十冋-2zJ%(Qtlh2鬲-6Z* + y/Zl + Zl cl】o、en to rniiHuuzt absoluH -aliK oF in呻mu吐知唯）atiniRtkrn土苑 2TR;sign opposite th机 cboscn in Icxus icrniaim%注意：在计算角度值（表中的Angle）的时候，如果分母 0 , angle angle + 180编辑版word波前差的RMS值和P-V值如果波前像差可以用 3级像差多项式来表示的话，我们就可以很方便地用每种3级像差所反映岀来的波数来描述波前像差。这种

13、方法在仅有一种3级像差存在的情况下尤为方便。而对于更加复杂一些的波前像差，我们也可以很方便地得岀其波前的P-V （ peak-to-valley）值（有时候也叫做 P-P值），其数值就是实际波前和理想波前之间的最大偏差，包括正向和负向两个方向。比如，如果正向的波前最大偏差10.+0.2波长，负向最大的波前偏差为-0.1波长，那么这个波前差的P-V值就是0.3波长。虽然用P-V值来描述一个波前质量很简单方便，但缺容易引起误导。一个波前差的P-V值比较大的光学系统的实际性能有可能会优于一个波前差的P-V值相对比较小的系统。这样，用RMS波前差来描述波前质量会更有意义一些。e ）是相对于理想球面波的偏差，其下述方程定义了圆形光阑下的波前差b及其方差（oA W （单位通常为波长数。AW 是平均波前光程差（OPD ）62)1=亠7T卩如血-（丽A0) pdpdO如果将波前像差用泽尼克多项式的形式表示的话,波前差的方差就可以很简单地通过泽尼克多项式的正