朱 軍,臧守濤,李 劍,李 汐,葉國(guó)駿
(1.安徽大學(xué) 電子信息工程學(xué)院,安徽 合肥 230601; 2. 中國(guó)科學(xué)院 上海微系統(tǒng)與信息技術(shù)研究所,上海 200050;3. 華為技術(shù)服務(wù)有限公司,上海 201206)
由于移動(dòng)設(shè)備和移動(dòng)流量的爆炸式增長(zhǎng),蜂窩系統(tǒng)正向5G無(wú)線系統(tǒng)發(fā)展.與4G網(wǎng)絡(luò)相比,5G網(wǎng)絡(luò)在頻譜、空口和網(wǎng)絡(luò)架構(gòu)上制訂了跨代的全新標(biāo)準(zhǔn),以滿足未來(lái)5G的應(yīng)用場(chǎng)景.這些新標(biāo)準(zhǔn)、新技術(shù),給5G無(wú)線網(wǎng)絡(luò)規(guī)劃帶來(lái)了挑戰(zhàn)[1-3].除此之外,大規(guī)模天線陣列增加了網(wǎng)絡(luò)優(yōu)化的難度,復(fù)雜的空間相關(guān)性使傳統(tǒng)的定點(diǎn)路測(cè)、定量覆蓋模型難以滿足快速、全面、準(zhǔn)確評(píng)估柵格級(jí)網(wǎng)絡(luò)性能的要求.傳統(tǒng)的網(wǎng)絡(luò)性能仿真雖能實(shí)現(xiàn)精確分析,但海量的計(jì)算量和時(shí)間開(kāi)銷使仿真無(wú)法在現(xiàn)場(chǎng)應(yīng)用[4].
隨著神經(jīng)網(wǎng)絡(luò)在信息、自動(dòng)化、醫(yī)學(xué)、經(jīng)濟(jì)等領(lǐng)域的廣泛應(yīng)用,已有使用神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)信道特征和網(wǎng)絡(luò)性能的報(bào)道[5-7].文獻(xiàn)[8]利用神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)無(wú)線正交頻分復(fù)用系統(tǒng)中的信道脈沖響應(yīng).文獻(xiàn)[9]使用神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)下行鏈路信道的質(zhì)量.文獻(xiàn)[10]利用神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)用戶信道中傳輸路徑的偏離角.這些方法均通過(guò)現(xiàn)網(wǎng)采集實(shí)測(cè)數(shù)據(jù),且將其作為訓(xùn)練樣本,由于樣本采集成本高,且網(wǎng)絡(luò)參數(shù)變更后需重新生成樣本數(shù)據(jù),因此這些方法很難真正應(yīng)用于實(shí)際網(wǎng)絡(luò).
文獻(xiàn)[11]利用神經(jīng)網(wǎng)絡(luò)在線預(yù)測(cè)5G信道狀態(tài)信息,但其預(yù)測(cè)的只有無(wú)線信道特性并沒(méi)有網(wǎng)絡(luò)性能,且存在預(yù)測(cè)參數(shù)有限、適用范圍小、數(shù)據(jù)采集周期長(zhǎng)等不足之處.文獻(xiàn)[12]基于11個(gè)基站發(fā)射機(jī)的路徑損耗數(shù)據(jù),利用前饋神經(jīng)網(wǎng)絡(luò)算法預(yù)測(cè)無(wú)線信道的路徑損耗,該方式在路損預(yù)測(cè)方面有較好的表現(xiàn),但不能在現(xiàn)網(wǎng)架設(shè)之前采集數(shù)據(jù),且通過(guò)發(fā)射機(jī)獲取數(shù)據(jù)成本較高.針對(duì)上述方法的不足,筆者提出利用時(shí)延仿真模型累積大量數(shù)據(jù),通過(guò)深度神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)數(shù)據(jù)特征,利用訓(xùn)練后的網(wǎng)絡(luò)模型預(yù)測(cè)柵格級(jí)用戶的網(wǎng)絡(luò)時(shí)延.
射線追蹤模型可描述為:利用傳輸射線模擬發(fā)射信號(hào),在傳輸過(guò)程中射線受到各種障礙,經(jīng)歷了發(fā)射、散射和衍射后,最終到達(dá)信號(hào)接收點(diǎn).通過(guò)射線追蹤模型,可計(jì)算信號(hào)的幅度、相位、到達(dá)角及傳輸時(shí)延等,進(jìn)而構(gòu)建出傳輸矩陣[13].射線追蹤模型中不包含小尺度信息,僅有射線幾何位置不能計(jì)算多徑效應(yīng),這些數(shù)據(jù)不足以構(gòu)建信道模型,因此使用相同場(chǎng)景下3GPP 38.901定義的統(tǒng)計(jì)模型,且補(bǔ)齊缺失的輸入?yún)?shù). 時(shí)延仿真模型的生成過(guò)程如圖1所示.
圖1 時(shí)延仿真模型的生成過(guò)程
5G無(wú)線仿真平臺(tái)以信道矩陣和工程參數(shù)作為輸入,計(jì)算信道衰落系數(shù)矩陣,需要在射線數(shù)據(jù)的基礎(chǔ)上補(bǔ)充天線布局信息、天線方向圖等.其中天線布局信息、天線方向圖采用3GPP 38.901標(biāo)準(zhǔn)的天線配置.信道矩陣按照3GPP 38.901定義的過(guò)程和方法計(jì)算,具體參數(shù)如表1所示.
表1 信道參數(shù)
信道系數(shù)的表達(dá)式為
(1)
其中:Frx,u,θ,Frx,u,φ分別是沿球面單位向量θ,φ方向接收天線的場(chǎng)方向圖;Ftx,s,θ,Ftx,s,φ分別是沿球面單位向量θ,φ方向發(fā)射天線的場(chǎng)方向圖;rrx,n,m,rtx,n,m分別是方位角為φn,m,AOA,φn,m,AOD的球面單位向量;drx,u,dtx,s分別為接收天線、發(fā)射天線的位置矢量.
假定消息包無(wú)須拆分、合并操作.此處的無(wú)線側(cè)時(shí)延主要考慮消息包從基站MAC層到用戶MAC層的下行單向時(shí)延.業(yè)務(wù)模型采用FULL BUFFER模型,單次傳輸?shù)臅r(shí)長(zhǎng)(無(wú)重傳)為4 ms, HARQ通道配置為8,消息包的最大重傳次數(shù)為4.消息包的傳輸時(shí)延為
T=t+N×Td,
(2)
其中:t為消息包的首次傳輸時(shí)延,N為消息包的重傳數(shù),Td為重傳的時(shí)間間隔.
統(tǒng)計(jì)用戶所有消息包的傳輸時(shí)延,可得到用戶的平均時(shí)延為
(3)
其中:TI表示用戶所有消息包的傳輸總時(shí)延,Pn為消息包總數(shù).
射線追蹤模型中,傳輸鏈路的輸出結(jié)果為射線起點(diǎn)(基站)和終點(diǎn)(用戶,用5 m×5 m柵格表示)間所有反射點(diǎn)的信息,包括基站和用戶間每條射線的起點(diǎn)、終點(diǎn)和中間反射點(diǎn)的3維坐標(biāo)、每條射線的路徑損耗和傳輸時(shí)長(zhǎng).
類比圖像識(shí)別中的特征提取, 筆者提出的三視圖特征模型將反射點(diǎn)分別投影到XOY、XOZ和YOZ平面上.3平面上反射點(diǎn)的投影位置均蘊(yùn)含了反射點(diǎn)3維空間信息的一部分,可從3投影平面中完整提取反射點(diǎn)3維空間信息.三視圖特征的提取過(guò)程如圖2所示.
圖2 三視圖特征提取過(guò)程
將反射點(diǎn)投影到3平面后,以基站為中心在X,Y軸的正負(fù)方向均延伸500 m,在Z軸方向以地面為起點(diǎn)向上延伸至基站上50 m.將3平面均劃分為32×32個(gè)網(wǎng)格,判定反射點(diǎn)的所在區(qū)間,構(gòu)建用于神經(jīng)網(wǎng)絡(luò)訓(xùn)練的數(shù)據(jù)集.該文使用卷積神經(jīng)網(wǎng)絡(luò)(convolutional neural network, 簡(jiǎn)稱CNN)訓(xùn)練模型,卷積神經(jīng)網(wǎng)絡(luò)由以下組成:
(1) 輸入層. 輸入層是整個(gè)神經(jīng)網(wǎng)絡(luò)的輸入,它描述了反射點(diǎn)在空間3平面的位置信息.
(2) 卷積層.卷積層對(duì)神經(jīng)網(wǎng)絡(luò)中的小塊進(jìn)行分析,得到抽象程度更高的特征.
(3) 池化層.池化層不改變3維矩陣的深度,但它能縮小矩陣的大小.在保持原有特征的基礎(chǔ)上,最大限度降低數(shù)組的維度.
(4) 全連接層.提取輸入特征后,進(jìn)行分類或回歸.
(5) 輸出層.對(duì)于回歸來(lái)說(shuō),輸出層只輸出一個(gè)節(jié)點(diǎn),這個(gè)節(jié)點(diǎn)是卷積神經(jīng)網(wǎng)絡(luò)對(duì)輸入數(shù)據(jù)的預(yù)測(cè)輸出.對(duì)于分類來(lái)說(shuō),輸出有多少類型,輸出層就輸出多少個(gè)節(jié)點(diǎn).
該文使用的卷積神經(jīng)網(wǎng)絡(luò)包含的卷積層、池化層、全連接層的參數(shù)如表2所示.
表2 卷積神經(jīng)網(wǎng)絡(luò)參數(shù)
時(shí)延相對(duì)誤差的表達(dá)式為
(4)
其中:Xi為網(wǎng)絡(luò)模型的預(yù)測(cè)值,Yi為實(shí)際值,n為用戶數(shù)量.
消息包重傳次數(shù)的相對(duì)誤差為
(5)
其中:Nmax為消息包重傳數(shù)之和的最大值.
誤差均值和誤差標(biāo)準(zhǔn)差的表達(dá)式分別為
(6)
(7)
信道路徑損耗(信道路損)是信道模型的一個(gè)重要指標(biāo),該文通過(guò)預(yù)測(cè)信道路損間接評(píng)估信道矩陣的準(zhǔn)確性.設(shè)置神經(jīng)網(wǎng)絡(luò)模型的訓(xùn)練樣本數(shù)為1.6×105.預(yù)測(cè)有以下2個(gè)場(chǎng)景:利用神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)同一基站下不同小區(qū)的信道路損(場(chǎng)景1);利用神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)不同基站下不同小區(qū)的信道路損(場(chǎng)景2).
3.2.1 利用神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)同一基站下不同小區(qū)的信道路損(場(chǎng)景1)
利用訓(xùn)練后的網(wǎng)絡(luò)模型預(yù)測(cè)信道路損,結(jié)果如圖3所示.由圖3可知,兩條曲線的趨勢(shì)大體吻合, 僅少數(shù)點(diǎn)存在差異.
圖3 場(chǎng)景1下預(yù)測(cè)的信道路損
場(chǎng)景1下預(yù)測(cè)的信道路損指標(biāo)的評(píng)估結(jié)果如表3所示,由表3可知,信道路損的相對(duì)誤差在5%以內(nèi),同時(shí)誤差均值為3.173 4 dB,可見(jiàn)預(yù)測(cè)值較準(zhǔn)確.
表3 場(chǎng)景1下預(yù)測(cè)的信道路損指標(biāo)的評(píng)估結(jié)果
3.2.2 利用神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)不同基站下不同小區(qū)的信道路損(場(chǎng)景2)
場(chǎng)景2下預(yù)測(cè)的信道路損如圖4所示. 由圖4可知,預(yù)測(cè)值與實(shí)際值的變化趨勢(shì)基本吻合, 僅少數(shù)點(diǎn)存在差異.
圖4 場(chǎng)景2下預(yù)測(cè)的信道路損
場(chǎng)景2下預(yù)測(cè)的信道路損指標(biāo)的評(píng)估結(jié)果如表4所示.由表4可知,信道路損的相對(duì)誤差為6.435 9%,誤差均值為7.670 1 dB,可見(jiàn)預(yù)測(cè)值較準(zhǔn)確.
表4 場(chǎng)景2下預(yù)測(cè)的信道路損指標(biāo)的評(píng)估結(jié)果
3.2.3 信道路損的預(yù)測(cè)分析
兩場(chǎng)景的信道路損相對(duì)誤差的概率分布如圖5所示.由圖5可知,場(chǎng)景1的預(yù)測(cè)結(jié)果明顯優(yōu)于場(chǎng)景2.兩場(chǎng)景預(yù)測(cè)結(jié)果差異較大的原因是場(chǎng)景1下的小區(qū)相似度較高.
圖5 兩場(chǎng)景的信道路損相對(duì)誤差的概率分布
無(wú)線側(cè)時(shí)延預(yù)測(cè)有以下3個(gè)場(chǎng)景: 利用神經(jīng)網(wǎng)絡(luò)模型直接預(yù)測(cè)網(wǎng)絡(luò)時(shí)延(場(chǎng)景3);將神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)的消息包重傳數(shù)作為回歸任務(wù)處理,間接預(yù)測(cè)網(wǎng)絡(luò)時(shí)延(場(chǎng)景4);將神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)的消息包重傳數(shù)作為分類任務(wù)處理,間接預(yù)測(cè)網(wǎng)絡(luò)時(shí)延(場(chǎng)景5).
3.3.1 網(wǎng)絡(luò)時(shí)延的直接預(yù)測(cè)(場(chǎng)景3)
利用神經(jīng)網(wǎng)絡(luò)模型直接預(yù)測(cè)網(wǎng)絡(luò)時(shí)延,神經(jīng)網(wǎng)絡(luò)模型的訓(xùn)練樣本數(shù)量為6×105,預(yù)測(cè)樣本為3個(gè)小區(qū)的部分用戶. 場(chǎng)景3下預(yù)測(cè)的網(wǎng)絡(luò)時(shí)延如圖6所示.由圖6可得, 預(yù)測(cè)值與實(shí)際值的變化趨勢(shì)基本吻合.
圖6 場(chǎng)景3下預(yù)測(cè)的網(wǎng)絡(luò)時(shí)延
場(chǎng)景3下預(yù)測(cè)的網(wǎng)絡(luò)時(shí)延指標(biāo)的評(píng)估結(jié)果如表5所示.由表5可知,相對(duì)誤差和誤差均值均較小.
表5 場(chǎng)景3下預(yù)測(cè)的網(wǎng)絡(luò)時(shí)延指標(biāo)的評(píng)估結(jié)果
3.3.2 基于回歸任務(wù)的網(wǎng)絡(luò)時(shí)延間接預(yù)測(cè)(場(chǎng)景4)
場(chǎng)景4下預(yù)測(cè)的消息包重傳數(shù)如圖7所示.由圖7可知,預(yù)測(cè)值與實(shí)際值的變化趨勢(shì)基本一致.
圖7 場(chǎng)景4下預(yù)測(cè)的消息包重傳數(shù)
場(chǎng)景4下預(yù)測(cè)的消息包重傳數(shù)指標(biāo)的評(píng)估結(jié)果如表6所示.由表6可知,與場(chǎng)景3相比,場(chǎng)景4下預(yù)測(cè)的相對(duì)誤差及誤差均值均偏大.
表6 場(chǎng)景4下預(yù)測(cè)的消息包重傳數(shù)指標(biāo)的評(píng)估結(jié)果
3.3.3 基于分類任務(wù)的用戶網(wǎng)絡(luò)時(shí)延間接預(yù)測(cè)(場(chǎng)景5)
場(chǎng)景5下預(yù)測(cè)的消息包重傳數(shù)如圖8所示.由圖8可知,預(yù)測(cè)值與實(shí)際值的變化趨勢(shì)基本一致,僅少數(shù)點(diǎn)存在差異.
圖8 場(chǎng)景5下預(yù)測(cè)的消息包重傳數(shù)
場(chǎng)景5下預(yù)測(cè)的消息包重傳數(shù)指標(biāo)的評(píng)估結(jié)果如表7所示.由表7可知,與場(chǎng)景4相比,場(chǎng)景5下預(yù)測(cè)的相對(duì)誤差及誤差均值均略有提高.
表7 場(chǎng)景5下預(yù)測(cè)的消息包重傳數(shù)指標(biāo)的評(píng)估結(jié)果
直接預(yù)測(cè)網(wǎng)絡(luò)時(shí)延、間接預(yù)測(cè)消息包重傳數(shù)的相對(duì)誤差概率分布分別如圖9,10所示.由圖9可得直接預(yù)測(cè)網(wǎng)絡(luò)時(shí)延的相對(duì)誤差集中在1%. 由圖10可知,場(chǎng)景4與場(chǎng)景5的相對(duì)誤差差異不大,二者消息包重傳數(shù)的相對(duì)誤差均集中在10%.
圖9 直接預(yù)測(cè)網(wǎng)絡(luò)時(shí)延的相對(duì)誤差概率分布 圖10 間接預(yù)測(cè)消息包重傳數(shù)的相對(duì)誤差概率分布
筆者提出的網(wǎng)絡(luò)時(shí)延預(yù)測(cè)方法可快速準(zhǔn)確預(yù)測(cè)柵格級(jí)用戶的網(wǎng)絡(luò)時(shí)延,不需采集海量路測(cè)數(shù)據(jù),節(jié)省了分析評(píng)估成本.筆者使用該文方法對(duì)上海市(占地6 340km2,柵格大小為5 m×5 m)進(jìn)行了網(wǎng)絡(luò)時(shí)延預(yù)測(cè), 使用4臺(tái)GPU服務(wù)器只需10 h左右即可完成網(wǎng)絡(luò)時(shí)延預(yù)測(cè),而采用時(shí)延仿真模型方法需5 d左右.該文三視圖特征模型僅提取射線反射點(diǎn)的特征,沒(méi)有完整提取用戶的時(shí)延特征,因此進(jìn)一步優(yōu)化三視圖特征模型為后續(xù)研究方向.