发明名称 基于位置相关色散特性的快速低相干干涉解调方法
摘要 基于位置相关色散特性的快速低相干干涉解调方法,属于光纤传感领域,具体涉及大气压力的测量,也适用于三维形貌检测以及光学相干断层扫描技术等测量领域。本方法基于低相干干涉技术,使用法布里-珀罗(F-P)传感器作为传感干涉仪,通过腔长变化感受大气压力,使用起偏器、双折射光楔和检偏器构成的解调干涉仪,对F-P传感器调制的光信号进行解调,在零光程差的局部区域形成空间低相干干涉条纹。该方法采用低精度的质心法得到包络峰值位置,根据系统位置相关色散特性及相移法计算得到准确的中心波长干涉条纹中心条纹的峰值位置,以及该中心条纹的级次,从而恢复出准确的中心波长干涉条纹的零级条纹峰值位置,进而实现距离的快速高精度解调。
申请公布号 CN103759748B 申请公布日期 2015.12.30
申请号 CN201410007974.2 申请日期 2014.01.08
申请人 天津大学 发明人 刘铁根;江俊峰;王双;刘琨;尹金德;石俊锋;吴振海;张学智
分类号 G01D5/26(2006.01)I;G02F2/00(2006.01)I 主分类号 G01D5/26(2006.01)I
代理机构 天津佳盟知识产权代理有限公司 12002 代理人 侯力
主权项 一种基于位置相关色散特性的快速低相干干涉解调方法,其特征在于该方法包括:第1、确定低相干干涉条纹包络峰值位置,并解调出包络峰值位置对应的低精度距离;第1.1、通过质心法确定低精度的低相干干涉条纹包络峰值位置s′<sub>m</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>s</mi><mi>m</mi><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><mrow><mi>&Sigma;</mi><mo>&lsqb;</mo><mi>s</mi><mo>&times;</mo><mo>|</mo><mi>I</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>|</mo><mo>&rsqb;</mo></mrow><mrow><mi>&Sigma;</mi><mo>|</mo><mi>I</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>|</mo></mrow></mfrac></mrow>]]></math><img file="FDA0000813508210000011.GIF" wi="391" he="172" /></maths>其中,s为低相干干涉信号对应的位置信息,I(s)为s处的低相干干涉信号强度;第1.2、根据上面的包络峰值位置s′<sub>m</sub>计算低精度距离h<sub>m</sub>:h<sub>m</sub>=(1+αk<sub>0</sub>/n<sub>0</sub>)s′<sub>m</sub>/2其中,k<sub>0</sub>=2π/λ<sub>0</sub>是光源光谱的中心波数,λ<sub>0</sub>为光源光谱的中心波长,n<sub>0</sub>是中心波长对应的折射率,α是光源光谱范围内位置相关色散引起的折射率变化率,是一种色散特征属性;第2、根据位置相关色散特性判定中心波长干涉条纹的中心条纹级次m:m=int[s′<sub>m</sub>αk<sub>0</sub>/(n<sub>0</sub>λ<sub>0</sub>)];第3、根据位置相关色散特性确定中心波长干涉条纹的中心条纹峰值位置;第3.1、确定低相干干涉条纹中心条纹峰值位置s′<sub>cm</sub>:s′<sub>cm</sub>=(1+αk<sub>0</sub>/n<sub>0</sub>)s′<sub>m</sub>‑[m+ΔΦ(s′<sub>m</sub>)/(2π)]λ<sub>0</sub>其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>&Phi;</mi><mrow><mo>(</mo><msubsup><mi>s</mi><mi>m</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>a</mi><mi>r</mi><mi>c</mi><mi>t</mi><mi>a</mi><mi>n</mi><mo>&lsqb;</mo><msup><mi>&alpha;&Delta;k</mi><mn>2</mn></msup><msubsup><mi>s</mi><mi>m</mi><mo>&prime;</mo></msubsup><mo>/</mo><mrow><mo>(</mo><mn>4</mn><msub><mi>n</mi><mn>0</mn></msub><mi>l</mi><mi>n</mi><mn>2</mn><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000813508210000012.GIF" wi="759" he="155" /></maths>Δk为光源光谱波数半峰全宽;第3.2、以s′<sub>cm</sub>为初始位置,利用五步相移公式得到低相干干涉条纹中心条纹对应的相对相位的正切函数tan[φ(s)]:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>t</mi><mi>a</mi><mi>n</mi><mo>&lsqb;</mo><mi>&phi;</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>=</mo><mfrac><mrow><mn>2</mn><mo>&lsqb;</mo><mi>I</mi><mrow><mo>(</mo><mi>s</mi><mo>-</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>/</mo><mn>4</mn><mo>)</mo></mrow><mo>-</mo><mi>I</mi><mrow><mo>(</mo><mi>s</mi><mo>+</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>/</mo><mn>4</mn><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mrow><mn>2</mn><mi>I</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>-</mo><mi>I</mi><mrow><mo>(</mo><mi>s</mi><mo>+</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>-</mo><mi>I</mi><mrow><mo>(</mo><mi>s</mi><mo>-</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>/</mo><mn>2</mn><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000813508210000013.GIF" wi="896" he="155" /></maths>其中,φ(s)∈(‑π,π)为低相干干涉条纹中心条纹对应的相对相位,可以通过对tan[φ(s)]进行四象限反正切运算获得;第3.3、对已经确定的低相干干涉条纹中心条纹对应的相对相位φ(s)进行最小二乘线性拟合,拟合直线纵坐标为零的点所对应的横坐标值即为准确的低相干干涉条纹的中心条纹峰值位置s′<sub>c</sub>;第3.4、根据位置相关色散特性及已获得的s′<sub>c</sub>,计算得到准确的中心波长干涉条纹中心条纹的峰值位置s<sub>c</sub>:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>s</mi><mi>c</mi></msub><mo>=</mo><msubsup><mi>s</mi><mi>c</mi><mo>&prime;</mo></msubsup><mo>+</mo><mfrac><mrow><mi>a</mi><mi>r</mi><mi>c</mi><mi>t</mi><mi>a</mi><mi>n</mi><mo>&lsqb;</mo><msup><mi>&alpha;k</mi><mn>2</mn></msup><msubsup><mi>s</mi><mi>c</mi><mo>&prime;</mo></msubsup><mo>/</mo><mrow><mo>(</mo><mn>4</mn><msub><mi>n</mi><mn>0</mn></msub><mi>l</mi><mi>n</mi><mn>2</mn><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mrow><mn>2</mn><msub><mi>k</mi><mn>0</mn></msub></mrow></mfrac><mo>-</mo><mfrac><mrow><msub><mi>&eta;</mi><mi>s</mi></msub><msup><mrow><mo>&lsqb;</mo><mi>&Delta;</mi><mi>k</mi><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>&alpha;k</mi><mn>0</mn></msub><mo>/</mo><msub><mi>n</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msubsup><mi>s</mi><mi>c</mi><mo>&prime;</mo></msubsup><mo>-</mo><msubsup><mi>s</mi><mi>m</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mrow><mrow><mn>16</mn><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msubsup><mi>&eta;</mi><mi>s</mi><mn>2</mn></msubsup><mo>)</mo></mrow><msub><mi>k</mi><mn>0</mn></msub><mi>l</mi><mi>n</mi><mn>2</mn></mrow></mfrac></mrow>]]></math><img file="FDA0000813508210000014.GIF" wi="1373" he="187" /></maths>其中,η<sub>s</sub>=αΔk<sup>2</sup>s′<sub>c</sub>/(4n<sub>0</sub>ln2);第4、利用准确的中心波长干涉条纹的中心条纹峰值位置s<sub>c</sub>和级次m,恢复出准确的中心波长干涉条纹的零级条纹峰值位置s<sub>0</sub>:s<sub>0</sub>=s<sub>c</sub>+mλ<sub>0</sub>从而确定准确的距离信息s<sub>0</sub>/2,实现距离的高精度解调。
地址 300072 天津市南开区卫津路92号