摘 要:以遠(yuǎn)距離無線電(Long Range Radio,LoRa)體制為代表的啁啾(Chirp)波形在低功耗廣覆蓋的物聯(lián)網(wǎng)(Internet of Things,IoT)技術(shù)中具有廣泛的應(yīng)用。針對低地球軌道(Low Earth Orbit,LEO)衛(wèi)星IoT(Satellite-based IoT,SIoT)場景下,Chirp 波形在LEO 高動態(tài)信道上的適應(yīng)性問題,在傳統(tǒng)LoRa 體制的基礎(chǔ)上提出了一種參數(shù)可配置和可擴(kuò)展的物聯(lián)波形設(shè)計(jì)方法,基于ZYNQ+AD9361 的軟件無線電(Software Defined Radio,SDR)平臺,終端側(cè)實(shí)現(xiàn)了帶寬(Band Width,BW)、擴(kuò)頻因子(Spread Factor,SF)等參數(shù)可配置、可擴(kuò)展的LoRa 信號調(diào)制功能,網(wǎng)關(guān)側(cè)實(shí)現(xiàn)了接收信號時(shí)頻同步以及快速傅里葉變換(Fast Fourier Transform,FFT)解調(diào)算法。研究分析了LoRa 調(diào)制的各項(xiàng)參數(shù)配置,在LEO衛(wèi)星信道中的適應(yīng)性;通過誤比特率(Bit Error Ratio,BER)測試驗(yàn)證了LoRa 調(diào)制的通信性能,并證明了參數(shù)可配置、可擴(kuò)展的LoRa 調(diào)制比傳統(tǒng)LoRa 調(diào)制,在LEO 衛(wèi)星信道下有更好的適應(yīng)性。說明以LoRa 體制為代表且基于Chirp 波形的這一類物聯(lián)信號在LEO SIoT 場景下具有應(yīng)用前景。
關(guān)鍵詞:低軌衛(wèi)星物聯(lián)網(wǎng);低軌高動態(tài)信道;啁啾;參數(shù)配置;適應(yīng)性分析
中圖分類號:TN919. 23 文獻(xiàn)標(biāo)志碼:A 開放科學(xué)(資源服務(wù))標(biāo)識碼(OSID):
文章編號:1003-3114(2024)04-0739-11
0 引言
衛(wèi)星物聯(lián)網(wǎng)(Satellite-based Internet of Things,SIoT)是指通過衛(wèi)星技術(shù)實(shí)現(xiàn)物聯(lián)網(wǎng)(Internet ofThings,IoT)設(shè)備之間的連接,利用衛(wèi)星網(wǎng)絡(luò)提供全球范圍內(nèi)的覆蓋,為IoT 設(shè)備提供穩(wěn)定、廣泛的連接能力。由于地球表面超80% 的陸地和95% 以上的海洋,不適合建設(shè)地面基站提供IoT 終端接入網(wǎng)關(guān),因此將SIoT 作為地面IoT 的補(bǔ)充和延伸,作為實(shí)現(xiàn)全球IoT 服務(wù)的一種途徑。由于和中高地球軌衛(wèi)星相比,低地球軌道(Low Earth Orbit,LEO)衛(wèi)星具有傳播時(shí)延短、傳播損耗低等優(yōu)點(diǎn),使地面IoT 終端更容易滿足小型化、低功耗和低成本等要求,所以LEO 衛(wèi)星是實(shí)現(xiàn)SIoT 系統(tǒng)的最佳選擇[1-2]。
在地面IoT 中,應(yīng)用比較廣泛的通信協(xié)議有遠(yuǎn)距離無線電(Long Range Radio,LoRa)和窄帶物聯(lián)網(wǎng)(Narrow Band-Internet of Things,NB-IoT)。LoRa是Semtech 公司擁有物理層專利的通信技術(shù)[3-5],具有覆蓋范圍廣、穿透性強(qiáng)、抗干擾性強(qiáng)的特點(diǎn),其物理層技術(shù)目前尚未公開。2013 年,Semtech 公司依據(jù)IEEE 802. 15. 4a 標(biāo)準(zhǔn)發(fā)布LoRa 芯片后就得到學(xué)術(shù)界極大的關(guān)注。關(guān)于LoRa 技術(shù)的研究只是對其通信能力的測試,包括對商用LoRa 芯片開展了鏈路性能和網(wǎng)絡(luò)容量等方面的性能測試和研究[6-7]。由于Semtech 公司未公開其物理層實(shí)現(xiàn)技術(shù),國外一些研究工作利用軟件無線電(Software Defined Ra-dio,SDR)平臺復(fù)現(xiàn)了LoRa 物理層的算法[8-10],并研究LoRa 調(diào)制在LEO 衛(wèi)星信道環(huán)境下的性能[11-12],包括動態(tài)多普勒頻偏對LoRa 調(diào)制的影響[13]。國外已有科研團(tuán)隊(duì)借助Semtech 公司LoRa芯片在真實(shí)的LEO 衛(wèi)星場景下測試了LoRa 的通信性能[14],研究了LoRa 芯片的不同帶寬(BandWidth,BW)和擴(kuò)頻因子(Spread Factor,SF)等參數(shù)配置,在衛(wèi)星不同軌道高度下的適應(yīng)性。以上是學(xué)術(shù)界對于LoRa 信號的采樣、同步、吞吐量、媒體介入控制層協(xié)議等方面在LEO SIoT 場景下的適應(yīng)性改造研究。在工業(yè)界,EchoStar 將部署由28 顆LEO衛(wèi)星構(gòu)成的S-Band 衛(wèi)星星座,該28 顆衛(wèi)星通過LoRa協(xié)議連接地面終端設(shè)備。中國近年來也在建設(shè)由38 顆衛(wèi)星組成天啟SIoT 系統(tǒng),使用了LoRa 技術(shù)。
LEO 衛(wèi)星信道實(shí)質(zhì)是一個(gè)高動態(tài)信道,由于不同SF 配置的LoRa 信號具有不同的覆蓋性能,這會導(dǎo)致在某些SF 配置的LoRa 信號,在一定的衛(wèi)星俯仰角范圍內(nèi),信號中的多普勒頻偏值出現(xiàn)動態(tài)變化,當(dāng)多普勒頻偏的動態(tài)變化量超過LoRa 信號的判決門限,會導(dǎo)致通信性能急劇惡化。為了適應(yīng)LEO 高動態(tài)信道,有研究提出了一種將不同調(diào)制信息映射為不同調(diào)頻率的折疊調(diào)頻斜率鍵控調(diào)制(FCrSK)方式[15],其對抗動態(tài)多普勒頻偏的能力優(yōu)于LoRa 調(diào)制,但是抗噪聲性能與頻譜效率有所下降。此外還有研究從波形設(shè)計(jì)的角度提出了對稱啁啾(Chirp)信號,對稱Chirp 信號具有較好的互相關(guān)性,更容易被捕獲[16-17]。為了解決在大頻偏情況下,對稱Chirp 信號互相關(guān)性較差以及在LEO SIoT 場景下海量終端連接的問題,在對稱Chirp 信號的基礎(chǔ)上進(jìn)一步提出了非對稱Chirp 信號,在多址干擾情況下具有更好的捕獲性能[18]。本文研究了基于Chirp 波形的LoRa 調(diào)制信號在不同參數(shù)(BW、SF 以及中心頻率)配置下,在LEO 高動態(tài)信道中的傳輸性能,首先在終端和網(wǎng)關(guān)側(cè)分別基于SDR 平臺實(shí)現(xiàn)了一種物理層波形BW、SF 和幀長等參數(shù)可擴(kuò)展配置的LoRa 調(diào)制信號,以及LoRa 調(diào)制信號時(shí)頻同步算法和信號解調(diào)方案[19]。在此基礎(chǔ)上利用衛(wèi)星信道模擬器測試了不同參數(shù)配置的LoRa 信號在不同仰角下的通信性能,為LoRa 調(diào)制在LEO 衛(wèi)星場景下的應(yīng)用奠定基礎(chǔ)。
1 系統(tǒng)模型與分析
1. 1 LEO SIoT 場景及信道模型
基于LoRa 體制的LEO SIoT,包括IoT 終端、LEO 衛(wèi)星星座、網(wǎng)關(guān)和網(wǎng)絡(luò)服務(wù)器。IoT 終端通過上行鏈路傳輸采集的物聯(lián)信息,星上接收機(jī)解調(diào)接收到的數(shù)據(jù)包,并且通過星間鏈路間采集信息傳輸?shù)叫抨P(guān)站可見的LEO 衛(wèi)星,最后LEO 衛(wèi)星通過饋電鏈路,將數(shù)據(jù)包傳輸給信關(guān)站到核心網(wǎng)提供物聯(lián)服務(wù)?;冢蹋铮遥?體制的LEO SIoT 場景圖如圖1 所示。
地面IoT 終端到衛(wèi)星的大尺度衰落模型為典型的自由空間損耗模型,表達(dá)式為:
g(l)= (λc / (4πl)) η , (1)
式中:l 為IoT 終端和LEO 衛(wèi)星之間的距離,η 為自由空間路徑損耗系數(shù),λc 為LoRa 載波信號的波長。
本文采用基于陰影萊斯分布的簡化陸地移動衛(wèi)星信道模型[13],其小尺度衰落系數(shù)包絡(luò)的概率密度函數(shù)(Probability Density Function)表達(dá)式為:
式中:Ω 為直射路徑的平均功率,2b 為散射路徑的平均功率,m 為衰落嚴(yán)重性參數(shù),F1()為合流超幾何函數(shù)。
在LEO SIoT 場景下,不同仰角下信號從終端到衛(wèi)星經(jīng)歷的小尺度衰落不同,所以仰角在該場景下也視為影響通信性能的一個(gè)重要因素。針對基于陰影萊斯分布的簡化陸地移動衛(wèi)星信道模型,其小尺度衰落系數(shù)的參數(shù)b、m 和Ω 與仰角θ 可參考如下的經(jīng)驗(yàn)式[20]:
LEO 衛(wèi)星接收信號功率表達(dá)式為:
Pr =Pt Gr g(l)|h|2, (4)
式中:Pt 為地面IoT 終端的信號發(fā)射功率,Gr 為LEO 衛(wèi)星天線的接收增益,| h |2 為陰影萊斯衰落的信道功率增益。LEO 衛(wèi)星天線的接收增益Gr 表達(dá)式為:
Gr =ξ(πD/ λc) 2 p J1(μ)/2μ+36J3(μ)/ μ3 y 2, (5)
式中:ξ 為LEO 衛(wèi)星天線效率,D 為LEO 衛(wèi)星天線等效口徑,λc 為載波波長,J1 和J3 為第一類一階和三階貝塞爾函數(shù),μ = 2. 071 23sin(θ)/ sin(θ3 dB )(θ為信號到達(dá)仰角,θ3 dB 為LEO 衛(wèi)星接收天線波束的3 dB 帶寬)。
LEO 衛(wèi)星接收信號信噪比(Signal to NoiseRatio,SNR)表達(dá)式為:
SNR=Pr / N, (6)
式中:N=kTBW 表示LEO 衛(wèi)星系統(tǒng)中噪聲平均功率(k為玻爾茲曼常數(shù),T 為LEO 衛(wèi)星等效噪聲溫度,BW 為LEO SIoT 上行鏈路中LoRa 信號的傳輸帶寬)。
按照LoRa 標(biāo)準(zhǔn)化文件規(guī)定,SF 的取值為{7,8,9,10,11,12},不同SF 的解調(diào)SNR 門限如表1 所示。
本文將統(tǒng)一使用接收機(jī)輸入信號的SNR 和在該SNR 下解調(diào)信號的誤比特率(Bit Error Ratio,BER),來分析LEO SIoT 場景下LoRa 信號的傳輸性能。由于LEO 衛(wèi)星運(yùn)行在各個(gè)位置時(shí)仰角是不同的,對應(yīng)的SNR 也不同,因此需要分析LEO 衛(wèi)星在不同仰角下LoRa 信號的通信性能。
1. 2 LEO 衛(wèi)星信道的動態(tài)多普勒頻偏
和地面IoT 幾乎零時(shí)延和富散射傳播的通信場景不同,LEO 衛(wèi)星信道具有高時(shí)延、高動態(tài)和欠散射傳播的特性。由于LEO 衛(wèi)星的軌道高度較低,運(yùn)動速度較快,不僅會產(chǎn)生較大的多普勒頻偏,而且在某些仰角下多普勒頻譜的變化速度也相對較快,將對LoRa 信號的解調(diào)和同步造成影響:一是多普勒頻偏會導(dǎo)致接收到的LoRa 信號起始頻率出現(xiàn)偏移,從而影響正確解調(diào);二是動態(tài)的多普勒頻偏會導(dǎo)致一個(gè)數(shù)據(jù)包內(nèi),后面一段信號的實(shí)際頻偏和時(shí)頻同步估算得到的頻偏差距過大,導(dǎo)致該段信號解調(diào)結(jié)果出現(xiàn)偏移。因此LEO 衛(wèi)星信道下的多普勒效應(yīng)會導(dǎo)致接收機(jī)的解調(diào)性能下降以及BER 上升。
LEO 衛(wèi)星的軌道模型,如圖2 所示。
圖2 中箭頭所指的兩個(gè)點(diǎn)分別為LEO 衛(wèi)星和地面終端節(jié)點(diǎn),r 為地球半徑,H 為LEO 衛(wèi)星軌道高度,δ 為衛(wèi)星天線張角,θ 為仰角,φ 為LEO 衛(wèi)星和終端節(jié)點(diǎn)之間的夾角,O 為地心,v 為衛(wèi)星運(yùn)動速度。
設(shè)IoT 終端發(fā)送的信號載波頻率為fc,信號傳輸時(shí)間為t,則信號從地面終端節(jié)點(diǎn)發(fā)送到LEO 衛(wèi)星后產(chǎn)生的多普勒頻偏可表示為:
式中:G 為引力常量,M 為地球質(zhì)量,GM =3. 986 005×1014 m3 / s2 為地心引力常量。
對fd(t)求導(dǎo)后可得多普勒頻偏一階變化率:
設(shè)定衛(wèi)星軌道高度h = 500 km,地球半徑r =6 371 km,天線張角δ = 60°,載波頻率分別設(shè)定為fc =480 MHz 和fc =1 668 MHz 時(shí)LEO 衛(wèi)星場景下的多普勒頻偏(fd(t))及多普勒變化率(dfd(t)/ dt)的曲線如圖3 所示,其中虛線和實(shí)線分別表示fc =480 MHz 和fc = 1 668 MHz 下的fd(t)和dfd(t)/ dt。從圖中可以看出,當(dāng)衛(wèi)星運(yùn)行在不同時(shí)刻所對應(yīng)的fd(t)和dfd(t)/ dt 不同,且在低仰角處,fd(t)值較大,而dfd(t)/ dt 較小,在高仰角處則相反。
1. 3 LoRa 信號的多普勒適應(yīng)性分析
LoRa 標(biāo)準(zhǔn)化文件中規(guī)定的不同SF 和BW 的單個(gè)LoRa 符號周期,如表2 所示。
在LEO 衛(wèi)星信道中,多普勒頻偏變化率最高可達(dá)到-3. 5×10-7 Hz/ s,按照LoRa 標(biāo)準(zhǔn)化文件規(guī)定,一幀數(shù)據(jù)包中最多可包含的載荷符號數(shù)量為200,則在不同SF 和BW 的參數(shù)配置下,地面終端發(fā)送200 載荷符號傳輸?shù)剑蹋牛?衛(wèi)星時(shí),一幀LoRa 數(shù)據(jù)包內(nèi)產(chǎn)生的最大多普勒頻偏變化量|fdmax| 是不同的。由于一幀數(shù)據(jù)包總包長加上前導(dǎo)碼部分最長可達(dá)212. 25 個(gè)符號,若要求在此LEO 衛(wèi)星信道下實(shí)現(xiàn)正確傳輸數(shù)據(jù),則在幀末尾處必須滿足| fdmax | <BW/2(SF+1),故BW/2(SF+1)就是在動態(tài)多普勒頻移環(huán)境下LoRa 接收信號的頻率判決門限。根據(jù)該門限值可以判斷出LoRa 標(biāo)準(zhǔn)中具體哪些SF 和BW 配置下,可以在LEO 衛(wèi)星信道中適應(yīng)動態(tài)多普勒頻偏。當(dāng)載荷符號數(shù)量一定時(shí),不同SF 和BW 配置會產(chǎn)生不同總包長的數(shù)據(jù)包,而不同總包長也對應(yīng)不同的BW/2(SF+1)和|fdmax |。當(dāng)載荷長度固定為200 個(gè)符號時(shí),載波頻率分別為fc = 480 MHz 和fc = 1 668 MHz,BW/2(SF+1)和|fdmax |隨LoRa 數(shù)據(jù)包總包長變化的曲線如圖4 所示。從圖中可以看出,在fc =480 MHz 的載波頻率下,當(dāng)數(shù)據(jù)包總包長大于434. 688 ms 時(shí),就會出現(xiàn)|fdmax |>BW/2(SF+1);而在fc =1 668 MHz 的載波頻率下,當(dāng)數(shù)據(jù)包總包長大于217. 344 ms 時(shí),會出現(xiàn)|fdmax | >BW/2(SF+1)。若| fdmax | >BW/2(SF+1),則一幀LoRa 數(shù)據(jù)包的最大多普勒頻偏將超過門限值。
由圖4 可知,在LoRa 標(biāo)準(zhǔn)中,當(dāng)fc = 480 MHz/fc =1 668 MHz 時(shí),只有使得LoRa 數(shù)據(jù)包總包長小于869. 376 ms/434. 688 ms 的SF 和BW 配置,才能適應(yīng)LEO 衛(wèi)星信道環(huán)境。
2 LoRa 物理層的FPGA 設(shè)計(jì)與實(shí)現(xiàn)
本文設(shè)計(jì)了一種基于Chirp 波形且參數(shù)可配置和擴(kuò)展的LoRa 調(diào)制,其物理層系統(tǒng)是基于ZYNQ+AD9361 的SDR 平臺,通過擴(kuò)展傳統(tǒng)LoRa 調(diào)制的參數(shù)集,實(shí)現(xiàn)發(fā)射信號的參數(shù)可配置。使用Xilinx 公司型號為XC7Z030 的SoC 芯片,在PS 端通過ARM處理器配置AD9361 的寄存器參數(shù);在PL 端通過Verilog HDL 和VHDL 語言設(shè)計(jì)了LoRa 調(diào)制信號的發(fā)射機(jī)和接收機(jī)。
2. 1 發(fā)射機(jī)物理層設(shè)計(jì)與實(shí)現(xiàn)
啁啾擴(kuò)頻(Chirp Spread Spectrum,CSS)是一種基于頻率掃描的擴(kuò)頻調(diào)制技術(shù),信號的頻率會隨著時(shí)間線性地變化,形成了Chirp 信號,頻率變化可以是向上或向下的線性斜坡,不同調(diào)制信息對應(yīng)于Chirp 信號不同的初始頻率。
標(biāo)準(zhǔn)Chirp 信號的表達(dá)式為:
s(t)= ej(2πfct+πμt2),0≤t≤Ts (9)
式中:fc 為載波頻率,Ts 為Chirp 信號周期,μ =BW/ Ts 為調(diào)頻率。
根據(jù)頻率隨時(shí)間的線性遞增(μ>0)或線性遞減(μ < 0),分為標(biāo)準(zhǔn)UpChirp 信號sup (t)和標(biāo)準(zhǔn)DownChirp 信號sdown(t)。
LoRa 調(diào)制則是將傳輸信息每SF 位比特組成十進(jìn)制調(diào)制信息Ki∈{1,2,…,2SF -1},將調(diào)制信息映射為LoRa 信號的初始頻率f0 =KiBW/2SF,且SF∈{7,8,9,10,11,12}、BW∈{125 kHz,250 kHz,500 kHz},傳統(tǒng)LoRa 調(diào)制信號的表達(dá)式為:
式中:T0 = (2SF -Ki )/ BW 為已調(diào)信號的頻率跳變時(shí)刻。
本文介紹的LoRa 調(diào)制,是在基于Chirp 波形的傳統(tǒng)LoRa 調(diào)制基礎(chǔ)上,擴(kuò)展了SF 和BW 的參數(shù)配置范圍。SF 增加了6,BW 增加了1 MHz。FPGA 中設(shè)計(jì)的LoRa 信號調(diào)制器結(jié)構(gòu)如圖5 所示,調(diào)制器輸入?yún)?shù)有BW、SF、每幀發(fā)送的載荷符號數(shù)sym_num和調(diào)制數(shù)據(jù)data_in,從而實(shí)現(xiàn)對初始頻率和累加頻率對應(yīng)的“控制字fw0”和“控制字dfw”的控制。頻率控制字累加器首先將fw0 的值賦給fw,再對fw 在每個(gè)時(shí)鐘周期都累加dfw 并將其輸出,當(dāng)fw 累加到BW 對應(yīng)的頻率控制字fwBW 時(shí),直接令fw 置零。輸出的fw 進(jìn)入直接數(shù)字頻率合成器(Direct DigitalFrequency Synthesis,DDS),用來生成LoRa 調(diào)制信號的I 路波形和Q 路波形。
LoRa 調(diào)制數(shù)據(jù)包的幀格式,包括前導(dǎo)碼(Preamble)、報(bào)頭(Header)、數(shù)據(jù)載荷(Payload),前導(dǎo)碼由可變前導(dǎo)(Variable Preamble)、幀同步字(Sync Word)和幀起始分隔符(Start-of-Frame Delimiter,SFD)三部分組成,如圖6 所示。數(shù)據(jù)包主要用于接收機(jī)的時(shí)頻同步。在前導(dǎo)碼中,可變前導(dǎo)用8 個(gè)UpChirp 信號表示,幀同步字用兩個(gè)LoRa 調(diào)制信號表示,幀起始分隔符用2. 25 個(gè)標(biāo)準(zhǔn)DownChirp 信號表示。
為了在FPGA 中實(shí)現(xiàn)LoRa 調(diào)制的數(shù)據(jù)包發(fā)送功能,在發(fā)射機(jī)中設(shè)計(jì)了一套狀態(tài)機(jī)控制程序,控制LoRa 調(diào)制的參數(shù)配置和數(shù)據(jù)包發(fā)送過程。發(fā)射機(jī)的狀態(tài)轉(zhuǎn)移圖如圖7 所示。其中S0 為上電復(fù)位狀態(tài),S1 為空閑狀態(tài),S2 為配置狀態(tài),S3 為前導(dǎo)碼發(fā)送狀態(tài),S4 為報(bào)頭發(fā)送狀態(tài),S5 為載荷發(fā)送狀態(tài),S6 為結(jié)束狀態(tài)。
2. 2 接收機(jī)物理層設(shè)計(jì)與實(shí)現(xiàn)
在無時(shí)延且頻率同步的理想條件下,直接通過接收信號和本地DownChirp 信號相乘,進(jìn)行快速傅里葉變換(Fast Fourier Transfrom,FFT)后,通過索引FFT 輸出的峰值對應(yīng)的頻率值就可以獲得正確的解調(diào)值。但是在現(xiàn)實(shí)場景下,發(fā)射機(jī)和接收機(jī)之間存在通信時(shí)延,且信道中存在多普勒效應(yīng),因此現(xiàn)實(shí)中無法做到時(shí)間和頻率都同步,上述條件將導(dǎo)致載波頻率偏移和符號定時(shí)偏移,載波頻率偏移分為整數(shù)倍載波頻率偏移CFOint 和小數(shù)倍載波偏移CFOfrac,符號定時(shí)偏移分為整數(shù)倍符號定時(shí)偏移STOint 和小數(shù)倍符號定時(shí)偏移STOfrac,。此外,由于發(fā)射機(jī)和接收機(jī)的晶振頻率不完全一致將導(dǎo)致采樣頻率偏移SFO,隨著時(shí)間的增加,SFO 產(chǎn)生的影響會越來越明顯,不利于接收機(jī)對長幀的解調(diào)。
為了得到全部的CFOint、STOint、CFOfrac、STOfrac和SFO 的估計(jì)值,使用了如圖8 所示展示的LoRa幀同步算法[8]。
圖8 中輸入信號的樣本按照從左到右的順序進(jìn)行處理,而算法的不同步驟對應(yīng)于幀同步的不同階段,執(zhí)行順序是從上到下,其中一些步驟同時(shí)發(fā)生。
參考圖8,在FPGA 開發(fā)環(huán)境中實(shí)現(xiàn)了參數(shù)可擴(kuò)展LoRa 信號的接收機(jī)時(shí)頻同步算法,即CFOint、STOint、CFOfrac 和STOfrac 的估計(jì)和補(bǔ)償,FPGA 中的算法流程如圖9 所示。
為了控制接收機(jī)的LoRa 信號解調(diào)和同步過程,還設(shè)計(jì)了狀態(tài)機(jī)程序,如圖10 所示。其中,S7為上電復(fù)位狀態(tài),S8 為參數(shù)配置狀態(tài),S9 為門限檢測狀態(tài),S10 為時(shí)頻同步狀態(tài),S11 為符號解調(diào)狀態(tài),S12 為結(jié)束狀態(tài)。
以下將介紹在FPGA 中實(shí)現(xiàn)估計(jì)CFOint、STOint、CFOfrac 和STOfrac 的具體步驟。
(1)估計(jì)CFOfrac 的實(shí)現(xiàn)步驟
① 求得并保存前兩個(gè)UpChirp 前導(dǎo)碼的FFT輸出最大幅值索引的相位φi(i=1,2);
② 對這兩個(gè)相位值前后分別作差,得Δφ =φ2 -φ1;
③ 對Δφ 進(jìn)行mod(2π)運(yùn)算后,判斷輸出結(jié)果是否大于π,若大于則再減π,得到Δφ′;
④ 令Δφ′/2π,求得CFOfrac 的估計(jì)值。(2)估計(jì)STOf rac 的實(shí)現(xiàn)步驟
① 累加第4 ~ 6 個(gè)UpChirp 前導(dǎo)碼的FFT 幅值,得到累加頻譜Yk;
② 設(shè)Yk 的最大幅值的索引為i,則計(jì)算TEraw =(|Yi+1 | -|Yi-1 |)/ |Yi |;
③ 將TEraw 與存放在ROM 表中的19 個(gè)TE 值做比較,找到滿足TE<TEraw 的最后一個(gè)索引值j(索引的范圍為0 ~18),STOfrac 值按下式計(jì)算得到:
式中:ROM 表中的TE 值源于文獻(xiàn)[7]中介紹的STOfrac 估計(jì)算法,在FPGA 中量化為19 個(gè)離散值存儲于ROM 中,如圖11 所示。
(3)聯(lián)合估計(jì)CFOint 和STOint 實(shí)現(xiàn)步驟
① 累加1 ~ 8 個(gè)全部的UpChirp 前導(dǎo)碼的FFT幅值并找到最大幅值的索引sup 和第1 個(gè)DownChirp的FFT 最大幅值的索引sdown;
② 求得sup +sdown 的值,再?。恚铮洌ǎ玻樱?),結(jié)果再除以2 后,幅值賦值給中間變量CFOtemp;
③ 若CFOtemp≤2SF -2,則CFOtemp 就作為整數(shù)倍頻偏的計(jì)算結(jié)果賦給CFOint;否則將CFOtemp ≤2SF -1賦給CFOint;
④ STOint 通過sdown -CFOint 再?。恚铮洌ǎ玻樱疲┯?jì)算得出。
接收機(jī)的FPGA 模塊設(shè)計(jì)如圖12 所示。
① I、Q 兩路數(shù)據(jù)i_data_i 和i_data_q 進(jìn)入頂層模塊后,將進(jìn)入下采樣模塊進(jìn)行奈奎斯特采樣;
② 采樣信號samp_i 和samp_q 進(jìn)入去線性調(diào)頻模塊和本地DownChirp 信號做復(fù)數(shù)乘法,執(zhí)行去線性調(diào)頻操作,輸出的兩路信號為o_data_i 和o_data_q;
③ 再輸入FFT 模塊進(jìn)行FFT,得到FFT 數(shù)據(jù)fft_data;
④ 將fft_data 輸出數(shù)據(jù)解調(diào)模塊和時(shí)頻同步模塊,時(shí)頻同步模塊將分別計(jì)算出各個(gè)估計(jì)值CFOfrac、STOfrac、CFOint 和STOint;
⑤ 將以上估計(jì)值輸送給下采樣模塊、去線性調(diào)頻模塊和FFT 模塊以實(shí)現(xiàn)對時(shí)頻偏的補(bǔ)償;
⑥ 時(shí)頻偏補(bǔ)償完成后,數(shù)據(jù)解調(diào)模塊會根據(jù)輸入的fft_data,解調(diào)得到正確的載荷數(shù)據(jù),并通過端口o_demod_data 輸出。
3 測試與結(jié)果分析
系統(tǒng)射頻連接如圖13 所示。
PC 機(jī)用于控制發(fā)射端發(fā)送參數(shù)可配置的LoRa信號,接收端PC 機(jī)則用于觀察解調(diào)數(shù)據(jù),頻譜儀用于觀察發(fā)射端輸出信號的頻譜,信號源用于產(chǎn)生加性高斯白噪聲,信道模擬器用于產(chǎn)生LEO 衛(wèi)星信道下的動態(tài)多普勒頻偏。在完成參數(shù)可擴(kuò)展LoRa 信號物理層的FPGA 系統(tǒng)設(shè)計(jì)后,開展了一系列測試工作,包括LoRa 信號收發(fā)測試及LEO 衛(wèi)星信道下的BER 性能測試。搭建的測試平臺主要設(shè)備包含發(fā)射端、接收端、信號源、頻譜儀和信道模擬器等,發(fā)射端和接收端是兩塊相同的ZYNQ+AD9361 開發(fā)板。測試平臺如圖14 所示。
3. 1 收發(fā)測試
采用型號為Keysight N9010B 的頻譜分析儀,觀察發(fā)射端發(fā)送LoRa 信號的頻譜。圖15 分別展示了BW 為125 kHz 和1 MHz 時(shí)LoRa 信號的頻譜,由頻譜圖可驗(yàn)證FPGA 成功產(chǎn)生對應(yīng)BW 的LoRa 信號。
配置FPGA 中的發(fā)射機(jī)和接收機(jī)的BW 為125 kHz,SF 為8,可變前導(dǎo)長度為8,載荷符號數(shù)為4,載荷內(nèi)容為十進(jìn)制的{1,10,20,30}。在接收機(jī)的PC 端,觀察接收到的LoRa 調(diào)制信號的時(shí)頻同步結(jié)果和解調(diào)數(shù)據(jù),測試結(jié)果如圖16 和圖17 所示。圖16 中上面兩行是一段收到的LoRa 信號,而下面4 行為CFOfrac、CFOint、STOfrac 和STOint 的估計(jì)值;圖17 為接收到的4 個(gè)載荷符號的解調(diào)數(shù)據(jù),可以看出該接收機(jī)可以正確同步解調(diào)LoRa 信號。
3. 2 LEO 衛(wèi)星信道下BER 測試
為了驗(yàn)證信道衰落和多普勒頻偏對LoRa 信號傳輸?shù)挠绊懛治龅玫降慕Y(jié)論,對SF∈{8,10},BW∈{125 kHz,1 MHz}的4 組參數(shù)配置在fc =1 668 MHz的L 頻段、20° 和80° 仰角、SNR 為-8 ~ 5 dB 的多普勒頻移條件下,開展誤碼性能測試,各項(xiàng)參數(shù)如表3所示。
根據(jù)測試結(jié)果,得到不同參數(shù)配置下的LoRa調(diào)制信號在仰角下的BER 曲線,如圖18 所示。其中實(shí)線和虛線分別表示20°和80°仰角下的BER。實(shí)測時(shí)由于在SF = 8,BW = 1 MHz 配置下接收機(jī)的SNR 只有-6 dB,因此只能測得-6 ~ 5 dB 的BER數(shù)據(jù)。
由圖18 可知,在低仰角條件下,所有參數(shù)配置的LoRa 調(diào)制的誤碼性能都較好,BER 和SNR 之間呈現(xiàn)明顯的相關(guān)性;在高仰角條件下,存在部分參數(shù)配置的LoRa 信號不能適應(yīng)LEO 衛(wèi)星信道環(huán)境,具體表現(xiàn)為其BER 在所有SNR 下都非常高,而且由于動態(tài)的多普勒引起的誤碼地板效應(yīng),BER 和SNR之間沒有表現(xiàn)出明顯的相關(guān)性。
3. 3 結(jié)果分析
設(shè)定衛(wèi)星軌道高度h = 500 km,地球半徑r =6 371 km,地面終端發(fā)射功率Pt =25 mW 時(shí),計(jì)算不同仰角θ 下SNR 數(shù)值,如表4 所示。
由上文的測試和分析可知,當(dāng)衛(wèi)星和地面終端在fc = 1 668 MHz 的L 頻段下通信時(shí),低仰角條件下,由于信道衰減較大,多普勒頻偏的變化率較小,LoRa 標(biāo)準(zhǔn)中SF 較小的那一部分配置(SF=7,8),由于SNR 過低而不能適應(yīng)低仰角(20°)條件下的信道環(huán)境;在高仰角條件下,由于信道衰減較小,多普勒頻偏的變化率較大,LoRa 標(biāo)準(zhǔn)中SF 較大的那一部分配置(SF>9),由于一幀數(shù)據(jù)包較長,其最大多普勒頻偏值將超過頻率判決門限,因此不能適應(yīng)高仰角(80°)條件下的信道環(huán)境。關(guān)于fc = 480 MHz,雖然在該頻段下動態(tài)多普勒頻偏的干擾較小,但是該頻段為非授權(quán)頻段,同頻干擾較為嚴(yán)重,因此選取未來SIoT 規(guī)劃的L 頻段進(jìn)行BER 測試。本文提出參數(shù)可擴(kuò)展的LoRa 調(diào)制在原有LoRa 標(biāo)準(zhǔn)的基礎(chǔ)上,將BW 配置增加至1 MHz,不僅可以在LEO 衛(wèi)星信道下以更高的信息速率通信,而且其最大的優(yōu)勢在于與傳統(tǒng)LoRa 調(diào)制相比,對LEO 衛(wèi)星信道的適應(yīng)能力有顯著提高。例如在SF = 10 的配置時(shí),80°高仰角信道環(huán)境下,傳統(tǒng)LoRa 調(diào)制無論如何選取BW都會由于數(shù)據(jù)包持續(xù)時(shí)間過長,而使得一幀內(nèi)的最大多普勒頻偏值超過頻率判決門限,導(dǎo)致解調(diào)出錯,因此不能適應(yīng)LEO 衛(wèi)星信道環(huán)境,而在所提參數(shù)可擴(kuò)展的LoRa 調(diào)制中,將BW 擴(kuò)展至1 MHz 時(shí),一幀內(nèi)的最大多普勒頻偏值不會超過頻率判決門限,因此能夠適應(yīng)高仰角的動態(tài)多普勒頻偏,并且在低仰角的信道環(huán)境SF = 10 的配置下,也能夠完全適應(yīng)信道中的衰減。上述優(yōu)點(diǎn)在BER 測試中均得到驗(yàn)證。
4 結(jié)束語
在傳統(tǒng)LoRa 調(diào)制的基礎(chǔ)上,提出了一種基于Chirp 波形且參數(shù)可配置和擴(kuò)展的LoRa 調(diào)制方式,且在FPGA 中設(shè)計(jì)并實(shí)現(xiàn)了其物理層系統(tǒng),在20°和80°衛(wèi)星仰角的信道環(huán)境下,進(jìn)行了誤碼性能測試。得出了傳統(tǒng)LoRa 信號的調(diào)制解調(diào)方法部分參數(shù)配置不能適應(yīng)高動態(tài)LEO 衛(wèi)星信道的結(jié)論,同時(shí)驗(yàn)證了參數(shù)可擴(kuò)展的LoRa 調(diào)制對比傳統(tǒng)LoRa 調(diào)制對LEO 衛(wèi)星信道有更好適應(yīng)性的這一優(yōu)勢。未來將研究接收機(jī)時(shí)頻同步算法改造和差分LoRa 調(diào)制的方法,以適應(yīng)LEO 高動態(tài)信道,使以LoRa 體制為代表的Chirp 波形能在LEO SIoT 場景中得到更好的應(yīng)用。
參考文獻(xiàn)
[1] 茍亮,左云鵬,萬揚(yáng)洋,等. 低軌衛(wèi)星物聯(lián)網(wǎng)綜述[J].信息化研究,2022,48(5):1-9.
[2] 張千,劉詠. 國外低軌物聯(lián)網(wǎng)星座系統(tǒng)發(fā)展分析[J].空間電子技術(shù),2023,20(4):20-26.
[3] KAMOONA Z J R,ILYAS M. Investigating thePerformance of LoRa Communication for Nominal LoRaand Interleaved Chirp Spreading LoRa[C]∥2022 International Conference on Artificial Intelligence of Things(ICAIoT). Istanbul:IEEE,2022:1-7.
[4] VANGELISTA L. Frequency Shift Chirp Modulation:TheLoRa Modulation [J]. IEEE Signal Processing Letters,2017,24(12):1818-1821.
[5] SAVAUX V,FEER? G. Simple Asymptotic BER Expressions for LoRa System over Rice and Rayleigh Channels[C]∥ 2021 Wireless Telecommunications Symposium(WTS). San Framsisco:IEEE,2021:1-4.
[6] CALLEBAUT G,VAN DER PEERE L. Characterization ofLoRa PointtoPoint Path Loss:Measurement Campaignsand Modeling Considering Censored Data [J ]. IEEEInternet of Things Journal,2020,7(3):1910-1918.
[7] FERNANDERS R,OLIVEIRA R,LU?S M,et al. On theReal Capacity of LoRa Networks:The Impact of Nondestructive Communications[J]. IEEE Communications Letters,2019,23(12):2437-2441.
[8] BERNIER C,DEHMAS F,DEPAIRS N. Low ComplexityLoRa Frame Synchronization for Ultralow Power Softwaredefined Radios [J ]. IEEE Transactions onCommunications,2020,68(5):3140-3152.
[9] XHONNEUS M,AFISIADIS O,BOL D,et al. A Lowcomplexity LoRa Synchronization Algorithm Robust toSampling Time Offsets[J]. IEEE Internet of Things Journal,2022,9(5):3756-3769.
[10]TAPPAREL J,AFISIADIS O,MAYORAZ P,et al. AnOpensource LoRa Physical Layer Prototype on GNURadio[C]∥IEEE 21st International Workshop on SignalProcessing Advances in Wireless Communications(SPAWC). Atlanta:IEEE,2020:1-5.
[11]WU T W,QU D X,ZHANG G X,Research on LoRaAdaptability in the LEO Satellites Internet of Things[C]∥2019 15th International Wireless Communications & MobileComputing Conference (IWCMC). Tangier:IEEE,2019:131-135.
[12]FERNANDEZ L,RUIZDEAZUA J A,CALVERASA,et al.Assessing LoRa for SatellitetoEarth CommunicationsConsidering the Impact of Ionospheric Scintillation[J].IEEE Access,2020,8:165570-165582.
[13]DOROSHKIN A A,ZADOROZHNY A M,KUS O N,et al. Experimental Study of LoRa Modulation Immunity toDoppler Effect in CubeSat Radio Communications [J].IEEE Access,2019,7:75721-75731.
[14]DOROSHKIN A A,ZADOROZHNY A M,GOREV V N,etal. First Flighttesting of LoRa Modulation in SatelliteRadio Communications in Lowearth Orbit[J]. IEEE Access,2022,10:100006-100023.
[15]YANG C,WANG M,ZHENG L,et al. Folded ChirprateShift Keying Modulation for LEO Satellite IoT[J]. IEEEAccess,2019,7:99451-99461.
[16]QIAN Y B,MA L,LIANG X W. Symmetry Chirp SpreadSpectrum Modulation Used in LEO Satellite Internet ofThings [J ]. IEEE Communications Letters,2018,22(11):2230-2233.
[17]QIAN Y B,MA L,LIANG X W. The Performance of ChirpSignal Used in LEO Satellite Internet of Things [J]. IEEECommunications Letters,2019,23(8):1319-1322.
[18]QIAN Y B,MA L,LIANG X W. The Acquisition Methodof Symmetry Chirp Signal Used in LEO Satellite Internetof Things [J]. IEEE Communications Letters,2019,23(9):1572-1575.
[19]樓才義,徐建良,楊小牛. 軟件無線電原理與應(yīng)用(第2版)[M]. 北京:電子工業(yè)出版社,2014:118-205.
[20]CORAZZA G E,VATALARO F. A Statistical Model forLand Mobile Satellite Channels and Its Application toNongeostationary Orbit Systems[J]. IEEE Transactions onVehicular Technology,1994,43(3):738-742.
作者簡介:
王楚鳴 男,(1996—),碩士研究生。主要研究方向:衛(wèi)星物聯(lián)網(wǎng)。
洪 濤 男,(1982—),博士,副教授。主要研究方向:衛(wèi)星物聯(lián)網(wǎng)、衛(wèi)星通信等。
鐘志偉 男,(2000—),碩士研究生。主要研究方向:衛(wèi)星物聯(lián)網(wǎng)。
丁曉進(jìn) 男,(1981—),博士,副教授。主要研究方向:空間信息網(wǎng)絡(luò)、衛(wèi)星物聯(lián)網(wǎng)、頻偏智能認(rèn)知等。
張更新 男,(1967—),博士,教授。主要研究方向:衛(wèi)星通信、深空通信、空間信息網(wǎng)絡(luò)。
基金項(xiàng)目:國家自然科學(xué)基金(U21A20450,62171234,61971440);江蘇省前沿引領(lǐng)技術(shù)基礎(chǔ)研究專項(xiàng)(BK20192002);國家重點(diǎn)研發(fā)計(jì)劃(2022YFB2902600)