基于PSO小波神經(jīng)網(wǎng)絡(luò)輔助卡爾曼濾波的BDS/INS定位

李衛(wèi)東，黃晨陽(yáng)，劉 楊，王新屏

(大連交通大學(xué)電氣信息學(xué)院，遼寧 大連 116023)

0 引言

在北斗定位系統(tǒng)(Beidou navigation satellite system，BDS)中，卡爾曼濾波是一種通過(guò)軟件消除定位信號(hào)噪聲的重要方法[1]。假如對(duì)列車(chē)運(yùn)行模型構(gòu)造不夠精確，會(huì)導(dǎo)致卡爾曼濾波后信號(hào)的發(fā)散、失真等不良結(jié)果[2]。通過(guò)構(gòu)建卡爾曼濾波的輸入/輸出小波神經(jīng)網(wǎng)絡(luò)，可以解決這些問(wèn)題[3]。

本文結(jié)合文獻(xiàn)[1]～文獻(xiàn)[3]的優(yōu)點(diǎn)，改進(jìn)了小波神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)率。在此基礎(chǔ)上，采用粒子群優(yōu)化(particle swarm optimization，PSO)算法結(jié)合遺傳算法(genetic algorithm，GA)優(yōu)化的小波神經(jīng)網(wǎng)絡(luò)對(duì)實(shí)際問(wèn)題進(jìn)行建模。仿真試驗(yàn)證明了該方法的有效性。

1 基于卡爾曼濾波的BDS/INS組合定位

目前，我國(guó)BDS總體性能略低于全球定位系統(tǒng)(global positioning system,GPS),它的主要問(wèn)題可能在于:定位精度低、易受到無(wú)線電的干擾、保密性和安全性弱、 數(shù)據(jù)傳輸率低。這些缺陷可以通過(guò)與慣性導(dǎo)航系統(tǒng)(inertial navigation system，INS)結(jié)合來(lái)補(bǔ)償[4]。BDS/INS可以使北斗信號(hào)在由于遮擋或高強(qiáng)度干擾的情況下，能單獨(dú)進(jìn)行導(dǎo)航定位，從而保證導(dǎo)航的可靠性。北斗系統(tǒng)可以修正INS的輸出信息，減少它隨時(shí)間積累的誤差，保證導(dǎo)航的精確性和工作的穩(wěn)定性。

卡爾曼濾波是一種通過(guò)算法估計(jì)被提取信號(hào)的最優(yōu)理論。自1980年以來(lái)，它被廣泛應(yīng)用于運(yùn)載體的導(dǎo)航和定位系統(tǒng)中[5]。

它的原理可以簡(jiǎn)述為：假設(shè)Xk為tk時(shí)刻的被估計(jì)狀態(tài)，Wk-1系統(tǒng)的驅(qū)動(dòng)噪聲序列驅(qū)動(dòng)方程如下：

Xk=φk,k-1Xk-1+Γk-1Wk-1

(1)

對(duì)Xk的量測(cè)方程滿足線性關(guān)系：

Zk=HxXk+Vk

(2)

狀態(tài)一步預(yù)測(cè)：

(3)

狀態(tài)估計(jì)：

(4)

濾波增益：

(5)

一步預(yù)測(cè)均方誤差：

(6)

估計(jì)均方誤差：

Pk=(I-KkHk)Pk/k-1

(7)

北斗定位系統(tǒng)與列車(chē)和地面設(shè)備時(shí)鐘一致，并向其上的接收機(jī)發(fā)送用于導(dǎo)航定位的調(diào)制波。導(dǎo)航電文發(fā)送隨機(jī)偽噪聲。理想情況下，接收機(jī)接收到電文后，將它與衛(wèi)星信號(hào)的原電文進(jìn)行對(duì)齊處理，求得衛(wèi)星與接收機(jī)的距離和列車(chē)的相關(guān)參數(shù)[6]。INS器件通過(guò)測(cè)量移動(dòng)載體的運(yùn)動(dòng)加速度及位移，再通過(guò)積分運(yùn)算，推導(dǎo)出列車(chē)相關(guān)參數(shù)。隨后，BDS系統(tǒng)中的接收機(jī)將收到的相關(guān)列車(chē)參數(shù)，經(jīng)電平轉(zhuǎn)換后，通過(guò)通信電路傳輸?shù)叫盘?hào)處理器進(jìn)行處理。同時(shí)在數(shù)字信號(hào)的裝置中將其中的慣性導(dǎo)航器件處理的信號(hào)導(dǎo)入，兩者信息相融合。

圖1 BDS/INS結(jié)構(gòu)圖

2 小波神經(jīng)網(wǎng)絡(luò)輔助的卡爾曼濾波

2.1 卡爾曼濾波的不足及彌補(bǔ)方法

實(shí)際運(yùn)行中，由于列車(chē)運(yùn)動(dòng)模型不能準(zhǔn)確構(gòu)造以及系統(tǒng)中存在硬件誤差等原因，會(huì)導(dǎo)致經(jīng)過(guò)卡爾曼濾波，獲得的速度波形與實(shí)際速度波形有偏差，進(jìn)而使得到的定位信息也不準(zhǔn)確[8]。雖然這些誤差的精確模型很難計(jì)算，但可以通過(guò)現(xiàn)有的卡爾曼濾波的輸入/輸出數(shù)據(jù)構(gòu)造神經(jīng)網(wǎng)絡(luò)模型，再對(duì)這個(gè)神經(jīng)網(wǎng)絡(luò)模型進(jìn)行校正。通過(guò)這種方式讓北斗定位系統(tǒng)獲取到較為準(zhǔn)確的速度波形，從而提高北斗定位的精度。

采用圖2所示的神經(jīng)網(wǎng)絡(luò)監(jiān)督控制替換卡爾曼濾波。通過(guò)在線比較誤差的方法對(duì)小波神經(jīng)網(wǎng)絡(luò)進(jìn)行修正，逐步替換反饋控制器。在受到外界噪聲或其他因素干擾的情況下，控制器的反饋將會(huì)再次發(fā)揮作用。神經(jīng)網(wǎng)絡(luò)監(jiān)督控制方法，不但確保了控制系統(tǒng)能夠比較穩(wěn)定地運(yùn)行，還有效提高了系統(tǒng)的精度及其自適應(yīng)的能力[9]。

圖2 神經(jīng)網(wǎng)絡(luò)監(jiān)督控制圖

2.2 小波神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

小波理論本身比較適合處理[10-11]小波神經(jīng)網(wǎng)絡(luò)(wavelet neural network，WNN)。與神經(jīng)網(wǎng)絡(luò)理論結(jié)合后，小波神經(jīng)網(wǎng)絡(luò)算法具有更加靈活有效的函數(shù)逼近能力以及較強(qiáng)的容錯(cuò)能力。其核心是將小波變換加入神經(jīng)網(wǎng)絡(luò)模型中,構(gòu)造以“Sigmoid函數(shù)”為基本構(gòu)架的小波神經(jīng)網(wǎng)絡(luò)。在構(gòu)造小波神經(jīng)網(wǎng)絡(luò)時(shí)，因?yàn)樾〔ǚ纸饽軌蜻M(jìn)行平移，理論上能夠逼近任意函數(shù)。

本文根據(jù)信號(hào)的特點(diǎn),采用了融合型的小波神經(jīng)網(wǎng)絡(luò)。該網(wǎng)絡(luò)主體為3層網(wǎng)絡(luò)[12],具有單輸入、單輸出的結(jié)構(gòu)，輸入層為并行輸入，隱含層、輸出層包含了m個(gè)神經(jīng)元。

小波神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)如圖3所示小波函數(shù)中的伸縮，平移參數(shù)a、b和神經(jīng)網(wǎng)絡(luò)參數(shù)都可以通過(guò)網(wǎng)絡(luò)學(xué)習(xí)進(jìn)行訓(xùn)練。圖3中，隱含層神經(jīng)元的激勵(lì)函數(shù)f1是 Morlet 小波函數(shù)：

圖3 小波神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)圖

(8)

隱含層輸出為：

(9)

(10)

式中：O為隱含層的輸出；x為輸入層的輸入；w為權(quán)值；h為Morlet小波函數(shù)。

輸出層輸出為：

(11)

2.3 改進(jìn)的小波神經(jīng)網(wǎng)絡(luò)訓(xùn)練算法

本文采用了改進(jìn)的梯度下降法[13]，對(duì)小波神經(jīng)網(wǎng)絡(luò)權(quán)值和參數(shù)進(jìn)行訓(xùn)練，誤差函數(shù)為：

(12)

式中:d(k)為期望輸出；y(k)為網(wǎng)絡(luò)實(shí)際輸出；p為樣本數(shù)。

首先對(duì)輸入數(shù)據(jù)進(jìn)行分割，并行輸入，得到權(quán)值平均值。這樣權(quán)值的反復(fù)變化次數(shù)減少，不但節(jié)省了時(shí)間，也提高了系統(tǒng)的穩(wěn)定性。神經(jīng)網(wǎng)絡(luò)主要基于誤差修正來(lái)進(jìn)行結(jié)構(gòu)的調(diào)整和學(xué)習(xí),學(xué)習(xí)率的改變會(huì)影響收斂的速度和訓(xùn)練的效果。如果學(xué)習(xí)率過(guò)小,系統(tǒng)的調(diào)節(jié)時(shí)間會(huì)太長(zhǎng);如果學(xué)習(xí)率太大,可能導(dǎo)致振蕩或發(fā)散。隨著時(shí)間增加，系統(tǒng)的調(diào)節(jié)幅度應(yīng)相應(yīng)減少。調(diào)整后的學(xué)習(xí)率單調(diào)增加且增加率逐漸變緩，log(a+x)無(wú)疑滿足這個(gè)條件。

學(xué)習(xí)率的調(diào)整式：

ηk+1=ηklog(1+k)e(k+1)

(13)

ηk+1=ηklog(2+1/k)e(k+1)≥e(k)

(14)

式中：η為學(xué)習(xí)因子；λ為動(dòng)量因子。

改進(jìn)的小波神經(jīng)網(wǎng)絡(luò)訓(xùn)練算法，一定程度上加快了小波神經(jīng)網(wǎng)絡(luò)的收斂速度。但它仍然較易陷入局部最小，這無(wú)疑會(huì)對(duì)系統(tǒng)的收斂速度和精確度造成影響[13]。下面采用PSO與GA的組合算法來(lái)對(duì)小波神經(jīng)網(wǎng)絡(luò)的權(quán)值和閾值進(jìn)行修正，以進(jìn)一步提升其性能。

3 粒子群優(yōu)化算法優(yōu)化的小波神經(jīng)網(wǎng)絡(luò)

1995年，美國(guó)的Kenney和Eeberhart對(duì)鳥(niǎo)類(lèi)的群體行為進(jìn)行探究，提出了PSO算法[11]。每個(gè)粒子都可以看作潛在的要解決的優(yōu)化問(wèn)題的解決方案。粒子的速度取決于它的飛行的距離及其方向，它的適應(yīng)值由優(yōu)化函數(shù)決定。在飛行過(guò)程中，粒子經(jīng)過(guò)多次迭代后，利用自身和飛行經(jīng)驗(yàn)組進(jìn)行動(dòng)態(tài)調(diào)整，最后取得最優(yōu)解。粒子本身所得到的最佳值可以稱(chēng)作個(gè)體的極值，而整個(gè)種群的最優(yōu)解應(yīng)是當(dāng)前的極值。這一極值會(huì)在每個(gè)粒子尋找這兩個(gè)最佳值時(shí)更新：

vk+1=c0+vkc1(pbestk-xk)+c2(gbestk-xk)

(15)

xk+1=xk+vk+1

(16)

式中：vk為粒子的速度矢量；xk為目前粒子所在的位置；c0一般取介于(0，1)之間的隨機(jī)常數(shù)，通常可以稱(chēng)為慣性系數(shù)或收縮因子。采用c1和c2的“自感知因子”和“社會(huì)認(rèn)知因子”來(lái)調(diào)節(jié)pbestk和gbestk對(duì)粒子的影響。

GA的本質(zhì)是構(gòu)建一個(gè)針對(duì)優(yōu)化問(wèn)題的目標(biāo)函數(shù)的適合度函數(shù)[12]。然后，根據(jù)一定的規(guī)則，生成遺傳密碼的初始組，并對(duì)該組進(jìn)行評(píng)估，進(jìn)行遺傳算法的交叉、變異、選擇等。

遺傳算法可以在不使算法變復(fù)雜的同時(shí)，提高算法的全局極值，從而提高算法的收斂速度和精度。在這種情況下，兩者優(yōu)勢(shì)互補(bǔ)。它們的融合可以在提升學(xué)習(xí)速度的同時(shí)，增強(qiáng)系統(tǒng)的泛化能力。

PSO優(yōu)化小波神經(jīng)網(wǎng)絡(luò)流程如圖4所示。

圖4 PSO優(yōu)化小波神經(jīng)網(wǎng)絡(luò)流程圖

4 仿真分析

對(duì)所提算法進(jìn)行仿真,試驗(yàn)數(shù)據(jù)來(lái)自哈爾濱西至長(zhǎng)春運(yùn)營(yíng)線路的實(shí)測(cè)數(shù)據(jù)，列車(chē)最高時(shí)速為200 km/h，平均車(chē)速為99.90 km/h[14]。設(shè)置北斗衛(wèi)星接收機(jī)數(shù)據(jù)采樣頻率為1 s/次，哈爾濱西站的坐標(biāo)作為初始位置的經(jīng)緯度坐標(biāo)，共收集了400個(gè)北斗Ⅱ衛(wèi)星的測(cè)量數(shù)據(jù)(列車(chē)實(shí)際運(yùn)行13.286 km)。為了便于分析，只采用了部分?jǐn)?shù)據(jù)點(diǎn)進(jìn)行仿真分析,并經(jīng)由傳統(tǒng)的卡爾曼濾波獲得輸出波形及數(shù)據(jù)。

在試驗(yàn)1中，小波神經(jīng)網(wǎng)絡(luò)采用并行輸入[15]，輸入和輸出的節(jié)點(diǎn)數(shù)都為1，是3層神經(jīng)網(wǎng)絡(luò)。其中隱層節(jié)點(diǎn)數(shù)為12，輸出層節(jié)點(diǎn)數(shù)為5。如果對(duì)數(shù)據(jù)分塊太多，則運(yùn)行時(shí)間過(guò)長(zhǎng)；如果數(shù)據(jù)分塊太少，則會(huì)使誤差偏大，所以將輸入數(shù)據(jù)分成了4塊。運(yùn)行時(shí)間及均方差對(duì)比如表1所示。學(xué)習(xí)率改變前后誤差曲線如圖5所示。

表1 運(yùn)行時(shí)間及均方差對(duì)比

圖5 學(xué)習(xí)率改變前后誤差曲線

在試驗(yàn)2中，遺傳代數(shù)選擇了300代，種群規(guī)模選擇了50，變異率為0.95。算法適應(yīng)度曲線如圖6所示。從圖6可以看出，組合算法在收斂速度上優(yōu)于粒子群算法，且遠(yuǎn)遠(yuǎn)優(yōu)于遺傳算法。

圖6 算法適應(yīng)度曲線

在試驗(yàn)3中，采用了北向的數(shù)據(jù)集，對(duì)基于GA改進(jìn)的PSO權(quán)值部分進(jìn)行訓(xùn)練。當(dāng)精度達(dá)到要求后結(jié)束訓(xùn)練,進(jìn)入在線狀態(tài)。在MATLAB 2015B中運(yùn)行程序，對(duì)WNN輔助的卡爾曼濾波器所獲得的速度曲線和原來(lái)實(shí)際運(yùn)行的曲線,以及它們的誤差曲線進(jìn)行比較。優(yōu)化前后誤差曲線、總體誤差對(duì)比曲線如圖7、圖8所示。

圖7 優(yōu)化前后誤差對(duì)比曲線

從圖7、圖8可以看出，與實(shí)際運(yùn)動(dòng)曲線相比，經(jīng)過(guò)改進(jìn)的小波神經(jīng)網(wǎng)絡(luò)輔助的卡爾曼濾波，對(duì)于速度的估計(jì)誤差更小,而且經(jīng)過(guò)長(zhǎng)時(shí)間的運(yùn)行，偏差仍然不大,表明了通過(guò)小波神經(jīng)網(wǎng)絡(luò)對(duì)卡爾曼濾波器輔助的方法是可行的。經(jīng)過(guò)PSO優(yōu)化后的神經(jīng)網(wǎng)絡(luò)具有運(yùn)行穩(wěn)定且誤差小和收斂快等優(yōu)點(diǎn)。在PSO優(yōu)化數(shù)據(jù)后，相比優(yōu)化前WNN 中100個(gè)輸出數(shù)據(jù)與理想輸出之差有所減少，整體的均方差也由原來(lái)的0.150 5降低到了0.019 4。

5 結(jié)束語(yǔ)

本文提出了一種采用WNN對(duì)實(shí)際問(wèn)題中的卡爾曼濾波進(jìn)行建模的方案。在此基礎(chǔ)上，對(duì)GA和PSO進(jìn)行了組合優(yōu)化。學(xué)習(xí)率則是通過(guò)加入對(duì)數(shù)函數(shù)進(jìn)行改進(jìn)。仿真表明，改進(jìn)后的系統(tǒng)不僅加快了運(yùn)行速度，而且減小了均方差。小波神經(jīng)網(wǎng)絡(luò)優(yōu)化不但提高了收斂速度,還有效提高了其定位的準(zhǔn)確性，同時(shí)在保持原速度波形的整體趨勢(shì)前提下，減少了濾波導(dǎo)致的信號(hào)誤差。在今后工作中，考慮在BDS/INS慣性導(dǎo)航系統(tǒng)失效的情況下，利用神經(jīng)網(wǎng)絡(luò)和其他算法融合進(jìn)行運(yùn)行路線的預(yù)測(cè)。

[1] 秦永元，張洪鉞，汪叔華.卡爾曼濾波與組合導(dǎo)航原理 [M].2版.西安：西北工業(yè)大學(xué)出版社，2012.

[2] 韓曉捷,陳德巧.基于神經(jīng)網(wǎng)絡(luò)酌高速列車(chē)位置計(jì)算模型在線學(xué)習(xí)算法設(shè)計(jì)[D].北京：北京交通大學(xué),2016.

[3] 周璐，郭超，鐘穎,等.小波網(wǎng)絡(luò)輔助卡爾曼濾波的捷聯(lián)慣導(dǎo)傳遞對(duì)準(zhǔn)[J].計(jì)算機(jī)測(cè)量與控制,2015，24(32)：2518-2520.

[4] 張宇陽(yáng)，蔡伯根.基于北斗/SINS融合的列車(chē)定位方法研究[D].北京：北京交通大學(xué)，2016.

[5] 羅建軍.組合導(dǎo)航原理及應(yīng)用[M].西安：西北工業(yè)大學(xué)出版社，2012.

[6] 李笠,王萬(wàn)良,徐新黎,等.基于網(wǎng)格排序的多目標(biāo)粒子群優(yōu)化算法[J],計(jì)算機(jī)研究與發(fā)展,2017，18(5)：1012-1023.

[7] 孟軍,史貫麗.融合粒子群優(yōu)化和遺傳算法的基因調(diào)控網(wǎng)絡(luò)構(gòu)建[J].計(jì)算機(jī)應(yīng)用,2016，811(11):2969-2973.

[8] 徐澤遠(yuǎn),伊國(guó)興,魏振楠.MEMS - SINS / GPS 組合導(dǎo)航系統(tǒng)設(shè)計(jì)[J].自動(dòng)化儀表,2016,37(12):67-70.

[9] 范廣坡,余學(xué)飛,盧廣文，等.改進(jìn)PSO-BP算法的壓力導(dǎo)絲溫度及非線性補(bǔ)償研究[J].自動(dòng)化儀表,2016，37(10)：12-16.

[10]周敏,劉健,楊林，等.一種自適應(yīng)卡爾曼濾波在動(dòng)態(tài)導(dǎo)航定位中的應(yīng)用[J].自動(dòng)化儀表,2016，37(8):1-4.

[11]黃杰，李軍，郭翔.遞推SOM神經(jīng)網(wǎng)絡(luò)在短時(shí)交通流預(yù)測(cè)中的應(yīng)用[J].自動(dòng)化儀表,2015，36(4):1-5.

[12]HU Z T,YUAN G Y,HU Y M.Training method of neural network based on cubature kalman filter[J].Control and Decision,2016，31(2):355-360.

[13]SISWANTORO J,PRABUWONO A S,ABDULLAH A.Linear model based on kalman filter for improving neural network classification performance[J].Expert System with Application,2016,49(15):112-122.

[14]李衛(wèi)東,賈洪明，馮祥雨.基于擴(kuò)展卡爾曼濾波的列車(chē)定位[J].大連交通大學(xué)學(xué)報(bào),2015(6):102-103.

[15]GU R H,BANERJEE A.Self-sensing SMA actuator using extended kalman filter and artificial neural network[J].Procedia,Engineering,2016，144(20):629-634.