发明名称 用于准确测量布里渊谱的频移的最大相似匹配分析方法
摘要 本发明公开了一种用于准确测量布里渊谱频移的最大相似匹配分析方法,预先测得光纤在未受应力和温度状态下的原始布里渊谱作为参考谱,将实际测量得到的布里渊谱作为测量谱,逐渐移动测量谱并比较其与截取的部分参考谱的形状,利用相关系数描述两者的相似程度,当计算出的相关系数最大时,此时测量谱的移动值即为此时光纤布里渊频移的改变值。使用本发明分析方法可更加准确测量最终的频移值,使测量的结果更加精确;克服了洛伦兹拟合只能在布里渊谱为单峰时有效的缺点,同时该方法运算简单,简化了数据分析的难度,测量时扫频较窄,频谱采样步进较小,减少了测量的时间。
申请公布号 CN103968864B 申请公布日期 2016.04.06
申请号 CN201410167156.9 申请日期 2014.04.23
申请人 南京大学 发明人 王峰;张旭苹;詹伟伟
分类号 G01D5/32(2006.01)I;G01D5/353(2006.01)I 主分类号 G01D5/32(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 朱小兵
主权项 用于准确测量布里渊谱的频移的最大相似匹配分析方法,其特征在于,预先测量光纤在未受应力和温度状态下的原始布里渊谱作为参考谱,将实际测量得到的布里渊谱作为测量谱,逐渐移动测量谱并比较其与截取的部分参考谱的形状,利用相关系数描述两者的形状相似程度,当相关系数取最大值时,测量谱的移动值即为实际的频移改变值;具体包括以下步骤:步骤一、在光纤未受应力和温度改变的情况下,布里渊光纤传感器测量光纤的布里渊信号,得到的布里渊谱作为参考谱为{G(f<sub>0</sub>+mF<sub>0</sub>)|m=0,1,...,M‑1},其中f<sub>0</sub>为参考谱的扫频的初始频率,M为参考谱的扫频的总数,F<sub>0</sub>为参考谱的扫频间隔;步骤二、将实际测量得到的布里渊谱作为测量谱,测量谱为{g(f+n·F(n)|n=0,1,...,N‑1},其中f为实际测量时扫频的初始频率,N为实际测量时扫描的频率总数,F(n)为实际测量的扫频间隔,应为参考谱的扫频间隔F<sub>0</sub>的整数倍,对于不同的扫频次数,它可以取相同或不同的值,测量谱的扫频范围在参考谱的扫频范围之内;步骤三、移动测量谱的位置,使该测量谱的初始频率对应于参考谱的某一频率f<sub>j</sub>,从参考谱中按照测量谱的扫频点截取部分参考谱,将f<sub>j</sub>作为截取的部分参考谱的初始频率,截取的部分参考谱为{G(f<sub>j</sub>+n·F(n))|n=0,1,...,N‑1};步骤四、将步骤三中截取的部分参考谱与测量谱按照如下公式进行计算得到相关系数C(f<sub>j</sub>):<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>C</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mrow><mo>(</mo><msub><mi>G</mi><mi>j</mi></msub><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo><mo>-</mo><mover><msub><mi>G</mi><mi>j</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mrow><mo>(</mo><mi>g</mi><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo><mo>-</mo><mover><mi>g</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow><mrow><mo>{</mo><msup><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>(</mo><mrow><msub><mi>G</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mover><msub><mi>G</mi><mi>j</mi></msub><mo>&OverBar;</mo></mover></mrow><mo>)</mo></mrow><mi>2</mi></msup><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>(</mo><mrow><mi>g</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mover><mi>g</mi><mo>&OverBar;</mo></mover></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo><mo>}</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup></mrow></mfrac></mrow>]]></math><img file="FDA0000858428580000011.GIF" wi="1030" he="319" /></maths>其中f<sub>i</sub>为截取的部分参考谱中的某一频率,G<sub>j</sub>(f<sub>i</sub>)为初始频率为f<sub>j</sub>的截取的部分参考谱在频率f<sub>i</sub>位置处的值,g(f<sub>i</sub>)为测量谱在频率f<sub>i</sub>位置处的值,其中<img file="FDA0000858428580000012.GIF" wi="71" he="84" />和<img file="FDA0000858428580000013.GIF" wi="46" he="77" />分别为截取的部分参考谱和测量谱的平均值;步骤五、改变截取的部分参考谱的初始频率f<sub>j</sub>,截取不同的参考谱部分,并依次计算其与测量谱的相关系数,当相关系数为最大时,此时测量谱的移动值即为此时光纤布里渊频移的改变值f<sub>BS</sub>,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>f</mi><mrow><mi>B</mi><mi>S</mi></mrow></msub><mo>=</mo><mi>f</mi><mo>-</mo><munder><mi>arg</mi><msub><mi>f</mi><mi>j</mi></msub></munder><mi>max</mi><mi> </mi><mi>C</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000858428580000021.GIF" wi="517" he="107" /></maths>
地址 210093 江苏省南京市鼓楼区汉口路22号