李翠然，季一凡，呂安琪，段寶峰

(蘭州交通大學(xué) 電子與信息工程學(xué)院,蘭州 730070)

0 引 言

高速列車運(yùn)行速度快、經(jīng)歷的地形(如高架橋、路塹、隧道和山地等)復(fù)雜多變。山地作為一種常見的高鐵地形，大量分布于我國(guó)西南和西北地區(qū)，例如，成貴高鐵沿線就穿越了大量的山地、隧道等。與其他地形相比，山地場(chǎng)景下的無線信號(hào)傳播有其自身的特點(diǎn)，該場(chǎng)景中除了列車周圍存在建筑物或樹木等散射體外，遠(yuǎn)端的山體同樣會(huì)對(duì)信號(hào)產(chǎn)生反射和散射。因此，在中國(guó)高鐵迅速發(fā)展的研究背景下，建立合理、準(zhǔn)確的山地地形信道模型是十分必要的。

目前，只有少量文獻(xiàn)針對(duì)山地地形的無線信號(hào)傳輸進(jìn)行了測(cè)量與建模研究[1-3]。信道模型的構(gòu)建主要基于部分測(cè)量數(shù)據(jù)，缺乏一定的普適性。文獻(xiàn)[4]采用射線跟蹤法對(duì)戶外場(chǎng)景和隧道內(nèi)的無線信號(hào)傳播特性進(jìn)行了建模研究，該模型計(jì)算復(fù)雜度較高，實(shí)現(xiàn)困難。文獻(xiàn)[5]基于隨機(jī)傳播圖理論構(gòu)建了髙鐵無線信道仿真模型，雖然模型簡(jiǎn)單，但其預(yù)測(cè)精度較低。文獻(xiàn)[6-8]針對(duì)開放空間、高架橋和路塹場(chǎng)景，構(gòu)建了幾何隨機(jī)分布模型(geometry-based stochastic model, GBSM)，假設(shè)散射體分布于以基站和移動(dòng)臺(tái)為焦點(diǎn)的多重共焦橢圓上。文獻(xiàn)[9-10]引入了單環(huán)模型和橢圓模型，對(duì)高速鐵路場(chǎng)景中的近、遠(yuǎn)端散射體進(jìn)行了有效區(qū)分。然而，上述信道模型均基于如下假設(shè)，即來自基站的無線電波僅從單一平面?zhèn)鞑ブ烈苿?dòng)臺(tái)，該假設(shè)條件對(duì)于地勢(shì)起伏較大的山地地形來說，顯然是不合理的。由此，學(xué)者們針對(duì)高鐵場(chǎng)景下的三維(3D)信道建模問題展開了研究。

文獻(xiàn)[11]以郊區(qū)和城市環(huán)境為研究場(chǎng)景，構(gòu)建了球形多輸入多輸出(multiple-input multiple-output, MIMO)信道模型，并對(duì)比分析了不同場(chǎng)景的空間相關(guān)性。文獻(xiàn)[12-13]分別針對(duì)深溝塹場(chǎng)景和隧道場(chǎng)景提出了一種基于3D圓柱的幾何隨機(jī)模型，研究了收發(fā)端天線距離和仰角對(duì)空時(shí)相關(guān)特性的影響。文獻(xiàn)[14]考慮了角度參數(shù)的時(shí)變性，將高架橋、路塹和開闊地場(chǎng)景下的2D多橢圓GBSM模型拓展到了3D。與2D模型不同的是，3D模型中考慮了散射體相對(duì)于基站和移動(dòng)臺(tái)的仰角。然而，上述信道建模方法并未區(qū)分近、遠(yuǎn)端散射體。在高鐵山地場(chǎng)景中，近、遠(yuǎn)端散射體對(duì)無線信號(hào)的反射和散射作用是不盡相同的，會(huì)影響3D MIMO信道建模的精確性和通用性。

針對(duì)上述缺陷，本文研究了高鐵山地場(chǎng)景下的3D MIMO無線信道建模問題。

1 山地場(chǎng)景的3D GBSM模型

高鐵山地場(chǎng)景下，由于鐵路沿線的建筑物、樹木以及遠(yuǎn)端的山體對(duì)無線信號(hào)的反射及散射作用，移動(dòng)臺(tái)接收信號(hào)是來自直射路徑和近、遠(yuǎn)端散射體產(chǎn)生的大量非直射路徑的合成信號(hào)。典型的高鐵山地地形見圖1。

MIMO技術(shù)作為4G LTE的關(guān)鍵技術(shù)之一，能夠在不增加系統(tǒng)帶寬需求的條件下改善傳輸可靠性[15]。本節(jié)以山地場(chǎng)景中下行鏈路2×2MIMO信道為研究對(duì)象，展開3D GBSM建模分析。圖2為考慮了視距(line of sight, LoS)和單次反射(single bounced, SB)路徑的3D GBSM模型。圖2中，由于基站作為發(fā)射端、移動(dòng)臺(tái)作為接收端，以移動(dòng)臺(tái)為中心、R為半徑的球面上分布著的散射體稱為近端散射體，以基站和移動(dòng)臺(tái)為焦點(diǎn)的橢圓體上分布著的散射體稱為遠(yuǎn)端散射體。近端散射體S(n1)(n1=1,2, …,N1)和遠(yuǎn)端散射體S(n2)(n2=1,2, …,N2)的數(shù)目分別為N1,N2。表1為3D GBSM模型的相關(guān)參數(shù)含義。

圖2 山地場(chǎng)景的3D GBSM模型Fig.2 3D GBSM model for mountain scenario

表1 3D GBSM模型中的參數(shù)Tab.1 Parameters in 3D GBSM model

考慮到信號(hào)經(jīng)過多次反射后功率會(huì)大大降低，且不同反射路徑反映了相同的信道統(tǒng)計(jì)特性[11]。因此，為降低GBSM模型的復(fù)雜度，只考慮LoS路徑和單次反射路徑，忽略了其他的多次反射路徑。圖2中點(diǎn)畫線為近端散射體的單次反射路徑，虛線為遠(yuǎn)端散射體的單次反射路徑。根據(jù)圖2的山地場(chǎng)景3D GBSM模型，基站天線p和移動(dòng)臺(tái)天線q之間的信道沖擊響應(yīng)(channel impulse response, CIR)是由LoS分量和近、遠(yuǎn)端散射體的SB分量疊加而成。CIR表達(dá)式為

hpq(t)=hLoS(t)+hSB1(t)+hSB2(t)

(1)

(1)式中，hLoS(t)，hSB1(t)和hSB2(t)可分別表示為[16]

(2)

(3)

(4)

(2)—(4)式中：K為萊斯因子；λ為載波波長(zhǎng)；fmax為移動(dòng)臺(tái)最大多普勒頻移；相位θS，θD分別為服從[-π,π)上均勻分布的隨機(jī)變量。

(5)

(6)

(7)

(8)

(9)

(5)—(9)式各參數(shù)的表達(dá)式見(10)—(11)式。

(10)

(11)

(12)

(12)式中：φ，φ∈[-π,π]，φ0，φ0∈[-π,π]分別為φ和φ的均值；k(k≥0)值的大小由散射體分布密度決定，散射體分布密度越大，k值越小[14]。由于近、遠(yuǎn)端散射體分布密度不同，可用k1，k2值加以區(qū)分。

圖3 3D信道模型的時(shí)變距離參數(shù)Fig.3 Time-varying distance parameters of 3D channel model

(13)

萊斯K因子與D(t)之間的關(guān)系可表示為[1]

(14)

2 3D GBSM理論模型統(tǒng)計(jì)特性

2.1 空時(shí)相關(guān)函數(shù)STCF

任意2個(gè)子信道hpq(t)和hp′q′(t)的空時(shí)相關(guān)函數(shù)(space-time correlation function, STCF)可表示為[19]

(15)

(15)式中：E{·}表示數(shù)學(xué)期望運(yùn)算；*表示復(fù)共軛運(yùn)算。由于LoS分量和近、遠(yuǎn)端散射體的SB分量相互獨(dú)立，(15)式可由不同分量的STCF之和表示為

(16)

LoS的STCF可寫為

(17)

近端散射體SB的STCF為

(18)

(18)式中，

(19)

遠(yuǎn)端散射體SB的STCF為

(20)

(20)式中，

(21)

2.2 空間互相關(guān)函數(shù)CCF

在 (15) 式，(17) 式，(18) 式和(20) 式中令τ=0，可以得到空間互相關(guān)函數(shù)(cross correlation function, CCF)的表達(dá)式為

(22)

LoS的CCF可寫為

(23)

近端散射體SB的CCF為

(24)

遠(yuǎn)端散射體SB的CCF為

(25)

2.3 時(shí)間自相關(guān)函數(shù)ACF

在(15) 式，(17) 式，(18) 式和(20) 式中令δT，δR=0，于是得到時(shí)間自相關(guān)函數(shù)(auto correlation function, ACF)的表達(dá)式

(26)

LoS的ACF可寫為

(27)

近端散射體SB的ACF為

(28)

遠(yuǎn)端散射體SB的ACF為

(29)

3 3D GBSM仿真模型

在3D GBSM理論模型中，考慮有限個(gè)散射體的情形，可以得到對(duì)應(yīng)的仿真模型。由 (16) 式，(23)—(25)式可得仿真模型的空間CCF，表示為

(30)

LoS的CCF為

(31)

近端散射體SB的CCF為

(32)

(32)式中，

(33)

遠(yuǎn)端散射體SB的CCF為

(34)

(34)式中，

(35)

由(16)式，(27)—(29)式可得仿真模型的時(shí)間ACF，可寫為

(36)

LoS的ACF為

(37)

近端散射體SB的CCF為

(38)

遠(yuǎn)端散射體SB的CCF為

(39)

4 仿真結(jié)果

圖4給出了不同的發(fā)射天線間隔δT和接收天線間隔δR對(duì)應(yīng)的空間CCF值。CCF是由直射分量、近端散射體反射分量和遠(yuǎn)端散射體反射分量相加所得，其中直射分量占主導(dǎo)?？梢钥闯?，隨著δR的增大，空間CCF的絕對(duì)值呈減小趨勢(shì)，而δT對(duì)空間CCF值的影響甚微。當(dāng)δR進(jìn)一步增大時(shí)，反射分量對(duì)CCF的影響增加，導(dǎo)致CCF值隨δT，δR的變化發(fā)生波動(dòng)。還可看出，當(dāng)δR值較大時(shí)，隨著δT的增大，CCF值波動(dòng)的幅度降低，波動(dòng)的頻次增加。因此，通過調(diào)整接收端的天線間隔即可獲得較低的空間CCF值。

圖4 不同天線間隔δT，δR的空間CCF絕對(duì)值Fig.4 Spatial CCF for different antenna spacing δT , δR

圖5 當(dāng)δT=λ, t=0 s時(shí)，2D及3D GBSM模型的 空間CCF值Fig.5 Spatial CCF for 2D and 3D GBSM model (when δT=λ and t=0 s)

圖6 當(dāng)t=0 s時(shí)，2D及3D GBSM模型的 時(shí)間ACF值Fig.6 Time ACF for 2D and 3D GBSM model (when t=0 s)

由圖5和圖6可看出：仿真模型和理論模型得到的空間CCF、時(shí)間ACF的曲線貼合度較高，驗(yàn)證了3D GBSM模型公式推導(dǎo)的正確性以及仿真過程的合理性。

由圖5可知，當(dāng)δR較小時(shí)(≤2λ)空間CCF受δR的影響顯著，但隨著δR的進(jìn)一步增大，空間CCF值的增加趨于緩慢。由圖6可知，當(dāng)時(shí)間間隔τ較小時(shí)(≤0.002 s)時(shí)間ACF受τ的影響顯著，但隨著τ的進(jìn)一步增大，時(shí)間ACF曲線將趨于平緩。由圖5、圖6還可看出：3D模型的空間CCF和時(shí)間ACF值均小于2D模型的值。這是由于，2D模型無法捕獲垂直平面上無線電波的空間多樣性，與3D模型相比，它高估了信號(hào)的空間相關(guān)性、低估了空間分集增益。

圖7、圖8分別為不同時(shí)刻(由移動(dòng)臺(tái)的移動(dòng)引起)的3D GBSM理論模型與仿真模型的空間CCF和時(shí)間ACF曲線。由圖7、圖8可知，不同時(shí)刻(t分別取值為0 s和2 s)的空間CCF和時(shí)間ACF曲線差異顯著。這是因?yàn)?，不同時(shí)刻對(duì)應(yīng)的D(t)值的變化引起了萊斯K因子值的變化，見(14)式。

圖7 當(dāng)δT=λ時(shí)，不同時(shí)刻的3D GBSM模型的 空間CCF值Fig.7 Spatial CCF for 3D GBSM model at different times (when δT=λ)

圖8 不同時(shí)刻的3D GBSM模型的時(shí)間ACF值Fig.8 Time ACF for 3D GBSM model at different times

圖9、圖10分別為基站和移動(dòng)臺(tái)之間的距離D(t)變化對(duì)3D GBSM模型的空間CCF和時(shí)間ACF的影響。由圖9和圖10可以看出：移動(dòng)臺(tái)無論處于近端區(qū)域還是遠(yuǎn)端區(qū)域，空間CCF值和時(shí)間ACF值都會(huì)隨著距離D(t)的增大呈下降趨勢(shì)，這是由于D(t)的增大引起了萊斯K因子值的減小。

圖9 距離D(t)與3D GBSM模型空間CCF值的關(guān)系Fig.9 Relationship between distance D(t) and spatial CCF

圖10 距離D(t)與3D GBSM模型時(shí)間ACF的關(guān)系Fig.10 Relationship between distance D(t) and time ACF

近、遠(yuǎn)端散射體分布密度(以參數(shù)k1，k2表征)對(duì)3D GBSM模型空間CCF性能的影響見圖11。由圖11可知，隨著k1，k2值的增大，空間CCF值呈減小趨勢(shì)。這說明，散射體分布密度較低的環(huán)境(較高的ki值)，3D模型的空間CCF值較小。此外，對(duì)比圖11a和圖11b，可看出，遠(yuǎn)端散射體分布密度k2的變化對(duì)空間CCF的性能影響更為顯著。

圖11 近、遠(yuǎn)端散射體分布密度與空間CCF值之間的關(guān)系Fig.11 Relationship between distribution density of near and remote scatters and spatial CCF value

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

本文研究了高鐵山地場(chǎng)景的3D MIMO無線信道建模問題，推導(dǎo)了信道沖激響應(yīng)函數(shù)CIR以及空時(shí)STCF、空間CCF和時(shí)間ACF的數(shù)學(xué)表達(dá)式，仿真得到結(jié)論：1)相比2D模型，提出的3D GBSM信道模型考慮了散射體相對(duì)于基站和移動(dòng)臺(tái)的仰角，并對(duì)近、遠(yuǎn)端散射體進(jìn)行了區(qū)分，能夠更真實(shí)地反映無線電波的空間多樣性；2)通過調(diào)整收發(fā)天線間隔δR，δT，可以獲得較低的空間CCF值；3)移動(dòng)臺(tái)在運(yùn)動(dòng)過程中會(huì)引起萊斯K因子值的變化，從而導(dǎo)致空間CCF和時(shí)間ACF值發(fā)生波動(dòng)；4)近、遠(yuǎn)端散射體分布密度較低信道環(huán)境下的空間CCF值較小。

未來的研究工作將圍繞以下內(nèi)容展開：模型仿真數(shù)據(jù)與實(shí)際測(cè)量數(shù)據(jù)的對(duì)比研究；分析反射路徑的數(shù)目、山地傳播環(huán)境的差異性對(duì)GBSM 3D信道模型的通用性影響；定性分析多維聯(lián)合多種參數(shù)下的GBSM 3D信道特性，研究信道模型隨影響參數(shù)的變化規(guī)律并提取參數(shù)值。