发明名称 施瓦茨光学系统波像差测量方法
摘要 一种施瓦茨光学系统波像差测量方法,该方法基于光栅横向剪切干涉测量装置,测量时将待测施瓦茨光学系统置于所述光栅横向剪切干涉测量装置的针孔掩模和剪切光栅之间,测量步骤如下:①获取x方向和y方向的横向剪切干涉图;②获取x方向和y方向的差分波前;③获取待测波前的环Zernike多项式系数。本发明的优点是将差分波前拟合得到差分环Zernike多项式,直接得到待测环形波前的环Zernike多项式系数,通过所述的前J项环Zernike多项式系数<img file="DDA0000567151120000011.GIF" wi="28" he="40" />与环Zernike多项式相乘再求和,即得到待测施瓦茨光学系统的波像差。本发明简化了波前重建过程,提高了波前测量精度和速度。
申请公布号 CN104236855A 申请公布日期 2014.12.24
申请号 CN201410455642.0 申请日期 2014.09.09
申请人 中国科学院上海光学精密机械研究所 发明人 戴凤钊;王向朝;蔡燕民
分类号 G01M11/02(2006.01)I 主分类号 G01M11/02(2006.01)I
代理机构 上海新天专利代理有限公司 31213 代理人 张泽纯;张宁展
主权项 一种施瓦茨光学系统波像差测量方法,该方法基于光栅横向剪切干涉测量装置,用于施瓦茨光学系统波像差测量,所述的光栅横向剪切干涉测量装置包括光源(1),沿该光源(1)输出的光束方向依次为针孔掩模(2)、剪切光栅(4)、双窗口掩模(5)和光电探测器(6),该光电探测器(6)的输出端与计算机(7)的输入端相连,测量时将待测施瓦茨光学系统(3)置于所述的针孔掩模(2)和剪切光栅(4)之间,特征在于该方法的测量步骤如下:①获取x方向和y方向的横向剪切干涉图:调整双窗口掩模(5),使所述的光栅横向剪切干涉测量装置的光轴通过该双窗口掩模的任意一个窗口的中心,并使两窗口中心连线平行于x轴,且该双窗口掩模的表面与所述的光栅横向剪切干涉测量装置的光轴垂直;启动光源(1),通过光电探测器(6)采集待测施瓦茨光学系统输出波前x方向的横向剪切干涉图,并输入至计算机保存;将剪切光栅(4)及双窗口掩模(5)旋转90°,通过光电探测器(6)采集y方向的横向剪切干涉图,并输入计算机保存;②获取x方向和y方向的差分波前:将所述的x方向的横向剪切干涉图通过傅里叶变换相位提取方法获取x方向的差分波前;将该x方向的差分波前采用补零法拓展为矩形区域,将所述的矩形区域离散点的值定义为M<sub>x</sub>矩阵,将M<sub>x</sub>矩阵的行数记为N,列数记为m;将所述的y方向的横向剪切干涉图通过傅里叶变换相位提取方法获取y方向的差分波前;将所述的y方向的差分波前采用补零法拓展为矩形区域,将该矩形区域离散点的值定义为M<sub>y</sub>矩阵,将M<sub>y</sub>矩阵的行数记为k,列数记为N;③获取待测波前的环Zernike多项式系数:在计算机中设置N×N的第一全零矩阵M<sub>1</sub>,用所述的M<sub>x</sub>矩阵所有m列的元素逐个替换第一全零矩阵M<sub>1</sub>的第1列至第m列的所有元素,获得x方向差分波前矩阵;将所述的x方向的差分波前矩阵采取逐行从左至右排序的方式变为列向量,将该列向量记为ΔW<sub>x</sub>;在计算机中设置一个N×N的第二全零矩阵M<sub>2</sub>,用所述的M<sub>y</sub>矩阵的所有行的元素逐个替换第二全零矩阵M<sub>2</sub>的第1行至第k行的所有元素,获得y方向的差分波前矩阵,将所述的y方向的差分波前矩阵采取逐行从左至右排序的方式变为列向量,将该列向量记为ΔW<sub>y</sub>;由所述的列向量ΔW<sub>x</sub>和ΔW<sub>y</sub>,利用最小二乘拟合法计算待测施瓦茨光学系统的环形波前的环Zernike多项式系数<img file="FDA0000567151090000011.GIF" wi="64" he="54" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>a</mi><mo>^</mo></mover><msup><mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mfenced open='(' close=')'><mtable><mtr><mtd><mi>&Delta;</mi><msub><mi>H</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>H</mi><mi>y</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mtd><mtd><mfenced open='(' close=')'><mtable><mtr><mtd><mi>&Delta;</mi><msub><mi>H</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>H</mi><mi>y</mi></msub></mtd></mtr></mtable></mfenced></mtd></mtr></mtable></mfenced></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mfenced open='(' close=')'><mtable><mtr><mtd><mi>&Delta;</mi><msub><mi>H</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>H</mi><mi>y</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mfenced open='(' close=')'><mtable><mtr><mtd><mi>&Delta;</mi><msub><mi>W</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>W</mi><mi>y</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000567151090000021.GIF" wi="701" he="159" /></maths>式中:ΔH<sub>x</sub>和ΔH<sub>y</sub>均为N<sup>2</sup>×(J‑1)矩阵,分别表示x方向和y方向除第1项外的前J项差分环Zernike多项式在所有N<sup>2</sup>点上的值,通过将x、y的坐标值、由装置决定的x和y方向的剪切率以及待测施瓦茨光学系统光瞳的中心遮拦比代入差分环形Zernike多项式计算得到;<img file="FDA0000567151090000022.GIF" wi="41" he="53" />为J‑1列向量,即表示待测波前除第1项外的前J项环Zernike多项式系数,<img file="FDA0000567151090000023.GIF" wi="349" he="78" />T表示转置,通过所述的前J项环Zernike多项式系数<img file="FDA0000567151090000024.GIF" wi="39" he="53" />与环Zernike多项式相乘再求和,即得到待测施瓦茨光学系统的波像差。
地址 201800 上海市嘉定区上海市800-211邮政信箱