CN106404354A - 非球面补偿器透射波前方程的测量装置和方法 - Google Patents
非球面补偿器透射波前方程的测量装置和方法 Download PDFInfo
- Publication number
- CN106404354A CN106404354A CN201610886907.1A CN201610886907A CN106404354A CN 106404354 A CN106404354 A CN 106404354A CN 201610886907 A CN201610886907 A CN 201610886907A CN 106404354 A CN106404354 A CN 106404354A
- Authority
- CN
- China
- Prior art keywords
- aspheric surface
- transmission wavefront
- sphere
- compensator
- surface compensator
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0242—Testing optical properties by measuring geometrical properties or aberrations
- G01M11/0271—Testing optical properties by measuring geometrical properties or aberrations by using interferometric methods
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
本发明提供了一种非球面补偿器透射波前方程的测量装置和方法,测量装置包括移相干涉仪、测距干涉仪和参考球面组件,非球面补偿器置于移相干涉仪和参考球面组件之间,且三者同轴设置,参考球面组件沿光轴方向移动,参考球面组件与非球面补偿器形成环形面形图,测量参考所述球面组件的相对位移,利用环形面形图和相对位移计算非球面补偿器透射波前方程。本发明可实现对非球面补偿器透射波前进行测量,从而在进行非球面加工之前确保非球面补偿器正确无误,避免非球面补偿器加工错误对最终光学***造成严重影响;测量方法的测量精度高,测量装置成本低,便于操作。
Description
技术领域
本发明涉及光学检测技术领域,尤其涉及一种非球面补偿器透射波前方程的测量装置和方法。
背景技术
在光学设计中,单个球面可以供优化的自由度只有曲率半径;而非球面除了顶点曲率半径之外,还有二次曲面常数和高阶项系数。由于非球面比球面拥有更多设计的自由度,所以在很多光学***中,都普遍采用非球面元件来减小***的复杂度,并提高***的成像质量。
当前,对于高精度非球面的检测,大都采用零位补偿法,即采用补偿镜或者计算全息图(Computer-Generated Hologram,CGH)作为补偿器,将干涉仪发出的球面波转化为与非球面相匹配的非球面波,从而实现零位检测。在检测过程中,非球面补偿器是限制检测精度的关键性因素;如果补偿器存在问题(如透镜的曲率半径加工出错或者透镜间隔装配错误),则直接导致加工出的非球面偏离设计值,从而给光学***造成严重的影响。为此,当非球面补偿器加工完成之后,非常有必要对其透射波前进行测量,从而保证加工出的非球面补偿器满足设计精度的要求。
因此,研发一种非球面补偿器透射波前方程的测量装置和方法,成为人们亟待解决的问题。
发明内容
有鉴于此,本发明的主要目的在于提供一种非球面补偿器透射波前方程的测量装置和方法,以实现对非球面补偿器的校验。
本发明提供了一种非球面补偿器透射波前方程的测量装置,包括:移相干涉仪、测距干涉仪和参考球面组件;其中,作为被测对象的非球面补偿器置于所述移相干涉仪和参考球面组件之间,且所述非球面补偿器、移相干涉仪和参考球面组件同轴设置,所述参考球面组件可沿光轴方向移动,所述移相干涉仪发出的光经过所述非球面补偿器,并由所述参考球面组件反射,反射光经过所述非球面补偿器后入射至所述移相干涉仪,所述参考球面组件沿光轴方向移动,在所述移相干涉仪内形成参考球面组件与非球面补偿器之间的环形面形图,所述测距干涉仪测量参考所述球面组件的相对位移,所述移相干涉仪利用所述环形面形图和相对位移计算非球面补偿器透射波前方程。
本发明还提供了一种利用上述非球面补偿器透射波前方程的测量装置的测量方法,所述测量方法包括:步骤A:移动参考球面,利用移相干涉仪测量非球面补偿器透射波前与参考球面形成的环形面形图,利用测距干涉仪测量所述参考球面的与所述环形面形图对应的相对位移;步骤B:依据所述移相干涉仪测量的环形面形图,计算所述环形面形图的“零条纹”半径;步骤C:基于所述参考球面与非球面补偿器透射波前的位置关系,得到非球面补偿器透射波前矢高的导数,进而得到非球面补偿器透射波前方程的第一表达式;步骤D:设定非球面补偿器透射波前方程的导数表达式,得到非球面补偿器透射波前方程第二表达式,进而得到非球面补偿器透射波前方程的迭代表达式;以及步骤E:利用所述参考球面的与环形面形图对应的相对位移和所述环形面形图的“零条纹”半径,对非球面补偿器透射波前方程迭代表达式进行迭代运算,得到非球面补偿器透射波前方程。
从上述技术方案可以看出,本发明的非球面补偿器透射波前方程的测量装置和方法具有以下有益效果:可实现对非球面补偿器透射波前进行测量,从而在进行非球面加工之前确保非球面补偿器正确无误,避免非球面补偿器加工错误对最终光学***造成严重影响;该测量方法的测量精度高,测量装置成本低,便于操作。
附图说明
图1为本发明实施例的非球面补偿器透射波前方程的测量装置示意图;
图2为本发明实施例的非球面补偿器透射波前方程的测量方法流程图;
图3为当干涉条纹与中心不连通时,满足设定阈值的“零条纹”位置分布图;
图4为利用图3中“零条纹”位置计算出的最佳拟合圆;
图5为当干涉条纹与中心连通时,仿真生成的环形干涉图;
图6为图5对应的环形面形图;
图7为仿真生成的环形干涉图;
图8为非球面补偿器透射波前传播示意图;
图9为非球面补偿器光学结构图;
图10为采用本发明方法计算出的非球面补偿器透射波前的矢高与设计值的偏差。
【符号说明】
S1-稳频激光器;S2-准直扩束***;S3-分光棱镜;S4-汇聚镜组;S5-移相参考面;S6-测距干涉仪;S7-参考球面;S8-夹持机构;S9-多维调整台;S10-成像镜;S11-探测器;S12-图像采集单元;S13-六维调整控制单元。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
参见图1,图1示出了本发明实施例的非球面补偿器透射波前方程的测量装置,该测量装置包括:移相干涉仪、测距干涉仪和参考球面组件;其中,
移相干涉仪包括:稳频激光器S1、准直扩束***S2、分光棱镜S3、汇聚镜组S4、移相参考面S5、成像镜S10、探测器S11和图像采集单元S12;
参考球面组件包括:参考球面S7、夹持机构S8、多维调整台S9以及六维调整控制单元S13,参考球面S7由夹持机构S8固定于多维调整台S9,六维调整控制单元S13可以控制多维调整台S9移动。
非球面补偿器S13置于移相干涉仪和参考球面组件之间,且非球面补偿器S13、移相干涉仪和参考球面S7同轴设置。
稳频激光器S1发出的光经准直扩束***S2后进入分光棱镜S3,分光棱镜S3将入射光分为两束,第一束光沿移相干涉仪光轴方向传播,第二束光沿垂直于移相干涉仪光轴方向传播,其中,第一束光经汇聚镜组S4入射至非球面补偿器S13,非球面补偿器S13的透射波前入射至参考球面S7,经参考球面S7反射的光依次经过非球面补偿器S13、汇聚镜组S4,并由分光棱镜S3反射至成像镜S10和探测器S11;第二束光入射至移相参考面S5,经移相参考面S5反射的光由分光棱镜S3透射至成像镜S10和探测器S11。
其中,六维调整控制单元S13控制多维调整台S9沿移相干涉仪光轴方向移动,从而带动参考球面S7由初始位置沿移相干涉仪光轴方向移动,经参考球面S7反射的光与经移相参考面S5反射的光发生干涉,参考球面S7与非球面补偿器S13在探测器S11上形成一系列环形面形图,与此同时,测距干涉仪S6测量参考球面S7相对于初始位置的移动距离,作为与一系列环形面形图相对应的相对位移,图像采集单元S12利用得到的环形面形图和相对位移计算出非球面补偿器透射波前方程。
本发明提供的非球面补偿器透射波前方程的测量装置,可实现对非球面补偿器透射波前的测量,从而在进行非球面加工之前得到非球面补偿器的参数,该测量装置成本低、结构简单,通过移动参考球面即可得到一系列环形面形图,从而得到非球面补偿器透射波前方程,操作方便灵活。
本发明的另一实施例还提供了一种测量方法,参见图2,该测量方法利用上述非球面补偿器透射波前方程的测量装置,对非球面补偿器透射波前方程进行测量,该测量方法包括:
步骤A:沿移相干涉仪光轴方向移动参考球面,利用移相干涉仪测量非球面补偿器透射波前与参考球面形成的环形面形图,利用测距干涉仪测量参考球面的与环形面形图对应的相对位移。
执行步骤A时,随着参考球面沿移相干涉仪光轴方向移动,非球面补偿器透射波前与参考球面会形成一系列环形面形图,与此同时,测距干涉仪测量参考球面相对于初始位置的移动距离,该移动距离作为参考球面的相对位移,并记录环形面形图相对应的相对位移,即对于一系列环形面形图中的第i个环形面形图对应一参考球面的相对位移si。
步骤B:依据所述移相干涉仪测量的环形面形图,计算环形面形图的“零条纹”半径。
由步骤A可知,随着参考球面的移动会形成一系列环形面形图,对于这些环形面形图来说,部分环形面形图的干涉条纹与环形面形图中心相连通,部分环形面形图的干涉条纹不与环形面形图中心连通,但这两种环形面形图的“零条纹”位置均满足:
其中,F(x,y)表示环形面形图,r表示环形面形图的“零条纹”半径。
对于步骤B来说,分以下两种情况计算环形面形图的“零条纹”半径:
情况一:当环形面形图的干涉条纹与环形面形图中心不连通时,为了便于计算,将公式(1)改写为:
其中,Fi(x,y)表示一系列环形面形图中的第i个环形面形图,ri表示第i个环形面形图的“零条纹”半径,θ表示极坐标下的角度坐标,i=1,2,…N。
对公式(2)采用数值微分法,设置合适的阈值,得到环形面形图的“零条纹”位置分布图,如图3所示,然后采用最小二乘法进行拟合,即可得到环形面形图的“零条纹”半径与中心,如图4所示。
情况二:当环形面形图的干涉条纹与环形面形图中心连通时,如图5、图6所示,环形面形图主要由离焦和球差项组成,即:
Fi(x,y)=a1 i(2ri 2-1)+a2 i(6ri 4-6ri+1) (3)
其中,Fi(x,y)表示一系列环形面形图中的第i个环形面形图,ri表示第i个环形面形图的“零条纹”半径,a1 i和a2 i分别表示第i个环形面形图的离焦和球差系数,其可以通过对环形面形图进行Zernike拟合获得。
然后将公式(3)带入公式(1),可得到环形面形图的“零条纹”半径:
步骤C:基于参考球面与非球面补偿器透射波前的位置关系,得到非球面补偿器透射波前矢高的导数,进而得到非球面补偿器透射波前方程的第一表达式。
在步骤C中,当非球面补偿器透射波前与参考球面径向任一位置的法线垂直时,会出现如图7所示的环形干涉图,参考球面与非球面补偿器透射波前的位置关系如图8所示。其中,P点为第i个环形面形图“零条纹”位置(Null Position)对应的参考球面上的点,F点为光线(参考球面上过P点的法线)与移相干涉仪光轴的交点,N点为光线(参考球面上过P点的法线)与非球面补偿器透射波前的交点,P点与F点的坐标可以分别表示为(ri,Zi)和(O,si),其中,ri为第i个环形面形图的“零条纹”半径,Zi为参考球面矢高,si为参考球面的与第i个环形面形图对应的相对位移,P点与N点的距离PN为非球面补偿器透射波前的传播距离,其表示为di,则N点坐标(ρi,z(ρi))为:
其中,θi表示P点和F点的连线与移相干涉仪光轴的夹角,即光线与光轴的夹角,ρi为光线与非球面补偿器透射波前的交点,即非球面补偿器透射波前径坐标,z(ρi)表示非球面补偿器透射波前的矢高,z′(ρi)表示非球面补偿器透射波前矢高的导数。
为了便于计算,将公式(6)改写为:
由公式(7)得到非球面补偿器透射波前方程第一表达式:
步骤D:设定非球面补偿器透射波前方程的导数表达式,得到非球面补偿器透射波前方程第二表达式,进而得到非球面补偿器透射波前方程的迭代表达式。
在步骤D中,为了提高测量精度,非球面补偿器透射波前方程的导数表达式为:
其中,c为最佳拟合球曲率,ρ为非球面补偿器透射波前径坐标,ρmax为非球面补偿器透射波前径坐标最大值,M由人工设定,其可以根据实际需要进行整数取值,优选M取6,bn表示多项式{Pn(x)}的系数,且多项式{Pn(x)}满足:
其中,当m=n时,δmn=1;当m≠n时,δmn=0。
根据公式(9),得到非球面补偿器透射波前方程第二表达式:
其中,Tn(x)=∫Pn(x)dx,n≥0。
令非球面补偿器透射波前方程第一表达式与第二表达式相等,得到非球面补偿器透射波前方程迭代表达式:
步骤E:利用参考球面的与环形面形图对应的相对位移(由步骤A得到)和环形面形图的“零条纹”半径(由步骤B得到),对非球面补偿器透射波前方程迭代表达式进行迭代运算,得到非球面补偿器透射波前方程。
在步骤E中,包括:
子步骤E1:令公式(5)的非球面补偿器透射波前的传播距离di等于参考球面的与环形面形图对应的相对位移si。
子步骤E2:由环形面形图的“零条纹”半径ri,利用参考球面方程可以得到光线与光轴的夹角θi,以及参考球面矢高Zi,将ri、θi、Zi、di代入公式(5),得到光线与非球面补偿器透射波前的交点、以及非球面补偿器透射波前矢高{ρi}、{z(ρi)}。
子步骤E3:设置最佳拟合球曲率c的初始值,将参考球面的与环形面形图对应的相对位移si、光线与非球面补偿器透射波前的交点ρi、以及非球面补偿器透射波前矢高z(ρi)代入公式(12)进行第一迭代运算,得到非球面补偿器透射波前方程第二表达式,该第二表达式的形式参见公式(11)。
子步骤E4:将光线与非球面补偿器透射波前的交点ρi代入子步骤E3得到的非球面补偿器透射波前方程第二表达式,重新计算非球面补偿器透射波前矢高z′(ρi),并利用公式(13)重新计算非球面补偿器透射波前的传播距离:
其中,di′为重新计算的非球面补偿器透射波前的传播距离,z′(ρi)为重新计算的非球面补偿器透射波前矢高。
子步骤E5:计算两次迭代之间非球面补偿器透射波前的传播距离差值di′-di,将传播距离差值与传感距离阈值比较;
若传播距离差值大于传播距离阈值,将子步骤E4重新计算的非球面补偿器透射波前的传播距离di′作为子步骤E2的di,重复子步骤E2至子步骤E5;
否则,将子步骤E3得到的非球面补偿器透射波前方程第二表达式,作为非球面补偿器透射波前方程。
其中,子步骤E3的第一迭代运算具体包括:
子分步骤E3a:将初始值赋予最佳拟合球曲率c。
子分步骤E3b:将最佳拟合球曲率c、参考球面的与环形面形图对应的相对位移si、光线与非球面补偿器透射波前的交点ρi、以及非球面补偿器透射波前矢高z(ρi)代入公式(12),采用最小二乘法拟合出多项式{Pn(x)}的系数{bn};最佳拟合球曲率c的初始值可由人工设定。
子分步骤E3c:得到多项式{Pn(x)}的系数{bn}后,可以得到非球面补偿器透射波前方程第二表达式,该第二表达式的形式参见公式(11)。
子分步骤E3d:将光线与非球面补偿器透射波前的交点ρi、以及非球面补偿器透射波前矢高z(ρi)代入子分步骤E3c得到的非球面补偿器透射波前方程第二表达式,重新计算最佳拟合球曲率,将重新计算的最佳拟合球曲率与最佳拟合球曲率的初始值相减,获得最佳拟合球曲率差值。
子分步骤E3e:将最佳拟合球曲率差值与最佳拟合球曲率阈值比较,若最佳拟合球曲率差值大于最佳拟合球曲率阈值,将子分步骤E3d重新计算的最佳拟合球曲率作为子分步骤E3b最佳拟合球曲率,重复执行子分步骤E3b至子分步骤E3e;否则,将子分步骤E3c得到的非球面补偿器透射波前方程第二表达式作为迭代运算结果输出。
通过上述方法针对一高次偶次非球面补偿器透射波前的方程进行测量,此非球面补偿器所补偿的非球面参数如表1所示,非球面补偿器本身的结构如图9所示。
表1
利用步骤A至步骤E完成了对非球面补偿器透射波前方程的测量,利用最终获得非球面透射波前方程计算出矢高,此矢高与设计值的偏差如图10所示。
由此可见,本发明的非球面补偿器透射波前方程的测量方法,可实现对非球面补偿器透射波前进行测量,从而在进行非球面加工之前确保非球面补偿器正确无误,避免非球面补偿器加工错误对最终光学***造严重的影响;该测量方法具有测量方法简单,成本低,便于操作、精度高等优点。
至此,已经结合附图对本实施例进行了详细描述。依据以上描述,本领域技术人员应当对本发明的非球面补偿器透射波前方程的测量装置和方法有了清楚的认识。
需要说明的是,在附图或说明书正文中,未绘示或描述的实现方式,均为所属技术领域中普通技术人员所知的形式,并未进行详细说明。此外,上述对各元件的定义并不仅限于实施例中提到的各种具体结构和形状,本领域普通技术人员可对其进行简单地更改或替换。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (11)
1.一种非球面补偿器透射波前方程的测量装置,其特征在于,包括:移相干涉仪、测距干涉仪和参考球面组件;其中,
作为被测对象的非球面补偿器置于所述移相干涉仪和参考球面组件之间,且所述非球面补偿器、移相干涉仪和参考球面组件同轴设置,所述参考球面组件可沿光轴方向移动;
所述移相干涉仪发出的光经过所述非球面补偿器,并由所述参考球面组件反射,反射光经过所述非球面补偿器后入射至所述移相干涉仪,所述参考球面组件沿光轴方向移动,在所述移相干涉仪内形成参考球面组件与非球面补偿器之间的环形面形图,所述测距干涉仪测量参考所述球面组件的相对位移,所述移相干涉仪利用所述环形面形图和相对位移计算非球面补偿器透射波前方程。
2.如权利要求1所述的测量装置,其特征在于,所述参考球面组件包括:参考球面、夹持机构、多维调整台以及六维调整控制单元,所述参考球面由所述夹持机构固定于所述多维调整台,所述六维调整控制单元控制所述多维调整台的移动,并带动所述参考球面移动。
3.如权利要求2所述的测量装置,其特征在于,所述移相干涉仪包括:稳频激光器、准直扩束***、分光棱镜、汇聚镜组、移相参考面、成像镜、探测器和图像采集单元;
所述稳频激光器发出的光经所述准直扩束***后进入所述分光棱镜,所述分光棱镜分出的第一束光经所述汇聚镜组入射至所述非球面补偿器,所述非球面补偿器的透射波前入射至所述参考球面,经所述参考球面反射的光依次经过所述非球面补偿器、所述汇聚镜组,并由所述分光棱镜反射至所述成像镜和探测器;分光棱镜分出的第二束光入射至所述移相参考面,经所述移相参考面反射的光由所述分光棱镜透射至所述成像镜和探测器。
4.如权利要求3所述的测量装置,其特征在于,所述六维调整控制单元控制所述多维调整台沿光轴方向移动,带动所述参考球面由初始位置沿光轴方向移动,所述参考球面与非球面补偿器在所述探测器上形成一系列环形面形图,所述测距干涉仪测量参考球面相对于初始位置的移动距离,作为与一系列环形面形图相对应的相对位移,所述图像采集单元利用环形面形图和相对位移计算出非球面补偿器透射波前方程。
5.一种利用权利要求1所述的非球面补偿器透射波前方程的测量装置的测量方法,其特征在于,所述测量方法包括:
步骤A:移动参考球面,利用移相干涉仪测量非球面补偿器透射波前与参考球面形成的环形面形图,利用测距干涉仪测量所述参考球面的与所述环形面形图对应的相对位移;
步骤B:依据所述移相干涉仪测量的环形面形图,计算所述环形面形图的“零条纹”半径;
步骤C:基于所述参考球面与非球面补偿器透射波前的位置关系,得到非球面补偿器透射波前矢高的导数,进而得到非球面补偿器透射波前方程的第一表达式;
步骤D:设定非球面补偿器透射波前方程的导数表达式,得到非球面补偿器透射波前方程第二表达式,进而得到非球面补偿器透射波前方程的迭代表达式;以及
步骤E:利用所述参考球面的与环形面形图对应的相对位移和所述环形面形图的“零条纹”半径,对非球面补偿器透射波前方程迭代表达式进行迭代运算,得到非球面补偿器透射波前方程。
6.如权利要求5所述的测量方法,其特征在于,对于干涉条纹与中心不连通的环形面形图,所述步骤B包括:
子步骤B1:得到环形面形图的“零条纹”半径的计算公式:
其中,Fi(x,y)表示第i个环形面形图,ri表示第i个环形面形图的“零条纹”半径,θ表示极坐标下的角度坐标;
子步骤B2:对子步骤B1得到的计算公式采用数值微分法,得到环形面形图的“零条纹”位置分布图,再采用最小二乘法进行拟合,得到环形面形图的“零条纹”半径。
7.如权利要求5所述的测量方法,其特征在于,对于干涉条纹与中心连通的环形面形图,所述步骤B包括:
子步骤B1:得到环形面形图的“零条纹”半径的计算公式:
Fi(x,y)=a1 i(2ri 2-1)+a2 i(6ri 4-6ri+1)
其中,Fi(x,y)表示第i个环形面形图,ri表示第i个环形面形图的“零条纹”半径,a1 i和a2 i分别表示第i个环形面形图的离焦和球差系数;
子步骤B2:计算得到环形面形图的“零条纹”半径:
8.如权利要求5所述的测量方法,其特征在于,所述步骤C包括:
子步骤C1:得到非球面补偿器透射波前矢高的导数:
其中,ρi为光线与非球面补偿器透射波前的交点,z(ρi)表示非球面补偿器透射波前的矢高,z'(ρi)表示非球面补偿器透射波前矢高的导数,且
其中,ri为第i个环形面形图的“零条纹”半径,Zi为参考球面矢高,si为参考球面的与第i个环形面形图对应的相对位移,di为非球面补偿器透射波前的传播距离,θi为光线与光轴的夹角;
子步骤C2:由非球面补偿器透射波前矢高的导数,得到非球面补偿器透射波前方程第一表达式:
9.如权利要求5所述的测量方法,其特征在于,所述步骤D包括:
子步骤D1:非球面补偿器透射波前方程的导数表达式为:
其中,c为最佳拟合球曲率,ρ为非球面补偿器透射波前径坐标,ρmax为非球面补偿器透射波前径坐标最大值,M为常数,bn表示多项式{Pn(x)}的系数,且多项式{Pn(x)}满足:
子步骤D2:根据非球面补偿器透射波前方程的导数表达式,得到非球面补偿器透射波前方程第二表达式:
其中,Tn(x)=∫Pn(x)dx,n≥0。
子步骤D3:令非球面补偿器透射波前方程第一表达式与第二表达式相等,得到非球面补偿器透射波前方程迭代表达式:
10.如权利要求5所述的测量方法,其特征在于,所述步骤E包括:
子步骤E1:令非球面补偿器透射波前的传播距离di等于参考球面的与环形面形图对应的相对位移si;
子步骤E2:由环形面形图的“零条纹”半径ri得到光线与光轴的夹角θi以及参考球面矢高Zi,由所述“零条纹”半径ri、所述光线与光轴的夹角θi、所述参考球面矢高Zi、所述传播距离di得到光线与非球面补偿器透射波前的交点{ρi}、以及非球面补偿器透射波前矢高{z(ρi)};
子步骤E3:设置最佳拟合球曲率c的初始值,将所述相对位移si、所述交点ρi、以及所述非球面补偿器透射波前矢高z(ρi)代入非球面补偿器透射波前方程迭代表达式进行第一迭代运算,得到非球面补偿器透射波前方程第二表达式;
子步骤E4:将所述交点ρi代入子步骤E3得到的非球面补偿器透射波前方程第二表达式,重新计算非球面补偿器透射波前矢高z'(ρi),并利用下式重新计算非球面补偿器透射波前的传播距离:
其中,di'为重新计算的非球面补偿器透射波前的传播距离,z'(ρi)为重新计算的非球面补偿器透射波前矢高;
子步骤E5:计算两次迭代之间非球面补偿器透射波前的传播距离差值di'-di,将传播距离差值与传播距离阈值比较;
若传播距离差值大于传播距离阈值,将子步骤E4重新计算的传播距离di'作为子步骤E2的传播距离di,重复子步骤E2至子步骤E5;
否则,将子步骤E3得到的非球面补偿器透射波前方程第二表达式,作为非球面补偿器透射波前方程。
11.如权利要求10所述的测量方法,其特征在于,所述子步骤E3的第一迭代运算包括:
子分步骤E3a:将初始值赋予最佳拟合球曲率c;
子分步骤E3b:将所述最佳拟合球曲率c、所述相对位移si、所述交点ρi、以及所述非球面补偿器透射波前矢高z(ρi)代入非球面补偿器透射波前方程迭代表达式,采用最小二乘法拟合出多项式{Pn(x)}的系数{bn};
子分步骤E3c:得到非球面补偿器透射波前方程第二表达式;
子分步骤E3d:将所述交点ρi、以及所述非球面补偿器透射波前矢高z(ρi)代入非球面补偿器透射波前方程第二表达式,重新计算最佳拟合球曲率,将重新计算的最佳拟合球曲率与最佳拟合球曲率的初始值相减,获得最佳拟合球曲率差值;
子分步骤E3e:将最佳拟合球曲率差值与最佳拟合球曲率阈值比较,
若最佳拟合球曲率差值大于最佳拟合球曲率阈值,将子分步骤E3d重新计算的最佳拟合球曲率作为子分步骤E3b最佳拟合球曲率,重复执行子分步骤E3b至子分步骤E3e;
否则,将子分步骤E3c得到的非球面补偿器透射波前方程第二表达式作为迭代运算结果输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610886907.1A CN106404354A (zh) | 2016-10-11 | 2016-10-11 | 非球面补偿器透射波前方程的测量装置和方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610886907.1A CN106404354A (zh) | 2016-10-11 | 2016-10-11 | 非球面补偿器透射波前方程的测量装置和方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106404354A true CN106404354A (zh) | 2017-02-15 |
Family
ID=59229347
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610886907.1A Pending CN106404354A (zh) | 2016-10-11 | 2016-10-11 | 非球面补偿器透射波前方程的测量装置和方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106404354A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018068199A1 (zh) * | 2016-10-11 | 2018-04-19 | 中国科学院长春光学精密机械与物理研究所 | 非球面补偿器透射波前方程的测量装置和方法 |
CN108507488A (zh) * | 2018-03-05 | 2018-09-07 | 中国科学院上海光学精密机械研究所 | 基于轴向扫描的锥镜面形检测***及检测方法 |
CN112902875A (zh) * | 2021-03-31 | 2021-06-04 | 中国科学院长春光学精密机械与物理研究所 | 一种非球面反射镜曲率半径检测装置及方法 |
CN113624456A (zh) * | 2021-08-05 | 2021-11-09 | 苏州维纳仪器有限责任公司 | 多波长激光干涉装置 |
CN113776455A (zh) * | 2021-08-30 | 2021-12-10 | 中国科学院西安光学精密机械研究所 | 一种离轴非球面反射镜零位补偿检测非线性误差校正方法 |
CN114252023A (zh) * | 2021-12-28 | 2022-03-29 | 中国科学院光电技术研究所 | 用于非球面计算全息检测的计算机辅助调整装置及方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002310624A (ja) * | 2001-04-09 | 2002-10-23 | Canon Inc | 面形状計測法 |
CN101545760A (zh) * | 2008-03-26 | 2009-09-30 | 南京理工大学 | 光学透射球面检测装置 |
CN102043175A (zh) * | 2009-10-15 | 2011-05-04 | 夏普株式会社 | 透镜、透镜阵列、透镜评价装置 |
EP2375217A2 (en) * | 2010-03-31 | 2011-10-12 | Fujifilm Corporation | Aspheric surface measuring apparatus |
CN102661719A (zh) * | 2012-04-16 | 2012-09-12 | 中国人民解放军国防科学技术大学 | 用于非球面子孔径拼接测量的近零位补偿器及面形测量仪和测量方法 |
CN105547179A (zh) * | 2015-12-01 | 2016-05-04 | 中国科学院长春光学精密机械与物理研究所 | 一种非球面方程的测量方法 |
-
2016
- 2016-10-11 CN CN201610886907.1A patent/CN106404354A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002310624A (ja) * | 2001-04-09 | 2002-10-23 | Canon Inc | 面形状計測法 |
CN101545760A (zh) * | 2008-03-26 | 2009-09-30 | 南京理工大学 | 光学透射球面检测装置 |
CN102043175A (zh) * | 2009-10-15 | 2011-05-04 | 夏普株式会社 | 透镜、透镜阵列、透镜评价装置 |
EP2375217A2 (en) * | 2010-03-31 | 2011-10-12 | Fujifilm Corporation | Aspheric surface measuring apparatus |
CN102661719A (zh) * | 2012-04-16 | 2012-09-12 | 中国人民解放军国防科学技术大学 | 用于非球面子孔径拼接测量的近零位补偿器及面形测量仪和测量方法 |
CN105547179A (zh) * | 2015-12-01 | 2016-05-04 | 中国科学院长春光学精密机械与物理研究所 | 一种非球面方程的测量方法 |
Non-Patent Citations (2)
Title |
---|
高松涛等: "用计算全息图对非球面的高精度检测与误差评估", 《光学学报》 * |
高松涛等: "用计算全息图校正非球面的畸变", 《光学精密工程》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018068199A1 (zh) * | 2016-10-11 | 2018-04-19 | 中国科学院长春光学精密机械与物理研究所 | 非球面补偿器透射波前方程的测量装置和方法 |
CN108507488A (zh) * | 2018-03-05 | 2018-09-07 | 中国科学院上海光学精密机械研究所 | 基于轴向扫描的锥镜面形检测***及检测方法 |
CN108507488B (zh) * | 2018-03-05 | 2019-12-20 | 中国科学院上海光学精密机械研究所 | 基于轴向扫描的锥镜面形检测***及检测方法 |
CN112902875A (zh) * | 2021-03-31 | 2021-06-04 | 中国科学院长春光学精密机械与物理研究所 | 一种非球面反射镜曲率半径检测装置及方法 |
CN113624456A (zh) * | 2021-08-05 | 2021-11-09 | 苏州维纳仪器有限责任公司 | 多波长激光干涉装置 |
CN113776455A (zh) * | 2021-08-30 | 2021-12-10 | 中国科学院西安光学精密机械研究所 | 一种离轴非球面反射镜零位补偿检测非线性误差校正方法 |
CN114252023A (zh) * | 2021-12-28 | 2022-03-29 | 中国科学院光电技术研究所 | 用于非球面计算全息检测的计算机辅助调整装置及方法 |
CN114252023B (zh) * | 2021-12-28 | 2022-10-21 | 中国科学院光电技术研究所 | 用于非球面计算全息检测的计算机辅助调整装置及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106404354A (zh) | 非球面补偿器透射波前方程的测量装置和方法 | |
CN105423948B (zh) | 采用变形镜的拼接干涉检测非球面面形的装置 | |
CN103591888B (zh) | 大口径离轴非球面光学元件几何参数的测算方法 | |
CN105547179B (zh) | 一种非球面方程的测量方法 | |
CN102997863B (zh) | 一种全口径光学非球面面形误差直接检测*** | |
CN101290218B (zh) | 用于校正非球面非零位检测时的原理误差的方法 | |
CN103335610B (zh) | 大口径凸高次非球面的检测*** | |
CN103134660B (zh) | 基于像散分解获得望远镜主次镜对准误差的方法 | |
CN110188321B (zh) | 一种基于神经网络算法的主次镜校准方法 | |
CN103776389A (zh) | 一种高精度非球面组合干涉检测装置与方法 | |
CN103471522B (zh) | 检测范围广的凹非球面的实时检测方法 | |
CN102620683B (zh) | 子孔径拼接检测非球面调整误差补偿方法 | |
CN103969787A (zh) | 一种离轴四反镜头的初装定位方法 | |
CN103234480A (zh) | 一种环形凸非球面的快速面形检测方法 | |
CN110715619B (zh) | 基于自适应非零干涉法的光学自由曲面面形误差测量方法 | |
CN103471521B (zh) | 快速、准确的光学非球面的实时检测方法 | |
CN104075668B (zh) | 凸双曲面Hindle检测中的高精度几何量测量方法 | |
CN103412391A (zh) | 一种基于激光跟踪仪实现光学***穿轴对心方法 | |
CN102645202A (zh) | 大口径非球面工件轮廓的测量方法 | |
CN105318847A (zh) | 基于***建模的非球面非零位环形子孔径拼接方法 | |
CN105352453B (zh) | 非零位干涉检测***中非球面顶点球曲率半径测量方法 | |
CN105157572A (zh) | 一种用于非球面环形子孔径拼接的中心偏移误差补偿方法 | |
CN105352451B (zh) | 一种基于可变形镜的准万能补偿镜及设计方法 | |
CN106705888A (zh) | 干涉检测中的ccd坐标系与镜面坐标系非线性关系标定方法 | |
CN206362310U (zh) | 一种非球面检测光路中光学间隔测量*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170215 |