发明名称 一种主被动结合的移相干涉波面测量方法
摘要 本发明为一种主被动结合的移相干涉波面测量方法,第一步:对被测口径内任意点引入环境干扰造成的被动移相量以及主动移相量;第二步:从光强-时间曲线中找到极大值和极小值,求得背景项和调制项;第三步:解得第i帧干涉图中点相位;第四步:得到该点在不同干涉图中的对应的相位信息;第五步:叠加各帧相位信息图求平均,得到被测的波面相位分布的测量结果,本发明的测量方法简单,对环境要求不高,无需特别的减振或隔绝气流措施;主被动结合,测量效率较高;有效消除环境干扰的影响,实现波面的高精度测量。
申请公布号 CN101639383B 申请公布日期 2010.12.29
申请号 CN200810188354.8 申请日期 2008.12.25
申请人 北京理工大学 发明人 郝群;朱秋东;汤磊;胡摇
分类号 G01J9/02(2006.01)I 主分类号 G01J9/02(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 张利萍
主权项 一种主被动结合的移相干涉波面测量方法,其特征在于:包括以下步骤:第一步:采集光强遍历干涉图序列,对被测口径内任意P(x,y)点引入主动移相量δa(x,y,t),同时包含环境干扰造成的被动移相量δp(x,y,t),对应总移相量为2 10个周期,主动移相条件为移相量单调增大且步长超过被动移相量,并且同一周期内主被动移相结合采样帧数超过100帧,则第i帧干涉图光强分布为I[x,y,δi(x,y,t)]=ia(x,y)+ib(x,y)·cos[Φ(x,y)+δi(x,y,t)] (1)其中,ia(x,y)是背景项,ib(x,y)为调制项,Φ(x,y)为被测波面的初始相位分布,δi(x,y,t)为该点在t时刻的移相量,即采集第i帧干涉图时的移相量,包括主被动两部分,δi(x,y,t)=δp(x,y,t)+δa(x,y,t);第二步:寻找干涉图各点光强在时间轴上的极值,满足第一步主动移相条件的同一点P(x,y)的光强值将随时间近似余弦规律变化,从该点光强 时间曲线中找到极大值Imax[x,y,δi(x,y,t)]和极小值Imin[x,y,δi(x,y,t)];第三步:求取干涉图各点背景项和调制项,第二步得到的极大极小值与公式(1)中的背景项和调制项满足公式(2)的关系,因此进一步利用公式(3)求得公式(1)中的背景项和调制项; <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>I</mi> <mi>max</mi> </msub> <mo>[</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <msub> <mi>&delta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> <mo>=</mo> <msub> <mi>i</mi> <mi>a</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>i</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mi>min</mi> </msub> <mo>[</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <msub> <mi>&delta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> <mo>=</mo> <msub> <mi>i</mi> <mi>a</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>i</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>a</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>I</mi> <mi>max</mi> </msub> <mo>[</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <msub> <mi>&delta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> <mo>+</mo> <msub> <mi>I</mi> <mi>min</mi> </msub> <mo>[</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <msub> <mi>&delta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> </mrow> <mn>2</mn> </mfrac> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>I</mi> <mi>max</mi> </msub> <mo>[</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <msub> <mi>&delta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> <mo>-</mo> <msub> <mi>I</mi> <mi>min</mi> </msub> <mo>[</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <msub> <mi>&delta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> </mrow> <mn>2</mn> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>第四步:求解各点相位,将公式(3)代入公式(1)中,利用反余弦函数可以解得第i帧干涉图中P(x,y)点相位如公式(4)所示 <mrow> <msub> <mi>&Phi;</mi> <mi>i</mi> </msub> <mo>[</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>&delta;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> <mo>=</mo> <mi>&Phi;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&delta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>=</mo> <mi>arccos</mi> <mo>{</mo> <mfrac> <mrow> <mn>2</mn> <mi>I</mi> <mo>[</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <msub> <mi>&delta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> <mo>-</mo> <msub> <mi>I</mi> <mi>max</mi> </msub> <mo>[</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>&delta;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> <mo>+</mo> <msub> <mi>I</mi> <mi>min</mi> </msub> <mo>[</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>&delta;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> </mrow> <mrow> <msub> <mi>I</mi> <mi>max</mi> </msub> <mo>[</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>&delta;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> <mo>-</mo> <msub> <mi>I</mi> <mi>min</mi> </msub> <mo>[</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>&delta;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> </mrow> </mfrac> <mo>}</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>将时间顺序记录下的各帧干涉图中同一点的光强值,分别带入到公式(4)所示的解相公式中,得到该点在不同干涉图中的对应的相位信息;第五步:叠加各帧波面分布信息求平均,得到被测的波面测量结果,如公式(5)所示: <mrow> <mi>&Phi;</mi> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>&prime;</mo> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mi>W</mi> <mo>{</mo> <msub> <mi>&Phi;</mi> <mi>i</mi> </msub> <mo>[</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>&delta;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> </mrow> <mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mi>W</mi> <mo>{</mo> <mi>&Phi;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&delta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mo>=</mo> <mi>&Phi;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mo>[</mo> <msub> <mi>&delta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>]</mo> </mrow> <mrow> <mo>=</mo> <mi>&Phi;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>L</mi> </mrow>其中,W为空间解包裹算符,Φ(x,y)′为被测波面相位分布的测量结果;Φ(x,y)为被测波面的实际相位分布;L为主动移相和环境振动带来的刚性移相量,对波面分布无影响。
地址 100081 北京市海淀区中关村南大街5号