用于大口径非球面的波前功率谱密度检测

　　引 言

　　在航天、天文、军事等领域中，大口径非球面光学元件在各种精密光学系统中的应用日益广泛，随着军事光学系统需求的增加和特殊要求，大口径高精度非球面光学元件的制造急需与之精度相适应的评价指标、检验方法和测试手段得到解决，因此如何实现对大口径非球面全频段面形质量检验与评价具有重要的研究意义[1]。

　　在现有的镜面质量评价方法中，传统干涉检验给出的波前峰谷值和均方根值反映的是镜面全口径的空间频率质量，受 CCD 分辨率影响，主镜镜面中高频差往往被低频误差所淹没。波前功率谱密度(PowerSpectral Density-PSD)是一种描述波前信息的新方法[2]，它与以往各种表面评价参数相比最大的优点在于，PSD 可以综合评价表面起伏的垂直量和水平量，并体现两者间的关系。波前PSD 目前还仅仅用于惯性约束聚变系统光学元件的中频差分析，对大口径非球面的中频差测量则未见相关报道，因此尚有更进一步研究的需要。

　　1 系统组成及工作原理

　　波前PSD 的检测，实质上是对光学元件透射或反射波前畸变进行检测，可以使用激光数字干涉仪或哈特曼检测装置作为检测仪器。图1 是干涉仪和哈特曼对非球面实施检验的光路布局简图。用干涉仪作为检测仪器时，在右边光源外需要球面反射镜自准;用哈特曼作为检测仪器时，需要提供激光光源。该装置布局展示了测试大口径非球面中高频误差的可行思路。利用局部测量方法是考虑到非球面加工的对称性，并有利于提高分辨率。例如对于 Φ100mm 口径，干涉仪 CCD 像素为 1024×1024，其采样周期可达到 0.09mm，将有效提升仪器可达到的高频限。每种面形检测系统都具有一定的测量带宽，这个带宽也就决定了检测系统的精度和测量范围。J MBennett 和E L Church 等人对此进行广泛的研究和讨论[3-4]。如按照E L Church 方法，对 PSD 有效频率宽度做保守估计，则有效频率宽度为[5]

　　2 数据处理

　　光学面形轮廓函数 z(x)的一维傅里叶变换如下[6-7]

　　PSD 的计算流程如图3 所示。首先读入检测仪器所采集到的原始数据，对原始数据做波前预处理，从频谱分析的角度，去除测试因数所带来的误差影响，由于使用计算机处理数据时，只能处理有限的离散形式的数据，因此“频率泄漏”是不可避免的，需要通过加窗处理来减小误差，修正原始数据，然后作 FFT，采用适当的PSD 算法作初步的PSD 估计。由于任何系统本身都不可能是理想的，其系统误差将对空间频率分量的测量产生影响，因此必须用系统传递函数对PSD 的估计值做修正，以便获得真实的PSD 结果。系统的传递函数计算公式为