本發(fā)明屬于無線射頻識(shí)別技術(shù)領(lǐng)域，具體涉及一種基于線性插值的射頻識(shí)別防碰撞標(biāo)簽數(shù)估計(jì)方法。

背景技術(shù)：

動(dòng)態(tài)幀時(shí)隙Aloha(DFSA,Dynamic Frame Slot Aloha)算法在RFID領(lǐng)域得到了廣泛的應(yīng)用,該算法根據(jù)剩余標(biāo)簽數(shù)動(dòng)態(tài)調(diào)整幀長(zhǎng)度，與固定幀時(shí)隙Aloha算法相比顯著提高了識(shí)別效率。理論研究和實(shí)驗(yàn)均證實(shí)，當(dāng)幀長(zhǎng)度與待識(shí)別標(biāo)簽數(shù)大致相等時(shí)，識(shí)別效率最高。為了將下一幀的幀長(zhǎng)度調(diào)整為與剩余標(biāo)簽數(shù)大致相等，以獲得最佳識(shí)別效率，需要對(duì)剩余標(biāo)簽進(jìn)行估計(jì)，其缺點(diǎn)在于，初始標(biāo)簽數(shù)估計(jì)會(huì)嚴(yán)重影響標(biāo)簽識(shí)別效率，如果估計(jì)標(biāo)簽數(shù)遠(yuǎn)大于實(shí)際標(biāo)簽數(shù)，則后續(xù)每幀調(diào)整的Q參數(shù)將大于實(shí)際待識(shí)別標(biāo)簽數(shù)，導(dǎo)致空閑時(shí)隙增加，浪費(fèi)時(shí)隙資源，進(jìn)而降低識(shí)別效率。反之，如果估計(jì)標(biāo)簽數(shù)遠(yuǎn)小于實(shí)際標(biāo)簽數(shù)，則后續(xù)幀長(zhǎng)調(diào)整時(shí)Q參數(shù)將會(huì)小于待識(shí)別標(biāo)簽數(shù)，導(dǎo)致碰撞概率增加，同樣會(huì)降低識(shí)別效率。

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

為了克服現(xiàn)有技術(shù)的不足,本發(fā)明提供一種基于線性插值的射頻識(shí)別防碰撞標(biāo)簽數(shù)估計(jì)方法，根據(jù)第一幀有效識(shí)別標(biāo)簽數(shù)，利用線性插值的方法估計(jì)初始標(biāo)簽數(shù)，進(jìn)而獲得剩余待識(shí)別的標(biāo)簽數(shù)，并按照幀長(zhǎng)度等于待識(shí)別標(biāo)簽數(shù)的最優(yōu)幀識(shí)別效率原則自適應(yīng)動(dòng)態(tài)調(diào)整后續(xù)識(shí)別幀的長(zhǎng)度，以達(dá)到最佳的識(shí)別效率，實(shí)現(xiàn)全局吞吐量最大化。

為了實(shí)現(xiàn)上述目的，本發(fā)明采用的技術(shù)方案是：

一種基于線性插值的射頻識(shí)別防碰撞標(biāo)簽數(shù)估計(jì)方法，其步驟如下：

步驟一：利用一個(gè)固定的幀長(zhǎng)度Q的查詢命令Query(Q)開啟第一幀識(shí)別；

步驟二：利用第一幀成功識(shí)別的標(biāo)簽數(shù)估計(jì)初始標(biāo)簽數(shù)；

步驟三：從第二幀開始，利用估計(jì)的原始標(biāo)簽數(shù)減去本幀之前已被成功識(shí)別的標(biāo)簽數(shù)作為本幀幀長(zhǎng)度，如此循環(huán)開啟新的識(shí)別幀，直至完成所有標(biāo)簽的識(shí)別。

所述的步驟一包括：接收到Query(Q)命令后，標(biāo)簽生成隨機(jī)數(shù)r_i，假設(shè)標(biāo)簽選擇時(shí)隙的概率服從均勻分布，則標(biāo)簽選擇Q個(gè)時(shí)隙中任意一個(gè)時(shí)隙的概率為：

p＝1/Q (1)

有且僅有一個(gè)標(biāo)簽選擇某個(gè)特定時(shí)隙r_i(i＝1,2,....Q)的概率(即成功識(shí)別時(shí)隙的概率)為：

那么，對(duì)于全部n個(gè)標(biāo)簽，Q個(gè)時(shí)隙中有Q_e個(gè)成功識(shí)別時(shí)隙的期望值為：

所述的步驟二包括：包括基于TPELI的初始標(biāo)簽數(shù)估計(jì)和“偽解去除”；

所述的基于TPELI的初始標(biāo)簽數(shù)估計(jì)包括如下步驟：

1：利用統(tǒng)計(jì)得到的有效識(shí)別標(biāo)簽數(shù)IdN代替Q_e，并解方程(3)可以求出初始標(biāo)簽數(shù)n的估值設(shè)

2：利用線性插值法求解式(4)得到x的近似解，并對(duì)x取整得到初始標(biāo)簽數(shù)的估計(jì)值Q>1時(shí)式(4)對(duì)x連續(xù)，如果在大于0的區(qū)間存在x_l和x_u，使得[f(x_l)-IdN][f(x_u)-IdN]＜0，則必然存在x∈[x_l,x_u]，使f(x)＝IdN，除最大值以外，每個(gè)IdN對(duì)應(yīng)兩個(gè)解，小于最優(yōu)值的解定義為n_S，大于最優(yōu)值的定義為n_B。

所述的“偽解去除”包括：

3：利用相似三角形公式，可以得到直線與f(x)＝IdN的交點(diǎn)的估計(jì)值

經(jīng)過變形得到

4：確定下界x_S,l和上界x_S,u，對(duì)于小于最優(yōu)值的解x_S，f(1)<f(x_S)<f(Q)，所以選擇x_S,l＝1，x_S,u＝Q；

5：計(jì)算估計(jì)值，令

6：計(jì)算并比較函數(shù)值，如果[f(x_S,l)-IdN]·[f(x_S)-IdN]<0，則x_S,u＝x_S；

7:反之，如果[f(x_S,u)-IdN]·[f(x_S)-IdN]<0，則x_S,l＝x_S，轉(zhuǎn)第5步；

8：終止條件，如果[f(x_S,u)-IdN]·[f(x_S,l)-IdN]＝0，或者|x_S,l-x_s,u|<ε_Th，則迭代終止，其中ε_Th為迭代終止判斷門限，由所需的精度確定；

結(jié)束迭代后，令得到初始標(biāo)簽數(shù)的估計(jì)。

本發(fā)明的有益效果：

1、根據(jù)第一幀有效識(shí)別標(biāo)簽數(shù)，利用線性插值的方法估計(jì)初始標(biāo)簽數(shù)，進(jìn)而獲得剩余待識(shí)別的標(biāo)簽數(shù)，并按照幀長(zhǎng)度等于待識(shí)別標(biāo)簽數(shù)的最優(yōu)幀識(shí)別效率原則，自適應(yīng)動(dòng)態(tài)調(diào)整后續(xù)識(shí)別幀的長(zhǎng)度，以達(dá)到最佳的識(shí)別效率，實(shí)現(xiàn)全局吞吐量最大化，本發(fā)明提供的方法標(biāo)簽數(shù)估計(jì)誤差小于4％，迭代次數(shù)小于10次，優(yōu)于現(xiàn)有標(biāo)簽估計(jì)方法。與目前射頻識(shí)別領(lǐng)域廣泛采用的經(jīng)典動(dòng)態(tài)幀時(shí)隙Aloha算法相比，標(biāo)簽識(shí)別全局吞吐量性能提高到了1.5倍以上，最高達(dá)3倍，顯著提高了批量識(shí)別的識(shí)別效率，縮短了識(shí)別時(shí)間，降低了動(dòng)態(tài)環(huán)境下標(biāo)簽漏讀概率。

2、本發(fā)明提供的標(biāo)簽估計(jì)算法具有標(biāo)簽估計(jì)精度高、算法迭代次數(shù)少、適用于大規(guī)模批量識(shí)別標(biāo)簽應(yīng)用環(huán)境等優(yōu)點(diǎn)。

3、本發(fā)明提供的TPELI標(biāo)簽防碰撞方法全局吞吐量性能遠(yuǎn)高于目前國(guó)際標(biāo)準(zhǔn)廣泛采用的經(jīng)典動(dòng)態(tài)幀時(shí)隙ALOHA防碰撞方法。

附圖說明

圖1為本發(fā)明所述基于線性插值的射頻識(shí)別防碰撞標(biāo)簽數(shù)估計(jì)方法處理流程。

圖2為本發(fā)明所述實(shí)驗(yàn)和理論計(jì)算的識(shí)別效率曲線。

圖3為按照本發(fā)明所述線性插值方法所實(shí)施例的實(shí)驗(yàn)結(jié)果。

圖4為按照本發(fā)明所述基于TPELI算法的初始標(biāo)簽數(shù)估計(jì)所實(shí)施例的實(shí)驗(yàn)結(jié)果。

圖5為按照本發(fā)明所述基于TPELI算法的標(biāo)簽數(shù)估計(jì)方法所實(shí)施例的估計(jì)精度實(shí)驗(yàn)結(jié)果。

圖6為按照本發(fā)明所述基于TPELI算法的標(biāo)簽數(shù)估計(jì)方法所實(shí)施例的估計(jì)迭代次數(shù)實(shí)驗(yàn)結(jié)果。

圖7為按照本發(fā)明所述基于線性插值的射頻識(shí)別防碰撞標(biāo)簽數(shù)估計(jì)方法全局吞吐量性能對(duì)比實(shí)驗(yàn)結(jié)果。

具體實(shí)施方式

以下結(jié)合附圖對(duì)本發(fā)明進(jìn)一步說明。

參照?qǐng)D1所示，一種基于線性插值的射頻識(shí)別防碰撞標(biāo)簽數(shù)估計(jì)方法，其步驟如下：

步驟一：利用一個(gè)固定的幀長(zhǎng)度Q的查詢命令Query(Q)開啟第一幀識(shí)別；

步驟二：利用第一幀成功識(shí)別的標(biāo)簽數(shù)估計(jì)初始標(biāo)簽數(shù)；

步驟三：從第二幀開始，利用估計(jì)的原始標(biāo)簽數(shù)減去本幀之前已被成功識(shí)別的標(biāo)簽數(shù)作為本幀幀長(zhǎng)度，如此循環(huán)開啟新的識(shí)別幀，直至完成所有標(biāo)簽的識(shí)別。

所述的步驟一包括：接收到Query(Q)命令后，標(biāo)簽生成隨機(jī)數(shù)r_i，假設(shè)標(biāo)簽選擇時(shí)隙的概率服從均勻分布，則標(biāo)簽選擇Q個(gè)時(shí)隙中任意一個(gè)時(shí)隙的概率為：

p＝1/Q (1)

有且僅有一個(gè)標(biāo)簽選擇某個(gè)特定時(shí)隙r_i(i＝1,2,....Q)的概率(即成功識(shí)別時(shí)隙的概率)為：

那么，對(duì)于全部n個(gè)標(biāo)簽，Q個(gè)時(shí)隙中有Q_e個(gè)成功識(shí)別時(shí)隙的期望值為：

參照?qǐng)D2給出實(shí)驗(yàn)和理論計(jì)算的識(shí)別效率曲線，標(biāo)簽數(shù)從10到3000變化，間隔為10，Q＝1000。實(shí)驗(yàn)采用蒙特卡羅方法模擬標(biāo)簽識(shí)別過程得到的成功識(shí)別標(biāo)簽數(shù)IdN與待識(shí)別標(biāo)簽數(shù)n之間的關(guān)系，理論曲線為利用式(3)計(jì)算的結(jié)果；由圖2看出，實(shí)驗(yàn)結(jié)果與理論值變化趨勢(shì)是一致的，所以，第一幀識(shí)別結(jié)束后，利用統(tǒng)計(jì)得到的有效識(shí)別標(biāo)簽數(shù)IdN代替Q_e，并解方程(3)可以求出初始標(biāo)簽數(shù)n的估值

所述的步驟二包括：包括基于TPELI的初始標(biāo)簽數(shù)估計(jì)和“偽解去除”；

所述的基于TPELI的初始標(biāo)簽數(shù)估計(jì)包括如下步驟：

1：由圖2看出，實(shí)驗(yàn)結(jié)果與理論值變化趨勢(shì)是一致的，所以，第一幀識(shí)別結(jié)束后，利用統(tǒng)計(jì)得到的有效識(shí)別標(biāo)簽數(shù)IdN代替Q_e，并解方程(3)可以求出初始標(biāo)簽數(shù)n的估值設(shè)

2：利用線性插值法求解式(4)得到x的近似解，并對(duì)x取整，得到初始標(biāo)簽數(shù)的估計(jì)值Q>1時(shí)式(4)對(duì)x連續(xù)，如果在大于0的區(qū)間存在x_l和x_u，使得[f(x_l)-IdN][f(x_u)-IdN]＜0，則必然存在x∈[x_l,x_u]，使f(x)＝IdN，除最大值以外，每個(gè)IdN對(duì)應(yīng)兩個(gè)解，小于最優(yōu)值的解定義為n_S，大于最優(yōu)值的定義為n_B。

所述的“偽解去除”包括：

3：利用相似三角形公式，可以得到直線與f(x)＝IdN的交點(diǎn)的估計(jì)值：

經(jīng)過變形得到：

4：確定下界x_S,l和上界x_S,u，對(duì)于小于最優(yōu)值的解x_S，f(1)<f(x_S)<f(Q)，所以選擇x_S,l＝1，x_S,u＝Q；

5：計(jì)算估計(jì)值，令

6：計(jì)算并比較函數(shù)值，如果[f(x_S,l)-IdN]·[f(x_S)-IdN]<0，則x_S,u＝x_S；7；

反之，如果[f(x_S,u)-IdN]·[f(x_S)-IdN]<0，則x_S,l＝x_S，轉(zhuǎn)第5步；

8：終止條件，如果[f(x_S,u)-IdN]·[f(x_S,l)-IdN]＝0，或者|x_S,l-x_s,u|<ε_Th，則迭代終止，其中ε_Th為迭代終止判斷門限，由所需的精度確定；

結(jié)束迭代后，令得到初始標(biāo)簽數(shù)的估計(jì)。

參照?qǐng)D2給出實(shí)驗(yàn)和理論計(jì)算的識(shí)別效率曲線，標(biāo)簽數(shù)從10到3000變化，間隔為10，Q＝1000，n分別取100、400、7000、1000，1300，1600等6種情況進(jìn)行實(shí)驗(yàn)，其中n_S,l＝1，n_S,u＝1000，n_B,l＝1000，n_B,u＝10000，ε_Th＝0.01；由理論計(jì)算得Q_e分別為91、268、238、368、354、323，結(jié)果如附圖3所示。

由圖2可知，當(dāng)n＝1000時(shí)曲線取最大值達(dá)到最優(yōu)識(shí)別效率，此時(shí)根據(jù)Q_e解方程可以得到唯一解，除此之外，對(duì)于每個(gè)Q_e值方程有兩個(gè)解，一個(gè)是偽解，不符合實(shí)際，需要去除。由圖3可以看出，在正確去除偽解的前提下，相對(duì)誤差最大為1.29％，迭代次數(shù)最多為12次。當(dāng)n＝1000時(shí)兩個(gè)解相差很小，分別為1000和1002，都是真解，不存在偽解，并且迭代次數(shù)也非常少，分別為2次和1次。

圖3給出數(shù)值解法的實(shí)驗(yàn)結(jié)果，是利用Q_e理論值作為函數(shù)值實(shí)驗(yàn)的結(jié)果，圖4給出了利用TPELI算法對(duì)初始標(biāo)簽數(shù)進(jìn)行估計(jì)的實(shí)驗(yàn)結(jié)果，即利用第一幀識(shí)別結(jié)束時(shí)統(tǒng)計(jì)得到的IdN代替Q_e進(jìn)行迭代實(shí)驗(yàn)。IdN為利用蒙特卡羅方法模擬第一幀成功識(shí)別的標(biāo)簽數(shù)，不是確定值，每次實(shí)驗(yàn)會(huì)有所差別，所以對(duì)原始標(biāo)簽數(shù)n的每個(gè)取值進(jìn)行三次實(shí)驗(yàn)進(jìn)行對(duì)比。

由圖4可知，除了n＝1000時(shí)第一次實(shí)驗(yàn)相對(duì)誤差較大為7.2％，其余情況相對(duì)誤差均小于4％。迭代次數(shù)，除了n＝1000時(shí)第3次實(shí)驗(yàn)出現(xiàn)錯(cuò)誤以外，其余情況迭代次數(shù)均小于等于11次。n＝1000時(shí)實(shí)驗(yàn)出現(xiàn)一次相對(duì)誤差較大，一次程序異常未能正確結(jié)束的情況。從附圖1可知，由于n＝1000時(shí)曲線取最大值，有可能出現(xiàn)f(x_S,l)和f(x_S,u)相差太小甚至相等(即線性插值斜率太小，甚至為0)，導(dǎo)致線性插值式(6)的分母值過小或?yàn)?，出現(xiàn)誤差太大或插值出現(xiàn)異常。這種情況在實(shí)際工程實(shí)現(xiàn)時(shí)可以采取特殊的處理方法加以解決，譬如可以依據(jù)判斷IdN的值，如果接近式(6)的最大值，則無需解方程直接用Q作為待識(shí)別標(biāo)簽數(shù)的近似。

當(dāng)n與Q相差較大時(shí)，對(duì)應(yīng)每個(gè)IdN會(huì)出現(xiàn)兩個(gè)解，有一個(gè)是偽解，需要加以去除。根據(jù)圖3和圖4實(shí)驗(yàn)結(jié)果，在第一幀識(shí)別完成后，獲得兩個(gè)解n_S和n_B，n_S為小于最優(yōu)值的解，n_B為大于最優(yōu)值的解。根據(jù)圖1所示，首先令Q＝n_S，閱讀器廣播Query(Q_S)命令，接收命令后標(biāo)簽Tag_i(其中，i＝1,2…..n)，根據(jù)Q_S產(chǎn)生隨機(jī)數(shù)RS_i∈[1,Q_S]作為期望向閱讀器發(fā)送信息占用的時(shí)隙，以及長(zhǎng)度為u的隨機(jī)數(shù)PS_i作為識(shí)別號(hào)ID的替代碼，u遠(yuǎn)遠(yuǎn)小于識(shí)別號(hào)ID的長(zhǎng)度，以減少時(shí)間開銷。在發(fā)送Query(Q_S)命令之后間隔一定時(shí)間，閱讀器發(fā)送NextSlot命令開始按時(shí)隙逐一讀標(biāo)簽數(shù)據(jù)，標(biāo)簽收到NextSlot命令將自己的期望時(shí)隙RS_i減1，結(jié)果為零的標(biāo)簽響應(yīng)閱讀器，回傳PS_i，結(jié)果不為零則等待下一個(gè)NextSlot命令繼續(xù)減1，直到結(jié)果為0響應(yīng)標(biāo)簽，閱讀器根據(jù)標(biāo)簽回傳信息統(tǒng)計(jì)有效識(shí)別標(biāo)簽數(shù)NS。用同樣方法處理Query(Q_B)命令得到有效標(biāo)簽識(shí)別數(shù)NB，如果NS>NB則n_B為偽解，令Q₂＝NS，反之，如果NS<NB則n_S為偽解，令Q₂＝NB，下標(biāo)2表示第2幀的參數(shù)。

去除偽解后，基于TPELI算法的標(biāo)簽數(shù)估計(jì)精度實(shí)驗(yàn)結(jié)果如圖5所示，迭代次數(shù)實(shí)驗(yàn)結(jié)果如圖6所示。實(shí)驗(yàn)參數(shù)：Q＝700，n從20到2000變化，步進(jìn)為20，n_S,l＝1，n_S,u＝700，n_B,l＝700，n_B,u＝7000，圖5所示的估計(jì)精度仿真ε_Th分別取0.1和0.01兩種，附圖6所示的迭代次數(shù)仿真ε_Th分別取0.1、0.01和0.001。

由圖5可知迭代終止判斷門限ε_Th對(duì)估計(jì)精度并沒有明顯影響，ε_Th＝0.1和ε_Th＝0.01兩條曲線沒有表現(xiàn)出明顯差異性。兩條曲線均在Q值附近，即大約在n＝600到800之間，估計(jì)精度有較大波動(dòng)，其余情況估計(jì)標(biāo)簽數(shù)均能很好地跟蹤待識(shí)別標(biāo)簽數(shù)的變化。

由圖6看出，門限ε_Th對(duì)迭代次數(shù)有明顯影響，迭代次數(shù)隨門限的減小而增加，這是由于估計(jì)精度提高，直接導(dǎo)致線性插值循環(huán)迭代次數(shù)增加。ε_Th＝0.1時(shí)除了在n＝40時(shí)迭代次數(shù)為6次，對(duì)于其余n值迭代次數(shù)均未超過4次，ε_Th＝0.01時(shí)迭代次數(shù)小于等于10次，而＝0.001時(shí)迭代次數(shù)最高達(dá)到25次。附圖6還可以看出待識(shí)別標(biāo)簽數(shù)在Q值附近時(shí)迭代次數(shù)波動(dòng)較大，標(biāo)簽數(shù)大于800時(shí)迭代次數(shù)比較平穩(wěn)，基本保持在5次以內(nèi)。

對(duì)比圖5和圖6，綜合考慮估計(jì)精度和迭代次數(shù)，終止判斷門限ε_Th取0.1和0.01時(shí)均具有較好的性能，能滿足大多數(shù)批量識(shí)別標(biāo)簽的實(shí)時(shí)性和估計(jì)精度要求，而ε_Th取0.001會(huì)導(dǎo)致迭代次數(shù)大幅增多，給系統(tǒng)帶來大的處理延遲。

全局吞吐量是綜合評(píng)價(jià)的一個(gè)射頻識(shí)別系統(tǒng)的理想指標(biāo)，其表示完全識(shí)別所有待識(shí)別標(biāo)簽與所消耗的總時(shí)間之比，定義為：

為了更清楚的評(píng)估TPELI算法的性能，給出TPELI算法與EPC C1 G2等標(biāo)準(zhǔn)所采用的經(jīng)典DFSA防碰撞算法的性能仿真結(jié)果，DFSA采用Q算法實(shí)現(xiàn)幀長(zhǎng)度的動(dòng)態(tài)調(diào)整，仿真結(jié)果如附圖7所示。仿真參數(shù)為：TPELI中Q＝512，ε_Th＝0.01，待識(shí)別標(biāo)簽數(shù)從10到1000變化，步進(jìn)為20；Q算法中c＝0.1，為使初始幀長(zhǎng)與TPELI算法相等，Q算法中幀長(zhǎng)的對(duì)數(shù)(以2為底)設(shè)為9。

由圖7可以看出，除了待識(shí)別標(biāo)簽規(guī)模很小的情況(標(biāo)簽數(shù)小于20)，TPELI全局吞吐量性能明顯高于DFSA。從比值曲線清楚地看出，當(dāng)標(biāo)簽數(shù)大于20時(shí)，比值迅速攀升，當(dāng)待識(shí)別標(biāo)簽數(shù)大于100時(shí)比值大于1.5，即TPELI的全局吞吐量性能是DFSA的1.5倍，最好時(shí)超過3倍。