基于勒让德多项式各项系数拟合的柱面镜面形检测研究

杨醒++王占山

摘要：

超薄玻璃热成形技术是目前国际上研制X射线掠入射天文成像望远镜的重要技术之一，超薄玻璃圆柱面镜片面形的高精度检测是实现超薄玻璃热成形的关键。围绕我国高能X射线掠入射天文望远镜研制对高面形精度的超薄玻璃柱面镜的需求，采用勒让德多项式对柱面镜面形的测量数据进行拟合，利用勒让德多项式各项系数与特定的几何像差的对应性，通过拟合勒让德多项式各项的系数，获得了待测柱面镜面形的特定几何像差，为全面评价柱面镜面形质量提供了参考。

关键词：

柱面镜； 勒让德多项式； 最小二乘法； 像差

中图分类号： TN 253 文献标志码： A doi： 10.3969/j.issn.10055630.2016.01.003

Figure error test of the cylindrical lens surface based on the fitting of

coefficients of Legendre polynomials

YANG Xing， WANG Zhanshan

（Key Laboratory of Advanced Microstructured Materials of MOE，School of

Physics Science and Engineering， Tongji University， Shanghai 200092， China）

Abstract：

The nominally cylindrical glass segments are forced to a conical form during the assembly process of an Xray imaging telescope as the ultrathin glass thermal slumping technique being fully successful. The surface should be measured in order to make highquality glass with cylindrical surface. Each coefficient of the Legendre polynomials corresponds to the specific geometrical aberration. In this paper， figure error of cylindrical lens is obtained based on the fitting coefficients of Legendre polynomials， which provides an exact method for cylindrical lens surface measurement.

Keywords： cylindrical lens； Legendre polynomials； least square method； aberrations

引 言

X射线天文观测是人类认识和了解宇宙，尤其是宇宙起源、超新星遗迹、黑体的重要手段。自20世纪60年代，人们开始运用WolterI型望远镜对各种天体进行观测。WolterI型望远镜是由共焦的掠入射抛物面和双曲面反射镜组成，为了获得大的有效面积，采用了嵌套结构。从1999年起，先后发展了基于常规光学方法制作镜片的Chandra望远镜[1]，基于Ni电镀方法制作镜片的Newton望远镜[2]，基于铝薄片成型分离方法制作镜片的ASCA[3]和Suzaku[4]X射线望远镜。铝薄片镜片的使用虽然大大降低了望远镜的造价，但是其自身性质限制了这种方法制作望远镜分辨率的进一步提升。随着技术的发展，人们可以采用超薄玻璃制作望远镜的反射镜片替代原有的铝片。2012年，美国发射的NuSTAR卫星[5]是采用超薄玻璃镜片制作望远镜的典范，其分辨率是58″。这种方法制作的反射镜片是先制作出圆柱面形状，在装配时将其压成圆锥面以达到设计要求。因此，在进行这种望远镜研制过程中，如何精确测量圆柱面面形是改进工艺、提高镜片制作质量的关键。

圆柱面面形检测方法有很多，常用的有三坐标测量法和采用计算全息图（CGH）的干涉测量法。一般的三坐标测量是用触针式探头接触待测圆柱面表面，通过探头位置读出并记录待测点的三坐标值，得到被测面的面形。当对厚度只有0.21 mm的超薄玻璃镜片测量时，这种方法会引入较大误差，不能满足面形测量精度的要求[6]。采用计算全息图的方法是用计算全息图产生柱面波，该柱波面与待测圆柱面表面反射的柱面波进行干涉，根据干涉结果计算获得待测圆柱面面形的偏差。

在分析干涉测量数据时，一般干涉仪会给出面形的均方根值（RMS）、最大差别值（PV）和等高线图等评价指标，更好的结果是用泽尼克（Zernike）多项式系数给出的面形测试结果，但其不适用于圆柱面的检测。在非圆域光学面形测量中，Paul研究了勒让德多项式应用于近柱面面形误差的表征，给出了一组包含勒让德多项式的方程，并成功地评价了掠入射光学系统的波前误差和成像误差[7]。Gao等将勒让德多项式用于方形域面形的重建，获得了相关的面形数据[8]。本文基于勒让德多项式的特点研究了其用于圆柱面面形检测，给出了类似于泽尼克多项式系数的圆柱面面形评价结果。

根据式（9），对十块不同柱面镜面形进行了拟合，求得的RMS和PV值如表2所示，其中柱面镜的面形测量数据采用干涉检测的计算全息方法获得[16]。

从表2可知，镜片3的面形RMS值和PV值最小，说明在十块镜片中它的面形质量最好。但仅从RMS和PV这两个评价指标来判定面形的优劣还不够明确。比如，镜片4的RMS值小于镜片5，但是它的PV值却大于镜片5。为了更全面评价柱面镜的面形精度，我们根据拟合出的勒让德多项式各项的系数所对应的像差给出另一种评价面形质量的方法。表3给出了镜片4和镜片5的勒让德多项式各项的拟合系数。

从表3可以得到如下的结论：

（1） a0=0，说明两块镜片的面形中都不存在平移像差；

（2） 勒让德多项式各项拟合系数除第6项外，镜片5对应系数的绝对值均大于镜片4，说明倾斜、离焦以及y轴初级彗差对镜片5面形的影响大于镜片4；

（3） 勒让德多项式各项拟合系数的第6项，镜片4对应系数的绝对值大于镜片5，说明x轴初级彗差在镜片4面形中的影响较大；

（4） 综上，镜片4的整体面形优于镜片5。

镜片4和镜片5的剩余表面误差图如图2和图3所示。

图2 镜片4剩余表面误差

Fig.2 Residual error of the surface of lens four

图3 镜片5剩余表面误差

Fig.3 Residual error of the surface of lens five

3 结 论

本文利用勒让德多项式拟合了被测的十块柱面镜面形测量结果，给出了勒让德多项式与初级像差之间的关系，获得了勒让德多项式系数的拟合值，确立了勒让德多项式拟合系数与初级像差在柱面镜面形测量中的影响信息。该方法比常用的评价指标RMS值和PV值更能准确地判定面形质量的优劣。

参考文献：

[1] WEISSKOPF M C.Chandra xray optics[J].Optical Engineering，2012，51（1）：011013.

[2] LUMB D H，SCHARTEL N，JANSEN F A.Xray multimirror mission （XMMNewton） observatory[J].Optical Engineering，2012，51（1）：011009.

[3] SERLEMITSOS P J，JALOTA L，SOONG Y，et al.The Xray telescope on board ASCA[J].PASJ：Publications of the Astronomical Society of Japan，1995，47（1）：105114.

[4] SERLEMITSOS P J，SOONG Y，CHAN K W，et al.The Xray telescope onboard Suzaku[J].PASJ：Publications of the Astronomical Society of Japan，2007，59（S1）：921.

[5] CRAIG W W，AN H J，BLAEDEL K L，et al.Fabrication of the NuSTAR flight optics[C]//O′DELL S L，PARESCHI G.Proceedings of SPIE8147，optics for EUV，Xray，and gammaray astronomy V.San Diego，California：SPIE，2011，8147：81470H181470H14.

[6] LEHAN J P，HADJIMICHAEL T，SKOCIK C.Testing of the mirrors for the constellationx spectroscopy xray telescope with a refractive null[C]//O′DELL S L，PARESCHI G.Proceedings of SPIE6688，optics for EUV，XRay，and gammaray astronomy III.San Diego，CA：SPIE，2007，6688：66881916688198.

[7] PAUL G.Set of orthonormal surface error descriptors for nearcylindrical optics[J].Optical Engineering，1984，23（4）：384390.

[8] GAO B，LI Q，HE Y H，et al.Absolute measurement of square flats with Legendre polynomial fitting method[C]//CZARSKE J，ZHANG S L，SAMPSON D， et al.Proceedings of SPIE9297， international symposium on optoelectronic technology and application 2014：laser and optical measurement technology； and fiber optic sensors.Beijing：SPIE，2014，9297：92970E192970E7.

[9] DAI G M，MAHAJAN V N.Orthonormal polynomials in wavefront analysis：error analysis[J].Applied Optics，2008，47（19）：34333445.

[10] SUN W Q，CHEN L，TUYA W L，et al.Analysis of the impacts of horizontal translation and scaling on wavefront approximation coefficients with rectangular pupils for Chebyshev and Legendre polynomials[J].Journal of the Optical Society of America A，2013，30（12）：25392546.

[11] MAHAJAN V N.Orthonormal aberration polynomials for anamorphic optical imaging systems with rectangular pupils[J].Applied Optics，2010，49（36）：69246929.

[12] BORN M，WOLF E.Principles of Optics：electromagnetic theory of propagation，interference and diffraction of light[M].6th ed.New York：Pergamon，1980.

[13] YUAN S，SASIAN J.Aberrations of anamorphic optical systems.II.Primary aberration theory for cylindrical anamorphic systems[J].Applied Optics，2009，48（15）：28362841.

[14] MAHAJAN V N，DAI G M.Orthonormal polynomials in wavefront analysis：analytical solution[J].Journal of the Optical Society of America A，2007，24（9）：29943016.

[15] RAYCES J L.Leastsquares fitting of orthogonal polynomials to the waveaberration function[J].Applied Optics，1992，31（13）：22232228.

[16] 马顺翔，王占山.用于柱面镜面形检测的CGH的设计与误差分析[J].光学仪器，2015，37（5）：419425.

（编辑：刘铁英）