侯春枝 劉永勝 吳振森 李清亮 盧昌勝 林樂(lè)科 陳安濤 張玉強(qiáng)
(1. 西安電子科技大學(xué)物理與光電工程學(xué)院,西安 710071;2. 中國(guó)電波傳播研究所 電波環(huán)境特性及模化技術(shù)重點(diǎn)實(shí)驗(yàn)室,青島266107;3. 中電科儀器儀表有限公司,青島 266555)
下一代無(wú)線通信技術(shù)要具備更高的數(shù)據(jù)傳輸速率、更低的時(shí)延和更好的移動(dòng)性,要達(dá)到相應(yīng)的技術(shù)指標(biāo)必須增加傳輸帶寬[1]. 但6 GHz以下的頻段已經(jīng)非常擁擠,無(wú)法為下一代移動(dòng)通信提供更大帶寬. 為了應(yīng)對(duì)這一挑戰(zhàn),未來(lái)移動(dòng)通信將使用毫米波頻段.由于可用帶寬范圍廣,毫米波可以提供高達(dá)每秒幾千兆比特的數(shù)據(jù)速率,但毫米波頻段同時(shí)具有更大的傳輸損耗. 由于毫米波頻段的波長(zhǎng)短,天線尺寸可以做的更小,因此可以采用大規(guī)模天線陣列技術(shù),利用高增益波束來(lái)補(bǔ)償其路徑損耗.
美國(guó)、日本和韓國(guó)等國(guó)家采用的毫米波頻段是28 GHz. 2016年7月第三代合作伙伴計(jì)劃(the 3rd Generation Partnership Project, 3GPP)制定并發(fā)布了5G信道標(biāo)準(zhǔn)化報(bào)告. 2019年11月在國(guó)際電聯(lián)世界無(wú)線通信會(huì)議(World Radiocommunication Conference 2019, WRC-19)上,國(guó) 際 電 信 聯(lián) 盟(International Telecommunication Union, ITU)為5G確定了更多的頻段,包括24.25~27.5 GHz、37~43.5 GHz、45.5~47 GHz、47.2~48.2 GHz和66~71 GHz[2]. 國(guó)內(nèi)外對(duì)28 GHz、32 GHz和39 GHz等毫米波段的信道測(cè)量和建模工作已經(jīng)很多[3-7]. 2015年,Maccartney G R等人進(jìn)行了28 GHz和73 GHz典型室內(nèi)辦公室環(huán)境的信道測(cè)量,研究了同極化和交叉極化方向的路徑損耗模型[8]. 2016年,趙雄文等人進(jìn)行了在室外微蜂窩環(huán)境中32 GHz的信道測(cè)量、建模、仿真和驗(yàn)證,并使用參數(shù)化方法和非參數(shù)化方法分析了大尺度信道參數(shù)[9]. 2013年,紐約大學(xué)的Shuai Nie等人使用滑動(dòng)相關(guān)信道測(cè)量系統(tǒng)研究了72 GHz室內(nèi)辦公室的傳輸和穿透特性[10]. 2014年,紐約大學(xué)George R等人研究了路徑損耗模型和信道統(tǒng)計(jì)數(shù)據(jù),提出在毫米波頻段使用定向天線能更好地應(yīng)用于未來(lái)的蜂窩通信[11].2017年,Xianyue Wu等人在辦公室環(huán)境中進(jìn)行了60 GHz毫米波信道測(cè)量和建模,研究三維空間中的60 GHz信道[12]. 2017年,張蕊等人研究了地物附加損耗模型等[13]. 總體上對(duì)于66~71 GHz頻段的研究仍然相對(duì)較少.
本文在會(huì)議室環(huán)境下進(jìn)行了65.5 GHz頻段的信道測(cè)量活動(dòng). 首先基于方向性路損模型,研究路徑損耗指數(shù)和陰影衰落特性;然后基于空間交替廣義期望最大化(space-alternating generalized expectationmaximization, SAGE)參數(shù)化分析算法,提取水平到達(dá)角、時(shí)延和幅度信息;最后根據(jù)多徑分量距離方法,對(duì)多徑分量進(jìn)行分簇,分析均方根時(shí)延擴(kuò)展和均方根角度擴(kuò)展特性. 本文旨在為5G毫米波60 GHz頻段鏈路和系統(tǒng)設(shè)計(jì)及其模擬仿真提供數(shù)據(jù)和模型支撐.
本文使用滑動(dòng)相關(guān)原理來(lái)進(jìn)行信道測(cè)試,如圖1所示. 本系統(tǒng)發(fā)射端主要由銣鐘、M8190任意波形信號(hào)發(fā)生器、E8267D矢量信號(hào)源、模擬信號(hào)源和N5152A上變頻器組成. 銣鐘提供時(shí)間參考信號(hào),M8190生成模擬測(cè)量信號(hào),E8267D將其變?yōu)? GHz中頻信號(hào),模擬信號(hào)源提供本振(local oscillation, LO)信號(hào),在N5152A上變頻器中將5 GHz中頻信號(hào)上變頻變?yōu)?5.5 GHz射頻信號(hào). N5152A上變頻范圍為57~66 GHz.
圖1 毫米波信道測(cè)量系統(tǒng)Fig. 1 Millimeter wave channel measurement system
接收端主要由N9029A下變頻器、N5183B信號(hào)源和M9703A數(shù)字化儀構(gòu)成. N9029A下變頻器射頻輸入范圍為50~75 GHz. N5183B為下變頻器提供本振信號(hào). N9029A輸出中頻信號(hào)頻率范圍為16 kHz~2.5 GHz,之后M9703A數(shù)字化儀對(duì)中頻信號(hào)進(jìn)行采樣. 本系統(tǒng)利用GPS信號(hào)對(duì)銣鐘進(jìn)行馴服,有非常好的同步精度,能夠支持系統(tǒng)進(jìn)行遠(yuǎn)距離測(cè)量.
發(fā)射天線和接收天線均采用喇叭天線,工作頻率為49.8~75.8 GHz,標(biāo)稱增益值為20 dBi. 65.5 GHz時(shí)喇叭天線的增益為21.68 dBi,E面和H面的半功率波束寬度分別為13.94°和13.83°,發(fā)射天線和接收天線均是垂直極化.
本文在室內(nèi)會(huì)議室環(huán)境下開展65.5 GHz的視距(line-of-sights, LoS)信道測(cè)量,測(cè)量場(chǎng)景如圖2所示. 測(cè)量帶寬為1 GHz,測(cè)量時(shí)E8267D中頻輸出的功率為?3 dBm,通過(guò)收發(fā)端直連校準(zhǔn)把上變頻器N5152A以及下變頻器N9029A對(duì)信號(hào)功率的影響(放大或衰減)消除掉. 發(fā)射端采用喇叭天線,固定在一個(gè)位置,使其為垂直極化天線. 接收端采用喇叭天線利用轉(zhuǎn)臺(tái)在水平方向上移動(dòng),形成虛擬陣列天線,水平方向移動(dòng)8次,每次移動(dòng)的水平距離為λ/2,測(cè)量了6個(gè)接收點(diǎn). 相鄰兩個(gè)點(diǎn)之間的距離為3 m. 發(fā)射天線和接收天線相對(duì)地面的高度均為1.8 m,墻的高度為2.75 m. 由于未來(lái)毫米波頻段會(huì)采用波束成形技術(shù),收發(fā)天線將會(huì)采用最大增益方向收發(fā)信號(hào),因此接收天線在任意測(cè)量點(diǎn)始終朝向發(fā)射天線. 具體測(cè)試參數(shù)信息如表1所示.
圖2 會(huì)議室測(cè)量場(chǎng)景Fig. 2 Meeting room measurement scene
表1 測(cè)試參數(shù)Tab. 1 Measurement parameters
路徑損耗模型是無(wú)線電系統(tǒng)設(shè)計(jì)的重要模型[14-16],分為全向路徑損耗模型和方向性路徑損耗模型. 通過(guò)合成所有發(fā)射接收天線波束的接收信號(hào),就能計(jì)算出每個(gè)測(cè)量位置的全向路徑損耗. 方向性路損模型是指取接收天線在所有旋轉(zhuǎn)方向上的最大功率或最小路徑損耗. 全向路徑損耗模型和方向性路損模型又都分為固定截距(close-in,CI)路徑損耗模型和浮動(dòng)截距(floating intercept,F(xiàn)I)路徑損耗模型. 對(duì)于LoS情況,按如下方式使用CI路徑損耗模型:
式中:LP(d)的單位為dB;d是距離,單位為m;LP,F(d0)是參考距離d0處的自由空間路徑損耗;線性斜率n是通過(guò)對(duì)數(shù)據(jù)進(jìn)行最小二乘法擬合得到的路徑損耗指數(shù);陰影效應(yīng)Xσ是均值為零和標(biāo)準(zhǔn)偏差為σ的高斯隨機(jī)變量,以dB為單位.
FI路徑損耗模型用于WINNER Ⅱ和3GPP標(biāo)準(zhǔn)中. 該模型需要兩個(gè)參數(shù),與式(1)有相似的形式:
式中:α是以dB為單位的FI;β是直線的斜率;Xσ是均值為零的陰影衰落隨機(jī)變量.
均方根時(shí)延擴(kuò)展[17]是無(wú)線電信道質(zhì)量的重要信道參數(shù),其是功率時(shí)延分布的二階矩,用于表征無(wú)線電信道的頻率選擇性. 為了提高參數(shù)分析準(zhǔn)確性,在進(jìn)行功率時(shí)延譜(power delay profile, PDP)處理時(shí),對(duì)20個(gè)快拍的數(shù)據(jù)進(jìn)行平均,均方根時(shí)延擴(kuò)展定義如下:
式中:Pl是第l條路徑的功率;τl是第l條路徑的時(shí)延;στ是平均附加時(shí)延. 當(dāng)計(jì)算均方根時(shí)延擴(kuò)展的時(shí)候,將低于閾值的所有路徑的功率設(shè)置為零,以減小噪聲的影響,此時(shí)有
根據(jù)3GPP空間信道模型規(guī)范所述,計(jì)算每個(gè)測(cè)量時(shí)刻的圓角度擴(kuò)展[18]. 利用角度估計(jì)值(例如:水平到達(dá)角、水平離開角、垂直到達(dá)角和垂直離開角)和每個(gè)多徑的接收功率計(jì)算角度擴(kuò)展,如下式所示:
式中,ψl是第l條路徑的到達(dá)角(水平到達(dá)角或者垂直到達(dá)角)或者離開角(水平離開角或者垂直離開角).
SAGE算法已經(jīng)廣泛應(yīng)用于無(wú)線信道參數(shù)估計(jì).SAGE算法是期望最大化的一種升級(jí),其由兩個(gè)步驟組成:E步和M步. SAGE算法的更多詳細(xì)理論和實(shí)現(xiàn)過(guò)程可以參考文獻(xiàn)[19-21]. 本文利用SAGE算法對(duì)多徑參數(shù)進(jìn)行提取,包括多徑相對(duì)時(shí)延、到達(dá)角、離開角、多普勒頻移和幅度信息. 在單輸入多輸出的場(chǎng)景中,接收天線由一個(gè)天線陣列組成,該天線陣列由位于不同位置的M個(gè)天線組成. 接收信號(hào)可表示為
式中:N(t)是高斯白噪聲;s(t;θl)是接收到的第l條路徑的信號(hào),可以表示為
式中:λ、e(?)和fm(?)分別為波長(zhǎng)、由?確定的單位向量和第m個(gè)天線元件的復(fù)電場(chǎng)方向圖.
多徑分量距離最初由文獻(xiàn)[22]引入,用于度量無(wú)線信道任意兩個(gè)多徑之間的關(guān)系[23]. 多徑分量距離定義為
由于多徑分量距離對(duì)角度域和時(shí)延域使用不同的處理方法,對(duì)于角度距離,給出了以下公式:
此式,適用于水平到達(dá)角和垂直到達(dá)角.
時(shí)延距離為
式中:?是合適的時(shí)延域和角度域的比例因子;τstd為時(shí)延的標(biāo)準(zhǔn)差;. 同時(shí)用時(shí)延比例因子表示時(shí)延域與角度域的比值. 本文測(cè)量的是二 維空間中的信道特性,因此可以簡(jiǎn)化垂直到達(dá)角.
根據(jù)2.1節(jié)介紹的方向性路損模型對(duì)會(huì)議室場(chǎng)景進(jìn)行建模,圖3是路徑損耗隨收發(fā)端距離變化曲線. 由圖3可知:CI路徑損耗模型的路徑損耗指數(shù)為2.090 34,陰影衰落標(biāo)準(zhǔn)差為0.385 14;FI路徑損耗模型路徑損耗指數(shù)為2.018 56,陰影衰落標(biāo)準(zhǔn)差為0.373 43. 圖3還給出了自由空間中隨著收發(fā)距離變化的路徑損耗,可以看出:CI路損模型和FI路損模型均與自由空間損耗符合得很好;路徑損耗指數(shù)與自由空間中的路徑損耗指數(shù)非常接近,說(shuō)明了本次測(cè)試的準(zhǔn)確性;CI路徑損耗模型的陰影衰落比FI路損模型的陰影衰落大,差別在1 dB內(nèi).
圖3 定向路徑損耗變化曲線Fig. 3 Directional path loss curve
利用SAGE算法提取信道沖激響應(yīng)中的時(shí)延、水平到達(dá)角和幅度信息. 運(yùn)行SAGE算法時(shí)參數(shù)設(shè)置如下:路徑條數(shù)為200,迭代次數(shù)為50,多徑功率閾值為?22 dB,可以從200條徑中選取有用的多徑信號(hào). 本文選取第1個(gè)測(cè)量點(diǎn)分析SAGE算法的提取結(jié)果,從圖4中可以看出方位角在85°附近、附加時(shí)延為2 ns時(shí),信號(hào)最強(qiáng). 由于接收天線的E面和H面半功率波束寬度約為14°,因此到達(dá)方位角大約為 83°~96°,相鄰兩個(gè)點(diǎn)的時(shí)延相差很小. 由于實(shí)際測(cè)量環(huán)境中會(huì)議室兩邊都是玻璃窗,因此本次測(cè)量中提取到的多徑信號(hào)主要來(lái)自鏡面反射.角度域的比例因子 ?=1.5,LoS徑的閾值為0.12,常規(guī)閾值為0.25. 本文對(duì)第一個(gè)測(cè)量點(diǎn)進(jìn)行分簇,如圖5所示,分為了兩個(gè)簇,由于主要是LoS和玻璃墻壁的鏡面反射引起的,因此符合客觀事實(shí). 包括紅色圓點(diǎn)
根據(jù)2.4節(jié)介紹的分簇算法(即多徑分量距離)來(lái)進(jìn)行分簇. 分簇時(shí),閾值的選擇要使分簇結(jié)果符合實(shí)際情況,并且簇的邊界在時(shí)延域和角度域上清晰可見. 本文進(jìn)行大量的試驗(yàn)可以得到時(shí)延域和的簇主要是由LoS引起的,另一個(gè)簇是由玻璃墻引起的,即體現(xiàn)出65.5 GHz頻率下傳播距離短,環(huán)境中的簇具有稀疏性.
圖4 第1個(gè)測(cè)量點(diǎn)提取到的多徑成分Fig. 4 Multipath component extracted from the first measurement point
圖5 第1個(gè)測(cè)量點(diǎn)的多徑分簇Fig. 5 Multipath clustering at the first measurement point
根據(jù)2.2節(jié)介紹的時(shí)延擴(kuò)展和角度擴(kuò)展對(duì)SAGE算法提取的時(shí)延、水平到達(dá)角和幅度進(jìn)行分析.SAGE算法提取的多徑是根據(jù)相同的功率值作為標(biāo)準(zhǔn)來(lái)篩選的,最終結(jié)果如圖6和7所示. 由圖6可知,根據(jù)所有徑的時(shí)延擴(kuò)展擬合出的均方根時(shí)延擴(kuò)展均值為?9.010 3 dBsec,標(biāo)準(zhǔn)差為0.148 65 dBsec,這些值與表2中3GPP TR38.900的值相近. 由圖7可知,擬合出所有徑的均方根角度擴(kuò)展均值為0.809 57 dB°,標(biāo)準(zhǔn)差為0.034 396 dB°,與3GPP TR38.900[24]中室內(nèi)65.5 GHz的時(shí)延擴(kuò)展與角度擴(kuò)展的均值和標(biāo)準(zhǔn)差十分接近.
圖6 均方根時(shí)延擴(kuò)展概率累積分布Fig. 6 Probability cumulative distribution of root mean square delay spread
圖7 均方根角度擴(kuò)展概率累積分布Fig. 7 Probability cumulative distribution of root mean square angle spread
表2 關(guān)于時(shí)延擴(kuò)展和角度擴(kuò)展參數(shù)3GPP TR38.900和測(cè)試結(jié)果的比較Tab. 2 A comparison between 3GPP TR38.900 and measurement in the delay spread and angle spread
為了研究5G 60 GHz毫米波頻段的信道特性,使用時(shí)域信道探測(cè)儀系統(tǒng)進(jìn)行65.5 GHz室內(nèi)LoS信道測(cè)量活動(dòng),獲得信道沖激響應(yīng). 基于路徑損耗模型,分析了路徑損耗指數(shù)和陰影衰落特性. 基于SAGE參數(shù)化分析方法和多徑分量距離方法,提取出水平到達(dá)角、時(shí)延和幅度,并對(duì)這些參數(shù)進(jìn)行分析,獲得多徑的分簇情況,研究了均方根時(shí)延擴(kuò)展和均方根角度擴(kuò)展特性. 結(jié)果表明:LoS情況下,方向性路損模型在CI和FI兩種方式的擬合下與自由空間路徑損耗模型符合得很好,并且路徑損耗指數(shù)與自由空間中的路徑損耗指數(shù)非常接近,CI路徑損耗模型的陰影衰落比FI路損模型的陰影衰落大,差別在1 dB內(nèi);電波發(fā)生的多徑信號(hào)主要是由鏡面反射引起的,鏡面反射在本文會(huì)議室場(chǎng)景中成為主要的傳播方式;在65.5 GHz下,傳播距離短,簇具有稀疏性. 時(shí)延擴(kuò)展和角度擴(kuò)展的結(jié)果與3GPP的結(jié)果相近. 作者將進(jìn)一步進(jìn)行其他會(huì)議室以及室內(nèi)場(chǎng)景的測(cè)量以及分析,這些分析將為5G毫米波會(huì)議室等室內(nèi)場(chǎng)景基站布局提供更可靠的統(tǒng)計(jì)數(shù)據(jù).