宋 宇，齊望東，趙衛(wèi)偉

(1. 陸軍工程大學(xué) 指揮控制工程學(xué)院，江蘇 南京 210007; 2. 國(guó)防科技大學(xué) 試驗(yàn)訓(xùn)練基地，陜西 西安 710106)

0 引言

在無(wú)線自組織網(wǎng)絡(luò)(Wireless Ad Hoc Networks,下簡(jiǎn)稱網(wǎng)絡(luò))中，當(dāng)一個(gè)節(jié)點(diǎn)所發(fā)送數(shù)據(jù)幀的目的節(jié)點(diǎn)超出了其傳輸范圍(Transmission Ranges，TR)時(shí)，它需要經(jīng)過(guò)中繼節(jié)點(diǎn)轉(zhuǎn)發(fā)。近年來(lái)，隨著模擬域、數(shù)字域自干擾(Self-Interference，SI)消除技術(shù)的進(jìn)步，節(jié)點(diǎn)可以實(shí)現(xiàn)同時(shí)同頻全雙工(Full-Duplex，F(xiàn)D)通信[1]。當(dāng)中繼節(jié)點(diǎn)以該模式工作時(shí)，其兩側(cè)的鏈路可以同時(shí)工作。但是，中繼節(jié)點(diǎn)以FD模式工作也帶來(lái)了新的干擾樣式。為此，有必要在媒體訪問(wèn)控制層(Media Access Control，MAC)設(shè)計(jì)新的機(jī)制。

基于載波偵聽(tīng)多路訪問(wèn)(Carrier Sensing Multiple Access，CSMA)機(jī)制[2]廣泛應(yīng)用于半雙工(Half-Duplex，HD)網(wǎng)絡(luò)的MAC層[3-4]，在發(fā)送端發(fā)送信號(hào)期間，位于其載波偵聽(tīng)范圍(Carrier Sensing Range，CSR)內(nèi)的節(jié)點(diǎn)都會(huì)偵聽(tīng)到信道忙，不會(huì)發(fā)送信號(hào)[5]。為了避免隱蔽終端問(wèn)題，典型的方法是要求位于接收端干擾范圍(Interference Range，IR)內(nèi)的節(jié)點(diǎn)在它接收幀期間始終能偵聽(tīng)到發(fā)送端的發(fā)送信號(hào)。該方法同樣適用于FD網(wǎng)絡(luò)。但是需要解決以下幾個(gè)新問(wèn)題。

發(fā)送端的CSR大小由其發(fā)送功率決定。假設(shè)鏈路速率固定，網(wǎng)絡(luò)容量大小取決于每條鏈路的有效吞吐率和網(wǎng)絡(luò)中可同時(shí)工作的鏈路數(shù)量。以節(jié)點(diǎn)A、B、C建立的兩跳路徑為例。為了增加網(wǎng)絡(luò)的空間復(fù)用，B不管是單純地向A發(fā)送反向信號(hào)(B未產(chǎn)生有效吞吐)，還是作為A到C的幀傳輸中繼(后續(xù)為表達(dá)方便，將這兩種模式分別簡(jiǎn)稱為FD模式1和FD模式2)，F(xiàn)D MAC機(jī)制都需要盡可能地減小各發(fā)送端的發(fā)送功率以減小它們的CSR。當(dāng)B以FD模式1工作時(shí)，A只能將Data幀發(fā)送至B，B需要再次發(fā)起到C的Data幀傳輸。當(dāng)B以FD模式2工作時(shí)，A可以直接將Data幀發(fā)送至C。如果不考慮控制幀開(kāi)銷，后者吞吐率是前者的兩倍。但是，當(dāng)B以FD模式1工作時(shí)隱蔽終端只在其周圍存在。當(dāng)B以FD模式2工作時(shí)隱蔽終端同時(shí)在B和C周圍存在，為了避免傳輸干擾，可能需要兩個(gè)發(fā)送端使用較大發(fā)送功率。

為此，F(xiàn)D MAC機(jī)制需要在吞吐增益和占用空間大小問(wèn)題上權(quán)衡B的模式選擇。具體地，分別建模和求解在FD模式1和2下兩個(gè)發(fā)送端發(fā)送功率的優(yōu)化模型。隨后，通過(guò)對(duì)比在兩種不同模式下兩跳路徑的吞吐率和占用空間大小以決定它們采用的通信模式。

1 系統(tǒng)模型

以節(jié)點(diǎn)A、B、C建立的兩跳路徑為例。將各發(fā)送端最大發(fā)送功率標(biāo)記為Ptmax。將A的發(fā)送功率標(biāo)記為PtA。DataA→B和TDataA→B分別標(biāo)記由A發(fā)送至B的數(shù)據(jù)幀以及其傳輸時(shí)間。B接收來(lái)自A的信號(hào)強(qiáng)度以及信噪比分別以PrA→B和SINRrA→B標(biāo)記。類似地，有PtB，DataB→C，TDataB→C，PrA→C，PrB→C以及SINRrB→C等標(biāo)記。此外，B還可以向A發(fā)送反向信號(hào)，以PrB→A標(biāo)記A接收來(lái)自B的信號(hào)強(qiáng)度。注意該信號(hào)對(duì)于當(dāng)前路徑?jīng)]有產(chǎn)生有效吞吐。PrA→B由式(1)給定：

(1)

其中，DAB標(biāo)記A和B之間的距離，c為常數(shù)。類似地，可以獲得PrA→C和PrB→C的表達(dá)式。SINRrA→B由式(2)給定：

(2)

其中，PnB標(biāo)記B所受到的總干擾功率，有：

(3)

Psthold代表節(jié)點(diǎn)偵聽(tīng)到信號(hào)的最小功率界。SINRthold代表最小節(jié)點(diǎn)成功接收信號(hào)的最小SINR界。將CSRA(PtA)標(biāo)記為A的CSR。位于CSRA(PtA)內(nèi)節(jié)點(diǎn)的接收信號(hào)需要滿足Psthold約束。由CSRA(PtA)所圍成的區(qū)域被認(rèn)為是A發(fā)送信號(hào)占用的空間大小，將其大小標(biāo)記為S(CSRA)。類似地，有CSRB(PtB)以及S(CSRB)。將CSRAB(PtA,PtB)標(biāo)記為A和B同時(shí)發(fā)送信號(hào)所形成的CSR，將其大小標(biāo)記為S(CSRAB)。IRB(FD)代表B以FD模式工作時(shí)的IR。由式(1)～(3)可得：

(4)

其中，c′由B以及與它最近且沒(méi)有干擾到B接收的節(jié)點(diǎn)天線的高度和增益決定。盡管C是以HD模式接收幀，但是，它會(huì)受到來(lái)自A的發(fā)送信號(hào)干擾。以IRC(HD)標(biāo)記C的IR，簡(jiǎn)寫(xiě)為IRC(HD)，其表達(dá)式如下：

(5)

為簡(jiǎn)化討論，假設(shè)各發(fā)送端天線高度和增益相同，c′和c均為1。同時(shí)令Pn=0。

2 針對(duì)FD中繼的MAC層高效空間復(fù)用機(jī)制

首先，分別針對(duì)B的兩種FD模式，計(jì)算各發(fā)送端的最優(yōu)發(fā)送功率。隨后，分別求解在FD模式1和2下兩跳路徑所占用的空間大小。

2.1 針對(duì)B發(fā)送補(bǔ)充信號(hào)時(shí)的功率控制方法

在接收DataA→B幀期間，B可以發(fā)送補(bǔ)充信號(hào)，以避免其接收數(shù)據(jù)幀時(shí)的隱蔽終端問(wèn)題[6]。此時(shí)，A和B之間的通信鏈路占用空間的大小由二者的發(fā)送信號(hào)共同決定。必須聯(lián)合求解PtA和PtB。

在A、B同時(shí)發(fā)送信號(hào)期間，為確保B對(duì)DataA→B幀的接收不被干擾，不允許位于IRB(FD)內(nèi)的節(jié)點(diǎn)發(fā)送信號(hào)。具體地，要求CSRAB(PtA,PtB)覆蓋IRB(FD)。將CSRAB(PtA,PtB)和IRB(FD)在極坐標(biāo)中建模，如圖1所示。將A的位置作為原點(diǎn)(0,0)，向量AB的方向作為θ=0方向，Q是該方向上的一個(gè)位置。將B與Q的距離標(biāo)記為DBQ。令DBQ=IRB(FD)。在某個(gè)θ方向上，假設(shè)M是可以偵聽(tīng)到A和B發(fā)送信號(hào)的最遠(yuǎn)位置。令：

DAM=rmax(θ)

(6)

圖1 極坐標(biāo)系中CSRAB(PtA,PtB)和IRB(FD)的示意圖

由DAM和DAB，根據(jù)余弦定理，可得DBM。進(jìn)一步地，由M位置的定義，根據(jù)Psthold定義，有式(7)成立：

(7)

若給定(PtA,PtB)值，結(jié)合DBM，可以通過(guò)啟發(fā)式方法求解rmax(θ)。該rmax(θ)即為在某個(gè)θ方向上CSRAB(PtA,PtB)的值。

注意到求解PtA和PtB的目標(biāo)是使得S(CSRAB)盡量小并且要求IRB(FD)必須被CSRAB(PtA,PtB)覆蓋。首先，討論P(yáng)tA和PtB取值的可行域。由于IRB(FD)邊界上的各個(gè)位置與B的距離相同而與A的距離不同，所以在這些位置上由B的發(fā)送信號(hào)產(chǎn)生的接收信號(hào)強(qiáng)度相同而由A的發(fā)送信號(hào)產(chǎn)生的接收信號(hào)強(qiáng)度不同。注意到IRB(FD)邊界上的Q位置距離A最遠(yuǎn)，所以A的發(fā)送信號(hào)在該位置產(chǎn)生的接收信號(hào)強(qiáng)度最小。因此，給定(PtA,PtB)，如果有下面的條件成立：

rmax(θ)|θ=0≥DAB+IRB(FD)

(8)

即節(jié)點(diǎn)在Q位置上可以偵聽(tīng)到來(lái)自A和B的發(fā)送信號(hào)，則整個(gè)IRB(FD)可以被CSRAB(PtA,PtB)覆蓋。此時(shí)，M位置也就一定位于IRB(FD)之外。進(jìn)一步地，式(9)給出了PtA和PtB需要滿足的約束：

(9)

滿足式(9)的PtA和PtB取值即為二者的可行解。下面根據(jù)(PtA,PtB)解的可行域中二者具體取值利用積分的思想近似求解S(CSRAB)。在圖1中，以Δθ標(biāo)記θ的增量。假設(shè)N是在θ+Δθ方向上可以偵聽(tīng)到A和B發(fā)送信號(hào)的最遠(yuǎn)位置。令：

DAN=rmax(θ+Δθ)

(10)

給定PtA和PtB取值，同樣采用啟發(fā)式方法求解rmax(θ+Δθ)值。從M和N位置向AB垂線，令交點(diǎn)分別為M′和N′。DMM′、DNN′以及DM′N′的值可以由基本的三角函數(shù)關(guān)系獲得。根據(jù)它們的值，可以得到由M、M′、N以及N′所圍成的梯形區(qū)域面積，標(biāo)記為SMM′N′N。

為了求解S(CSRAB)，可以將[0,π]劃分為若干Δθ。隨后，采用上述方法求解每個(gè)θ變化區(qū)間由A和B發(fā)送信號(hào)強(qiáng)度所圍成的梯形面積。將它們累加起來(lái)乘以2即為S(CSRAB)的近似值。在PtA和PtB取值的可行域中，使得S(CSRAB)取得最小值的PtA和PtB組合即為該模式下的最優(yōu)解。

2.2 針對(duì)B作為FD中繼時(shí)的功率控制方法

在接收DataA→B幀期間，B可以以FD模式將該幀直接中繼給C。在信號(hào)層面上，該模式下也是A和B同時(shí)發(fā)送信號(hào)。此時(shí)PtA和PtB除了需要確保B對(duì)DataA→B幀的接收不被干擾外，還需要避免當(dāng)C接收信號(hào)時(shí)其周圍可能存在的隱蔽終端問(wèn)題。也就是說(shuō)，在A和B同時(shí)發(fā)送幀期間，要求CSRAB(PtA,PtB)同時(shí)覆蓋IRB(FD)和IRC(HD)。前一個(gè)覆蓋條件與上一節(jié)討論一致，下面主要討論后一個(gè)覆蓋條件。

當(dāng)B作為A向C發(fā)送數(shù)據(jù)幀的中繼節(jié)點(diǎn)時(shí)，盡管C是以HD模式接收來(lái)自B的數(shù)據(jù)幀，但是由式(5)，IRC(HD)與PtA、PtB、DBC以及DAC均相關(guān)。圖2采用的極坐標(biāo)系與圖1一致。分別延著向量AC和BC方向作延長(zhǎng)線，將二者與IRC(HD)邊界的交點(diǎn)標(biāo)記為P和Q。此外，IRC(HD)還分別與AC和BC(或者是二者反方向的延長(zhǎng)線)相交。將這兩個(gè)交點(diǎn)標(biāo)記為A′和B′。圖2所示為前一種的情況。下面討論要求IRC(HD)必須被CSRAB(PtA,PtB)覆蓋條件對(duì)PtA和PtB可行域的約束。

圖2 極坐標(biāo)系中CSRAB(PtA,PtB)和IRC(HD)的示意圖

在PB′弧(該弧在IRC(HD)邊界上)的各個(gè)位置上，P位置對(duì)于A和B來(lái)說(shuō)都是最遠(yuǎn)的位置。所以，節(jié)點(diǎn)在該位置接收到的來(lái)自二者的發(fā)送信號(hào)強(qiáng)度最小。同理，在QA′弧的各個(gè)位置上，節(jié)點(diǎn)在Q位置接收到的來(lái)自A和B的發(fā)送信號(hào)強(qiáng)度最小。所以，在數(shù)據(jù)幀發(fā)送期間，為了讓PB′弧和QA′弧上面的節(jié)點(diǎn)不發(fā)送信號(hào)，只需讓處于P和Q位置的節(jié)點(diǎn)偵聽(tīng)到A和B的發(fā)送信號(hào)強(qiáng)度即可。當(dāng)節(jié)點(diǎn)沿著PQ弧(或QP弧)移動(dòng)時(shí)，節(jié)點(diǎn)與其中一個(gè)發(fā)送端的距離增大，而與另一個(gè)發(fā)送端的距離減小。此時(shí)，如果節(jié)點(diǎn)在端點(diǎn)P和Q位置上可以偵聽(tīng)到兩個(gè)發(fā)送端的發(fā)送信號(hào)強(qiáng)度，則近似認(rèn)為節(jié)點(diǎn)在位于PQ弧的各個(gè)位置也能偵聽(tīng)到。所以，將避免C接收幀被干擾的條件建模如下：

(11)

(12)

由余弦定理，可以獲得DAQ和DBP的值。將B作為FD中繼節(jié)點(diǎn)中繼A到B數(shù)據(jù)幀傳輸，要求PtA和PtB必須滿足這兩個(gè)公式和上述公式(9)。在(PtA,PtB)所有可行解中使得S(CSRAB)取得最小值的解即為該問(wèn)題的解。與上述方法一致，不再贅述。

以單位面積吞吐率測(cè)度判定B的FD模式在仿真分析中討論。

3 仿真分析

分別取SIB=1×10-9以及5×10-9。給出B分別以FD模式1和2工作時(shí)S(CSRAB)的數(shù)值仿真及分析。在此基礎(chǔ)上，討論B應(yīng)采用的FD模式。

3.1 B節(jié)點(diǎn)以FD模式1工作時(shí)S(CSRAB)的數(shù)值仿真及分析

該模式下B接收DataA→B幀期間僅發(fā)送反向信號(hào)。隨著DAB和SIB的變化，S(CSRAB)如圖3所示。

圖3 根據(jù)DAB和SIB取值的不同S(CSRAB)的變化示意圖

假設(shè)PtA不變，隨著DAB增大，PrA→B會(huì)減弱。同時(shí)，由于SIB≠0，當(dāng)B發(fā)送補(bǔ)充信號(hào)時(shí)它會(huì)有自干擾。這兩個(gè)因素均會(huì)降低SINRA→B。所以，在SIB兩種取值下，隨著DAB增大，PtA均會(huì)增大。保持DAB不變，當(dāng)SIB由1×10-9增大到5×10-9時(shí)，相同的PtB會(huì)產(chǎn)生更大的自干擾。為確保DataA→B幀的正確接收，PtA和PtB分別需要增大和減小。對(duì)應(yīng)地，S(CSRAB)有略微增大，如圖3所示。由于在該模式下其他節(jié)點(diǎn)并不需要接收B發(fā)送的信號(hào)，對(duì)PtB的約束只是讓S(CSRAB)盡量小，所以，當(dāng)PtA隨著DAB增大時(shí)，有可能與該P(yáng)tA只需要搭配一個(gè)較小的PtB值即可滿足約束(9)。仿真中發(fā)現(xiàn)，隨著DAB變化，PtB的變化并無(wú)明確的單調(diào)性。需要說(shuō)明的是，當(dāng)DAB=100 m時(shí)，PtA和PtB在SIB=1×10-9情況下有解(盡管此時(shí)PtB值很小)，而在SIB=5×10-9情況下無(wú)解。這是因?yàn)樵谠揇AB和SIB取值下，即使是很小的PtB也需要足夠大的PtA以確保DataA→B幀的正確接收，但是所要求的PtA值大于Ptmax。

3.2 B節(jié)點(diǎn)以FD模式2工作時(shí)S(CSRAB)的數(shù)值仿真及分析

該模式下B接收DataA→B幀的同時(shí)需要發(fā)送DataB→C幀。為簡(jiǎn)化討論，僅選取DAB=40 m和60 m兩種情況。假設(shè)圖2中φ=5π/6，隨著DAB、DBC、SIB的變化，S(CSRAB)分別如圖4和5所示(針對(duì)SIB取值的兩種情況分別展示)。

圖4 當(dāng)SIB=1×10-9時(shí)根據(jù)DAB和DBC取值的不同S(CSRAB)的變化

圖5 當(dāng)SIB=5×10-9時(shí)根據(jù)DAB和DBC取值的不同S(CSRAB)的變化

不論SIB=1×10-9還是5×10-9，若保持DAB和PtB不變，隨著DBC增大，PrB→C會(huì)減弱。為確保DataB→C幀的正確接收，PtB會(huì)增大。類似地，保持DBC不變，隨著DAB增大，為確保DataA→B幀的正確接收，PtA會(huì)增大。相應(yīng)地，隨著DAB和DBC增大，S(CSRAB)會(huì)增大，如圖4和5所示。仿真中發(fā)現(xiàn)，當(dāng)DBC取較大值時(shí)(圖4中DBC取80 m或90 m，圖5中DBC取60 m或70 m), 隨著DAB增大，除了PtA增大外，PtB也會(huì)增大。這是因?yàn)镻tA增大會(huì)擴(kuò)大IRC(HD)，盡管PtA增大使得CSRAB(PtA,PtB)擴(kuò)大，但是由于A距離C較遠(yuǎn)，該擴(kuò)大的CSRAB(PtA,PtB)不足以覆蓋擴(kuò)大后的IRC(HD)，所以需要適當(dāng)增大PtB。而當(dāng)SIB=5×10-9時(shí)，PtB增大會(huì)產(chǎn)生較大的自干擾，為確保DataA→B幀的正確接收，PtA需要進(jìn)一步增大。所以，當(dāng)DAB增大到一定程度時(shí)，所要求的PtB大于Ptmax。比如，當(dāng)DAB=60 m，DBC=60 m或90 m時(shí)，(PtA,PtB)無(wú)解。

3.3 節(jié)點(diǎn)之間通信模式的選擇

當(dāng)SIB=1×10-9時(shí)，對(duì)應(yīng)DAB=40 m和60 m，當(dāng)B以FD模式1工作時(shí)，S(CSRAB)分別為7.54×104m2和11.76×104m2。當(dāng)B以FD模式2工作時(shí)，對(duì)應(yīng)DAB=40 m，S(CSRAB)最小為9.05×104m2(當(dāng)DBC=40 m時(shí))，最大為20.48×104m2(當(dāng)DBC=90 m時(shí))。當(dāng)S(CSRAB)大于兩倍的7.54×104m2時(shí)(仿真中DBC>70 m)， B不適宜以FD模式2工作。否則，可以以該模式工作以增加單位面積的吞吐率。同樣，對(duì)應(yīng)DAB=60 m，S(CSRAB)最小為12.09×104m2(當(dāng)DBC=40 m時(shí))，最大為22.3×104m2(當(dāng)DBC=90 m時(shí))。與前面不同的是，該范圍內(nèi)所有S(CSRAB)均無(wú)法大于11.76×104m2的兩倍。此時(shí)，B不宜以FD模式2工作。當(dāng)SIB=5×10-9時(shí)，有類似討論，不再贅述。

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

針對(duì)無(wú)線網(wǎng)絡(luò)中的兩跳路徑，可以采用FD節(jié)點(diǎn)作為中繼節(jié)點(diǎn)。針對(duì)該節(jié)點(diǎn)以FD工作的兩種模式，以最小化兩個(gè)發(fā)送端的發(fā)送信號(hào)占用的空間大小為目標(biāo)，分別求解了它們的最小發(fā)送功率。并提出了單位面積吞吐率測(cè)度以讓FD中繼節(jié)點(diǎn)選擇所采用的FD模式。給出的數(shù)值仿真驗(yàn)證了該機(jī)制的有效性。下一步，需要設(shè)計(jì)相關(guān)的控制幀以及協(xié)議交互過(guò)程，以將該機(jī)制融入具體的FD MAC協(xié)議。