CN107121114A - 基于功率谱的大口径平面镜低阶像差估计方法 - Google Patents

基于功率谱的大口径平面镜低阶像差估计方法 Download PDF

Info

Publication number
CN107121114A
CN107121114A CN201710316123.XA CN201710316123A CN107121114A CN 107121114 A CN107121114 A CN 107121114A CN 201710316123 A CN201710316123 A CN 201710316123A CN 107121114 A CN107121114 A CN 107121114A
Authority
CN
China
Prior art keywords
mrow
msub
mover
msup
sub
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.)
Granted
Application number
CN201710316123.XA
Other languages
English (en)
Other versions
CN107121114B (zh
Inventor
安其昌
张景旭
杨飞
赵宏超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changchun Institute of Optics Fine Mechanics and Physics of CAS
Original Assignee
Changchun Institute of Optics Fine Mechanics and Physics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changchun Institute of Optics Fine Mechanics and Physics of CAS filed Critical Changchun Institute of Optics Fine Mechanics and Physics of CAS
Priority to CN201710316123.XA priority Critical patent/CN107121114B/zh
Publication of CN107121114A publication Critical patent/CN107121114A/zh
Application granted granted Critical
Publication of CN107121114B publication Critical patent/CN107121114B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/20Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring contours or curvatures, e.g. determining profile

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)

Abstract

基于功率谱的大口径平面镜低阶像差估计方法,涉及大口径平面镜面形评价领域,解决了现有估计方法存在的计算量大、全面性低的问题。该方法包括:将Zernike多项式在频域上的表达代入不同子孔径的不同阶Zernike多项式系数之间的互相关系数在频域上的表达,利用上述获得的公式以及全口径Zernike多项式系数βj的定义进行计算得到子孔径Zernike多项式系数αi与全口径Zernike多项式系数βj之间的关系: 为全口径Zernike多项式系数平方的期望;将上述公式进行展开讨论,当阶数i大于3时:本发明利用子孔径镜面功率谱的统计学特征估计大口径平面镜的低阶像差,节省了检测成本,提高了检测精度。

Description

基于功率谱的大口径平面镜低阶像差估计方法
技术领域
本发明涉及大口径平面镜面形评价技术领域,具体涉及一种基于功率谱的大口径平面镜低阶像差估计方法。
背景技术
子孔径拼接技术最早是由美国专家C.Kim和J.Wyant共同提出的,其基本原理是利用互相重叠的子区域,结合适合的拼接算法以获得全口径的面形数据。目前,国内外利用子孔径拼接技术检测大口径平面镜都取得了一定的进展。然而,随着大口径平面镜口径的增加,子孔径的数量也会随之增加,同时由于拼接算法的关系,将各个子孔径拼接完成之后提取的低阶像差(离焦、像散)测量敏感度也随之降低。如果为了增加其估计精度而使用更大口径的平面干涉仪,其成本也难以进行控制。另一方面,进行子孔径拼接时,往往需要进行多次测量,仅使用一次测量的数据对时间、人员成本也造成了极大的浪费。因此,为了精确估计大口径平面镜低阶像差,亟需对子孔径拼接算法加以改进。
随着大口径望远镜的发展,仅依靠波前均方根是不能有效的表征其全频域的特征。功率谱是功率谱密度函数的简称,定义为单位频带内的信号功率,表示信号功率随着频率的变化情况,即信号功率在频域的分布状况。功率谱评价方法是上个世纪末NIF所提出的,已经有了较为成熟的标准(ISO 10110)。但是,功率谱评价方法只能评价某一方向的波前(波在介质中传播时,某时刻刚刚开始位移的质点构成的面,称为波前,它代表某时刻波能量到达的空间位置,它是运动着的)起伏,对于超光滑表面来说,可以取某个方向作为评价的标准。但是对于大口径望远镜来说,由于成本限制以及大气的影响,其面形没有必要达到亚纳米级,在此尺度下,面形在较大尺度上的各项异性就变得十分明显。之后对于功率谱的研究引入了二维功率谱,最后坍陷到一维来进行面形评价。但是傅里叶变换必须是针对正交的数据进行处理,而对于其它角度上的面形的统计学特征就无法评价。另一方面,在低频的部分,功率谱评价方法由于缺乏平均来降低噪声,故其评价的效果也会受影响。
发明内容
为了解决现有估计方法存在的计算量大、全面性低的问题,本发明提供一种基于功率谱的大口径平面镜低阶像差估计方法。
本发明为解决技术问题所采用的技术方案如下:
本发明的基于功率谱的大口径平面镜低阶像差估计方法,包括以下步骤:
步骤一、将Zernike多项式在频域上的表达代入不同子孔径的不同阶Zernike多项式系数之间的互相关系数在频域上的表达,得到式(5):
式(5)中,为第一个子孔径第i阶Zernike多项式系数的共轭,αj'为第二个子孔径第j阶Zernike多项式系数,Φsub为子孔径面形数据,i为整数,j'为整数,m为第一个子孔径周向对称数,m′为第二个子孔周向对称数,n为第一个子孔径轴向对称数,n′为第二个子孔径轴向对称数,R为大口径平面镜的全口径半径,为第一个子孔径所对应的空间频率矢量,为n+1阶0型贝塞尔函数,为n′+1阶0型贝塞尔函数;
步骤二、利用式(5)以及全口径Zernike多项式系数βj的定义进行计算得到子孔径Zernike多项式系数αi与全口径Zernike多项式系数βj之间的关系,如式(9)所示:
式(9)中,为全口径Zernike多项式系数平方的期望;
步骤三、将式(9)进行展开得到式(10)、式(11)和式(12):
其中,为全孔径第1阶Zernike系数平方的期望,为全孔径第2阶Zernike系数平方的期望,为全孔径第3阶Zernike系数平方的期望,β1为全孔径第1阶Zernike系数,β2为全孔径第2阶Zernike系数,β3为全孔径第3阶Zernike系数,Rsub为大口径平面镜的子孔径半径,Rsub=μR,即α1为子孔径第1阶Zernike系数,α2为子孔径第2阶Zernike系数,α3为子孔径第3阶Zernike系数,α7为子孔径第7阶Zernike系数,α8为子孔径第8阶Zernike系数,ρO为子孔径偏心距离;
当阶数i大于3时:
进一步的,步骤一中,获得不同子孔径的不同阶Zernike多项式系数之间的互相关系数在频域上的表达的具体过程如下:
①不同子孔径的不同阶Zernike多项式系数之间的互相关系数的定义如式(1)所示:
式(1)中,为第一个子孔径位置矢量,为第二个子孔径位置矢量,为第一个子孔径的共轭,为第二个子孔径的波前,为第一个子孔径对应的i阶Zernike多项式,第二个子孔径对应的j'阶Zernike多项式;
②结合快速傅里叶变换以及功率谱对式(1)进行计算得到式(2):
式(2)中,PSDsub(f/R,f′/R)为各个子孔径平均所获得的功率谱,f为第一个子孔径所对应的空间频率,f′为第二个子孔径所对应的空间频率;
③将式(2)带入式(1)中,使式(1)在频域上重新表达,获得式(3):
式(3)中,Fsub为子孔径的频域数据,为第一个子孔径的第i阶Zernike多项式的频域数据,为第二个子孔径的第j'阶Zernike多项式的频域数据,为第一个子孔径所对应的空间频率矢量,为第二个子孔径所对应的空间频率矢量。
进一步的,步骤一中,所述Zernike多项式在频域上的表达如式(4)所示:
式(4)中,为Zernike多项式在频域上的表达式,为n+1阶0型贝塞尔函数,φ为波前整体的相位偏置。
进一步的,步骤二的具体过程如下:
①子孔径面形数据Φsub与全口径面形数据Φfull在对应位置的数据相同,如式(6)和式(7)所示:
其中,αi为子孔径Zernike多项式系数,βj为全口径Zernike多项式系数,i为整数,j'为整数,为子孔径Zernike多项式,为全口径Zernike多项式,O点为全口径覆盖区域的镜面圆心,O′点为子孔径覆盖区域的镜面圆心,A点为子孔径覆盖区域的镜面圆周上的点,为O点到A点的位置矢量,为O′到A点的位置矢量,为O点到O′点的位置矢量,Rsub为大口径平面镜的子孔径半径,Rsub=μR,即
②对式(5)、式(6)和式(7)进行计算得到式(8):
式(8)中,为全孔径面形数,为第j阶Zernike多项式,为第i阶Zernike多项式;
③对式(8)进行计算得到子孔径Zernike多项式系数αi与全口径Zernike多项式系数βj之间的关系,如式(9)所示。
本发明的有益效果是:本发明利用子孔径镜面功率谱的统计学特征,估计大口径平面镜的低阶像差,计算量大幅度降低,检测比较全面,不但节省了大口径平面镜拼接检测成本,而且提高了低阶像差的检测精度。
附图说明
图1为子孔径参数示意图。
图2为利用功率谱估计大口径平面镜低阶像差流程图。
具体实施方式
以下结合附图对本发明作进一步详细说明。
本发明的一种基于功率谱的大口径平面镜低阶像差估计方法,可以在不获得镜面全面形的情况下,实现对大口径平面镜低阶像差的估计。该方法如图2所示,其实现的具体过程如下:
一、不同子孔径的不同阶Zernike多项式系数之间的互相关系数的基本定义如式(1)所示:
式(1)中,为第一个子孔径第i阶Zernike多项式系数的共轭,αj'为第二个子孔径第j'阶Zernike多项式系数,Φsub为子孔径面形数据,为第一个子孔径位置矢量,为第二个子孔径位置矢量,i为整数,j'为整数,为第一个子孔径的共轭,一般来说,第一个子孔径的波前为纯实数,在此为了和频域计算相对应,故取其共轭;为第二个子孔径的波前,为第一个子孔径对应的i阶Zernike多项式,第二个子孔径对应的j'阶Zernike多项式。
二、通过观察式(1),结合快速傅里叶变换(FFT)以及功率谱(PSD)计算方法可得:
式(2)中,PSDsub(f/R,f′/R)为各个子孔径平均所获得的功率谱,R为大口径平面镜的全口径半径,f为第一个子孔径所对应的空间频率,f′为第二个子孔径所对应的空间频率,为第一个子孔径的共轭,为第二个子孔径的波前。实际计算的时候,首先分别计算各个子孔径的功率谱,之后进行平均,要求子孔径尽量覆盖全口径。
三、将式(2)带入式(1)中,使式(1)在频域上重新表达,获得式(3):
式(3)中,为第一个子孔径第i阶Zernike多项式系数的共轭,αj'为第二个子孔径第j阶Zernike多项式系数,Φsub为子孔径面形数据,i为整数,j'为整数,Fsub为子孔径的频域数据,为第一个子孔径的第i阶Zernike多项式的频域数据,为第二个子孔径的第j'阶Zernike多项式的频域数据,f为第一个子孔径所对应的空间频率,为第一个子孔径所对应的空间频率矢量,f′为第二个子孔径所对应的空间频率,为第二个子孔径所对应的空间频率矢量,R为大口径平面镜的全口径半径,而其中的Zernike多项式在频域上的表达如式(4)所示:
式(4)中,为Zernike多项式在频域上的表达式,为n+1阶0型贝塞尔函数,m为第一个子孔径周向对称数,n为第一个子孔径轴向对称数,为第一个子孔径所对应的空间频率矢量,φ为波前整体的相位偏置,i为整数,j'为整数。
四、将Zernike多项式在频域上的表达代入不同子孔径的不同阶Zernike多项式系数之间的互相关系数在频域上重新表达,即将式(4)代入式(3)可以得到基于功率谱(PSD)的相关系数表达,如式(5)所示:
式(5)中,为第一个子孔径第i阶Zernike多项式系数的共轭,αj'为第二个子孔径第j阶Zernike多项式系数,Φsub为子孔径面形数据,i为整数,j'为整数,m为第一个子孔径周向对称数,m′为第二个子孔周向对称数,n为第一个子孔径轴向对称数,n′为第二个子孔径轴向对称数,R为大口径平面镜的全口径半径,为第一个子孔径所对应的空间频率矢量,为n+1阶0型贝塞尔函数,为n′+1阶0型贝塞尔函数。
五、利用式(5)可以建立起功率谱和Zernike多项式。利用子孔径数据对大口径平面镜镜面面形估计的核心思想在于子孔径面形数据Φsub与全口径面形数据Φfull在对应位置的数据相同,如式(6)和式(7)所示:
即:
式(6)和式(7)中,αi为子孔径Zernike多项式系数,βj为全口径Zernike多项式系数,i为整数,j'为整数,为子孔径Zernike多项式,为全口径Zernike多项式。
如图1所示:O点为全口径覆盖区域的镜面圆心,O′点为子孔径覆盖区域的镜面圆心,A点为子孔径覆盖区域的镜面圆周上的点,为O点到A点的位置矢量,为O′到A点的位置矢量,为O点到O′点的位置矢量,R为大口径平面镜的全口径半径,Rsub为大口径平面镜的子孔径半径,Rsub=μR,即
六、利用式(5)以及全口径Zernike多项式系数βj的最基本定义式(6)和式(7)可得:
式(8)中,为全孔径面形数,为第j阶Zernike多项式,为第i阶Zernike多项式,Φsub为子孔径面形数据,αi为子孔径Zernike多项式系数,i为整数,j'为整数。
七、对式(8)进行计算可得子孔径Zernike多项式系数αi与全口径Zernike多项式系数βj之间的关系,如式(9)所示:
式(9)中,为全口径Zernike多项式系数平方的期望,βj为全口径Zernike多项式系数,αi为子孔径Zernike多项式系数,Φsub为子孔径面形数据。
八、将式(9)展开并进行分情况讨论可得:
式(10)、式(11)和式(12)中,为全孔径第1阶Zernike系数平方的期望,为全孔径第2阶Zernike系数平方的期望,为全孔径第3阶Zernike系数平方的期望,β1为全孔径第1阶Zernike系数,β2为全孔径第2阶Zernike系数,β3为全孔径第3阶Zernike系数,R为大口径平面镜的全口径半径,Rsub为大口径平面镜的子孔径半径,Rsub=μR,即α1为子孔径第1阶Zernike系数,α2为子孔径第2阶Zernike系数,α3为子孔径第3阶Zernike系数,α7为子孔径第7阶Zernike系数,α8为子孔径第8阶Zernike系数,ρO为子孔径偏心距离。
当阶数i大于3时:
式(13)中,αi为子孔径Zernike多项式系数,i为整数,R为大口径平面镜的全口径半径,Rsub为大口径平面镜的子孔径半径,Rsub=μR,即
以低阶像差中的离焦为例,对于离焦n=2,m=0,结合式(5)与式(10)、式(11)、式(12)可得离焦的表达式如(14)所示:
式(14)中,β4为全孔径第4阶Zernike系数,为第一个子孔径所对应的空间频率矢量,为3阶0型贝塞尔函数,R为大口径平面镜的全口径半径,Rsub为大口径平面镜的子孔径半径,Rsub=μR,即
需要说明的是:上述基于功率谱的大口径平面镜低阶像差估计方法中,涉及所有公式中的同一字母表示的含义是相同的。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (4)

1.基于功率谱的大口径平面镜低阶像差估计方法,其特征在于,包括以下步骤:
步骤一、将Zernike多项式在频域上的表达代入不同子孔径的不同阶Zernike多项式系数之间的互相关系数在频域上的表达,得到式(5):
式(5)中,为第一个子孔径第i阶Zernike多项式系数的共轭,αj'为第二个子孔径第j阶Zernike多项式系数,Φsub为子孔径面形数据,i为整数,j'为整数,m为第一个子孔径周向对称数,m′为第二个子孔周向对称数,n为第一个子孔径轴向对称数,n′为第二个子孔径轴向对称数,R为大口径平面镜的全口径半径,为第一个子孔径所对应的空间频率矢量,为n+1阶0型贝塞尔函数,为n′+1阶0型贝塞尔函数;
步骤二、利用式(5)以及全口径Zernike多项式系数βj的定义进行计算得到子孔径Zernike多项式系数αi与全口径Zernike多项式系数βj之间的关系,如式(9)所示:
<mrow> <mo>&lt;</mo> <msubsup> <mi>&amp;beta;</mi> <mi>j</mi> <mn>2</mn> </msubsup> <msub> <mo>&gt;</mo> <mrow> <mi>f</mi> <mi>u</mi> <mi>l</mi> <mi>l</mi> </mrow> </msub> <mo>=</mo> <mo>&lt;</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;infin;</mi> </munderover> <msub> <mi>&amp;alpha;</mi> <mi>i</mi> </msub> <munder> <mrow> <mo>&amp;Integral;</mo> <mo>&amp;Integral;</mo> </mrow> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </munder> <msub> <mi>Z</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>R</mi> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>f</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>Z</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mi>d</mi> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>f</mi> </msub> <msub> <mo>&gt;</mo> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
式(9)中,为全口径Zernike多项式系数平方的期望;
步骤三、将式(9)进行展开得到式(10)、式(11)和式(12):
<mrow> <mo>&lt;</mo> <msubsup> <mi>&amp;beta;</mi> <mn>1</mn> <mn>2</mn> </msubsup> <msub> <mo>&gt;</mo> <mrow> <mi>f</mi> <mi>u</mi> <mi>l</mi> <mi>l</mi> </mrow> </msub> <mo>=</mo> <msup> <mi>&amp;mu;</mi> <mrow> <mn>5</mn> <mo>/</mo> <mn>3</mn> </mrow> </msup> <mo>&lt;</mo> <msubsup> <mi>&amp;alpha;</mi> <mn>1</mn> <mn>2</mn> </msubsup> <msub> <mo>&gt;</mo> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <mo>-</mo> <mn>4</mn> <msubsup> <mi>&amp;rho;</mi> <mn>0</mn> <mn>2</mn> </msubsup> <mo>&lt;</mo> <msubsup> <mi>&amp;alpha;</mi> <mn>2</mn> <mn>2</mn> </msubsup> <msub> <mo>&gt;</mo> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <mo>-</mo> <mn>4</mn> <msqrt> <mn>2</mn> </msqrt> <msub> <mi>&amp;rho;</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mn>3</mn> <msubsup> <mi>&amp;rho;</mi> <mn>0</mn> <mn>3</mn> </msubsup> <mo>-</mo> <mn>2</mn> <msub> <mi>&amp;rho;</mi> <mn>0</mn> </msub> <mo>+</mo> <mn>3</mn> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>&lt;</mo> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <msub> <mi>&amp;alpha;</mi> <mn>8</mn> </msub> <msub> <mo>&gt;</mo> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mo>&lt;</mo> <msubsup> <mi>&amp;beta;</mi> <mn>2</mn> <mn>2</mn> </msubsup> <msub> <mo>&gt;</mo> <mrow> <mi>f</mi> <mi>u</mi> <mi>l</mi> <mi>l</mi> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msup> <mi>&amp;mu;</mi> <mrow> <mn>5</mn> <mo>/</mo> <mn>3</mn> </mrow> </msup> <mo>-</mo> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>&lt;</mo> <msubsup> <mi>&amp;alpha;</mi> <mn>2</mn> <mn>2</mn> </msubsup> <msub> <mo>&gt;</mo> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <mo>-</mo> <msqrt> <mn>2</mn> </msqrt> <mrow> <mo>(</mo> <mn>9</mn> <msubsup> <mi>&amp;rho;</mi> <mn>0</mn> <mn>2</mn> </msubsup> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>2</mn> <msup> <mi>&amp;mu;</mi> <mn>4</mn> </msup> <mo>+</mo> <mn>2</mn> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>&lt;</mo> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <msub> <mi>&amp;alpha;</mi> <mn>8</mn> </msub> <msub> <mo>&gt;</mo> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mo>&lt;</mo> <msubsup> <mi>&amp;beta;</mi> <mn>3</mn> <mn>2</mn> </msubsup> <msub> <mo>&gt;</mo> <mrow> <mi>f</mi> <mi>u</mi> <mi>l</mi> <mi>l</mi> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msup> <mi>&amp;mu;</mi> <mrow> <mn>5</mn> <mo>/</mo> <mn>3</mn> </mrow> </msup> <mo>-</mo> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>&lt;</mo> <msubsup> <mi>&amp;alpha;</mi> <mn>3</mn> <mn>2</mn> </msubsup> <msub> <mo>&gt;</mo> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <mo>-</mo> <msqrt> <mn>2</mn> </msqrt> <mrow> <mo>(</mo> <mn>3</mn> <msubsup> <mi>&amp;rho;</mi> <mn>0</mn> <mn>2</mn> </msubsup> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>2</mn> <msup> <mi>&amp;mu;</mi> <mn>4</mn> </msup> <mo>+</mo> <mn>2</mn> <msup> <mi>&amp;mu;</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>&lt;</mo> <msub> <mi>&amp;alpha;</mi> <mn>3</mn> </msub> <msub> <mi>&amp;alpha;</mi> <mn>7</mn> </msub> <msub> <mo>&gt;</mo> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
其中,为全孔径第1阶Zernike系数平方的期望,为全孔径第2阶Zernike系数平方的期望,为全孔径第3阶Zernike系数平方的期望,β1为全孔径第1阶Zernike系数,β2为全孔径第2阶Zernike系数,β3为全孔径第3阶Zernike系数,Rsub为大口径平面镜的子孔径半径,Rsub=μR,即α1为子孔径第1阶Zernike系数,α2为子孔径第2阶Zernike系数,α3为子孔径第3阶Zernike系数,α7为子孔径第7阶Zernike系数,α8为子孔径第8阶Zernike系数,ρO为子孔径偏心距离;
当阶数i大于3时:
<mrow> <mo>&lt;</mo> <msubsup> <mi>&amp;beta;</mi> <mrow> <mi>i</mi> <mo>&gt;</mo> <mn>3</mn> </mrow> <mn>2</mn> </msubsup> <msub> <mo>&gt;</mo> <mrow> <mi>f</mi> <mi>u</mi> <mi>l</mi> <mi>l</mi> </mrow> </msub> <mo>=</mo> <msup> <mi>&amp;mu;</mi> <mrow> <mn>5</mn> <mo>/</mo> <mn>3</mn> </mrow> </msup> <mo>&lt;</mo> <msubsup> <mi>&amp;alpha;</mi> <mi>i</mi> <mn>2</mn> </msubsup> <msub> <mo>&gt;</mo> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
2.根据权利要求1所述的基于功率谱的大口径平面镜低阶像差估计方法,其特征在于,步骤一中,获得不同子孔径的不同阶Zernike多项式系数之间的互相关系数在频域上的表达的具体过程如下:
①不同子孔径的不同阶Zernike多项式系数之间的互相关系数的定义如式(1)所示:
<mrow> <mo>&lt;</mo> <msub> <mover> <mi>&amp;alpha;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <msub> <mi>a</mi> <msup> <mi>j</mi> <mo>&amp;prime;</mo> </msup> </msub> <msub> <mo>&gt;</mo> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </msub> <mo>=</mo> <munder> <mrow> <mo>&amp;Integral;</mo> <mo>&amp;Integral;</mo> </mrow> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </munder> <mover> <mi>&amp;Phi;</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>)</mo> </mrow> <msup> <mi>&amp;Phi;</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msup> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> <msub> <mover> <mi>Z</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>)</mo> </mrow> <msub> <mi>Z</mi> <msup> <mi>j</mi> <mo>&amp;prime;</mo> </msup> </msub> <mrow> <mo>(</mo> <msup> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> <mi>d</mi> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>d</mi> <msup> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
式(1)中,为第一个子孔径位置矢量,为第二个子孔径位置矢量,为第一个子孔径的共轭,为第二个子孔径的波前,为第一个子孔径对应的i阶Zernike多项式,第二个子孔径对应的j'阶Zernike多项式;
②结合快速傅里叶变换以及功率谱对式(1)进行计算得到式(2):
<mrow> <mi>F</mi> <mi>F</mi> <mi>T</mi> <mo>&amp;lsqb;</mo> <mover> <mi>&amp;Phi;</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>)</mo> </mrow> <msup> <mi>&amp;Phi;</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msup> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>=</mo> <mfrac> <mn>1</mn> <msup> <mi>R</mi> <mn>2</mn> </msup> </mfrac> <msub> <mi>PSD</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>/</mo> <mi>R</mi> <mo>,</mo> <msup> <mi>f</mi> <mo>&amp;prime;</mo> </msup> <mo>/</mo> <mi>R</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
式(2)中,PSDsub(f/R,f′/R)为各个子孔径平均所获得的功率谱,f为第一个子孔径所对应的空间频率,f′为第二个子孔径所对应的空间频率;
③将式(2)带入式(1)中,使式(1)在频域上重新表达,获得式(3):
<mrow> <mo>&lt;</mo> <msub> <mover> <mi>&amp;alpha;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <msub> <mi>a</mi> <msup> <mi>j</mi> <mo>&amp;prime;</mo> </msup> </msub> <msub> <mo>&gt;</mo> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </msub> <mo>=</mo> <mo>&lt;</mo> <msub> <mover> <mi>&amp;alpha;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <msub> <mi>a</mi> <msup> <mi>j</mi> <mo>&amp;prime;</mo> </msup> </msub> <msub> <mo>&gt;</mo> <msub> <mi>F</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <msup> <mi>R</mi> <mn>2</mn> </msup> </mfrac> <munder> <mrow> <mo>&amp;Integral;</mo> <mo>&amp;Integral;</mo> </mrow> <msub> <mi>F</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </munder> <msub> <mover> <mi>Q</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mover> <mi>f</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>)</mo> </mrow> <msub> <mi>Q</mi> <msup> <mi>j</mi> <mo>&amp;prime;</mo> </msup> </msub> <mrow> <mo>(</mo> <msup> <mover> <mi>f</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> <msub> <mi>PSD</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>/</mo> <mi>R</mi> <mo>,</mo> <msup> <mi>f</mi> <mo>&amp;prime;</mo> </msup> <mo>/</mo> <mi>R</mi> <mo>)</mo> </mrow> <mi>d</mi> <mover> <mi>f</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>d</mi> <msup> <mover> <mi>f</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
式(3)中,Fsub为子孔径的频域数据,为第一个子孔径的第i阶Zernike多项式的频域数据,为第二个子孔径的第j'阶Zernike多项式的频域数据,为第一个子孔径所对应的空间频率矢量,为第二个子孔径所对应的空间频率矢量。
3.根据权利要求1所述的基于功率谱的大口径平面镜低阶像差估计方法,其特征在于,步骤一中,所述Zernike多项式在频域上的表达如式(4)所示:
<mrow> <msub> <mi>Q</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mover> <mi>f</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>,</mo> <mi>&amp;phi;</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <mi>n</mi> <mo>+</mo> <mn>1</mn> </mrow> </msqrt> <mfrac> <mrow> <msubsup> <mi>J</mi> <mrow> <mi>n</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mn>2</mn> <mi>&amp;pi;</mi> <mover> <mi>f</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&amp;pi;</mi> <mover> <mi>f</mi> <mo>&amp;RightArrow;</mo> </mover> </mrow> </mfrac> <mo>&amp;times;</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <mrow> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mi>m</mi> </mrow> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> </mrow> </msup> <msup> <mi>j</mi> <mi>m</mi> </msup> <msqrt> <mn>2</mn> </msqrt> <mi>cos</mi> <mrow> <mo>(</mo> <mrow> <mi>m</mi> <mi>&amp;phi;</mi> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>m</mi> <mo>&amp;NotEqual;</mo> <mn>0</mn> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mo>=</mo> <mi>e</mi> <mi>v</mi> <mi>e</mi> <mi>n</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <mrow> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mi>m</mi> </mrow> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> </mrow> </msup> <msup> <mi>j</mi> <mi>m</mi> </msup> <msqrt> <mn>2</mn> </msqrt> <mi>sin</mi> <mrow> <mo>(</mo> <mrow> <mi>m</mi> <mi>&amp;phi;</mi> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>m</mi> <mo>&amp;NotEqual;</mo> <mn>0</mn> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mo>=</mo> <mi>o</mi> <mi>d</mi> <mi>d</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <msup> <mrow> <mo>(</mo> <mrow> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <mrow> <mi>n</mi> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mtd> <mtd> <mrow> <mi>m</mi> <mo>=</mo> <mn>0</mn> </mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
式(4)中,为Zernike多项式在频域上的表达式,为n+1阶0型贝塞尔函数,φ为波前整体的相位偏置。
4.根据权利要求1所述的基于功率谱的大口径平面镜低阶像差估计方法,其特征在于,步骤二的具体过程如下:
①子孔径面形数据Φsub与全口径面形数据Φfull在对应位置的数据相同,如式(6)和式(7)所示:
<mrow> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>f</mi> <mi>u</mi> <mi>l</mi> <mi>l</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>R</mi> <mo>,</mo> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>f</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>&amp;mu;</mi> <mi>R</mi> <mo>,</mo> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;infin;</mi> </munderover> <msub> <mi>&amp;alpha;</mi> <mi>i</mi> </msub> <msub> <mi>Z</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;mu;</mi> <mi>R</mi> <mo>,</mo> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;infin;</mi> </munderover> <msub> <mi>&amp;beta;</mi> <mi>j</mi> </msub> <msub> <mi>Z</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>R</mi> <mo>,</mo> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>f</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
其中,αi为子孔径Zernike多项式系数,βj为全口径Zernike多项式系数,i为整数,j'为整数,为子孔径Zernike多项式,为全口径Zernike多项式,O点为全口径覆盖区域的镜面圆心,O′点为子孔径覆盖区域的镜面圆心,A点为子孔径覆盖区域的镜面圆周上的点,为O点到A点的位置矢量,为O′到A点的位置矢量,为O点到O′点的位置矢量,Rsub为大口径平面镜的子孔径半径,Rsub=μR,即
②对式(5)、式(6)和式(7)进行计算得到式(8):
<mrow> <msub> <mi>&amp;beta;</mi> <mi>j</mi> </msub> <mo>=</mo> <munder> <mrow> <mo>&amp;Integral;</mo> <mo>&amp;Integral;</mo> </mrow> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </munder> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>f</mi> <mi>u</mi> <mi>l</mi> <mi>l</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>R</mi> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>f</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>Z</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>R</mi> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>f</mi> </msub> <mo>)</mo> </mrow> <mi>d</mi> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>f</mi> </msub> <mo>=</mo> <munder> <mrow> <mo>&amp;Integral;</mo> <mo>&amp;Integral;</mo> </mrow> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </munder> <mrow> <mo>(</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;infin;</mi> </munderover> <msub> <mi>&amp;alpha;</mi> <mi>i</mi> </msub> <msub> <mi>Z</mi> <mi>i</mi> </msub> <mo>(</mo> <mrow> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>s</mi> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <msub> <mi>Z</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>R</mi> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>f</mi> </msub> <mo>)</mo> </mrow> <mi>d</mi> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>f</mi> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;infin;</mi> </munderover> <msub> <mi>&amp;alpha;</mi> <mi>i</mi> </msub> <munder> <mrow> <mo>&amp;Integral;</mo> <mo>&amp;Integral;</mo> </mrow> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </munder> <msub> <mi>Z</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>R</mi> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>f</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>Z</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mi>d</mi> <msub> <mover> <mi>&amp;rho;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>f</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
式(8)中,为全孔径面形数,为第j阶Zernike多项式,为第i阶Zernike多项式;
③对式(8)进行计算得到子孔径Zernike多项式系数αi与全口径Zernike多项式系数βj之间的关系,如式(9)所示。
CN201710316123.XA 2017-05-08 2017-05-08 基于功率谱的大口径平面镜低阶像差估计方法 Active CN107121114B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710316123.XA CN107121114B (zh) 2017-05-08 2017-05-08 基于功率谱的大口径平面镜低阶像差估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710316123.XA CN107121114B (zh) 2017-05-08 2017-05-08 基于功率谱的大口径平面镜低阶像差估计方法

Publications (2)

Publication Number Publication Date
CN107121114A true CN107121114A (zh) 2017-09-01
CN107121114B CN107121114B (zh) 2019-05-07

Family

ID=59727776

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710316123.XA Active CN107121114B (zh) 2017-05-08 2017-05-08 基于功率谱的大口径平面镜低阶像差估计方法

Country Status (1)

Country Link
CN (1) CN107121114B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111649701A (zh) * 2020-06-30 2020-09-11 长春博信光电子有限公司 一种复曲面镜偏心值检测方法和装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05296879A (ja) * 1992-04-17 1993-11-12 Olympus Optical Co Ltd 光学性能測定方法及び装置
JPH10325710A (ja) * 1997-03-28 1998-12-08 Fuji Photo Optical Co Ltd 平面形状測定解析方法
CN1746648A (zh) * 2005-10-31 2006-03-15 中国科学院光电技术研究所 一种大口径深型非球面镜检测系统
CN102243068A (zh) * 2011-04-27 2011-11-16 南京理工大学 一种子孔径拼接中系统误差的修正方法
CN103791854A (zh) * 2014-01-23 2014-05-14 中国科学院长春光学精密机械与物理研究所 高空间分辨率子孔径拼接方法
CN104776812A (zh) * 2015-04-14 2015-07-15 中国科学院上海光学精密机械研究所 一种提高平面面形子孔径拼接检测精度的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05296879A (ja) * 1992-04-17 1993-11-12 Olympus Optical Co Ltd 光学性能測定方法及び装置
JPH10325710A (ja) * 1997-03-28 1998-12-08 Fuji Photo Optical Co Ltd 平面形状測定解析方法
CN1746648A (zh) * 2005-10-31 2006-03-15 中国科学院光电技术研究所 一种大口径深型非球面镜检测系统
CN102243068A (zh) * 2011-04-27 2011-11-16 南京理工大学 一种子孔径拼接中系统误差的修正方法
CN103791854A (zh) * 2014-01-23 2014-05-14 中国科学院长春光学精密机械与物理研究所 高空间分辨率子孔径拼接方法
CN104776812A (zh) * 2015-04-14 2015-07-15 中国科学院上海光学精密机械研究所 一种提高平面面形子孔径拼接检测精度的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
安其昌 等: "基于功率谱的大口径望远镜面型评价", 《光子学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111649701A (zh) * 2020-06-30 2020-09-11 长春博信光电子有限公司 一种复曲面镜偏心值检测方法和装置
CN111649701B (zh) * 2020-06-30 2021-10-29 长春博信光电子有限公司 一种复曲面镜偏心值检测方法和装置

Also Published As

Publication number Publication date
CN107121114B (zh) 2019-05-07

Similar Documents

Publication Publication Date Title
CN102841344B (zh) 一种少阵元近场宽带信号源参数估计方法
CN111947045B (zh) 基于gvmd参数优化和奇异值分解的流体管道泄漏定位方法
Doroud et al. Gauge theory dynamics and Kähler potential for Calabi-Yau complex moduli
CN106123897A (zh) 基于多特征的室内融合定位方法
CN103596267A (zh) 一种基于欧氏距离的指纹图匹配方法
CN109407045A (zh) 一种非均匀传感器阵列宽带信号波达方向估计方法
CN104316049B (zh) 高精度低信噪比椭圆化星点光斑细分定位方法
CN105792356A (zh) 一种基于wifi的位置指纹定位方法
CN106501770A (zh) 基于幅相误差阵列的远近场宽带混合源中近场源定位方法
CN106199500A (zh) 指纹特征定位方法及装置
Kothandaraman et al. An adaptive ICA-based cross-correlation techniques for water pipeline leakage localization utilizing acousto-optic sensors
CN104616059B (zh) 一种基于量子粒子群的波达方向估计方法
CN103235888B (zh) 一种精确计算双各向异性介质球电磁散射的方法
CN102147458A (zh) 一种针对宽带声源的波达方向估计方法及其装置
CN103558498B (zh) 基于小波分析的绝缘子污闪泄漏电流信号稀疏表示方法
CN107341284A (zh) 高精度预测低频电波传播特性的双向抛物方程方法
CN107656240A (zh) 一种基于信号非圆特性的分步秩损远近场参数估计算法
CN113155973A (zh) 一种基于自适应奇异值分解的梁损伤识别方法
CN106295395A (zh) 基于图划分的不确定轨迹隐私保护方法
CN104780127A (zh) 一种时延-多普勒r-l解卷多路径信道估计方法
CN109725290A (zh) 一种误差提取方法、装置、电子设备及可读存储介质
CN107121114A (zh) 基于功率谱的大口径平面镜低阶像差估计方法
CN108120898B (zh) 基于vmd和sdeo的低采样率行波故障定位方法
CN108566256B (zh) 一种频谱地图的构建方法
Bao et al. Numerical reconstruction of electromagnetic inclusions in three dimensions

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant