牛和昊 林 志 張峻寧 王 勇 王 磊
1. 國防科技大學(xué)電子對抗學(xué)院 安徽 合肥 230037 2. 電子制約技術(shù)安徽省重點實驗室 安徽 合肥 230037
現(xiàn)有的通信能源通常用固定電網(wǎng)和太陽能等方式獲取, 然而, 隨著無線傳感網(wǎng)絡(luò)的飛速發(fā)展, 傳統(tǒng)的供電方式無法滿足大規(guī)模布局的無線傳感器對能源供給的需求[1]. 當(dāng)前, 有學(xué)者提出利用空間中電磁信號自帶的電磁能量, 通過無線攜能終端設(shè)備從射頻信號中獲取能量, 該技術(shù)也被稱為無線攜能傳輸技術(shù)(simultaneous wireless information and power transfer, SWIPT)[2]. 無線攜能傳輸技術(shù)被認(rèn)為是提升無線網(wǎng)絡(luò)能源利用率和能量效率的有效途徑之一.
盡管多天線技術(shù)可以有效地增強(qiáng)信息和能量的傳輸效率[3-6], 考慮到每根發(fā)射端天線均需連接一條射頻鏈路, 這使得硬件成本、計算復(fù)雜度和能量消耗均隨著天線數(shù)的增多而大幅增加[7]. 近年來, 可重構(gòu)智能反射面(reconfigurable intelligent surface, RIS)技術(shù)由于其能大幅提高網(wǎng)絡(luò)的吞吐量和頻譜效率, 受到了學(xué)者的廣泛關(guān)注并成為當(dāng)前的研究熱點[8]. RIS 配備大量的反射元素, 每個元素均可以獨立改變?nèi)肷湫盘柕南辔? 從而實現(xiàn)對發(fā)端信號波束方向的調(diào)整, 進(jìn)而在目的端增強(qiáng)對期望信號的接收, 同時抑制干擾信號[9]. 另一方面, 由于RIS 的金屬貼片單元結(jié)構(gòu)尺寸很小, 在較小的面積附著數(shù)以百計的金屬貼片單元, 從而顯著地提升信號傳輸?shù)臒o源增益[10].
基于上述優(yōu)點, 已有學(xué)者針對RIS 在無線通信中的應(yīng)用進(jìn)行了研究. 文獻(xiàn)[11]針對RIS 輔助的無線網(wǎng)絡(luò), 研究了聯(lián)合的主動和被動波束成形設(shè)計. 文獻(xiàn)[12]研究了離散情況下的反射系數(shù)(reflecting coefficient, RC)優(yōu)化設(shè)計. 此外, 基于RIS 的傳輸設(shè)計已在多組多播網(wǎng)絡(luò)[13]和多小區(qū)網(wǎng)絡(luò)[14-15]中得到了深入的研究. 除了可以增強(qiáng)網(wǎng)絡(luò)的頻譜效率, RIS 也用于提升無線攜能網(wǎng)絡(luò)的能量采集效率[16]和竊聽網(wǎng)絡(luò)的安全性能[17-18]. 然而, 上述文獻(xiàn)中的優(yōu)化方法復(fù)雜度較高,難以直接用于實時的傳輸設(shè)計. 因此, 如何基于RIS設(shè)計出低復(fù)雜度的傳輸方案具有重要的意義.
本文針對多輸入單輸出(multiple-input singleoutput, MISO)的無線攜能傳輸網(wǎng)絡(luò), 研究了基于RIS的混合波束形成設(shè)計方案, 其中, 數(shù)字波束成形(beam forming, BF)在基站發(fā)射端實現(xiàn), 模擬BF 在RIS 進(jìn)行處理. 本文假設(shè)接收端采用了非線性能量采集模型, 并使用功率分割技術(shù)來同時實現(xiàn)信息的提取和能量的采集. 通過聯(lián)合優(yōu)化設(shè)計數(shù)字BF、模擬BF 和功率分割(power splitting, PS)因子, 在多個信干噪比(signal-to-noise-plus-noise ratio, SINR)和能量采集的約束下, 研究系統(tǒng)所有接收機(jī)信息速率總和的最大化問題. 針對此非凸優(yōu)化問題, 采用交替利用優(yōu)化和罰凸方法將原優(yōu)化問題轉(zhuǎn)化為凸問題, 并采用迭代算法進(jìn)行求解.
在實際應(yīng)用中,實現(xiàn)連續(xù)RC 的硬件成本較高.因此, 將離散RC 應(yīng)用于反射元件更為現(xiàn)實. 基于連續(xù)RC 的優(yōu)化方法可拓展至離散RC 設(shè)計,更具有一般性.
本文研究一個RIS 輔助的MISO 無線攜能傳輸系統(tǒng), 如圖1 所示, 該系統(tǒng)由一個基站、一個RIS 和N 個接收機(jī)組成. 基站配備有Nt根天線, RIS 有M 個反射元素, 每個接收機(jī)為單天線節(jié)點. 將分別定義為基站至RIS、基站至第n 個接收機(jī)以及RIS 至第n 個接收機(jī)的信道系數(shù).
圖1 MISO 下行RIS 輔助無線攜能網(wǎng)絡(luò)Fig.1 The MISO downlink SWIPT system with RIS
假設(shè)所有的信道均服從平穩(wěn)衰落信道, 則第n個接收機(jī)收到的信息為
本文假設(shè)每個接收機(jī)使用PS 方式來將接收信號分成獨立的兩部分, 即每個接收機(jī)處的接收信號被功率分配器分割, 一部分進(jìn)行信息解碼, 另一部分進(jìn)行能量收集. 因此, 第n 個接收機(jī)用于信息解碼的信號可表示為
因此, 第n 個接收機(jī)的信息速率為
此外, 第n 個接收機(jī)用于能量采集的信號可表示為
則第n 個接收機(jī)收集到的能量為
本文采用文獻(xiàn)[9]中新提出的非線性能量采集模型來反映實際的電路特性. 因此, 第n 個接收機(jī)處獲得的實際功率可由下式給出
本文研究了聯(lián)合混合BF 和PS 因子設(shè)計, 以最大限度提高所有接收機(jī)信息速率總和. 優(yōu)化問題可表述為
其中, 式(12)表示第n 個接收機(jī)的最小信息速率閾值, 式(13)表示第n 個接收機(jī)的最小接收功率閾值.此外, 式(15)表示模擬BF 的單位模約束, 因為只改變傳輸信號的相位而不改變幅度.
由于耦合變量和非凸約束的存在, 式(11)~式(15)很難直接求解.然而,當(dāng)部分變量固定時,式(11)~式(15)可以轉(zhuǎn)化為其他變量的近似凸問題. 因此, 可將式(11)~式(15)分解成幾個子問題, 并以交替優(yōu)化的方式更新這些變量[19], 下面將詳細(xì)研究該方法.
由于式(16)~式(18)是非凸的, 難以直接求解式(19). 可利用半正定松弛(semi-definite relaxation,SDR)[20]方法, 并定義, 可將優(yōu)化問題轉(zhuǎn)換為
可知限制條件式(21)~式(23)使得式(20)仍是非凸的. 采用Taylor 展開式, 在給定點附近, 可近似為
下面將針對非凸約束式(24)進(jìn)行處理, 采用迭代懲罰方法, 首先將式(24)重寫為
因此, 式(32)~式(33)可等價表示為
因此, 式(32)~式(33)可轉(zhuǎn)換為如下
由于原優(yōu)化式(36)~式(37)是凸問題, 可利用凸優(yōu)化工具包CVX 求解[21].
本節(jié)利用懲罰凹凸過程法(penalty concave convex procedure, PCCP), 在給定w 和ρn的情況下優(yōu)化. 關(guān)于的優(yōu)化可表述為如下的可行性問題
然后, 需要將非凸約束式(43)和式(44)轉(zhuǎn)化為近似的凸約束. 具體而言, 通過一階展開式, 式(43)可近似為
同上述方法, 可將式(44)近似為
至此, 優(yōu)化問題中唯一的非凸約束是恒定模約束式(45). 接下來, 將采用PCCP 方法來處理此限制條件. 根據(jù)懲罰原則, 首先將式(45)等價轉(zhuǎn)換為. 然后, 引入下述引理來轉(zhuǎn)換非凸的部分.
引理1[22]:令a 為模值為1 的復(fù)變量, 則對任意的確定點, 有下式成立:
利用以上引理, 轉(zhuǎn)換為以下問題
本部分的仿真結(jié)果評估了所提出算法的性能.仿真場景如圖2 所示, 包含1 個基站、1 個RIS、3 個接收機(jī). 基站和RIS 的坐標(biāo)分別為(0 m, 0 m),(30 m,10 m). 接收機(jī)假設(shè)均勻分布在一個圓內(nèi), 圓心為(40 m, 0 m), 半徑為5 m. 此外, 基站、RIS 和接收機(jī)的高度分別為10 m、10 m 和2 m.
圖2 仿真場景圖Fig.2 The deployment diagram of the simulation scenario
1)隨機(jī)RIS 方案, 即隨機(jī)選擇RIS 的相移矩陣;
2)無RIS 輔助的方案;
3)固定PS 的方案, 即設(shè)置固定的ρn=0.5, 僅優(yōu)化w 和θ. 這些方法分別被標(biāo)記為“本文方法”“隨機(jī)RIS 方法”“無RIS 方法”和“固定PS 方法”.
圖3 展示了本文算法的收斂性. 顯示了迭代過程的每一步中該網(wǎng)絡(luò)所有接收機(jī)的信息速率總和與迭代次數(shù)的關(guān)系. 由圖3 可以看出, 本文所提的方法在經(jīng)過15 次的迭代后, 能收斂到優(yōu)化問題的最優(yōu)解.此外, 從圖3 可以看出, 接收機(jī)信息速率總和隨著Nt或者M(jìn) 的增加而增大, 但計算復(fù)雜度也會增加, 收斂速度會變慢. 因此, 更大的Nt或者M(jìn) 可以帶來更好的性能, 代價是系統(tǒng)的復(fù)雜度升高.
圖3 本文方法的收斂性Fig.3 The convergence feature of the proposed method
圖4 展示了接收機(jī)信息速率總和與基站的最大發(fā)射功率Ps的關(guān)系. 可以看出, 與其他方法相比, 本文所提方法的性能更好, 其次是固定PS 方法, 而沒有RIS 輔助的方法性能最差. 綜合對比這幾個方法的性能可知, RIS 的存在增加了額外的基站至接收機(jī)的通信鏈路, 通過合理設(shè)計RIS 的反射系數(shù), 能使得更多的無線信號能量匯聚在接收機(jī)處, 進(jìn)而提升接收機(jī)獲取無線能量和信息的可能性, 以提高系統(tǒng)性能.
圖4 接收機(jī)信息速率總和與發(fā)射功率關(guān)系圖Fig.4 The relationship diagram between sum rate of receiver information and the transmission power
圖5 展示了接收機(jī)信息速率總和與基站的天線數(shù)Nt之間的關(guān)系. 從圖中可知, 所有方法的接收機(jī)信息速率總和均隨著Nt的增加而增加, 這是因為當(dāng)基站配備更多的天線時, 可以獲得更高的空間自由度.在此基礎(chǔ)上通過合理的設(shè)計基站數(shù)字波束成型, 能使無線信號更準(zhǔn)確地匯聚至接收機(jī)信道方向, 進(jìn)而提升信息速率.
圖5 接收機(jī)信息速率總和與發(fā)射天線數(shù)關(guān)系圖Fig.5 The relationship diagram between sum rate of receiver information and the number of transmission antenna
圖6 描繪了接收機(jī)信息速率總和與RIS 反射元素個數(shù)M 間的關(guān)系. 可以看出, 所有方案所獲得的接收機(jī)信息速率總和均隨著M 的增加而增大. 這個結(jié)果主要取決于2 個因素:1)增大M 后, RIS 處可獲得更高的陣列增益, RIS 能接收到更多的來自基站的入射信號;2)當(dāng)RIS 的反射系數(shù)得到合理優(yōu)化時, 更大的M 可使RIS 反射的信號更準(zhǔn)確地匯聚至接收機(jī)信道, 從而提高接收機(jī)接收的信號質(zhì)量.
圖6 接收機(jī)信息速率總和與反射面元素個數(shù)關(guān)系圖Fig.6 The relationship diagram between sum rate of receiver information and the number of transmission surfaces
圖7 描述了接收機(jī)信息速率總和與能量采集門限Ψn,th的關(guān)系. 可以看出, 所有方案所獲得的接收機(jī)信息速率總和均隨著Ψn,th的增加而減小. 這是因為增大Ψn,th后, 接收信號中更多的比例需要用于能量采集, 因此, 用于信息譯碼的信號比例減小, 從而降低了信息速率. 該結(jié)果表明對于無線攜能通信系統(tǒng), 需要在能量采集和信息傳輸中做折衷選擇.
圖7 接收機(jī)信息速率總和與能量采集門限關(guān)系圖Fig.7 The relationship diagram between sum rate of receiver information and the EH threshold
本文針對RIS 輔助無線攜能網(wǎng)絡(luò)中的信息速率總和的最大化問題, 設(shè)計一種混合波束成形算法, 采用交替優(yōu)化、罰凹凸過程和迭代算法求解出波束成形權(quán)矢量、反射相位和功率分割比. 仿真結(jié)果表明:通過在無線攜能網(wǎng)絡(luò)中設(shè)置RIS 并聯(lián)合優(yōu)化其反射系數(shù)和接收機(jī)的功率分配, 能有效提升無線攜能網(wǎng)絡(luò)的信息和能量傳輸效率.