[發(fā)明專利]一種大風(fēng)區(qū)電氣化鐵路接觸網(wǎng)正饋線舞動(dòng)定位算法

在審

• 申請(qǐng)?zhí)枺?/td>: 202110578103.6
• 申請(qǐng)日: 2021-05-26
• 公開（公告）號(hào): CN113324646A
• 公開（公告）日: 2021-08-31
• 發(fā)明（設(shè)計(jì)）人: 趙珊鵬;趙少翔;張友鵬;王思華;王文豪;王鵬;何亞萍;葛威;張海喜;馮強(qiáng);張永豐;金相龍;岳永文;張宸瑞;王天宇;趙全齊
• 申請(qǐng)（專利權(quán)）人: 蘭州交通大學(xué)
• 主分類號(hào): G01H17/00
• 分類號(hào): G01H17/00
• 代理公司: 蘭州錦知源專利代理事務(wù)所(普通合伙) 62204
• 代理人: 鐘國(guó)

1.一種大風(fēng)區(qū)電氣化鐵路接觸網(wǎng)正饋線舞動(dòng)定位算法，包括：

S10、獲取監(jiān)測(cè)點(diǎn)舞動(dòng)信號(hào)，得到x、y、z三軸對(duì)應(yīng)的信號(hào)分量；

S20、對(duì)步驟S10獲得的信號(hào)分量進(jìn)行分解，得到IMF(Intrinsic Mode Function,本征模態(tài)函數(shù))序列；

S30、對(duì)步驟S20的IMF序列提取趨勢(shì)項(xiàng)；

S40、將信號(hào)分量的原始信號(hào)減去步驟S30得到的趨勢(shì)項(xiàng)，得到修正的信號(hào)分量的波形；

S50、根據(jù)三軸信號(hào)分量的波形，獲得監(jiān)測(cè)點(diǎn)真實(shí)的舞動(dòng)軌跡。

2.如權(quán)利要求1所述的方法，其特征在于，

步驟S20包括：

S21、設(shè)定往積分信號(hào)中添加的高斯白噪聲的幅值系數(shù)k和分解次數(shù)M，并執(zhí)行初始次數(shù)m＝1，然后m＝m+1；

S22、在原信號(hào)x(t)中加入高斯白噪聲信號(hào)g(t)，得到待分解的混合信號(hào)x_i(t)：

x_i(t)＝x(t)+kg(t) (1)

S23、對(duì)x_i(t)進(jìn)行EMD(Empirical Mode Decomposition,經(jīng)驗(yàn)?zāi)B(tài)分解)分解，得到n個(gè)IMF序列，第m次分解得到的第i個(gè)IMF序列為：

imf_i,m(i＝1,2,...,n),n＝n+1 (2)

S24、若mM，重復(fù)執(zhí)行S22和S23；

S25、對(duì)M次分解得到的IMF序列求均值：

S26、EEMD(Ensemble Empirical Mode Decomposition，集合經(jīng)驗(yàn)?zāi)B(tài)分解)的最終分解結(jié)果為：

其中，δ為EEMD分解后的余項(xiàng)。

3.如權(quán)利要求1所述的方法，其特征在于，

步驟S30包括：

S31、計(jì)算原始信號(hào)x(t)的均值：

S32、計(jì)算第j個(gè)IMF序列至余項(xiàng)δ后n個(gè)信號(hào)的均值和：

S33、判斷當(dāng)兩均值之比c在誤差允許的范圍之內(nèi)時(shí)，確定j的值：

S34、判斷c是否滿足設(shè)定的條件，如果為真，則判斷結(jié)束，輸出第j個(gè)IMF序列至余項(xiàng)δ和為該信號(hào)的趨勢(shì)項(xiàng)r(t)：

如果結(jié)果為假，令j＝j(luò)-1，并轉(zhuǎn)入步驟S32，繼續(xù)判斷，直到滿足條件，退出循環(huán)。

4.如權(quán)利要求3所述的方法，其特征在于，

步驟S34中設(shè)定的條件是判斷c是否在0.95和1.05之間，c表示的是原始信號(hào)平均值與提取的趨勢(shì)項(xiàng)平均值的比值，c的取值范圍是根據(jù)大量研究結(jié)果而設(shè)置的經(jīng)驗(yàn)值。當(dāng)c在設(shè)定的0.95～1.05范圍內(nèi)時(shí)，可以表明提取的趨勢(shì)項(xiàng)是正確的。