本發(fā)明涉及一種檢測無線傳感網(wǎng)絡(luò)中高維信號重建質(zhì)量的方法，特別涉及一種在應(yīng)用迫零估計(ZF)與線性最小均方誤差估計(LMMSE)方法后檢測被重建信號的均方誤差指標的計算方法，屬于無線傳感網(wǎng)絡(luò)中信號重建的技術(shù)領(lǐng)域。

背景技術(shù)：

在無線傳感網(wǎng)絡(luò)中的信號重建過程中，隨機范德蒙矩陣扮演著重要的角色。無線傳感器所采集到的信號可以用傅里葉級數(shù)展開的形式表示。考慮到每個傳感器所處的地理位置不同，采集到的信號可以表示為一個范德蒙矩陣與一個表示該信號頻譜的向量之間的乘積?？紤]到采集信號中的噪聲影響，需要利用估計器對含噪信號進行重建。

經(jīng)典的算法有迫零估計與最小均方誤差估計，在利用這兩種估計器重建原始信號后，需要對重建的信號計算其均方誤差，從而確定信號重建的準確度。均方誤差的計算需要求解之前提到的范德蒙矩陣的特征值的概率密度函數(shù)。

當采集信號的維數(shù)、傳感器的個數(shù)、以及采集到的信號的帶寬較低時，范德蒙矩陣的特征值的概率密度函數(shù)可以通過蒙特卡洛數(shù)值計算的方法予以估計。然而，當上述參數(shù)值較大時，蒙特卡洛數(shù)值計算所需的計算量呈指數(shù)增長，因此需要一種能夠簡便又準確地獲得此矩陣特征值的概率密度函數(shù)的方法。

技術(shù)實現(xiàn)要素：

針對上述問題，本發(fā)明所要解決的技術(shù)問題是提供一種檢測無線傳感網(wǎng)絡(luò)中高維信號重建質(zhì)量的方法，本發(fā)明提供一種能夠簡便又準確地獲得范德蒙矩陣特征值的概率密度函數(shù)的方法，從而準確便捷地可以計算出系統(tǒng)的檢測性能指標。

本發(fā)明為解決上述技術(shù)問題采用以下技術(shù)方案：

本發(fā)明提供一種檢測無線傳感網(wǎng)絡(luò)中高維信號重建質(zhì)量的方法，所述無線傳感網(wǎng)絡(luò)中由m個無線傳感器構(gòu)成，每個無線傳感器分別采集一個d維的信號，該信號的帶寬表示為n，被采集的信號表示為y＝V^Ha+n_w，其中，y表示信號的測量值，a表示需要被重建的信號，n_w表示白噪聲信號，V表示一個維度為n^d×m的范德蒙矩陣。

當采用迫零估計或最小均方誤差估計方法重建信號時，重建結(jié)果的均方誤差分別為和其中，為白噪聲信號的方差，β＝n^d/m，E_λ{·}表示求解關(guān)于λ的期望值，λ為矩陣βVV^H的特征值。

通過以下方法求解λ的概率密度函數(shù)f_λ(d,β,λ)，具體為：

1)構(gòu)造優(yōu)化方程：且滿足條件方程：其中，μ^p為λ的p階矩；

2)通過高斯求積規(guī)則把優(yōu)化方程和條件方程中的積分運算用固定點位置f_λ(d,β,λ)取值的加權(quán)和表示出來，即：

其中，λ_j為高斯求積規(guī)則中所定義的第j個采樣點的坐標，N為采樣點個數(shù)，w_j為對應(yīng)λ_j的權(quán)值；

3)求解優(yōu)化方程和條件方程，得到對應(yīng)于λ_j的f_λ(d,β,λ_j)最優(yōu)值，從而擬合出f_λ(d,β,λ)。

根據(jù)得到的f_λ(d,β,λ)即能夠計算重建結(jié)果的均方誤差，從而依據(jù)均方誤差來衡量采用迫零估計或最小均方誤差估計方法進行信號重建的質(zhì)量。

作為本發(fā)明的進一步技術(shù)方案，采用迫零估計方法重建的信號為：

作為本發(fā)明的進一步技術(shù)方案，采用最小均方誤差估計方法重建的信號為：其中，I表示單位矩陣。

作為本發(fā)明的進一步技術(shù)方案，范德蒙矩陣V的元素為：其中，V_s,t表示V的第s行第t列元素，x_t表示第t個無線傳感器的地理位置。

作為本發(fā)明的進一步技術(shù)方案，a是一個n^d維向量。

作為本發(fā)明的進一步技術(shù)方案，y是一個m維向量。

本發(fā)明采用以上技術(shù)方案與現(xiàn)有技術(shù)相比，具有以下技術(shù)效果：當采集信號的維數(shù)、傳感器的個數(shù)以及采集到的信號的帶寬值較大時，采用本發(fā)明的技術(shù)方案能夠簡便又準確地獲得范德蒙矩陣特征值的概率密度函數(shù)。

附圖說明

圖1是本發(fā)明具體實施例中計算得出的系統(tǒng)均方誤差的數(shù)值圖。

具體實施方式

下面結(jié)合附圖以及具體實施例對本發(fā)明的技術(shù)方案做進一步的詳細說明：

在一個由m個無線傳感器構(gòu)成的無線傳感網(wǎng)絡(luò)中，每個傳感器負責采集一個d維的信號，該信號的帶寬表示為n。被采集的信號可以用需要被重建的信號表示為：y＝V^Ha+n_w，其中，y表示信號的測量值，是一個m維向量；a表示需要被重建的信號，是一個n^d維向量；n_w表示白噪聲信號；V表示一個維度為n^d×m的范德蒙矩陣，此矩陣的第s行第t列元素可以表示為：其中，x_t表示第t個無線傳感器的地理位置。

當應(yīng)用迫零估計(ZF)時，重建的信號表示為：該重建結(jié)果的均方誤差因此可以寫為：

當應(yīng)用最小均方誤差估計(LMMSE)時，重建的信號表示為：該重建結(jié)果的均方誤差因此可以寫為：其中，I表示單位矩陣。

當m、n、d的值較大時，假設(shè)n^d與m的比值表示為β＝n^d/m，則應(yīng)用迫零估計或最小均方誤差估計重建信號的均方誤差可以進一步被準確估計為：

其中，λ為矩陣βVV^H的特征值。因此，λ的概率密度函數(shù)(PDF)f_λ(d,β,λ)需要一種方法予以給出，其中，d與β為其兩個參數(shù)。

在計算λ的概率密度函數(shù)(PDF)之前，λ的p階矩已知，記為μ^p。λ的p階矩在A.Nordio,C.-F.Chiasserini,E.Viterbo,“Performance of linear field reconstruction techniques with noise and uncertain sensor locations,”IEEE Trans.on Signal Processing,Vol.56,No.8,pp.3535–3547,Aug.2008.中有詳細介紹。

構(gòu)造優(yōu)化方程：且滿足條件方程：

通過高斯求積規(guī)則把優(yōu)化方程和條件方程中的積分運算用固定點位置f_λ(d,β,λ)取值的加權(quán)和表示出來，即：其中，λ_j為高斯求積規(guī)則中所定義的第j個采樣點的坐標，N為采樣點個數(shù)，w_j為對應(yīng)λ_j的權(quán)值，也通過高斯求積規(guī)則獲得。因此，通過求解優(yōu)化方程和條件方程，對應(yīng)于每個λ_j采樣點的f_λ(d,β,λ_j)最優(yōu)值可以被求出，從而擬合得到f_λ(d,β,λ)。

在仿真測試環(huán)境中，我們測試當d＝1或4、且m＝n^d/β的參數(shù)設(shè)置，對應(yīng)于特征值λ的前12階矩均已知。假設(shè)系統(tǒng)分別采用迫零估計與最小均方誤差估計方法重建信號，我們分別比較由蒙特卡洛方法得出的系統(tǒng)均方誤差與由本發(fā)明所提出的方法所計算出的系統(tǒng)均方誤差。

如圖1所示的系統(tǒng)均方誤差的數(shù)值圖，通過觀察圖1可以看出，最小均方誤差估計方法的表現(xiàn)強于迫零估計方法，因為最小均方誤差的目標即是最小化均方誤差。更重要的是，隨著β值得增大，由本發(fā)明計算出的均方誤差與蒙特卡洛方法得出的均方誤差結(jié)果非常相近。但是本發(fā)明使得均方誤差的計算量大大縮小，圖1中，“□”和“×”的點表示蒙特卡洛方法得出的仿真結(jié)果，虛線和實線的線條表示的是本發(fā)明方法的仿真結(jié)果。

以上所述，僅為本發(fā)明中的具體實施方式，但本發(fā)明的保護范圍并不局限于此，任何熟悉該技術(shù)的人在本發(fā)明所揭露的技術(shù)范圍內(nèi)，可理解想到的變換或替換，都應(yīng)涵蓋在本發(fā)明的包含范圍之內(nèi)，因此，本發(fā)明的保護范圍應(yīng)該以權(quán)利要求書的保護范圍為準。