夏田奇,劉瑞華

(中國(guó)民航大學(xué) 電子信息與自動(dòng)化學(xué)院,天津 300300)

0 引 言

北斗衛(wèi)星導(dǎo)航系統(tǒng)(BDS)由我國(guó)自主發(fā)展、獨(dú)立運(yùn)行,是國(guó)家正在建設(shè)的重要空間信息基礎(chǔ)設(shè)施. BDS采用衛(wèi)星無(wú)線電導(dǎo)航業(yè)務(wù)(RNSS)和衛(wèi)星無(wú)線電測(cè)定業(yè)務(wù)(RDSS)雙模結(jié)構(gòu)體制,在航空、氣象、地質(zhì)、海洋等領(lǐng)域有著廣泛的應(yīng)用. 例如,在通用航空器監(jiān)視、指揮、位置信息收集、氣象信息發(fā)布等應(yīng)用領(lǐng)域,北斗短報(bào)文通信都能發(fā)揮出重要的作用[1-2]. 而在對(duì)信息傳輸實(shí)時(shí)性要求較高的應(yīng)用場(chǎng)合,對(duì)北斗短報(bào)文的傳輸成功率和傳輸延時(shí)這兩項(xiàng)性能指標(biāo)要求相對(duì)苛刻,故需要具體分析北斗短報(bào)文傳輸延時(shí)分布情況.

1 北斗RDSS通信性能測(cè)試

1.1 測(cè)試方法

文中測(cè)試使用某公司用戶機(jī)類型為二類用戶機(jī)設(shè)備，做用戶端北斗硬件部分與地面指揮端硬件部分,分別通過(guò)串口與計(jì)算機(jī)連接,編寫程序?qū)Υ谶M(jìn)行讀寫,實(shí)現(xiàn)北斗短報(bào)文發(fā)送和接收. 北斗用戶卡通信等級(jí)為3級(jí)非密用戶,短報(bào)文服務(wù)頻度為60 s,通信長(zhǎng)度為628 bit[3].

測(cè)試分為兩組:1) 由用戶端向指揮端發(fā)送測(cè)試報(bào)文;2) 由指揮端向用戶端發(fā)送測(cè)試報(bào)文.

測(cè)試報(bào)文內(nèi)容包括:報(bào)文生成日期時(shí)間,精確到千分之一秒;報(bào)文計(jì)數(shù),用于統(tǒng)計(jì)報(bào)文發(fā)送個(gè)數(shù).

發(fā)送時(shí),指揮端或用戶端記錄時(shí)間和報(bào)文計(jì)數(shù)并生成報(bào)文,同時(shí)通過(guò)串口向北斗設(shè)備發(fā)送報(bào)文,用戶端或指揮端通過(guò)串口解析到報(bào)文的同時(shí)記錄時(shí)間. 發(fā)送時(shí)以固定的頻度發(fā)送測(cè)試報(bào)文,接收時(shí)進(jìn)行CRC校驗(yàn)數(shù)據(jù)是否出錯(cuò),將報(bào)文內(nèi)容、接收時(shí)間以及校驗(yàn)結(jié)果保存成TXT格式文件.

1.2 測(cè)試數(shù)據(jù)

測(cè)試報(bào)文的一條保存內(nèi)容如下:

2018-01-21 18:14:25.357 | 24 54 58 58 58 00 1D 06 FC 60 62 06 FC 60 00 00 00 48 01 21 18 14 24 03 61 00 01 00 74 CRC正確.

翻譯成明文為:接收時(shí)間2018年1月21日18時(shí)14分25.357秒;接收方ID:457***;發(fā)送方ID:235***;發(fā)送時(shí)間:1月21日18時(shí)14分24.361秒;第1條報(bào)文;CRC校驗(yàn)正確.

我們?cè)趦蓚€(gè)多月間斷測(cè)試了約1 169個(gè)小時(shí),共發(fā)送70 151條報(bào)文,接收67 648條有效報(bào)文. 測(cè)試樣本如表1所示.

表1 北斗短報(bào)文測(cè)試樣本

其中序號(hào)為2、4、5和7的測(cè)試樣本為用戶端發(fā)指揮端收,共測(cè)試約701個(gè)小時(shí),發(fā)送42 116條報(bào)文,接收41 067條有效報(bào)文;序號(hào)為1、3和6的測(cè)試樣本為指揮端發(fā)用戶端收,共測(cè)試約467個(gè)小時(shí),發(fā)送28 035條報(bào)文,接收26 581條有效報(bào)文.

1.3 測(cè)試數(shù)據(jù)統(tǒng)計(jì)分析

1.3.1 傳輸成功率統(tǒng)計(jì)

北斗短報(bào)文傳輸成功率是評(píng)價(jià)北斗短報(bào)文通信性能的一個(gè)重要參數(shù)[4-5],定義為

(1)

總計(jì)發(fā)送短報(bào)文個(gè)數(shù)為Tsum=70 151條,接收有效的短報(bào)文個(gè)數(shù)為Rsum=67 648條,即短報(bào)文傳輸成功率psum=96.43%.

1.3.2 平均傳輸延時(shí)統(tǒng)計(jì)

北斗短報(bào)文平均傳輸延時(shí)定義為[6]:

(2)

2 北斗RDSS傳輸延時(shí)核密度估計(jì)

在概率密度分布未知的條件下,通過(guò)已知的樣本數(shù)據(jù)對(duì)未知的概率密度分布進(jìn)行估計(jì),從而預(yù)測(cè)概率密度分布. 常用的概率密度估計(jì)方法有兩種:參數(shù)估計(jì)和非參數(shù)估計(jì)[7].

參數(shù)法根據(jù)樣本數(shù)據(jù)頻率直方圖的輪廓推斷總體的概率密度分布,使用樣本數(shù)據(jù)對(duì)推斷的概率分布參數(shù)進(jìn)行估計(jì). 參數(shù)法憑借對(duì)樣本數(shù)據(jù)分布的推斷,常采用復(fù)雜的多元函數(shù)對(duì)概率密度分布參數(shù)進(jìn)行估計(jì),需要在計(jì)算機(jī)上使用最大似然估計(jì)等方法得到概率密度參數(shù)估計(jì)值,過(guò)程較為復(fù)雜不便.

(3)

式中:n為樣本總數(shù);h為帶寬;K為核函數(shù),核函數(shù)常選用以0為中心的對(duì)稱單峰概率密度函數(shù)[8]. 某點(diǎn)x處的概率密度估計(jì)值與該點(diǎn)附近所含有的樣本個(gè)數(shù)相關(guān). 若x附近樣本個(gè)數(shù)較多,則概率密度估計(jì)值較大,反之較小. 核函數(shù)K的形式對(duì)概率密度估計(jì)值的影響較小,帶寬h的大小會(huì)直接影響概率密度估計(jì)的結(jié)果. 帶寬h取較小的值,估計(jì)結(jié)果的曲線會(huì)為不光滑的折線,但能說(shuō)明每個(gè)樣本數(shù)據(jù)所包含的信息;取較大的值則會(huì)使較多的樣本數(shù)據(jù)對(duì)概率密度估計(jì)結(jié)果產(chǎn)生影響,估計(jì)結(jié)果的曲線也越光滑,但丟失了部分樣本數(shù)據(jù)的信息. 使用核密度估計(jì),可以得到每個(gè)樣本點(diǎn)的概率密度估計(jì)值,但不能得到像參數(shù)法那樣總體的概率密度函數(shù).

以1.2節(jié)中北斗短報(bào)文實(shí)測(cè)傳輸時(shí)延樣本1、2、3、4、5、6和7的數(shù)據(jù)為基礎(chǔ),使用數(shù)學(xué)工具M(jìn)ATLAB軟件中ksdensity函數(shù)對(duì)北斗短報(bào)文實(shí)測(cè)延時(shí)數(shù)據(jù)進(jìn)行核密度估計(jì),獲得傳輸延時(shí)分布的概率密度估計(jì)值并繪制曲線,如圖1所示.

利用多個(gè)含權(quán)重系數(shù)的高斯函數(shù)對(duì)概率密度估計(jì)曲線進(jìn)行擬合[9],其表達(dá)式為

基地選擇：海拔200～2 000 m；土壤為沙壤土或紅壤土，中性至微酸性(pH值5.5～7.0)，肥沃；坡度10°～30°，便于排灌。

(4)

式中:n為所用高斯函數(shù)個(gè)數(shù);ai、bi、ci為分布參數(shù),i=1,2，…，n. 擬合方法通過(guò)MATLAB軟件中的擬合工具實(shí)現(xiàn),所用高斯函數(shù)個(gè)數(shù)由函數(shù)f(x)曲線與概率密度估計(jì)值曲線的符合程度決定. 當(dāng)使用的高斯函數(shù)的個(gè)數(shù)n=6時(shí),f(x)曲線與實(shí)測(cè)延時(shí)概率密度估計(jì)值曲線最為接近,擬合所獲得的函數(shù)f(x)的參數(shù)如表2所示.

同時(shí)利用威布爾函數(shù)進(jìn)行擬合并且與高斯函數(shù)擬合的結(jié)果進(jìn)行對(duì)比,如圖2所示,高斯函數(shù)擬合結(jié)果對(duì)延時(shí)概率密度估計(jì)值曲線描述更加準(zhǔn)確.

3 北斗RDSS延時(shí)偽隨機(jī)數(shù)生成

根據(jù)已知概率密度曲線生成偽隨機(jī)數(shù)的方法主要有反函數(shù)法、變換法和舍選法,本文應(yīng)用的舍選法是馮·諾曼最早提出來(lái)的. 它的基本思想是:按照給定的概率密度分布函數(shù)f(x),對(duì)均勻分布的隨機(jī)數(shù)序列{r1,r2,r3,…,ri}進(jìn)行舍選. 舍選的規(guī)則是:在密度函數(shù)f(x)較大的地方,保留較多的隨機(jī)數(shù)ri;在密度函數(shù)f(x)較小的地方,保留較少的隨機(jī)數(shù)ri,這樣得到的子樣本中ri的分布將滿足概率密度分布函數(shù)的要求[10].

生成概率密度函數(shù)分布為f(x)的偽隨機(jī)數(shù)X的步驟如下[11]:設(shè)隨機(jī)變量X的概率密度函數(shù)為f(x),又存在實(shí)數(shù)a

1) 選取常數(shù)λ,使得λf(x)≤1,x∈(a,b);

3) 比較r2與λf(y),若r2≤λf(y),則令x=y;否則剔除r1和r2,重返步驟2).

如此重復(fù)循環(huán),產(chǎn)生的隨機(jī)數(shù)序列x1,x2,…,xn的分布由概率密度f(wàn)(x)確定.

本文按照上述方法和步驟分別生成1 000、10 000和100 000個(gè)傳輸延時(shí)偽隨機(jī)數(shù),進(jìn)行核密度估計(jì)后繪制曲線與箱線圖并與實(shí)測(cè)延時(shí)數(shù)據(jù)進(jìn)行比較,結(jié)果如圖3～5所示.

從三組數(shù)據(jù)的統(tǒng)計(jì)對(duì)比結(jié)果可以看出,隨著實(shí)驗(yàn)生成的隨機(jī)數(shù)的個(gè)數(shù)增加,生成隨機(jī)數(shù)的概率密度曲線與實(shí)測(cè)延時(shí)數(shù)據(jù)概率密度曲線吻合程度越來(lái)越高,箱線圖中兩者分布形狀也越來(lái)越接近,幾乎無(wú)法區(qū)分,直觀上說(shuō)明生成的延時(shí)隨機(jī)數(shù)與實(shí)測(cè)的延時(shí)是同分布的,滿足模擬北斗短報(bào)文傳輸延時(shí)偽隨機(jī)數(shù)的要求.

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

本文對(duì)北斗短報(bào)文通信性能進(jìn)行了長(zhǎng)時(shí)間的測(cè)試,獲取了大量的短報(bào)文通信測(cè)試樣本數(shù)據(jù), 在此基礎(chǔ)上,對(duì)北斗短報(bào)文傳輸延遲和平均傳輸成功率進(jìn)行了統(tǒng)計(jì)分析,利用多個(gè)含權(quán)重系數(shù)的高斯函數(shù)對(duì)短報(bào)文傳輸延時(shí)的核密度估計(jì)值曲線進(jìn)行擬合,采用舍選法模擬生成符合北斗短報(bào)文傳輸延時(shí)分布的偽隨機(jī)數(shù).論文內(nèi)容可為今后北斗短報(bào)文應(yīng)用場(chǎng)景中工程仿真提供模擬延時(shí)數(shù)據(jù).