发明名称 一种基于白光干涉测试系统的中心波长误差补偿方法
摘要 本发明公开了一种基于白光干涉测试系统的中心波长误差补偿方法,属于光纤端面测量仪技术领域。所述的解算方法首先通过相移器来改变参考光束和测试光束之间的光程差,同时采集若干幅干涉图,然后对干涉图进行求解得到光纤端面的高度。基于Carré相移算法对提取到的光强极大值进行相位修正,由于光源的中心波长存在标定误差而影响计算精度,因此首先使用单刻线标准样块校准压电陶瓷步距,通过校对后的步距对波长进行在线的标定。可以减小由于光源特性及环境扰动所带来的误差,提高计算精度。
申请公布号 CN105606338A 申请公布日期 2016.05.25
申请号 CN201610034498.2 申请日期 2016.01.19
申请人 北京航空航天大学 发明人 李慧鹏;谭朦曦;郑晓;宋凝芳;高爽;李皎
分类号 G01M11/00(2006.01)I 主分类号 G01M11/00(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 赵文颖
主权项 一种基于白光干涉测试系统的中心波长误差补偿方法,包括以下几个步骤:第一步,标定压电陶瓷马达步距;首先,通过单刻线标准样块标定压电陶瓷马达的步距,压电陶瓷马达的步距为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>d</mi><mo>=</mo><mfrac><mi>H</mi><mrow><mo>|</mo><msub><mi>N</mi><mn>1</mn></msub><mo>-</mo><msub><mi>N</mi><mn>2</mn></msub><mo>|</mo></mrow></mfrac></mrow>]]></math><img file="FDA0000909803170000011.GIF" wi="317" he="150" /></maths>其中,N<sub>1</sub>、N<sub>2</sub>为对同一行刻线内的点以及刻线外的点分别提取其光强极值对应的图像帧数,H为刻线深度;第二步,以标定好的步距在线实时标定白光光源的中心波长;在对被测物进行扫描时,提取到零级条纹所包含的图片帧数P,实时计算中心波长<img file="FDA0000909803170000012.GIF" wi="70" he="71" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mo>=</mo><mn>2</mn><mi>P</mi><mo>&times;</mo><mi>d</mi></mrow>]]></math><img file="FDA0000909803170000013.GIF" wi="270" he="70" /></maths>第三步,基于Carré相移算法利用提取到的波长解算相对高度;获取两束光振动合成的光的强度值:<img file="FDA0000909803170000014.GIF" wi="708" he="93" />式中:I<sub>a</sub>表示由光纤表面反射的光的光强值,I<sub>b</sub>表示由参考镜反射的光的光强值,I表示合成的光的强度值,<img file="FDA0000909803170000015.GIF" wi="50" he="55" />表示经过分光镜到达被测物的反射后的光与经过分光镜到达参考经反射后的光的位相差,设:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>I</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>=</mo><msub><mi>I</mi><mi>a</mi></msub><mo>+</mo><msub><mi>I</mi><mi>b</mi></msub><mo>+</mo><mn>2</mn><msqrt><mrow><msub><mi>I</mi><mi>a</mi></msub><msub><mi>I</mi><mi>b</mi></msub></mrow></msqrt></mrow>]]></math><img file="FDA0000909803170000016.GIF" wi="500" he="82" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>I</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>=</mo><msub><mi>I</mi><mi>a</mi></msub><mo>+</mo><msub><mi>I</mi><mi>b</mi></msub><mo>-</mo><mn>2</mn><msqrt><mrow><msub><mi>I</mi><mi>a</mi></msub><msub><mi>I</mi><mi>b</mi></msub></mrow></msqrt></mrow>]]></math><img file="FDA0000909803170000017.GIF" wi="552" he="93" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>&gamma;</mi><mo>=</mo><mfrac><mrow><msub><mi>I</mi><mi>max</mi></msub><mo>-</mo><msub><mi>I</mi><mi>min</mi></msub></mrow><mrow><msub><mi>I</mi><mi>max</mi></msub><mo>+</mo><msub><mi>I</mi><mi>min</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000909803170000018.GIF" wi="349" he="143" /></maths>其中,I<sub>max</sub>为光强极大值,I<sub>min</sub>为光强极小值,γ为干涉条纹对比度;令I<sub>0</sub>=I<sub>a</sub>+I<sub>b</sub>,得到:<img file="FDA0000909803170000019.GIF" wi="453" he="79" />当光程差为零时,即<img file="FDA00009098031700000110.GIF" wi="38" he="55" />为零时,光强度为峰值;对相位提取误差φ,进行修正,压电陶瓷步距为d,两相邻图像之间的相位变化为<img file="FDA0000909803170000021.GIF" wi="150" he="125" />则相位变化为:<img file="FDA0000909803170000022.GIF" wi="390" he="125" />其中,<img file="FDA0000909803170000023.GIF" wi="81" he="61" />为相位变化值,ε表示压电陶瓷线性误差;基于Carré相移算法对提取到的光强极大值进行相位修正,其基本原理如下:<img file="FDA0000909803170000024.GIF" wi="707" he="343" />其中:I<sub>bg</sub>为背景光强值,I<sub>‑3</sub>、I<sub>‑1</sub>、I<sub>1</sub>、I<sub>3</sub>分别为相对于所定位的零级条纹光强值前3帧、前1帧、后1帧、后3帧对应的光强值;解得<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msup><mi>tan</mi><mn>2</mn></msup><mi>&phi;</mi><mo>=</mo><mfrac><mrow><mo>|</mo><mo>&lsqb;</mo><mrow><mo>(</mo><msub><mi>I</mi><mrow><mo>-</mo><mn>3</mn></mrow></msub><mo>-</mo><msub><mi>I</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msub><mi>I</mi><mrow><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>I</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>&rsqb;</mo><mo>&lsqb;</mo><mn>3</mn><mrow><mo>(</mo><msub><mi>I</mi><mrow><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>I</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msub><mi>I</mi><mrow><mo>-</mo><mn>3</mn></mrow></msub><mo>-</mo><msub><mi>I</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>&rsqb;</mo><mo>|</mo></mrow><msup><mrow><mo>&lsqb;</mo><mrow><mo>(</mo><msub><mi>I</mi><mrow><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>I</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msub><mi>I</mi><mrow><mo>-</mo><mn>3</mn></mrow></msub><mo>+</mo><msub><mi>I</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mfrac></mrow>]]></math><img file="FDA0000909803170000025.GIF" wi="1364" he="167" /></maths>通过对相位的解算,得到被测物的相对高度:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>h</mi><mo>=</mo><mi>N</mi><mo>&times;</mo><mi>d</mi><mo>-</mo><mfrac><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mrow><mn>4</mn><mi>&pi;</mi></mrow></mfrac><mi>&phi;</mi></mrow>]]></math><img file="FDA0000909803170000026.GIF" wi="430" he="143" /></maths>其中,h为相对高度,N为光强极大值对应的图片帧数。
地址 100191 北京市海淀区学院路37号