[發(fā)明專利]一種基于FBG光譜奇偶分解的精確尋峰方法

有效

• 申請?zhí)枺?/td>: 201811260209.6
• 申請日: 2018-10-26
• 公開（公告）號: CN109341855B
• 公開（公告）日: 2020-10-30
• 發(fā)明（設計）人: 黃俊斌;丁朋;唐勁松;顧宏燦;劉文
• 申請（專利權）人: 中國人民解放軍海軍工程大學
• 主分類號: G01J3/28
• 分類號: G01J3/28
• 代理公司: 武漢宇晨專利事務所 42001
• 代理人: 王敏鋒

1.一種基于FBG光譜奇偶分解的精確尋峰方法，其特征在于：在所述方法中，具體包括以下步驟：

1)把光譜信號波長所對應的幅值生成一元原始數(shù)列，記錄為x(n),判斷數(shù)列的長度N是奇數(shù)還是偶數(shù)；設置循環(huán)計算的控制參數(shù)i，i為正整數(shù)，賦予i初值為1；

若N是偶數(shù)：

2.1)在x(n)前添加N-2i+1個0元素，生成第一新數(shù)列，記為y(m)；翻轉(zhuǎn)y(m)，生成第二新數(shù)列，記為yf(m)；數(shù)列y(m)減去數(shù)列yf(m)，生成第三新數(shù)列，記為yof(m)；尋找yof(m)的最大值，最大值記為c(i)；令循環(huán)計算的控制參數(shù)i加上1，并把計算得到的值賦予i；

3.1)執(zhí)行步驟2.1)，直到i大于N/2；就獲得了N/2個c(i)，組成第四新數(shù)列記為c(n/2)；令循環(huán)計算的控制參數(shù)i等于1；

4.1)在x(n)后添加2i-1個0元素，生成第五新數(shù)列，記為z(m)；翻轉(zhuǎn)z(m)，生成第六新數(shù)列，記為zf(m)；數(shù)列z(m)減去數(shù)列zf(m)，生成第七新數(shù)列，記為zof(m)；尋找zof(m)的最大值，最大值記為d(i)；令循環(huán)計算的控制參數(shù)i加1，并把計算得到的值賦予i；

5.1)執(zhí)行步驟4.1)，直到i大于N/2；就獲得了N/2個d(i)，組成第八新數(shù)列記為d(n/2)；

6.1)把第四新數(shù)列記c(n/2)和第八新數(shù)列d(n/2)合成第九新數(shù)列，記為[c(n/2)d(n/2)]，尋找這個數(shù)列中的最小值，所對應的位置即為FBG的中心波長；

若N是奇數(shù)：

2.2)在x(n)前添加N-2i+2個0元素，生成第十新數(shù)列，記為y(m)；翻轉(zhuǎn)y(m)，生成第十一新數(shù)列，記為yf(m)；數(shù)列y(m)減去數(shù)列yf(m)，生成第十二新數(shù)列，記為yof(m)；尋找yof(m)的最大值，最大值記為c(i)；令循環(huán)計算的控制參數(shù)i加1，并把計算得到的值賦予i；

3.2)執(zhí)行步驟2.2)，直到i大于(N-1)/2；就獲得了(N-1)/2個c(i)，組成十三新數(shù)列，記為c((n-1)/2)；設置循環(huán)計算的控制參數(shù)i，賦予i初值為1；

4.2)在x(n)后添加2i個0元素，生成第十四新數(shù)列，記為z(m)；翻轉(zhuǎn)z(m)，生成第十五新數(shù)列，記為zf(m)；數(shù)列z(m)減去數(shù)列zf(m)，生成第十六新數(shù)列，記為zof(m)；尋找zof(m)的最大值，最大值記為d(i)；令循環(huán)計算的控制參數(shù)i加1，并把計算得到的值賦予i；

5.2)執(zhí)行步驟4.2)，直到i大于(N-1)/2；就獲得了(N-1)/2個d(i)，組成第十七新數(shù)列，記為d((n-1)/2)；

6.2)翻轉(zhuǎn)x(n)，生成第十八新數(shù)列，記為xf(n)；數(shù)列x(n)減去數(shù)列xf(n)，生成第十九新數(shù)列，記為xof(n)；尋找xof(n)的最大值，最大值記錄為xn_mid；

7.2)把數(shù)列c((n-1)/2)、xn_mid和數(shù)列d((n-1)/2)，合成第二十新數(shù)列，記為[c((n-1)/2)xn_mid d((n-1)/2)]，尋找這個數(shù)列中的最小值，所對應的位置即為FBG的中心波長。