薛翠薇 朱秋明,,3 楊穎 陳小敏 廖志忠

(1. 南京航空航天大學(xué) 江蘇省物聯(lián)網(wǎng)與控制技術(shù)重點實驗室,南京 211100;2. 中國空空導(dǎo)彈研究院,洛陽 471009;3. 英國赫瑞瓦特大學(xué)工程與物理科學(xué)學(xué)院,愛丁堡 EH144AS)

三維空間多天線相關(guān)性及互耦影響研究

薛翠薇1朱秋明1,2,3楊穎1陳小敏1廖志忠2

(1. 南京航空航天大學(xué) 江蘇省物聯(lián)網(wǎng)與控制技術(shù)重點實驗室,南京 211100;2. 中國空空導(dǎo)彈研究院,洛陽 471009;3. 英國赫瑞瓦特大學(xué)工程與物理科學(xué)學(xué)院,愛丁堡 EH144AS)

針對三維傳播環(huán)境下多天線系統(tǒng)空域相關(guān)性評估問題,首先將達(dá)波信號二維角度譜建模為更為通用的von Mises-Fisher(vMF)分布,并推導(dǎo)了不同陣列流型天線空域相關(guān)系數(shù)的解析式和近似表達(dá)式;然后分析互耦機(jī)理及其對接收信號矢量的影響,據(jù)此推導(dǎo)了互耦效應(yīng)下天線空域相關(guān)性的閉式表達(dá)式,并給出了存在互耦的相關(guān)性與無互耦相關(guān)性之間的內(nèi)在關(guān)系．仿真結(jié)果表明,本文推導(dǎo)的天線相關(guān)性表達(dá)式與數(shù)值計算結(jié)果非常吻合．另外,小角度擴(kuò)展時,考慮互耦效應(yīng)的天線空域相關(guān)曲線圍繞著無互耦時的相關(guān)性曲線上下波動,且面型陣列具有更好的抗互耦能力．

陣列流型;空間相關(guān)性;互耦效應(yīng);von Mises-Fisher分布

DOI 10.13443/j.cjors.2016040601

引 言

大規(guī)模多輸入多輸出(Massive Multiple-input Multiple-output,Massive MIMO)技術(shù)能夠成倍地提高通信系統(tǒng)容量和頻譜資源利用率,已成為未來移動通信系統(tǒng)的關(guān)鍵技術(shù)之一[1]．同時,智能終端設(shè)備體積越來越緊湊,天線間距進(jìn)一步減小,導(dǎo)致互耦效應(yīng)不可忽略,因而包括第四代移動通信技術(shù)(4thGeneration Long Term Evolution,4G LTE)和第五代 (5thGeneration,5G)移動通信技術(shù)在內(nèi)的各種信道模型均考慮了互耦因素[2-3]．

針對平面波入射的二維空間多天線信號相關(guān)性研究已較成熟,如文獻(xiàn)[4]假設(shè)達(dá)波信號水平角度譜服從均勻、高斯、拉普拉斯和von Mises(vM)分布,推導(dǎo)給出了多天線信號的空域相關(guān)性解析式．然而,近年來信道實測數(shù)據(jù)表明接收信號在垂直面上存在一定的角度擴(kuò)展(Azimuth Spread, AS),分布區(qū)間通常為[16°,39°] ,對于高樓密集的城區(qū)環(huán)境可達(dá)[0°,50°]．因而,隨著4G LTE標(biāo)準(zhǔn)化和5G移動通信標(biāo)準(zhǔn)的發(fā)展,三維(3-Dimensional, 3D)信道模型也成為目前通信領(lǐng)域的研究熱點．Mammasis等利用數(shù)學(xué)方法對大量實測數(shù)據(jù)分析后指出,接收信號的方位角(Azimuth of Arrival, AOA)和俯仰角(Elevation of Arrival, EOA)之間存在相關(guān)性,二維聯(lián)合角度譜可建模為von Mises-Fisher(vMF)分布,并利用數(shù)值積分方法分析了信號的空域相關(guān)性[5];文獻(xiàn)[6]則假設(shè)接收端附近散射體服從vMF分布,提出了一個3D-MIMO信道模型,并通過仿真研究了接收信號的相關(guān)特性．需要指出的是,若角度譜服從vMF分布,空域相關(guān)系數(shù)理論推導(dǎo)非常困難,目前研究結(jié)論僅限于數(shù)值積分或仿真方法,文獻(xiàn)[7]雖然推導(dǎo)給出了三種不同天線陣列的空域相關(guān)性表達(dá)式,但最終結(jié)果仍非常復(fù)雜．

針對互耦效應(yīng)影響的問題,文獻(xiàn)[8-9]研究了二維空間互耦因素對相關(guān)性及系統(tǒng)容量的影響,指出互耦效應(yīng)在一定程度上會降低MIMO信道容量;周杰等研究了三維均勻散射空間下互耦對不同陣列相關(guān)性及容量的影響,表明互耦效應(yīng)對面陣和圓陣的容量影響要比線陣小[10]．然而,針對更為通用的vMF角度譜(非均勻散射場景),多天線空域相關(guān)性理論及互耦效應(yīng)影響,尚有很多問題需要進(jìn)一步深入研究．

針對上述問題,本文基于三維空間信號空域相關(guān)性定義,推導(dǎo)簡化了vMF角度譜下不同天線陣列空域相關(guān)性的解析式;并針對常見的小角度擴(kuò)展傳播場景,推導(dǎo)給出了近似計算方法．在此基礎(chǔ)上,文中還推導(dǎo)了互耦因素對空域相關(guān)系數(shù)影響的表達(dá)式,并給出針對vMF角度譜的數(shù)值仿真結(jié)果．

1 天線陣列相關(guān)性理論

實際散射環(huán)境下,到達(dá)接收端的電波信號具有非全向性以及非均勻性的特點,從而導(dǎo)致不同位置的天線陣元接收信號具有不同程度的相關(guān)性． 圖1以線型天線陣列接收系統(tǒng)為例給出了電波信號入射模型,其中:φ,Δφ,φ0分別表示信號AOA,角度區(qū)間和方位角均值；θ,Δθ,θ0分別表示信號EOA,角度區(qū)間和俯仰角均值．

圖1 多天線系統(tǒng)模型

對于N元接收天線,接收信號矢量可表示為

y(φ,θ) =[y1(φ,θ),…,yN(φ,θ)]

=[a1(φ,θ)g1(φ,θ),…,aN(φ,θ)gN(φ,θ)].

(1)

式中:gk(φ,θ)(k=1,2,…,N)表示各天線增益;a(φ,θ)=[a1(φ,θ),a2(φ,θ),…,aN(φ,θ)]為天線陣列流型矢量．第k根天線接收信號的平均功率可表示為

k=1,2,…,N.

(2)

根據(jù)定義,第k根與第i根(k,i=1,2,…,N)天線之間的空域互相關(guān)系數(shù)可表示為

(3)

式中,P(φ,θ)表示達(dá)波AOA和EOA二維聯(lián)合角度譜．

考慮不同流型的天線陣列,假設(shè)各陣元均為增益歸一化的全向天線,不同陣列接收信號的功率歸一化(Pk=1)矢量可表示為:

(4)

(5)

(6)

式中:yL(φ,θ)為線陣信號矢量,d為相鄰天線間距,N為天線的總數(shù)目;yC(φ,θ)為圓陣信號矢量,φi(i=1,2,…,N)為圓陣天線的方向角;yR(φ,θ)為面陣的信號矢量,vec表示矢量化函數(shù)．

2 三維空間不同陣列空域相關(guān)性

2.1 三維空間信號二維角度譜

早期的標(biāo)準(zhǔn)化信道模型只考慮一維角度譜,如3GPP信道標(biāo)準(zhǔn)規(guī)定郊區(qū)宏蜂窩、市區(qū)宏蜂窩和市區(qū)微蜂窩三種場景下,角度譜服從不同均值與AS的拉普拉斯分布[11],WINNER模型和COST2100模型則規(guī)定角度譜服從正態(tài)或?qū)?shù)正態(tài)分布[12]．近年來,隨著三維傳播理論研究深入及信道實測數(shù)據(jù)積累,3D-MIMO信道建模成為研究熱點,比如假設(shè)AOA及EOA服從二維拉普拉斯分布的3D-MIMO模型[13-14];文獻(xiàn)[5]則進(jìn)一步指出二維角度譜用vMF分布逼近實測數(shù)據(jù)效果更好,該分布目前已在MIMO信道建模中得到了應(yīng)用[6]．vMF分布的概率密度函數(shù)可表示為

(7)

式中:0≤φ≤2π;0≤θ≤π;k因子表示功率的集中程度;Im/2-1(k)表示零階修正貝塞爾函數(shù)．需要指出的是,當(dāng)m=3時,式(7)表示二維譜分布;當(dāng)m=2時,該分布退化為一維譜分布(即vM分布)．

基于文獻(xiàn)[5]的實測數(shù)據(jù),圖2(a)給出了兩種k值的vMF聯(lián)合概率密度函數(shù)．為了直觀觀察兩種情況下入射信號的角度分布情況,筆者利用文獻(xiàn)[16]方法模擬產(chǎn)生了10 000支路的入射信號,對應(yīng)的二維角度在球面上的分布如圖2(b)所示．由圖2可看出:vMF分布與k的取值密切相關(guān),k值越大,分布越集中密集;k值越小,分布越趨向于均勻,直至趨向二維均勻分布(k=0)．

(a) 二維角度譜理論分布

(b) 入射信號球面分布(10 000支路)圖2 三維空間信號角度譜及其分布

2.2 vMF空域相關(guān)系數(shù)解析式

根據(jù)相關(guān)系數(shù)定義,假設(shè)天線接收信號的功率歸一化,任意兩個天線陣元k,i的空域相關(guān)系數(shù)可表示為

(8)

(9)

利用歐拉公式,公式(9)可化簡為實部和虛部的形式,二者可分別表示為：

cos(zcos(φ+τ)sinθ)sinθdφdθ;

(10)

sin(zcos(φ+τ)sinθ)sinθdφdθ．

(11)

根據(jù)三角函數(shù)的級數(shù)展開公式

(12)

可得:

sin (θ)ek(sin θsin θ0cos(φ-φ0)+cos θcos θ0)dφdθ;

(13)

(14)

式中,Jn(Dk,i)表示第一類n階貝塞爾函數(shù)．

2.3 空域相關(guān)系數(shù)近似表達(dá)式

實際傳播環(huán)境下,達(dá)波信號AS通常都比較小,比如3GPP標(biāo)準(zhǔn)公布的郊區(qū)宏蜂窩、市區(qū)宏蜂窩、市區(qū)微蜂窩三種場景下基站端信號AS約小于5°[11],WINNER模型提供的AS參數(shù)一般也比較小[12];并且三維場景下的實測數(shù)據(jù)表明每簇信號的AS范圍一般小于10°[15]．鑒于該傳播特點,空域相關(guān)系數(shù)式(12)和(13)可進(jìn)一步近似化簡．首先,令v=φ-φ0,q=θ-θ0,則

φ=v+φ0,θ=q+θ0

(15)

將式(15)代入式(9)中可得

ejzcos(v+φ0+τ)sin(q+θ0)dvdq．

(16)

由于v,q很小,采用近似式

cosv≈1,sinv≈v,cosq≈1,sinq≈q,

(17)

式(16)可化簡為

e-jz(qvsin(φ0+τ)cos θ0-vsin(φ0+τ)sin θ0)+kdvdq．

(18)

(19)

利用積分公式對變量v積分可得

(20)

鑒于q很小,對式(20)化簡并對變量q積分,最終空域相關(guān)系數(shù)的近似表達(dá)式可表示為

(21)

為了驗證本文空域相關(guān)系數(shù)近似算法的準(zhǔn)確性及實用性,又不失一般性,我們參考文獻(xiàn)[5]提供的城市環(huán)境下二維vMF角度譜實測數(shù)據(jù)k=3.95,θ0=84.479°,φ0=232.48°,仿真比較近似表達(dá)式的均方誤差隨天線間距和AS的變化情況，結(jié)果如圖3所示.圖中AS的變化范圍為[0° 10°],天線歸一化間距d/λ∈[0，3],均方誤差定義為

(22)

圖3 近似表達(dá)式的均方誤差

由圖3可以看出:均方誤差隨AS的增大而增大,原因在于近似式(20)的誤差與AS密切相關(guān);隨著d/λ的增大,均方誤差迅速減小,但呈現(xiàn)一定的起伏．另外,當(dāng)AS小于5°時,均方誤差約小于0.000 1;當(dāng)AS小于10°時,均方誤差約小于0.005．

圖4 不同AS下空域相關(guān)性模值

圖4仿真給出了AS分別為2°,5°,10°時,利用本文方法獲得的任意兩天線空域相關(guān)系數(shù)隨歸一化間距的變化情況．綜合圖3和圖4的結(jié)果可看出：當(dāng)AS小于等于5°時,近似表達(dá)式和解析式的相關(guān)曲線幾乎重合,誤差很小,可用近似表達(dá)式代替解析式或數(shù)值積分結(jié)果,從而簡化空域相關(guān)性的計算;當(dāng)AS大于5°時,兩者相關(guān)曲線誤差有所增大．

3 互耦效應(yīng)影響機(jī)理及分析

3.1 互耦機(jī)理及其對信號影響

當(dāng)天線距離很近時,每個陣元上接收的電壓信號感應(yīng)出的電流信號會激勵出一個新磁場,并影響相鄰陣元上的信號,導(dǎo)致天線陣列中的天線方向圖發(fā)生畸變,即互耦效應(yīng)．互耦效應(yīng)使得陣列天線信號之間發(fā)生串?dāng)_,考慮互耦因素的信號矢量可改寫為

yc=C*y．

(23)

式中,C為互耦矩陣．對于任意天線陣列,假設(shè)各天線特性相同,負(fù)載阻抗相等,互耦矩陣可由下式計算獲得:

C=(ZL+Z)-1(zL+zA)．

(24)

式中:ZL表示負(fù)載阻抗矩陣;Z為互阻抗矩陣．zL和zA分別為負(fù)載阻抗和天線自阻抗．

假設(shè)各陣元均為歸一化的全向半波偶極子天線,長度LA=λ/2,載波fc=2 GHz,對于任意陣列的任意兩個天線,設(shè)定其歸一化距離d/λ=0.1,0.6,1.1,2,俯仰角θ=90°,圖5給出了考慮互耦因素情況下,信號矢量的模值,圖中ync表示不考慮互耦的接收信號矢量模值．由圖可以看出：不考慮互耦時的信號矢量模值恒為1;考慮互耦時,d/λ值越小,互耦對接收信號矢量影響越顯著;當(dāng)d/λ值小于1時,信號矢量畸變較嚴(yán)重;而當(dāng)d/λ大于1時,信號矢量模值趨近于1,互耦影響很小．

圖5 互耦效應(yīng)對陣列信號矢量模值影響

3.2 互耦對空域相關(guān)性影響

互耦效應(yīng)使負(fù)載電阻所接收到的電壓發(fā)生改變,進(jìn)而影響接收信號功率及空域相關(guān)性．考慮互耦效應(yīng)情況下,任意天線接收信號功率為

(25)

根據(jù)式(23)可知

(26)

(27)

式中,ρn,m表示無互耦效應(yīng)下任意兩個天線的相關(guān)性．進(jìn)一步,根據(jù)式(3)及式(27),最終可推導(dǎo)獲得考慮互耦效應(yīng)時,任意天線陣列的天線相關(guān)性為

(28)

式中:

(29)

圖6 互耦對任意兩天線相關(guān)性影響

考慮三維傳播場景下,方位角和俯仰角服從vMF分布,根據(jù)文獻(xiàn)[5]對符合該分布的角度譜測量,并結(jié)合3GPP標(biāo)準(zhǔn)中的三種場景基站端的AS值(2°,2°,5°),設(shè)定兩組k,φ0,θ0分別為(17.71,331.89°,91.414°),(3.952 7,232.48°,84.479°)．圖6仿真比較了耦合與無耦合情況下,任意陣列任意兩天線相關(guān)性近似表達(dá)式隨歸一化距離d/λ∈[0,5]變化情況．從圖中可以看出:不考慮互耦情況時,k值越大,天線相關(guān)性越大,原因在于此時信號角度分布比較集中,對應(yīng)AS較小;考慮互耦情況時,兩天線相關(guān)性曲線在不計互耦效應(yīng)曲線的上下振蕩,k值越小,互耦效應(yīng)對天線相關(guān)性影響越大．

4 數(shù)值仿真及應(yīng)用

本節(jié)將針對3GPP中的市區(qū)微蜂窩傳播環(huán)境,利用本文近似算法,計算均勻線陣(UniformLinearArray,ULA)、均勻圓陣(UniformCircularArray,UCA)、均勻方陣(UniformRectangularArray,URA)三種陣列耦合和無耦合情況下的相關(guān)系數(shù),并進(jìn)一步分析研究互耦對三種天線陣列在小角度擴(kuò)展情況下的影響情況．

設(shè)定天線陣列為N=8的半波偶極子,ULA和URA相鄰天線間距d/λ為0.3,UCA的圓半徑R/λ為0.3,結(jié)合文獻(xiàn)[5]實測數(shù)據(jù),令k,φ0,θ0分別為3.952 7,232.48°, 84.479°．

(a) ULA無互耦情況 (b) ULA互耦情況

(c) UCA無互耦情況 (d) UCA互耦情況

(e) URA無互耦情況 (f) URA互耦情況圖7 互耦對不同陣列流型的相關(guān)性影響

圖7分別計算給出了是否考慮互耦兩種情況下的三種不同天線陣列的相關(guān)系數(shù)模值,其中,顏色越深表示相關(guān)性越小,反之，顏色越淺相關(guān)性越大．從圖中可以看出:對于不同的陣列,互耦效應(yīng)均在一定程度上降低了不同陣列天線的空間相關(guān)性;互耦效應(yīng)對面陣列的相關(guān)性影響最?。?/p>

5 結(jié) 論

本文考慮三維傳播場景并假設(shè)達(dá)波信號角度譜服從vMF分布,詳細(xì)推導(dǎo)給出不同流型陣列的天線空域相關(guān)系數(shù)通用解析式和近似表達(dá)式．在此基礎(chǔ)上,通過研究分析互耦效應(yīng)對陣列接收信號矢量的影響,推導(dǎo)獲得了綜合互耦因素的天線空域相關(guān)系數(shù)的閉合表達(dá)式,該結(jié)論適用于線型、圓型和面型等各種不同的天線陣列,后續(xù)筆者也將應(yīng)用上述結(jié)論研究天線陣列流型及陣元位置優(yōu)化設(shè)計．

[1] LARSSON E G, EDFORS O, TUFVESSON F, et al. Massive MIMO for next generation wireless systems[J]. IEEE communications magazine, 2014, 52(2): 186-195.

[2] GHOSH C K. A compact 4-channel microstrip MIMO antenna with reduced mutual coupling[J]. AEU-international journal of electronics and communications, 2016, 70(7): 873-879.

[3] AMIN M B, ZIRWAS W, HAARDT M. Advanced channel prediction concepts for 5G radio systems[C]//IEEE International Symposium on Wireless Communication Systems (ISWCS), 2015: 166-170.

[4] TSAI J A, BUEHRER R M, WOERNER B D. BER performance of a uniform circular array versus a uniform linear array in a mobile radio environment[J]. IEEE transactions on wireless communications, 2004, 3(3): 695-700.

[5] MAMMASIS K, STEWART R W, THOMPSON J S. Spatial fading correlation model using mixtures of Von Mises Fisher distributions[J]. IEEE transactions on wireless communications, 2009, 8(4): 2046-2055.

[6] YUAN Y, WANG C X, CHENG X, et al. Novel 3D geometry-based stochastic models for non-isotropic MIMO vehicle-to-vehicle channels[J]. IEEE transactions on wireless communications, 2014, 13(1):298-309.

[7] ZHOU J, CAO Z, KIKUCHI H. Analysis of MIMO antenna array based on 3D Von Mises Fisher distribution[J]. Journal of China Universities of Posts & Telecommunications, 2015, 22(2):15-23.

[8] CLERCKX B, CRAEYE C, VANHOENACKER-JANVIER D, et al. Impact of antenna coupling on 2×2 MIMO communications[J]. IEEE transactions on vehicular technology, 2007, 56(3): 1009-1018.

[9] CHEN K H, KIANG J F. Effect of mutual coupling on the channel capacity of MIMO systems[J]. IEEE transactions on vehicular technology, 2015.

[10] 周杰, 陳靖峰, 邱琳,等. 三維空間MIMO信道接收天線陣列互耦效應(yīng)及系統(tǒng)容量分析[J]. 通信學(xué)報, 2012, 33(6):1-10.

ZHOU J, CHEN J F, QIU L, et al. Effect of mutual coupling and antenna correlation on MIMO system in three-dimensional spatial channel models[J]. Journal on communications, 2012, 33(6):1-10.(in Chinese)

[11]UMTS. ETSI TR 125.996 V12.0.0 spatial channel model for multiple input multiple output[S]. European Telecommunications Standards Institute, 2014.

[12]KY?STI P,MEINILJ,HENTILL,et al.IST-4-027756 WINNER II D1.1.2 v1.2 WINNER II channel models[S/OL].[2016-04-06].https://www.researchgate.net/publication/259900906.

[13]CHENG X, YU B, YANG Y, ZHANG J H. Communicating in the real word: 3D MIMO[J]. IEEE wireless communication, 2014, 21(4):136-144.

[14]ZHANG J H, PAN C, PEI F, et al. Three-dimensional fading channel models: A survey of elevation angle research[J]. IEEE communications magazine, 2014, 52(6): 218-226.

[15]LUO Q, PEI F, ZHANG J, et al. 3D MIMO channel model based on field measurement campaign for UMa scenario[C]//IEEE Wireless Communications and Networking Conference, 2014:171-176.

[16]KURZ G, GILITSCHENSKI I, HANEBECK U D. Unscented von Mises-Fisher filtering[J]. IEEE signal processing letters, 2016, 23(4): 463-467.

朱秋明 (1979—),男,江蘇人,博士,南京航空航天大學(xué)副教授,現(xiàn)為英國赫瑞瓦特大學(xué)訪問學(xué)者,研究領(lǐng)域為無線信道勘測、建模和模擬以及航空數(shù)據(jù)鏈技術(shù)等．

楊穎 (1994—),女,江蘇人,碩士研究生,研究方向為無線信道建模及天線特性影響．

On spatial correlation and effect of mutual coupling for multiple antennas under 3D environment

XUE Cuiwei1ZHU Qiuming1,2,3YANG Ying1CHEN Xiaomin1LIAO Zhizhong2

(1.JiangsuKeyLaboratoryofInternetofThingsandControlTechnologies,NanjingUniversityofAeronauticsandAstronautics,Nanjing210016,China; 2.ChinaAirborneMissileAcademy,Luoyang471009,China；3.SchoolofEngineering&PhysicalSciences,Heriot-WattUniversity,EdinburghEH144AS,UK)

For evaluating the spatial correlation(SC) of multiple antennas under three-dimensional(3D) environment, the power angle spectrum of arrival signal is modeled as von Mises-Fisher (vMF) distributions on the sphere coordinate. On this basis, a closed-form expression and an approximate expression of SC for different multiple antenna arrays are derived. The mutual coupling (MC) and its impact on the received signal vector are analyzed, and the general expressions of SC with MC is derived. Simulation results validate that the proposed analytical expressions consistent with the numerical ones. In addition, spatial correlation coefficient with MC waves around the one without MC when the angle spread is small and the uniform rectangular array has the least effect caused by MC.

array manifold; spatial correlation; mutual coupling; von Mises-Fisher distribution

2016-04-06

國家自然科學(xué)基金(61571225)；江蘇省博士后基金(1601017C)；中央高校基本科研業(yè)務(wù)費(fèi)青年科技創(chuàng)新基金

10.13443/j.cjors.2016040601

TN98

A

1005-0388(2017)01-0050-08

薛翠薇 (1991—),女,河北人,碩士研究生,研究方向為無線信道建模及天線特性影響．

(NS2015046, NS2016044)；江蘇省物聯(lián)網(wǎng)與控制技術(shù)重點實驗室基金(NJ20160027)

聯(lián)系人： 朱秋明E-mail：zhuqiuming@nuaa.edu.cn

薛翠薇,朱秋明,楊穎, 等. 三維空間多天線相關(guān)性及互耦影響研究[J]. 電波科學(xué)學(xué)報,2017,32(1):50-57.

XUE C W, ZHU Q M, YANG Y, et al. On spatial correlation and effect of mutual coupling for multiple antennas under 3D environment[J]. Chinese journal of radio science,2017,32(1):50-57. (in Chinese). DOI: 10.13443/j.cjors.2016040601