通過(guò)少量數(shù)據(jù)包估計(jì)無(wú)線傳感網(wǎng)絡(luò)鏈路質(zhì)量的方法
【專利摘要】本發(fā)明公開了通過(guò)少量數(shù)據(jù)包估計(jì)無(wú)線傳感網(wǎng)絡(luò)鏈路質(zhì)量的方法。方法包括如下幾步：接收數(shù)據(jù)包的信噪比、鏈路質(zhì)量指數(shù)和數(shù)據(jù)包接收率的統(tǒng)計(jì)、回歸參數(shù)計(jì)算和模型選擇、鏈路質(zhì)量估計(jì)后修改估計(jì)鏈路質(zhì)量表達(dá)式幾步，并給出具體的實(shí)施方式，該方法可以在鏈路質(zhì)量動(dòng)態(tài)變化的環(huán)境中根據(jù)少量接收數(shù)據(jù)包的參數(shù)估計(jì)將來(lái)一段時(shí)間內(nèi)的數(shù)據(jù)包接收率，并且具有可區(qū)分PRR不同的過(guò)渡型鏈路、能更準(zhǔn)確度量鏈路的SNR和LQI、具備自我修復(fù)能力以及易于實(shí)施的優(yōu)點(diǎn)。
【專利說(shuō)明】通過(guò)少量數(shù)據(jù)包估計(jì)無(wú)線傳感網(wǎng)絡(luò)鏈路質(zhì)量的方法

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及一種根據(jù)少量數(shù)據(jù)包的信噪比、鏈路質(zhì)量指數(shù)和包接收率估計(jì)無(wú)線傳 感網(wǎng)絡(luò)鏈路質(zhì)量的方法，屬于無(wú)線傳感網(wǎng)絡(luò)【技術(shù)領(lǐng)域】。

【背景技術(shù)】
[0002] 真實(shí)環(huán)境下，受自由衰落、多徑衰落和陰影衰落等因素的影響，無(wú)線傳感網(wǎng)絡(luò)的數(shù) 據(jù)包接收率會(huì)有差異，一般根據(jù)鏈路質(zhì)量分為穩(wěn)定連接型（數(shù)據(jù)包接收率> 90%)、過(guò)渡型 (10 % <數(shù)據(jù)包接收率<90 %)和低接收率型（接收率固定較低），大量的節(jié)點(diǎn)的鏈路質(zhì)量 屬于過(guò)渡型，其數(shù)據(jù)包接收率（以下簡(jiǎn)稱PRR)在10% - 90%之間動(dòng)態(tài)變化。
[0003] PRR能最直接反應(yīng)鏈路質(zhì)量，獲得PRR需要對(duì)大量數(shù)據(jù)包的收發(fā)情況進(jìn)行統(tǒng)計(jì)，但 是收發(fā)大量數(shù)據(jù)包會(huì)引發(fā)能量消耗和通信延遲等一系列的問(wèn)題。
[0004] 信噪比（以下簡(jiǎn)稱SNR)和鏈路質(zhì)量指數(shù)（以下簡(jiǎn)稱LQI)是衡量鏈路質(zhì)量的另外 兩個(gè)指標(biāo)，SNR可由無(wú)線射頻芯片中獲得的參數(shù)通過(guò)一定計(jì)算得到，LQI可由無(wú)線射頻芯片 直接獲得。以CC2420為例，該芯片有一個(gè)內(nèi)置的接收信號(hào)強(qiáng)度指示寄存器RSSI_VAL，讀取 該值再減去45dB即可得到接收信號(hào)功率，空閑（沒有接收數(shù)據(jù)包）時(shí)的接收信號(hào)功率就是 噪聲功率，正在接收數(shù)據(jù)包時(shí)的接收信號(hào)功率減去噪聲功率可得SNR。SNR和PRR之間有 單調(diào)增加的對(duì)應(yīng)關(guān)系，鏈路質(zhì)量指數(shù)（LQI)指示接收到的數(shù)據(jù)包的差錯(cuò)率，CC2420芯片提 供了一個(gè)平均相關(guān)值Corr，LQI可由（Corr-a) ?b計(jì)算得到（a和b是常數(shù)，在0 - 255之 間），通常，計(jì)算得到的LQI在50 - 110之間，LQI的均值和PRR之間的單調(diào)遞增關(guān)系較好。
[0005]目前，已有的采用SNR和LQI估計(jì)鏈路質(zhì)量有如下缺點(diǎn)：采用單一的SNR估計(jì)鏈路 質(zhì)量，3順在3.4-7.3(18之間變化，即可導(dǎo)致？1^在10%-90%之間劇烈變化，因此，僅采 用信噪比這一參數(shù)估計(jì)鏈路質(zhì)量存在鏈路質(zhì)量的"分辨率"過(guò)低（SNR很小的變化就會(huì)導(dǎo)致 PRR劇烈變化）、難以區(qū)分PRR不同的過(guò)渡型鏈路的問(wèn)題；由于CC2420只記錄接收到的數(shù)據(jù) 包的接收信號(hào)強(qiáng)度和LQI，因此，采用SNR、LQI或者兩者的均值估計(jì)鏈路質(zhì)量時(shí)，記錄的SNR 和LQI比實(shí)際鏈路偏高，從而導(dǎo)致SNR-PRR曲線和LQI-PRR曲線有偏差。


【發(fā)明內(nèi)容】

[0006] 針對(duì)現(xiàn)有的采用PRR、SNR或LQI估計(jì)鏈路質(zhì)量存在的不足，本發(fā)明提出一種將 PRR、SNR和PRR融合，僅通過(guò)少量數(shù)據(jù)包就可進(jìn)行鏈路質(zhì)量估計(jì)的方法。
[0007] -種通過(guò)少量數(shù)據(jù)包估計(jì)無(wú)線傳感網(wǎng)絡(luò)鏈路質(zhì)量的方法，步驟包括：接收數(shù)據(jù)包 SNR、LQI和PRR的統(tǒng)計(jì)、回歸參數(shù)計(jì)算和模型選擇、鏈路質(zhì)量估計(jì)和修改估計(jì)鏈路質(zhì)量表達(dá) 式；
[0008] 步驟1):接收數(shù)據(jù)包SNR、LQI和PRR的統(tǒng)計(jì)，
[0009]發(fā)射機(jī)和接收機(jī)的位置固定后，發(fā)射機(jī)發(fā)送StartNum個(gè)StartTransmitMsg消息， 消息的時(shí)間間隔為tpStartNum彡100,StartTransmitMsg消息中包含消息序列號(hào)；接著， 發(fā)射機(jī)無(wú)時(shí)間空隙地發(fā)送m*w個(gè)長(zhǎng)度為L(zhǎng)ength的數(shù)據(jù)包，m、w分別為10 - 20的自然數(shù)， Length為20 -30bytes,數(shù)據(jù)包的時(shí)間間隔為t2，每m個(gè)數(shù)據(jù)包記為一批，接收機(jī)接收到第一 個(gè)StartTransmitMsg消息后，對(duì)于第1批數(shù)據(jù)包，根據(jù)接收到的StartTransmitMsg消息的 時(shí)間和序列號(hào)分別計(jì)算發(fā)射機(jī)發(fā)送的第一批數(shù)據(jù)包中第1個(gè)、第m個(gè)數(shù)據(jù)包的時(shí)間點(diǎn)，在第 1個(gè)時(shí)間點(diǎn)之前將SNRt(rtal和LQIt(rtal置為0,在第1個(gè)和第2個(gè)時(shí)間點(diǎn)之間接收發(fā)射機(jī)發(fā)送 的數(shù)據(jù)包，計(jì)算每個(gè)包的信噪比的方法如下：接收機(jī)以60ms為周期定時(shí)采樣信道空閑時(shí)的 接收信號(hào)功率，作為一次噪聲基底的樣本值并更新噪聲基底表中的條目，每8次采樣后由8 個(gè)條目的平均值計(jì)算此時(shí)的噪聲基底，接收到的數(shù)據(jù)包的信噪比由該數(shù)據(jù)包的接收信號(hào)功 率減去噪聲基底得到；把接收到的數(shù)據(jù)包的SNR和LQI的值累加，統(tǒng)計(jì)接收到數(shù)據(jù)包的數(shù)目 m'，在第2個(gè)時(shí)間點(diǎn)將SNR和LQI的累加值分別記為SNRt(rtal和LQIt(rtal，按下述公式計(jì)算第 一批數(shù)據(jù)包的平均信噪比SNRAvgl、第一批數(shù)據(jù)包的平均鏈路質(zhì)量指數(shù)LQIAvgl和第一批數(shù)據(jù) 包的數(shù)據(jù)包接收率PR&:

【權(quán)利要求】
1. 一種通過(guò)少量數(shù)據(jù)包估計(jì)無(wú)線傳感網(wǎng)絡(luò)鏈路質(zhì)量的方法，其特征在于，步驟包括： 接收數(shù)據(jù)包SNR、LQI和PRR的統(tǒng)計(jì)、回歸參數(shù)計(jì)算和模型選擇、鏈路質(zhì)量估計(jì)和修改估計(jì)鏈 路質(zhì)量表達(dá)式； 步驟1):接收數(shù)據(jù)包SNR、LQI和PRR的統(tǒng)計(jì)， 發(fā)射機(jī)和接收機(jī)的位置固定后，發(fā)射機(jī)發(fā)送StartNum個(gè)StartTransmitMsg消息，消息 的時(shí)間間隔為tpStartNum彡100,StartTransmitMsg消息中包含消息序列號(hào)；接著，發(fā)射 機(jī)無(wú)時(shí)間空隙地發(fā)送m*w個(gè)長(zhǎng)度為L(zhǎng)ength的數(shù)據(jù)包，m、w分別為10 - 20的自然數(shù)，Length 為20 - 30bytes,數(shù)據(jù)包的時(shí)間間隔為t2,每m個(gè)數(shù)據(jù)包記為一批,接收機(jī)接收到第一個(gè) StartTransmitMsg消息后，對(duì)于第1批數(shù)據(jù)包，根據(jù)接收到的StartTransmitMsg消息的時(shí) 間和序列號(hào)分別計(jì)算發(fā)射機(jī)發(fā)送的第一批數(shù)據(jù)包中第1個(gè)、第m個(gè)數(shù)據(jù)包的時(shí)間點(diǎn)，在第1 個(gè)時(shí)間點(diǎn)之前將SNRtrtal和LQIttrtal置為0,在第1個(gè)和第2個(gè)時(shí)間點(diǎn)之間接收發(fā)射機(jī)發(fā)送 的數(shù)據(jù)包，計(jì)算每個(gè)包的信噪比的方法如下：接收機(jī)以60ms為周期定時(shí)采樣信道空閑時(shí)的 接收信號(hào)功率，作為一次噪聲基底的樣本值并更新噪聲基底表中的條目，每8次采樣后由8 個(gè)條目的平均值計(jì)算此時(shí)的噪聲基底，接收到的數(shù)據(jù)包的信噪比由該數(shù)據(jù)包的接收信號(hào)功 率減去噪聲基底得到；把接收到的數(shù)據(jù)包的SNR和LQI的值累加，統(tǒng)計(jì)接收到數(shù)據(jù)包的數(shù)目 m'，在第2個(gè)時(shí)間點(diǎn)將SNR和LQI的累加值分別記為SNRttrtal和LQIttrtal，按下述公式計(jì)算第 一批數(shù)據(jù)包的平均信噪比SNRavs1、第一批數(shù)據(jù)包的平均鏈路質(zhì)量指數(shù)LQIavs1和第一批數(shù)據(jù) 包的數(shù)據(jù)包接收率PRR1:
對(duì)于發(fā)射機(jī)發(fā)送的第2, 3……w批數(shù)據(jù)包，接收機(jī)也進(jìn)行類似的統(tǒng)計(jì)和運(yùn)算，分別得到SNRavs2,LQIavs2,PRR2,……，SNRAvgw，LQIAvgw，PRRw，去除PRR= 1 的數(shù)據(jù)點(diǎn)，假設(shè)去除后只有p組 數(shù)據(jù)點(diǎn)SNRAvgl，LQIAvgl，PRR1,......，SNRAvgp，LQIAvgp，PRRp ; 步驟2):模型參數(shù)計(jì)算和模型選擇， 2. 1)使用的模型 鏈路質(zhì)量估計(jì)中，步驟1)中每一批PRR〈1的m個(gè)數(shù)據(jù)包的接收率和數(shù)據(jù)包的SNRAvgjP LQIAvgi符合下述單調(diào)遞增的模型表達(dá)式之一： 模型一， Pimi^ Ai-SNR1hvj+ Br LQI1lnv⑷， 模型二， PRRi =A2 ·SNRAvgi+B2 ·LQIAvgi+C2 (5)， 其中，i= 1，2,……參數(shù)； 2. 2)模型參數(shù)計(jì)算 分別按上述2. 1)中兩種模型的表達(dá)式對(duì)p個(gè)數(shù)據(jù)點(diǎn)（SNRAvgl，LQIAvgl，PRR1)，…… ，（SNRavsp,LQIAvgp，PRRp)進(jìn)行線性回歸，具體計(jì)算方法如下： 第一種模型， 對(duì)上述P個(gè)數(shù)據(jù)點(diǎn)采用二元線性回歸，得參數(shù)A1和B1的值計(jì)算公式如下：
由此得到公式(5)中的參數(shù)ApB1，然后確定模型一的表達(dá)式；再將p個(gè)數(shù)據(jù)點(diǎn)對(duì) (SNRAvgl,LQIavs1),……，（SNRavsp,LQIavsp)帶入模型一的表達(dá)式，可得p個(gè)計(jì)算得到的數(shù)據(jù)包 接收率PRReall,PRReal2,……,PRRealp，記^>士|>撤,氣麗=士！洲氣則皮爾遜相 Pi~! Pi~l 關(guān)系數(shù)R1由以下公式求出
第二種模型， 當(dāng)P個(gè)數(shù)據(jù)點(diǎn)符合公式(5)的模型時(shí)，采用二元線性回歸的方法，得參數(shù)A2、B2和C2的值 計(jì)算公式如下： ^Dl 2=￥ (9)， 2=~D (10)， 1=? 紙 其中D、DpD2、D3的計(jì)算方法如下：

由此得到公式(5)中的參數(shù)A2、B2、C2，進(jìn)而確定模型二的表達(dá)式，再將p個(gè)數(shù)據(jù)點(diǎn)對(duì) (SNRAvgl,LQIavs1),……，（SNRavsp,LQIavsp)帶入模型二的表達(dá)式，得到p個(gè)計(jì)算得到的數(shù)據(jù)包 接收率PRReall,PRReal2,……，PRRm1p，記= 士輯氣Z5撤=士玄漬,.，則皮爾遜 相關(guān)系數(shù)R2可由以下公式求出
步驟2.3):模型選擇， 當(dāng)Rl和R2中只有一個(gè)大于0. 7時(shí)，皮爾遜相關(guān)系數(shù)大于0. 7的模型計(jì)算值能準(zhǔn)確反 應(yīng)實(shí)際的PRR變化，用該模型表達(dá)式對(duì)鏈路質(zhì)量進(jìn)行估計(jì)； 當(dāng)Rl和R2都大于0. 7時(shí)，這兩個(gè)模型計(jì)算值都能準(zhǔn)確反應(yīng)實(shí)際的PRR變化，因此這兩 個(gè)模型都可對(duì)鏈路質(zhì)量進(jìn)行估計(jì)； 當(dāng)Rl和R2都小于0. 7時(shí)，這兩個(gè)模型計(jì)算值都能不能準(zhǔn)確反應(yīng)實(shí)際的PRR變化，因此 重復(fù)第一步和第二步，直到選出皮爾遜相關(guān)系數(shù)大于0.7的模型，用該模型表達(dá)式對(duì)鏈路 質(zhì)量進(jìn)行估計(jì)； 步驟3):鏈路質(zhì)量估計(jì)， 發(fā)射機(jī)發(fā)送q個(gè)長(zhǎng)度為L(zhǎng)ength的數(shù)據(jù)包，q為10 - 20,接收機(jī)將收到的數(shù)據(jù)包的SNR和LQI分別累加，得到SNRsum和LQIsum，通過(guò)下述公式計(jì)算出SNR'Avg和LQI'Avg :
將計(jì)算得到的SNR'Avg和LQI'Avg代入步驟2. 3)選擇的模型的表達(dá)式即可估計(jì)出這段時(shí) 間內(nèi)本條鏈路的包接收率PRRimsmt ; 發(fā)射機(jī)在緊接著5秒內(nèi)發(fā)送長(zhǎng)度為L(zhǎng)ength數(shù)據(jù)包時(shí)， DDD --pnn ΓIVIVprediction Γ J\J\presenl： j 步驟4):鏈路質(zhì)量估計(jì)后修改估計(jì)鏈路質(zhì)量表達(dá)式， 在步驟3)中，如果發(fā)射機(jī)第1次發(fā)送了q個(gè)長(zhǎng)度為L(zhǎng)ength的數(shù)據(jù)包，按表達(dá)式?)計(jì)算 出的PRRpralietitm記為PRRpralietitml，緊接著5秒內(nèi)發(fā)送了多個(gè)長(zhǎng)度為L(zhǎng)ength的數(shù)據(jù)包，統(tǒng)計(jì) 得到的數(shù)據(jù)包接收率為PRR1，鏈路質(zhì)量估計(jì)的相對(duì)誤差為
同理可得第2, 3, 4,......，r次發(fā)送了q個(gè)長(zhǎng)度為L(zhǎng)ength的數(shù)據(jù)包并實(shí)測(cè)數(shù)據(jù)包接收率 后，鏈路質(zhì)量估計(jì)的相對(duì)誤差為S2,δ3,δ4,......，δρ當(dāng)相對(duì)誤差的值大于10%的次數(shù) 有3次時(shí)，就開始修改估計(jì)鏈路質(zhì)量表達(dá)式的過(guò)程，即重復(fù)步驟1)、步驟2)、步驟3)。
2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述的步驟2. 3):當(dāng)Rl和R2都大于0. 7 時(shí)，優(yōu)選的，以皮爾遜相關(guān)系數(shù)更大的模型表達(dá)式對(duì)鏈路質(zhì)量進(jìn)行估計(jì)。
【文檔編號(hào)】H04W84/18GK104320805SQ201410584342
【公開日】2015年1月28日 申請(qǐng)日期:2014年10月27日 優(yōu)先權(quán)日:2014年10月27日
【發(fā)明者】魯琛 申請(qǐng)人:浙江科技學(xué)院