李國權(quán)，林金朝，徐勇軍，黃正文，劉挺

（1.重慶郵電大學通信與信息工程學院，重慶 400065；2.光電信息感測與傳輸技術(shù)重慶市重點實驗室，重慶 400065；3.西安郵電大學陜西省信息通信網(wǎng)絡及安全重點實驗室，陜西 西安 710121；4.布魯內(nèi)爾大學電子與計算機工程系，倫敦 UB8 3PH）

1 引言

面對重大自然災害和緊急突發(fā)事件，有效的應急通信對于提高救援效率和保障人民生命安全具有重要的意義。當發(fā)生事故或災難時，當?shù)氐墓潭ɑ就ǔ２荒苷Ｊ褂?，而下一代通信系統(tǒng)中的無人機（UAV,unmanned aerial vehicle）具有更好的靈活性[1]，特別是對于偏遠地區(qū)或缺乏基站設(shè)施的地區(qū)，應快速構(gòu)建無人機空中基站。因此，無人機為有序應對突發(fā)事件，盡可能地降低其危害程度做出了重要的貢獻[2]。針對應急通信場景下頻譜資源稀缺及基站毀壞的問題，提升無人機系統(tǒng)的頻譜效率具有重要意義。

非正交多址接入（NOMA,non-orthogonal multiple access）技術(shù)具有較高的頻譜效率與較好的公平性，被視為5G乃至下一代無線通信系統(tǒng)的關(guān)鍵技術(shù)[3-5]。在功率域NOMA系統(tǒng)中，信道條件好的用戶被分配較低的功率分配因子，而信道條件差的用戶被分配較高的功率分配因子。在接收端，利用串行干擾消除（SIC,successive interference cancellation）技術(shù)來消除部分用戶的干擾，實現(xiàn)更多用戶在同一時頻域復用，從而提升頻譜效率[6]。

將NOMA引入無人機通信系統(tǒng)，能夠充分利用功率域資源，從而保證空中基站為地面用戶更好地傳輸信號（下行鏈路）或者保證空中無人機為地面用戶更好地傳輸信號（上行鏈路）。無人機系統(tǒng)主要受到視線（LoS,line of sight）鏈路的影響，所以傳統(tǒng)的瑞利衰落不適合表示其信道特性[7]。

關(guān)于UAV-NOMA系統(tǒng)的研究，按照信道特性分類，主要有空對地（A2G,air to ground）信道、Nakagami-m衰落信道、路徑損耗信道和萊斯信道。針對考慮A2G信道的UAV-NOMA系統(tǒng)，文獻[8]研究了不同城市環(huán)境下的和速率最大化問題，同時對比了固定與動態(tài)UAV高度的影響，從而減少能耗。在相同的模型下，文獻[9]考慮了多用戶服務質(zhì)量約束，將能效最大化問題等效為非線性分式規(guī)劃問題，其中基于信道條件的用戶分組方案被考慮。針對考慮Nakagami-m衰落信道的UAV-NOMA系統(tǒng)，文獻[10]提出了以UAV為中心的卸載操作策略和以用戶為中心的應急通信策略，分別適用于密集網(wǎng)絡與需要同時服務所有用戶的場景，從而提升系統(tǒng)覆蓋概率。文獻[11]則考慮了LoS鏈路與非視線（NLoS,non-line of sight）鏈路的影響，利用隨機幾何模型來模擬用戶與UAV的位置，并且推導了系統(tǒng)中斷概率與遍歷速率的閉式表達式。針對LoS鏈路與NLoS鏈路場景，文獻[12]首先根據(jù)訪問優(yōu)先級確定用戶分組方案，再利用消息傳遞算法進行子信道分配，最后聯(lián)合優(yōu)化UAV-NOMA系統(tǒng)的發(fā)射功率。文獻[13]則研究了相同模型下的多個旋轉(zhuǎn)無人機輔助的上行NOMA系統(tǒng)，并提出了基于K均值的UAV部署方案和基于坐標的用戶分組方案。同時，考慮了最小服務質(zhì)量約束，聯(lián)合優(yōu)化基于總功率最小化的資源分配問題。結(jié)果表明，UAV-NOMA系統(tǒng)的頻譜效率和能量效率得到有效提升[12-13]。針對考慮路徑損耗模型的UAV-NOMA系統(tǒng)，文獻[14]考慮了最小可達速率需求，研究了聯(lián)合UAV位置與功率分配優(yōu)化問題來減小發(fā)射功率。文獻[15]研究了聯(lián)合UAV位置、接納控制與功率分配優(yōu)化問題，利用懲罰函數(shù)法與連續(xù)凸逼近法來求解卡羅需-庫恩-塔克（KKT,Karush-Kuhn-Tucher）收斂點，從而提高UAV-NOMA系統(tǒng)的用戶接入數(shù)。此外，文獻[16]將無人機引入車聯(lián)網(wǎng)，利用對偶分解法，通過聯(lián)合優(yōu)化鏈路切換次數(shù)與鏈路傳輸速率來提高UAV-NOMA系統(tǒng)吞吐量。針對考慮萊斯信道的UAV-NOMA系統(tǒng)，文獻[17]考慮了上行干擾消除約束，研究了聯(lián)合UAV位置、地面基站位置、預編碼優(yōu)化問題，從而提升地面用戶的和速率。在相同的信道模型下，文獻[18]研究了上行UAV-NOMA系統(tǒng)的預編碼優(yōu)化問題，考慮了SIC約束和速率需求約束，利用一階泰勒展開將上述優(yōu)化問題轉(zhuǎn)化為二階圓錐規(guī)劃問題，從而提升系統(tǒng)和速率。文獻[19]研究了總發(fā)射功率約束、總帶寬約束、UAV高度約束、天線波束寬度約束條件下的最大最小速率優(yōu)化問題，利用路徑追蹤算法求解，得到有效的UAV-NOMA系統(tǒng)和速率提升。此外，文獻[20]研究了UAV旋轉(zhuǎn)軌跡的優(yōu)化問題，從而降低了下行UAV-NOMA系統(tǒng)的中斷概率。

上述關(guān)于UAV-NOMA系統(tǒng)的工作大多只考慮UAV高度、位置、軌跡的影響，而缺乏對地面用戶分組與功率分配聯(lián)合優(yōu)化問題的研究。因此，本文主要研究UAV-NOMA系統(tǒng)中的動態(tài)用戶分組與功率分配算法，利用圖論中的最大割定理與輔助變量法，使UAV-NOMA系統(tǒng)的和速率性能得到較好的提升。本文主要貢獻如下。

1)在下行單小區(qū)多組場景下，建立了一個基于無人機輔助的NOMA通信網(wǎng)絡模型，其中一個無人機作為空中基站，服務多個地面用戶，而這些地面用戶被均分到多個組。傳統(tǒng)的UAV-NOMA系統(tǒng)通常假設(shè)每組用戶數(shù)為2，即每組只包含一個近用戶和一個遠用戶。為了更符合實際，本文所提的UAV-NOMA系統(tǒng)將考慮組內(nèi)用戶數(shù)更多的情形。

2)為了進一步提升UAV-NOMA系統(tǒng)的頻譜效率，考慮了基站最大發(fā)射功率約束與用戶分組約束，構(gòu)建了基于多用戶和速率最大化的用戶分組與功率分配聯(lián)合優(yōu)化問題。該問題是混合整數(shù)非線性規(guī)劃問題，故難以求其閉式解。

3)為了解決上述問題，提出了兩步策略。首先，基于圖論中的最大割定理，進行用戶分組，保證用戶與無人機的相對距離較小。然后，對每組用戶分別執(zhí)行基于多用戶和速率最大化的功率分配子問題，通過輔助變量法將其轉(zhuǎn)換為凸優(yōu)化問題，從而得到每組用戶的功率分配因子。

4)仿真結(jié)果表明，本文算法在不同的覆蓋半徑、無人機高度、發(fā)射功率和信噪比下具有較好的和速率性能。

2 系統(tǒng)模型及問題描述

考慮一個下行UAV-NOMA系統(tǒng)，如圖1所示。該系統(tǒng)包括一個含Nt根發(fā)射天線的UAV作為空中基站，工作在固定高度，且地面覆蓋范圍是半徑為R的圓。此外，該系統(tǒng)具有N個單天線地面用戶。假設(shè)N個用戶被分為T組，每組有K個用戶，并記第t組的第k個用戶為，定義組號集合為t∈{1,2,…,T}，用戶序號集合為n∈{1,2,…,N}，組內(nèi)用戶序號集合為k∈{1,2,…,K}。所有組的用戶1均勻分布在半徑為R1的小圓內(nèi)，而所有組的用戶i均勻分布在半徑為Ri-1和Ri的同心圓環(huán)內(nèi)。其中，Ri-1＜Ri，且RK=R，同時定義用戶覆蓋半徑集合為。此外，組內(nèi)用戶利用NOMA技術(shù)共享相同的時頻域資源，而各組之間保持正交多址接入（OMA,orthogonal multiple access），即忽略組間干擾[21]。

圖1 基于無人機的非正交多址接入系統(tǒng)模型

為了提升UAV-NOMA系統(tǒng)頻譜效率和傳輸質(zhì)量，考慮基于系統(tǒng)和速率最大的用戶分組與功率分配聯(lián)合優(yōu)化問題，該優(yōu)化問題可以表示為

其中，C1表示基站最大功率約束，C2表示用戶功率分子因子約束，C3表示一個用戶至多分到一個組，C4表示每組用戶數(shù)約束，=1表示用戶i屬于第t組，表示用戶i不屬于第t組[22]。

顯然，上述問題是混合整數(shù)非線性規(guī)劃問題，故難以求解。為了使和速率達到最大化，需要窮舉搜索才能得到全局最優(yōu)的用戶分組方案，而搜索次數(shù)高達次[23]。當用戶數(shù)很多時，該方式不合理。

為此，提出一種兩步策略來減少復雜度。首先，執(zhí)行用戶分組算法。然后，對每組用戶執(zhí)行功率分配算法來得到次優(yōu)解。

3 基于最大割定理的用戶分組算法

本節(jié)將基于圖論中的最大割定理提出低復雜度的用戶分組算法。傳統(tǒng)的用戶分組方案通常只考慮用戶的信道條件，將信道條件最好的分到一組[24]，或?qū)⑿诺啦町愖畲蟮姆值揭唤M[22]。但是，上述方案都只考慮瑞利信道的影響。在復合信道模型下，UAV-NOMA系統(tǒng)會受到Nakagami-m衰落信道的影響。因此，傳統(tǒng)的方案也許不適于無人機場景。

此外，資源分配可等價為圖論問題。文獻[25]將車聯(lián)網(wǎng)領(lǐng)域中的配對問題轉(zhuǎn)化成圖論中的最大割問題，將用戶間干擾最小的用戶分為一組，從而保證車輛連接的可靠性。文獻[26]將用戶分組問題轉(zhuǎn)換成圖論中的尋找負價環(huán)問題來滿足最小功率需求。然而，在UAV-NOMA系統(tǒng)中，用戶與UAV的距離會影響系統(tǒng)性能。雖然假設(shè)用戶與UAV的位置固定不變，但是可通過圖論中的最大割定理來減少組內(nèi)用戶到UAV的相對距離，從而改善組內(nèi)用戶的信道狀況。

假設(shè)分組前N個用戶具有個初始組，并且每個初始組包含個用戶。此外，需要滿足式(7)和式(8)中的約束條件C3和C4。UAV-NOMA系統(tǒng)中用戶與UAV的相對距離關(guān)系可表示為無向圖，如圖2所示。其中，一個節(jié)點表示一個用戶，每2個節(jié)點間由一條邊相連。邊的權(quán)值為即UAV與用戶a和用戶b的距離之和。在滿足分組約束的條件下，保證組內(nèi)用戶與UAV的總距離最近。這是因為用戶距離UAV越近，則信道狀況越好。

圖2 用戶與UAV的相對距離關(guān)系

圖2中優(yōu)化用戶與UAV的相對距離等價于圖論中的最大割問題。已知圖G=(V,E)，其中，V表示地面所有用戶的集合，Vi表示第i個初始組的用戶集合；E表示邊集合，表示每次可被選擇的邊集合，表示K的補集，并且當i=1時，否則，

對于無向圖，分組后所有用戶被分到T個組，因此顯然，而分組后每組用戶數(shù)。此外，根據(jù)所有已選邊與未選邊之和為邊集E，可以得到

下面，證明利用K維相對距離來實現(xiàn)最大割的可行性，有命題1成立。

命題1如果需要完成K個用戶分組，并且所有用戶則可以通過計算ωK，而不是判斷來自所有組合邊的權(quán)值之和來確保K個用戶與UAV的相對距離最小。

證明對于任意的組內(nèi)用戶數(shù)K與組數(shù)T，已知每個用戶均來自用戶集合V，則有

證畢。

本文所提基于最大割定理的用戶分組算法如算法1所示。

4 基于和速率最大化的功率分配算法

本節(jié)研究的功率分配算法是基于第2節(jié)中定義的信道增益排序假設(shè)。不失一般性，定義第3節(jié)的分組結(jié)果為t*。基于該結(jié)果，式(4)問題可以等價為如下功率分配子問題

其中，約束條件為每組用戶的總功率分配因子約束。

雖然該約束條件是凸條件，但是目標函數(shù)分子分母存在耦合項，不易求解。利用輔助變量法[27]定義非負輔助變量，可以將目標函數(shù)縮放成

其中，C5與式(11)相同，C6為輔助變量上界。

由于目標函數(shù)和約束條件是正項式，因此式(13)是幾何規(guī)劃問題。因此，通過改變變量，分別可轉(zhuǎn)化為。約束條件C2可轉(zhuǎn)換為

其中，約束條件C7為凸條件。至此，可以求得次優(yōu)功率分配因子[28]。

5 仿真結(jié)果

本節(jié)將通過仿真驗證本文算法的有效性。當Nakagami-m衰落信道參數(shù)m＞1時，表示LoS鏈路。不失一般性，假設(shè)m=2，路徑損耗參數(shù)α=3，且每組總發(fā)射功率P=1W。此外，假設(shè)UAV的發(fā)射天線數(shù)為4。為了方便地分析問題，定義如下幾類算法：本文算法，即聯(lián)合本文所提的用戶分組算法和功率分配算法；算法1，即聯(lián)合最強近用戶與最強遠用戶（BNBF,best-near best-far）分組算法[24]和固定功率分配算法[21]；算法2，即聯(lián)合最強近用戶與最弱遠用戶（BNWF,best-near worst-far）分組算法[22]和固定功率分配算法[21]，其中，地面用戶的功率分配因子固定為分別對應兩用戶、三用戶、四用戶和五用戶的情形；算法3，即聯(lián)合BNBF用戶分組算法[24]與基于最小速率需求的功率分配算法[29]。

圖3對比了不同覆蓋半徑下兩用戶正交與非正交UAV系統(tǒng)和速率性能。其中，非正交UAV系統(tǒng)應用本文算法。此外，假設(shè)R1=0.5km、H=0.5km、SNR=5dB。由圖3可知，與OMA方案相比，無人機在NOMA方案下具有更好的和速率性能。顯然，隨著組數(shù)T增大，非正交UAV系統(tǒng)的性能優(yōu)勢更明顯。這表明應用本文算法的非正交UAV系統(tǒng)在用戶數(shù)較多時仍具有較好的頻譜效率。

圖3 不同覆蓋半徑下兩用戶正交與非正交UAV系統(tǒng)和速率性能

圖4給出了不同高度下兩用戶非正交UAV系統(tǒng)和速率性能。其中，R1=0.5 km、R2=1.5 km、T=10、SNR=5 dB。由圖4可知，與其他動態(tài)算法相比，應用本文算法的UAV-NOMA系統(tǒng)具有更好的和速率性能。這表明本文算法的有效性，從而保證UAV-NOMA系統(tǒng)的工作高度更具靈活性。

圖4 不同高度下兩用戶非正交UAV系統(tǒng)和速率性能

圖5對比了不同功率下兩用戶非正交UAV系統(tǒng)和速率性能。其中，假設(shè)R1=0.5 km、R2=3 km、H=1 km、T=10、SNR=5 dB。顯然，隨著總發(fā)射功率的增加，UAV-NOMA系統(tǒng)和速率也增加。其中，本文算法具有最優(yōu)的和速率性能。因此，對比動態(tài)BNBF、BNWF算法、基于最小速率需求的功率分配算法，本文算法使UAV-NOMA系統(tǒng)的頻譜效率得到較好提升。

圖5 不同功率下兩用戶非正交UAV系統(tǒng)和速率性能

圖6給出了不同信噪比下兩用戶非正交UAV系統(tǒng)和速率性能。其中，R1=0.5km、R2=3km、H=1km、T=10、SNR=5dB。由圖6可知，本文算法具有最優(yōu)的和速率性能。對比傳統(tǒng)的動態(tài)用戶分組算法[22,24]與功率分配算法[5,29]，本文算法與基于和速率最大化的功率分配算法更適于Nakagami-m衰落信道下的兩用戶UAV-NOMA系統(tǒng)。

圖6 不同信噪比下兩用戶非正交UAV系統(tǒng)和速率性能

圖7對比了不同信噪比下三用戶非正交UAV系統(tǒng)和速率性能。其中，R1=0.5 km、R2=1 km、R3=3 km、H=1 km、T=10、SNR=5 dB。由圖7可知，本文算法在三用戶場景下也具有相對最優(yōu)的和速率性能。

圖7 不同信噪比下三用戶非正交UAV系統(tǒng)和速率性能

圖8給出了不同信噪比下四用戶和五用戶非正交UAV系統(tǒng)和速率性能比較。其中，=[0.5,1.5,2.5,3]和=[0.5,1,1.5,2,3]分別代表四用戶和五用戶UAV-NOMA系統(tǒng)中各用戶的覆蓋半徑。此外，H=1 km、T=10、SNR=5 dB。顯然，在相同的組內(nèi)用戶數(shù)K下，本文算法具有更優(yōu)的和速率性能。另外，五用戶系統(tǒng)和速率性能比四用戶系統(tǒng)和速率性能更好，這表明本文算法比其他傳統(tǒng)算法更適用于用戶數(shù)較多的情形。

圖8 不同信噪比下四用戶和五用戶非正交UAV系統(tǒng)和速率性能

6 結(jié)束語

本文研究了基于無人機輔助的NOMA網(wǎng)絡的聯(lián)合用戶分組與功率分配算法，從而進一步提升了無人機系統(tǒng)的頻譜利用率。本文算法能夠動態(tài)地進行用戶分組，并且得到每組用戶的功率分配因子。仿真結(jié)果表明，對比正交UAV系統(tǒng)，應用本文算法的非正交UAV系統(tǒng)在不同覆蓋半徑條件下具備更好的系統(tǒng)和速率性能。此外，針對不同UAV高度、不同發(fā)射功率、不同信噪比等情形，非正交UAV系統(tǒng)在應用本文算法時都具備較好的和速率性能。同時，對比隨機用戶分組，動態(tài)BNBF、BNWF用戶分組，以及固定功率分配、最小速率需求功率分配算法，本文所提的基于最大割定理的用戶分組算法與基于和速率最大化的功率分配算法具有更好的穩(wěn)健性，從而有效地提升了Nakagami-m衰落信道下的UAV-NOMA系統(tǒng)的頻譜效率。后續(xù)工作將重點考慮組間干擾與不完備信道狀態(tài)信息的影響，同時設(shè)計上行多小區(qū)多組場景下旋轉(zhuǎn)無人機輔助的NOMA網(wǎng)絡啟發(fā)式功率分配算法。