武加文，李光輝,3

（1. 江南大學 物聯(lián)網(wǎng)工程學院，江蘇 無錫 214122; 2. 物聯(lián)網(wǎng)技術(shù)應(yīng)用教育部工程技術(shù)研究中心，江蘇 無錫214122; 3. 江蘇省模式識別與計算智能工程實驗室，江蘇 無錫 214122）

近年來，隨著無線通信，計算機網(wǎng)絡(luò)和傳感器技術(shù)的巨大進步，無線傳感器網(wǎng)絡(luò)(wireless sensor networks，WSN)已成為重要的研究對象并被廣泛應(yīng)用于醫(yī)療、軍事、環(huán)境監(jiān)測等眾多領(lǐng)域[1]。在環(huán)境監(jiān)測領(lǐng)域，大規(guī)模的無線傳感器網(wǎng)絡(luò)通常被部署于無人監(jiān)管、氣候復(fù)雜的野外環(huán)境中。若長期使用，節(jié)點讀數(shù)容易產(chǎn)生數(shù)據(jù)漂移。這對于終端用戶來說是一個嚴重的問題，因為在無線傳感器網(wǎng)絡(luò)的大規(guī)模應(yīng)用場景中必須使用精確數(shù)據(jù)，才能制定科學決策。因此，對節(jié)點數(shù)據(jù)進行漂移校準至關(guān)重要。然而，在許多大規(guī)模無線傳感器網(wǎng)絡(luò)的應(yīng)用中，數(shù)百到數(shù)千個節(jié)點被分散部署到氣候復(fù)雜的野外環(huán)境中，手動校準這些節(jié)點是非常困難的[2]。因此，需要自動跟蹤和校準節(jié)點數(shù)據(jù)漂移的方法，現(xiàn)有的方法可以分為兩類：盲校準和非盲校準。

非盲校準方法是基于傳感器節(jié)點讀數(shù)和已知的參考信息作為輸入來調(diào)整參數(shù)[3]。一種簡單的非盲校準方法是將已知的激勵作用于傳感器網(wǎng)絡(luò)中并測量其響應(yīng)，將網(wǎng)絡(luò)響應(yīng)與參考信息(網(wǎng)絡(luò)預(yù)期值)進行比較，然后相應(yīng)地調(diào)整參數(shù)。參考信息一般為出廠標準值或用戶手動校準的結(jié)果。例如，Ramanathan等[4]采用額外的高精度節(jié)點讀數(shù)作為參考信息。另一種類型的方法是手動校準一組節(jié)點，然后未校準節(jié)點組基于該組來調(diào)整參數(shù)[5]。非盲校準方法適用于特定空間(例如室內(nèi))小規(guī)模網(wǎng)絡(luò)的場景下[6]，然而在大規(guī)模無線傳感器網(wǎng)絡(luò)中，使用上述方法需要在每個階段進行大量的工作，這是不切實際的且成本高昂。

鑒于非盲校準方法存在的局限性，盲校準方法得到了廣泛的關(guān)注。許多盲校準方法是基于節(jié)點被密集部署的假設(shè)，因此當網(wǎng)絡(luò)中未產(chǎn)生數(shù)據(jù)漂移時，節(jié)點間應(yīng)具有相關(guān)的讀數(shù)[7]。然而，在許多現(xiàn)有的大型網(wǎng)絡(luò)中，由于環(huán)境因素和成本因素，可能無法實現(xiàn)節(jié)點的密集部署。Balzano等[8]首先提出了一種無需密集部署的校準方法，假定節(jié)點的測量位于測量空間的較低維信號子空間中，然后通過信號子空間投影將校準過程轉(zhuǎn)換為求解線性方程。來自清華大學的楊華中等[9,10]對這種方法進行了利用和擴充，取得了非常良好的實用價值。

基于相鄰節(jié)點間具有相關(guān)測量值的假設(shè)，Takruri等[11]提出了一種分布式遞歸貝葉斯(distributed recursive Bayesian，DRB)盲校準算法，其使用鄰居節(jié)點讀數(shù)的平均值作為漂移節(jié)點的預(yù)測值。為了更好地擬合傳感器節(jié)點間的時空相關(guān)性，Takruri等[12]又提出了一種基于支持向量回歸(support vector regression，SVR)來預(yù)測真實值的方法。針對SVR可能存在的訓練精度較低的情況，Kumar等[13]采用克里金插值作為預(yù)測方法，以避免SVR的訓練階段。這些方法基于節(jié)點間的時空相關(guān)性以預(yù)測節(jié)點數(shù)據(jù)的漂移，并都使用了卡爾曼濾波器(Kalman filter，KF)以跟蹤漂移。然而，一旦預(yù)測值不精確，偏差較大的預(yù)測值也會被用來進一步校準漂移值，從而導(dǎo)致最終校準值的精度受累積誤差的限制。

針對上述問題，本文提出了一種基于遺傳算法(genetic algorithm，GA)優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)(BP neural network optimized by GA，GABP)與卡爾曼濾波器相結(jié)合的無線傳感器網(wǎng)絡(luò)數(shù)據(jù)漂移盲校準算法(GABP-KF)。該算法首先對節(jié)點讀數(shù)進行二次采集和降噪處理。然后使用基于遺傳算法優(yōu)化的神經(jīng)網(wǎng)絡(luò)對目標節(jié)點和其鄰居節(jié)點間的時空相關(guān)性進行數(shù)學建模，從而得到目標節(jié)點的預(yù)測值。最后將目標節(jié)點的預(yù)測值和測量值反饋入卡爾曼濾波器中以跟蹤并校準其數(shù)據(jù)漂移。仿真結(jié)果表明，GABP-KF具有良好的漂移校準性能。

1 預(yù)備知識

1.1 遺傳算法

遺傳算法是模擬自然界遺傳機制的生物進化模型，是一種啟發(fā)式隨機搜索最優(yōu)化方法。它將“物競天擇，適者生存”的進化原理和遺傳變異理論引入編碼串聯(lián)群體中，即在每一次的進化過程中，按照所構(gòu)造的適應(yīng)度函數(shù)，通過選擇操作、交叉操作和變異操作來對個體進行篩選，把適應(yīng)度值好的個體傳給下一代，反復(fù)循環(huán)，直到最終個體滿足特定的條件表示進化結(jié)束[14]。遺傳算法具備良好的全局搜索能力，不受求導(dǎo)和函數(shù)連續(xù)性的限定，能自動獲取待優(yōu)化的搜索空間，從而自適應(yīng)地控制搜索過程[15]。

1.2 卡爾曼濾波器

卡爾曼濾波器是由R.E.Kalman提出的線性最小方差遞推估算的方法，是一種離散數(shù)據(jù)線性濾波問題的遞歸解決方案[16]。它是利用當前時刻的觀測值和前一時刻的估計值來更新當前時刻的最佳估計，其所需數(shù)據(jù)存儲量較小，便于實時處理[17]。

2 基于GABP-KF的WSN數(shù)據(jù)漂移校準方法

2.1 定義測量環(huán)境

假設(shè)在某一區(qū)域內(nèi)的傳感器網(wǎng)絡(luò)中，隨機分布多個節(jié)點。其中每個節(jié)點采集該區(qū)域內(nèi)的環(huán)境參數(shù)(溫度、濕度或其他類型參數(shù))。隨著時間的推移，一些節(jié)點產(chǎn)生數(shù)據(jù)漂移，降低了數(shù)據(jù)質(zhì)量，從而影響了傳感器網(wǎng)絡(luò)的可靠性。為了更好地模擬實際場景，本文對傳感器網(wǎng)絡(luò)的測量環(huán)境做了如下定義。

定義1傳感器節(jié)點的數(shù)據(jù)漂移指的是節(jié)點測量值受其內(nèi)部固有偏差或外界環(huán)境影響所產(chǎn)生的緩慢、單向、長期的變化。在本文中，節(jié)點的數(shù)據(jù)漂移表示如下：

式中：X表示節(jié)點的測量值；T表示所測環(huán)境的真實值；d表示漂移值；W表示測量高斯白噪聲。其中，漂移d是平滑的，因為漂移通常是一個緩慢的線性或指數(shù)的變化過程，沒有突變，激增或尖峰現(xiàn)象[18]。同時，漂移的產(chǎn)生是隨機的，與其內(nèi)部構(gòu)造以及環(huán)境因素密切相關(guān)。因此，不同節(jié)點是否產(chǎn)生漂移以及漂移值的大小不具有相關(guān)性。由于節(jié)點在部署前都經(jīng)過預(yù)先校準，以確保它們處于工作狀態(tài)。因此，節(jié)點部署后的短時間內(nèi)，漂移值應(yīng)為零。在本文中，設(shè)定節(jié)點不會同時產(chǎn)生數(shù)據(jù)漂移，且產(chǎn)生的漂移值不同。

定義2在以目標節(jié)點為圓心的測量區(qū)域內(nèi)，該節(jié)點向其周圍節(jié)點傳輸信號，該信號的強度會隨距離的增加而減弱，當信號強度小于閾值時，接收節(jié)點無法接收到數(shù)據(jù)。因此，通信半徑指的是信號衰減至閾值的距離，而鄰居節(jié)點指目標節(jié)點通信半徑內(nèi)的其他節(jié)點。在本文中，假設(shè)目標節(jié)點i共有m-1個鄰居節(jié)點，分別標號為1,2,···,m,i-1, i+1,···,m。

定義3目標節(jié)點i與其鄰居節(jié)點的時空相關(guān)性指的是同一時段內(nèi)的節(jié)點測量數(shù)據(jù)的變化趨勢是一致的。在本文中，假設(shè)目標節(jié)點i與其鄰居節(jié)點間具有時空相關(guān)性，因此，可以用來預(yù)測目標節(jié)點的測量值，即：

2.2 算法框架

GABP-KF算法的設(shè)計目標是在隨機誤差(噪聲)和系統(tǒng)誤差(漂移)的干擾下精確地校準目標節(jié)點i的測量值 xi。該算法有兩個運行階段：訓練階段和測試階段。在訓練階段，首先對所有節(jié)點數(shù)據(jù)進行去噪處理。然后使用來自鄰居節(jié)點的去噪后測量值 { X }m-1作為GABP神經(jīng)網(wǎng)絡(luò)的輸入值，并輸出節(jié)點i的預(yù)測值。在測試階段，將預(yù)測值和節(jié)點i的去噪后測量值 y 輸入到卡爾曼i濾波器中以跟蹤其數(shù)據(jù)漂移 di。最后通過測量值yi減去估計漂移值 di以得到最終校準后的讀數(shù)。圖1描述了本文算法的具體框架。

圖1 GABP-KF算法的數(shù)據(jù)漂移校準框架Fig.1 Data drift calibration framework of GABP-KF algorithm

2.3 基于GABP神經(jīng)網(wǎng)絡(luò)的建模過程

2.3.1 數(shù)據(jù)預(yù)處理

在傳感器網(wǎng)絡(luò)的實際部署中，由于受到傳感器節(jié)點內(nèi)部固有偏差或其他環(huán)境因素，節(jié)點在測量過程中會受到噪聲的干擾。在節(jié)點校準之前，必須考慮消除噪聲的影響。本文采用小波去噪[19]方法來降低噪聲的影響，首先使用閾值去噪法對含噪數(shù)據(jù)進行小波分解，再對小波系數(shù)進行閾值處理，最后利用處理后結(jié)果重構(gòu)原信號。實驗證明該方法在信號保護和噪聲抑制之間達到了良好的平衡。

2.3.2 GABP神經(jīng)網(wǎng)絡(luò)的建模

傳統(tǒng)的BP神經(jīng)網(wǎng)絡(luò)具有良好的自學習、自適應(yīng)能力，但同時也存在容易陷入局部最優(yōu)解、收斂速度差等問題。為了克服這些缺點，本文使用遺傳算法來優(yōu)化BP神經(jīng)網(wǎng)絡(luò)，將遺傳算法的全局搜索能力與BP神經(jīng)網(wǎng)絡(luò)的局部搜索能力有效地結(jié)合起來，從而提高了網(wǎng)絡(luò)的預(yù)測精度和泛化能力。具體方法如下：首先，BP神經(jīng)網(wǎng)絡(luò)根據(jù)數(shù)據(jù)集的劃分情況確定網(wǎng)絡(luò)結(jié)構(gòu)，從而確定遺傳算法的個體長度。其次，算法通過選擇操作、交叉操作和變異操作來優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的參數(shù)。最后，BP神經(jīng)網(wǎng)絡(luò)使用最優(yōu)參數(shù)來預(yù)測目標節(jié)點的測量值。GABP神經(jīng)網(wǎng)絡(luò)的模型建立方法如圖2所示。

本文使用在傳感器網(wǎng)絡(luò)部署初期(一般為前五天)所采集到的節(jié)點數(shù)據(jù)作為訓練數(shù)據(jù)集，訓練數(shù)據(jù)集劃分為： T rain={TY;TX}。TY表示數(shù)據(jù)測試集，TX表示訓練集。 T Y={xi,t:t=1,2,···,n}，TX={xj,t:j=1,2,···,i-1,i+1,···,m,t=1,2,···,n}，其中，j表示節(jié)點序號，t表示時序。

圖2 GABP神經(jīng)網(wǎng)絡(luò)的建模算法Fig.2 Modeling algorithm of GABP neural network

2.4 基于卡爾曼濾波的漂移跟蹤算法

本文設(shè)定t時刻節(jié)點i的平滑漂移值為 di,t，t時刻節(jié)點i所測量的真實溫度值為 Ti,t，t時刻節(jié)點i的測量值為 ri,t。其中i表示節(jié)點序號，t表示時間序列。為了預(yù)測漂移值 di,t，首先需要對節(jié)點的數(shù)據(jù)漂移進行數(shù)學建模，見式(3)：

式中： wi,t表示高斯噪聲，而 Qi,t表示狀態(tài)噪聲協(xié)方差。

式(3)表示節(jié)點數(shù)據(jù)漂移的狀態(tài)跟蹤方程，用來跟蹤漂移值隨時間的變化情況。除了狀態(tài)方程，還需要建立一個觀測方程來模擬節(jié)點的觀測結(jié)果。假設(shè)現(xiàn)實中存在一種測量儀器，可以測量出節(jié)點的數(shù)據(jù)漂移值，同時一定存在與節(jié)點測量時相關(guān)聯(lián)的誤差，見式(4)：

式中： vi,t表示高斯噪聲，而 Ri,t表示觀測噪聲協(xié)方差。

實際上并不存在可以直接測量節(jié)點數(shù)據(jù)漂移的儀器。在這種情況下，理想的數(shù)據(jù)漂移值 di,t等于測量值減去實際值，見式(5)：

然而在實際場景中，實際值是無法得知的，只能使用其預(yù)測值 x 。因此：

綜上所述，式(3)和式(4)建立了一個傳感器節(jié)點的狀態(tài)-觀測模型，使用卡爾曼濾波器以分散迭代的方式進行跟蹤漂移，具體算法參見圖3。

圖3 針對目標節(jié)點的漂移跟蹤算法Fig.3 Drift tracking algorithm for target node

如圖3所示，其中 Zt、 Qt、 Rt分別表示t時刻目標節(jié)點的漂移測量值、過程噪聲和測量噪聲。由于目標節(jié)點i共測量了n個數(shù)據(jù)，因此算法需要迭代n次，即，算法結(jié)束條件為：t ＞ n。在進行第t次迭代時，首先基于t-1時刻的漂移估計值dt-1和均方誤差 pt-1預(yù)測t時刻的漂移估計值d_pre和均方誤差p_ pre，然后計算t時刻的卡爾曼濾波增益K，最后更新t時刻的漂移值 dt和均方誤差 Pt。算法迭代完成后，基于濾波后漂移值dt來校準節(jié)點數(shù)據(jù)，見式(7)：

3 實驗結(jié)果及分析

為了驗證GABP-KF的算法性能，本文在mat-lab 2013b 的實驗環(huán)境下針對國際通用的Intel Berkeley(LBRL)無線傳感器網(wǎng)絡(luò)數(shù)據(jù)集和LUCE無線傳感器網(wǎng)絡(luò)數(shù)據(jù)集分別實現(xiàn)了GABP-KF算法、未優(yōu)化的BP-KF算法、Takruri提出的SVRKF算法[12]以及DRB-KF算法[11]，并進行了實驗結(jié)果的比較。

3.1 數(shù)據(jù)集

3.1.1 IBRL數(shù)據(jù)集

IBRL數(shù)據(jù)集來自于部署在Intel Berkeley實驗室內(nèi)的無線傳感器網(wǎng)絡(luò)，包含54個節(jié)點，用于監(jiān)控實驗室環(huán)境(參見圖4)。

圖4 IBRL的節(jié)點分布Fig.4 Deployment scheme of IBRL sensor network

選取由IBRL部署中的兩組傳感器節(jié)點的數(shù)據(jù)子集進行實驗。第1組數(shù)據(jù)子集(IBRL_1)包含的節(jié)點 ID 分別是 1、2、3、4、6、7、10、33。第2組數(shù)據(jù)子集(IBRL_2)包含的節(jié)點ID分別是17、18、19、20、21、22、23。兩組數(shù)據(jù)集都對應(yīng)于2004-02-28～2004-03-07 9天內(nèi)所收集的數(shù)據(jù)。

3.1.2 LUCE數(shù)據(jù)集

LUCE數(shù)據(jù)集(洛桑城市冠層實驗)來自于2006年7月以來部署在洛桑聯(lián)邦理工學院內(nèi)的無線傳感器網(wǎng)絡(luò)。該網(wǎng)絡(luò)共包含97個節(jié)點，根據(jù)節(jié)點之間的時空相關(guān)性分為10組傳感器節(jié)點集。

選取LUCE數(shù)據(jù)集中的兩組傳感器節(jié)點的數(shù)據(jù)子集作為實驗對象，第1組數(shù)據(jù)子集(LUCE_1)包含的節(jié)點 ID 分別是 10、14、15、17、18、19。第2組數(shù)據(jù)子集(LUCE_2)包含的節(jié)點ID分別是 21、23、24、25、26、27、28。兩組數(shù)據(jù)子集都對應(yīng)于2016-10-10～2016-10-13 4天內(nèi)所收集的數(shù)據(jù)。

由于上述數(shù)據(jù)集中節(jié)點的測量數(shù)據(jù)量龐大，其中存在缺漏值，且采樣周期均為31 s。因此，對每個節(jié)點分別以5 min的時間間隔進行重新采樣，并均使用溫度作為評估數(shù)據(jù)。為了精確評估GABP-KF算法的性能，排除實驗數(shù)據(jù)由于缺漏現(xiàn)象造成的影響，人工建立了一組具有時空相關(guān)性的測試數(shù)據(jù)集(DR_d)。表1列出了所使用的所有數(shù)據(jù)集。

表1 實驗所用數(shù)據(jù)集Table1 Experimental data sets

3.2 去噪處理

針對上述數(shù)據(jù)集中的所有節(jié)點數(shù)據(jù)，采用硬閾值去噪的方法進行處理，以減少噪聲對GABPKF算法評估所造成的影響。由于本文所采用的數(shù)據(jù)集中噪聲干擾較小，為進行去噪算法的評估，故選擇在數(shù)據(jù)集Ⅲ中10號節(jié)點處對節(jié)點數(shù)據(jù)手動添加高斯白噪聲，去噪結(jié)果參見圖5。如圖5所示，原始數(shù)據(jù)與去噪后數(shù)據(jù)幾近重合，這表明本方法對噪聲有很好的抑制作用。

圖5 小波去噪實驗結(jié)果Fig.5 Experimental results of wavelet denoising

3.3 模型對比實驗

為比較GABP-KF算法與其他同類算法的模型擬合程度，基于表1中的5個數(shù)據(jù)集分別隨機選擇其中一個節(jié)點的數(shù)據(jù)針對4種算法進行了5次對比實驗，并選取平均值作為最終數(shù)據(jù)，表2給出了4種算法的模型擬合程度對比結(jié)果。本文選用均方誤差(式(8))與決定系數(shù)(式(9))作為判斷依據(jù)。

式中：n表示訓練樣本數(shù)目， xi與分別為第i個樣本的真實值與預(yù)測值。一般情況下，決定系數(shù)越高，均方誤差越小，表明模型的擬合程度越好。

表2 4種算法的模型擬合程度對比結(jié)果Table2 Comparison of the fitting degree among four algorithms

從表2可以看出，基于本文所采用的數(shù)據(jù)集，GABP算法的均方誤差比DRB算法、BP、SVR算法平均減少了 0.110 8、0.121 0、0.249 5；GABP 算法的決定系數(shù)比DRB、BP、SVR算法平均提高了1.74%、2.56%、4.23%。其中，4種算法在數(shù)據(jù)集V下的擬合程度均好于其余數(shù)據(jù)集，這是由于數(shù)據(jù)集V為人工建立，不受外界誤差影響，其余4個數(shù)據(jù)集選自真實數(shù)據(jù)集，均受到數(shù)據(jù)缺失，數(shù)據(jù)突變等誤差的影響。DRB、BP、SVR算法受這些誤差的影響較大，而GABP神經(jīng)網(wǎng)絡(luò)通過GA優(yōu)化有效地避免了局部誤差，因此，GABP神經(jīng)網(wǎng)絡(luò)的模型擬合程度要優(yōu)于其他算法。

圖6表示4種算法在數(shù)據(jù)集I下針對1號節(jié)點所建立的預(yù)測模型，其中real value描述的曲線表示真實值，其余4條曲線分別表示GABP、SVR、BP和DRB算法所獲得的預(yù)測結(jié)果。從圖6中可以直觀地看出，GABP神經(jīng)網(wǎng)絡(luò)的擬合程度更好，這與表2的結(jié)果是一致的。綜上所述，GABP算法具有更好的模型擬合程度，所輸出的預(yù)測值與真實值的誤差更小。

3.4 基于GABP-KF的漂移校準對比實驗

在漂移校準過程中，于不同節(jié)點測試數(shù)據(jù)的隨機位置處引入不同的指數(shù)漂移值與測量噪聲。其中，當觀測噪聲協(xié)方差R值或狀態(tài)噪聲協(xié)方差Q值設(shè)置過高時，預(yù)測值會產(chǎn)生較大誤差。因此，采用試湊法對參數(shù)值進行調(diào)整，在數(shù)據(jù)集I下，基于網(wǎng)格搜索的原理將卡爾曼濾波器的參數(shù)Q和R分別設(shè)置為0.1和0.01。

圖6 4種算法的模型對比結(jié)果Fig.6 Comparison among the training models of four algorithms

為比較GABP-KF算法與其他同類算法的漂移校準性能，基于表1中的5個數(shù)據(jù)集分別隨機選擇其中一個節(jié)點的數(shù)據(jù)針對4種算法進行了3次對比實驗，并選取平均值作為最終數(shù)據(jù)，表3給出了實驗結(jié)果。

表3 4種算法的漂移校準性能對比結(jié)果Table3 Comparison of the drift calibration performance among four algorithms

從表3中可以看出，基于本文所采用的數(shù)據(jù)集，GABP-KF算法的均方誤差比DRB-KF、BPKF、SVR-KF算法平均減少了0.295 2、0.202 0、0.370 2；GABP-KF算法的決定系數(shù)比DRB-KF、BP-KF、SVR-KF算法平均提高了1.43%、1.72%、2.40%。實驗表明，基于本文所采用的數(shù)據(jù)集，GABP-KF算法具有更好的漂移校準性能。

圖7表示在數(shù)據(jù)集Ⅲ下針對14號節(jié)點的數(shù)據(jù)使用GABP-KF算法進行漂移校準后的結(jié)果。從圖7中可以看出，節(jié)點真實值曲線與濾波后數(shù)值曲線基本保持一致。這表明本算法成功地消除了引入的漂移值和測量噪聲。

圖7 卡爾曼濾波后的校準示例Fig.7 Example of data drifts calibration after Kalman filter

3.5 整體評估

在進行整體評估時，考慮以下兩種情形。第1種情形為所有節(jié)點數(shù)據(jù)均未添加漂移值，稱為N-data，第2種情形為在指定的節(jié)點數(shù)據(jù)中加入漂移值，稱為D-data。兩組數(shù)據(jù)N-data和D-data在經(jīng)過GABP-KF算法的校準后分別稱為CN-data和CD-data。其中D-data中在不同節(jié)點測試數(shù)據(jù)的隨機位置處引入不同的指數(shù)漂移，即每個節(jié)點開始產(chǎn)生漂移的時刻不同，產(chǎn)生漂移的大小不同，是否產(chǎn)生漂移互不相關(guān)。

本文比較了5種情形下節(jié)點的平均絕對誤差MAE 值 (見式 (10))。如圖8所示，情形 1、2、3、4、5分別表示僅有一個節(jié)點漂移、兩個節(jié)點同時漂移、3個節(jié)點同時漂移、4個節(jié)點同時漂移、5個節(jié)點同時漂移等5種情形。具體計算方法如下：在每個時刻分別計算每個節(jié)點的測量數(shù)據(jù)N-data與D-data之間的MAE值和校準后數(shù)據(jù)CN-data和CD-data之間的MAE值，然后計算網(wǎng)絡(luò)中5種情形下的整體MAE值。本文選擇在數(shù)據(jù)集IV下21號節(jié)點的數(shù)據(jù)上進行上述對比實驗。

式中：n為訓練樣本數(shù)目，Xi與分別為第i個樣本的真實值和預(yù)測值。

圖8 5種情形下MAE值的變化情況Fig.8 Comparison of MAE among five cases

圖8 (a)表示5種情形下測量數(shù)據(jù)N-data和D-data之間的平均絕對誤差。圖8(b)表示5種情形下校準數(shù)據(jù)CN-data和CD-data之間的平均絕對誤差。比較圖8(a)與圖8(b)，顯然，應(yīng)用GABPKF算法導(dǎo)致所有情形下的MAE值較小。本文假設(shè)無線傳感器網(wǎng)絡(luò)中閾值MAE值為1.2，超過閾值則表示該網(wǎng)絡(luò)失效，閾值的選擇取決于實際網(wǎng)絡(luò)中允許的容錯誤差。在圖8(a)中，情景2～5的曲線在第一天時就跨越了閾值線。相比之下，圖8(b)中的所有情形的曲線在實驗的整個周期中均未跨越閾值線。這證明使用本算法可以允許多個節(jié)點同時產(chǎn)生數(shù)據(jù)漂移的情況，這更符合實際情況。

4 結(jié)束語

本文提出了一種使用基于遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)和卡爾曼濾波器相結(jié)合的無線傳感器網(wǎng)絡(luò)數(shù)據(jù)漂移盲校準算法。仿真實驗表明，與以往同類算法相比，GABP-KF算法對無線傳感器網(wǎng)絡(luò)節(jié)點數(shù)據(jù)流漂移有更好的跟蹤和校準性能，使用該算法有效地提高了傳感器網(wǎng)絡(luò)數(shù)據(jù)的可靠性。在大規(guī)模無線傳感器網(wǎng)絡(luò)的應(yīng)用中，節(jié)點可以同時測量多種屬性的數(shù)據(jù)，本算法只考慮了單一屬性(溫度)下的漂移校準情況。在今后的研究將考慮節(jié)點同時采集多種屬性數(shù)據(jù)時產(chǎn)生或未產(chǎn)生漂移的各種情況，以便更好地模擬實際情況。