• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于Chirp波形的物聯(lián)信號在低軌衛(wèi)星信道上的適應(yīng)性分析

    2024-09-14 00:00:00王楚鳴洪濤鐘志偉丁曉進張更新
    無線電通信技術(shù) 2024年4期

    摘 要:以遠(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ù)可配置和可擴展的物聯(lián)波形設(shè)計方法,基于ZYNQ+AD9361 的軟件無線電(Software Defined Radio,SDR)平臺,終端側(cè)實現(xiàn)了帶寬(Band Width,BW)、擴頻因子(Spread Factor,SF)等參數(shù)可配置、可擴展的LoRa 信號調(diào)制功能,網(wǎng)關(guān)側(cè)實現(xiàn)了接收信號時頻同步以及快速傅里葉變換(Fast Fourier Transform,FFT)解調(diào)算法。研究分析了LoRa 調(diào)制的各項參數(shù)配置,在LEO衛(wèi)星信道中的適應(yīng)性;通過誤比特率(Bit Error Ratio,BER)測試驗證了LoRa 調(diào)制的通信性能,并證明了參數(shù)可配置、可擴展的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 文獻標(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ù)實現(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 的補充和延伸,作為實現(xiàn)全球IoT 服務(wù)的一種途徑。由于和中高地球軌衛(wèi)星相比,低地球軌道(Low Earth Orbit,LEO)衛(wèi)星具有傳播時延短、傳播損耗低等優(yōu)點,使地面IoT 終端更容易滿足小型化、低功耗和低成本等要求,所以LEO 衛(wèi)星是實現(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],具有覆蓋范圍廣、穿透性強、抗干擾性強的特點,其物理層技術(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 公司未公開其物理層實現(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]。國外已有科研團隊借助Semtech 公司LoRa芯片在真實的LEO 衛(wèi)星場景下測試了LoRa 的通信性能[14],研究了LoRa 芯片的不同帶寬(BandWidth,BW)和擴頻因子(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)星信道實質(zhì)是一個高動態(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è)計的角度提出了對稱啁啾(Chirp)信號,對稱Chirp 信號具有較好的互相關(guān)性,更容易被捕獲[16-17]。為了解決在大頻偏情況下,對稱Chirp 信號互相關(guān)性較差以及在LEO SIoT 場景下海量終端連接的問題,在對稱Chirp 信號的基礎(chǔ)上進一步提出了非對稱Chirp 信號,在多址干擾情況下具有更好的捕獲性能[18]。本文研究了基于Chirp 波形的LoRa 調(diào)制信號在不同參數(shù)(BW、SF 以及中心頻率)配置下,在LEO 高動態(tài)信道中的傳輸性能,首先在終端和網(wǎng)關(guān)側(cè)分別基于SDR 平臺實現(xiàn)了一種物理層波形BW、SF 和幀長等參數(shù)可擴展配置的LoRa 調(diào)制信號,以及LoRa 調(diào)制信號時頻同步算法和信號解調(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)信息,星上接收機解調(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)歷的小尺度衰落不同,所以仰角在該場景下也視為影響通信性能的一個重要因素。針對基于陰影萊斯分布的簡化陸地移動衛(wèi)星信道模型,其小尺度衰落系數(shù)的參數(shù)b、m 和Ω 與仰角θ 可參考如下的經(jīng)驗式[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)一使用接收機輸入信號的SNR 和在該SNR 下解調(diào)信號的誤比特率(Bit Error Ratio,BER),來分析LEO SIoT 場景下LoRa 信號的傳輸性能。由于LEO 衛(wèi)星運行在各個位置時仰角是不同的,對應(yīng)的SNR 也不同,因此需要分析LEO 衛(wèi)星在不同仰角下LoRa 信號的通信性能。

    1. 2 LEO 衛(wèi)星信道的動態(tài)多普勒頻偏

    和地面IoT 幾乎零時延和富散射傳播的通信場景不同,LEO 衛(wèi)星信道具有高時延、高動態(tài)和欠散射傳播的特性。由于LEO 衛(wèi)星的軌道高度較低,運動速度較快,不僅會產(chǎn)生較大的多普勒頻偏,而且在某些仰角下多普勒頻譜的變化速度也相對較快,將對LoRa 信號的解調(diào)和同步造成影響:一是多普勒頻偏會導(dǎo)致接收到的LoRa 信號起始頻率出現(xiàn)偏移,從而影響正確解調(diào);二是動態(tài)的多普勒頻偏會導(dǎo)致一個數(shù)據(jù)包內(nèi),后面一段信號的實際頻偏和時頻同步估算得到的頻偏差距過大,導(dǎo)致該段信號解調(diào)結(jié)果出現(xiàn)偏移。因此LEO 衛(wèi)星信道下的多普勒效應(yīng)會導(dǎo)致接收機的解調(diào)性能下降以及BER 上升。

    LEO 衛(wèi)星的軌道模型,如圖2 所示。

    圖2 中箭頭所指的兩個點分別為LEO 衛(wèi)星和地面終端節(jié)點,r 為地球半徑,H 為LEO 衛(wèi)星軌道高度,δ 為衛(wèi)星天線張角,θ 為仰角,φ 為LEO 衛(wèi)星和終端節(jié)點之間的夾角,O 為地心,v 為衛(wèi)星運動速度。

    設(shè)IoT 終端發(fā)送的信號載波頻率為fc,信號傳輸時間為t,則信號從地面終端節(jié)點發(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 時LEO 衛(wèi)星場景下的多普勒頻偏(fd(t))及多普勒變化率(dfd(t)/ dt)的曲線如圖3 所示,其中虛線和實線分別表示fc =480 MHz 和fc = 1 668 MHz 下的fd(t)和dfd(t)/ dt。從圖中可以看出,當(dāng)衛(wèi)星運行在不同時刻所對應(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 的單個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)星時,一幀LoRa 數(shù)據(jù)包內(nèi)產(chǎn)生的最大多普勒頻偏變化量|fdmax| 是不同的。由于一幀數(shù)據(jù)包總包長加上前導(dǎo)碼部分最長可達(dá)212. 25 個符號,若要求在此LEO 衛(wèi)星信道下實現(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ù)量一定時,不同SF 和BW 配置會產(chǎn)生不同總包長的數(shù)據(jù)包,而不同總包長也對應(yīng)不同的BW/2(SF+1)和|fdmax |。當(dāng)載荷長度固定為200 個符號時,載波頻率分別為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 時,就會出現(xiàn)|fdmax |>BW/2(SF+1);而在fc =1 668 MHz 的載波頻率下,當(dāng)數(shù)據(jù)包總包長大于217. 344 ms 時,會出現(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 時,只有使得LoRa 數(shù)據(jù)包總包長小于869. 376 ms/434. 688 ms 的SF 和BW 配置,才能適應(yīng)LEO 衛(wèi)星信道環(huán)境。

    2 LoRa 物理層的FPGA 設(shè)計與實現(xiàn)

    本文設(shè)計了一種基于Chirp 波形且參數(shù)可配置和擴展的LoRa 調(diào)制,其物理層系統(tǒng)是基于ZYNQ+AD9361 的SDR 平臺,通過擴展傳統(tǒng)LoRa 調(diào)制的參數(shù)集,實現(xiàn)發(fā)射信號的參數(shù)可配置。使用Xilinx 公司型號為XC7Z030 的SoC 芯片,在PS 端通過ARM處理器配置AD9361 的寄存器參數(shù);在PL 端通過Verilog HDL 和VHDL 語言設(shè)計了LoRa 調(diào)制信號的發(fā)射機和接收機。

    2. 1 發(fā)射機物理層設(shè)計與實現(xiàn)

    啁啾擴頻(Chirp Spread Spectrum,CSS)是一種基于頻率掃描的擴頻調(diào)制技術(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ù)頻率隨時間的線性遞增(μ>0)或線性遞減(μ < 0),分為標(biāo)準(zhǔn)UpChirp 信號sup (t)和標(biāo)準(zhǔn)DownChirp 信號sdown(t)。

    LoRa 調(diào)制則是將傳輸信息每SF 位比特組成十進制調(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)信號的頻率跳變時刻。

    本文介紹的LoRa 調(diào)制,是在基于Chirp 波形的傳統(tǒng)LoRa 調(diào)制基礎(chǔ)上,擴展了SF 和BW 的參數(shù)配置范圍。SF 增加了6,BW 增加了1 MHz。FPGA 中設(shè)計的LoRa 信號調(diào)制器結(jié)構(gòu)如圖5 所示,調(diào)制器輸入?yún)?shù)有BW、SF、每幀發(fā)送的載荷符號數(shù)sym_num和調(diào)制數(shù)據(jù)data_in,從而實現(xiàn)對初始頻率和累加頻率對應(yīng)的“控制字fw0”和“控制字dfw”的控制。頻率控制字累加器首先將fw0 的值賦給fw,再對fw 在每個時鐘周期都累加dfw 并將其輸出,當(dāng)fw 累加到BW 對應(yīng)的頻率控制字fwBW 時,直接令fw 置零。輸出的fw 進入直接數(shù)字頻率合成器(Direct DigitalFrequency Synthesis,DDS),用來生成LoRa 調(diào)制信號的I 路波形和Q 路波形。

    LoRa 調(diào)制數(shù)據(jù)包的幀格式,包括前導(dǎo)碼(Preamble)、報頭(Header)、數(shù)據(jù)載荷(Payload),前導(dǎo)碼由可變前導(dǎo)(Variable Preamble)、幀同步字(Sync Word)和幀起始分隔符(Start-of-Frame Delimiter,SFD)三部分組成,如圖6 所示。數(shù)據(jù)包主要用于接收機的時頻同步。在前導(dǎo)碼中,可變前導(dǎo)用8 個UpChirp 信號表示,幀同步字用兩個LoRa 調(diào)制信號表示,幀起始分隔符用2. 25 個標(biāo)準(zhǔn)DownChirp 信號表示。

    為了在FPGA 中實現(xiàn)LoRa 調(diào)制的數(shù)據(jù)包發(fā)送功能,在發(fā)射機中設(shè)計了一套狀態(tài)機控制程序,控制LoRa 調(diào)制的參數(shù)配置和數(shù)據(jù)包發(fā)送過程。發(fā)射機的狀態(tài)轉(zhuǎn)移圖如圖7 所示。其中S0 為上電復(fù)位狀態(tài),S1 為空閑狀態(tài),S2 為配置狀態(tài),S3 為前導(dǎo)碼發(fā)送狀態(tài),S4 為報頭發(fā)送狀態(tài),S5 為載荷發(fā)送狀態(tài),S6 為結(jié)束狀態(tài)。

    2. 2 接收機物理層設(shè)計與實現(xiàn)

    在無時延且頻率同步的理想條件下,直接通過接收信號和本地DownChirp 信號相乘,進行快速傅里葉變換(Fast Fourier Transfrom,FFT)后,通過索引FFT 輸出的峰值對應(yīng)的頻率值就可以獲得正確的解調(diào)值。但是在現(xiàn)實場景下,發(fā)射機和接收機之間存在通信時延,且信道中存在多普勒效應(yīng),因此現(xiàn)實中無法做到時間和頻率都同步,上述條件將導(dǎo)致載波頻率偏移和符號定時偏移,載波頻率偏移分為整數(shù)倍載波頻率偏移CFOint 和小數(shù)倍載波偏移CFOfrac,符號定時偏移分為整數(shù)倍符號定時偏移STOint 和小數(shù)倍符號定時偏移STOfrac,。此外,由于發(fā)射機和接收機的晶振頻率不完全一致將導(dǎo)致采樣頻率偏移SFO,隨著時間的增加,SFO 產(chǎn)生的影響會越來越明顯,不利于接收機對長幀的解調(diào)。

    為了得到全部的CFOint、STOint、CFOfrac、STOfrac和SFO 的估計值,使用了如圖8 所示展示的LoRa幀同步算法[8]。

    圖8 中輸入信號的樣本按照從左到右的順序進行處理,而算法的不同步驟對應(yīng)于幀同步的不同階段,執(zhí)行順序是從上到下,其中一些步驟同時發(fā)生。

    參考圖8,在FPGA 開發(fā)環(huán)境中實現(xiàn)了參數(shù)可擴展LoRa 信號的接收機時頻同步算法,即CFOint、STOint、CFOfrac 和STOfrac 的估計和補償,FPGA 中的算法流程如圖9 所示。

    為了控制接收機的LoRa 信號解調(diào)和同步過程,還設(shè)計了狀態(tài)機程序,如圖10 所示。其中,S7為上電復(fù)位狀態(tài),S8 為參數(shù)配置狀態(tài),S9 為門限檢測狀態(tài),S10 為時頻同步狀態(tài),S11 為符號解調(diào)狀態(tài),S12 為結(jié)束狀態(tài)。

    以下將介紹在FPGA 中實現(xiàn)估計CFOint、STOint、CFOfrac 和STOfrac 的具體步驟。

    (1)估計CFOfrac 的實現(xiàn)步驟

    ① 求得并保存前兩個UpChirp 前導(dǎo)碼的FFT輸出最大幅值索引的相位φi(i=1,2);

    ② 對這兩個相位值前后分別作差,得Δφ =φ2 -φ1;

    ③ 對Δφ 進行mod(2π)運算后,判斷輸出結(jié)果是否大于π,若大于則再減π,得到Δφ′;

    ④ 令Δφ′/2π,求得CFOfrac 的估計值。(2)估計STOf rac 的實現(xiàn)步驟

    ① 累加第4 ~ 6 個UpChirp 前導(dǎo)碼的FFT 幅值,得到累加頻譜Yk;

    ② 設(shè)Yk 的最大幅值的索引為i,則計算TEraw =(|Yi+1 | -|Yi-1 |)/ |Yi |;

    ③ 將TEraw 與存放在ROM 表中的19 個TE 值做比較,找到滿足TE<TEraw 的最后一個索引值j(索引的范圍為0 ~18),STOfrac 值按下式計算得到:

    式中:ROM 表中的TE 值源于文獻[7]中介紹的STOfrac 估計算法,在FPGA 中量化為19 個離散值存儲于ROM 中,如圖11 所示。

    (3)聯(lián)合估計CFOint 和STOint 實現(xiàn)步驟

    ① 累加1 ~ 8 個全部的UpChirp 前導(dǎo)碼的FFT幅值并找到最大幅值的索引sup 和第1 個DownChirp的FFT 最大幅值的索引sdown;

    ② 求得sup +sdown 的值,再取mod(2SF ),結(jié)果再除以2 后,幅值賦值給中間變量CFOtemp;

    ③ 若CFOtemp≤2SF -2,則CFOtemp 就作為整數(shù)倍頻偏的計算結(jié)果賦給CFOint;否則將CFOtemp ≤2SF -1賦給CFOint;

    ④ STOint 通過sdown -CFOint 再?。恚铮洌ǎ玻樱疲┯嬎愕贸?。

    接收機的FPGA 模塊設(shè)計如圖12 所示。

    ① I、Q 兩路數(shù)據(jù)i_data_i 和i_data_q 進入頂層模塊后,將進入下采樣模塊進行奈奎斯特采樣;

    ② 采樣信號samp_i 和samp_q 進入去線性調(diào)頻模塊和本地DownChirp 信號做復(fù)數(shù)乘法,執(zhí)行去線性調(diào)頻操作,輸出的兩路信號為o_data_i 和o_data_q;

    ③ 再輸入FFT 模塊進行FFT,得到FFT 數(shù)據(jù)fft_data;

    ④ 將fft_data 輸出數(shù)據(jù)解調(diào)模塊和時頻同步模塊,時頻同步模塊將分別計算出各個估計值CFOfrac、STOfrac、CFOint 和STOint;

    ⑤ 將以上估計值輸送給下采樣模塊、去線性調(diào)頻模塊和FFT 模塊以實現(xiàn)對時頻偏的補償;

    ⑥ 時頻偏補償完成后,數(shù)據(jù)解調(diào)模塊會根據(jù)輸入的fft_data,解調(diào)得到正確的載荷數(shù)據(jù),并通過端口o_demod_data 輸出。

    3 測試與結(jié)果分析

    系統(tǒng)射頻連接如圖13 所示。

    PC 機用于控制發(fā)射端發(fā)送參數(shù)可配置的LoRa信號,接收端PC 機則用于觀察解調(diào)數(shù)據(jù),頻譜儀用于觀察發(fā)射端輸出信號的頻譜,信號源用于產(chǎn)生加性高斯白噪聲,信道模擬器用于產(chǎn)生LEO 衛(wèi)星信道下的動態(tài)多普勒頻偏。在完成參數(shù)可擴展LoRa 信號物理層的FPGA 系統(tǒng)設(shè)計后,開展了一系列測試工作,包括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 時LoRa 信號的頻譜,由頻譜圖可驗證FPGA 成功產(chǎn)生對應(yīng)BW 的LoRa 信號。

    配置FPGA 中的發(fā)射機和接收機的BW 為125 kHz,SF 為8,可變前導(dǎo)長度為8,載荷符號數(shù)為4,載荷內(nèi)容為十進制的{1,10,20,30}。在接收機的PC 端,觀察接收到的LoRa 調(diào)制信號的時頻同步結(jié)果和解調(diào)數(shù)據(jù),測試結(jié)果如圖16 和圖17 所示。圖16 中上面兩行是一段收到的LoRa 信號,而下面4 行為CFOfrac、CFOint、STOfrac 和STOint 的估計值;圖17 為接收到的4 個載荷符號的解調(diào)數(shù)據(jù),可以看出該接收機可以正確同步解調(diào)LoRa 信號。

    3. 2 LEO 衛(wèi)星信道下BER 測試

    為了驗證信道衰落和多普勒頻偏對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 的多普勒頻移條件下,開展誤碼性能測試,各項參數(shù)如表3所示。

    根據(jù)測試結(jié)果,得到不同參數(shù)配置下的LoRa調(diào)制信號在仰角下的BER 曲線,如圖18 所示。其中實線和虛線分別表示20°和80°仰角下的BER。實測時由于在SF = 8,BW = 1 MHz 配置下接收機的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 時,計算不同仰角θ 下SNR 數(shù)值,如表4 所示。

    由上文的測試和分析可知,當(dāng)衛(wèi)星和地面終端在fc = 1 668 MHz 的L 頻段下通信時,低仰角條件下,由于信道衰減較大,多普勒頻偏的變化率較小,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 頻段進行BER 測試。本文提出參數(shù)可擴展的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 的配置時,80°高仰角信道環(huán)境下,傳統(tǒng)LoRa 調(diào)制無論如何選?。拢锥紩捎跀?shù)據(jù)包持續(xù)時間過長,而使得一幀內(nèi)的最大多普勒頻偏值超過頻率判決門限,導(dǎo)致解調(diào)出錯,因此不能適應(yīng)LEO 衛(wèi)星信道環(huán)境,而在所提參數(shù)可擴展的LoRa 調(diào)制中,將BW 擴展至1 MHz 時,一幀內(nèi)的最大多普勒頻偏值不會超過頻率判決門限,因此能夠適應(yīng)高仰角的動態(tài)多普勒頻偏,并且在低仰角的信道環(huán)境SF = 10 的配置下,也能夠完全適應(yīng)信道中的衰減。上述優(yōu)點在BER 測試中均得到驗證。

    4 結(jié)束語

    在傳統(tǒng)LoRa 調(diào)制的基礎(chǔ)上,提出了一種基于Chirp 波形且參數(shù)可配置和擴展的LoRa 調(diào)制方式,且在FPGA 中設(shè)計并實現(xiàn)了其物理層系統(tǒng),在20°和80°衛(wèi)星仰角的信道環(huán)境下,進行了誤碼性能測試。得出了傳統(tǒng)LoRa 信號的調(diào)制解調(diào)方法部分參數(shù)配置不能適應(yīng)高動態(tài)LEO 衛(wèi)星信道的結(jié)論,同時驗證了參數(shù)可擴展的LoRa 調(diào)制對比傳統(tǒng)LoRa 調(diào)制對LEO 衛(wèi)星信道有更好適應(yīng)性的這一優(yōu)勢。未來將研究接收機時頻同步算法改造和差分LoRa 調(diào)制的方法,以適應(yīng)LEO 高動態(tài)信道,使以LoRa 體制為代表的Chirp 波形能在LEO SIoT 場景中得到更好的應(yīng)用。

    參考文獻

    [1] 茍亮,左云鵬,萬揚洋,等. 低軌衛(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)。

    丁曉進 男,(1981—),博士,副教授。主要研究方向:空間信息網(wǎng)絡(luò)、衛(wèi)星物聯(lián)網(wǎng)、頻偏智能認(rèn)知等。

    張更新 男,(1967—),博士,教授。主要研究方向:衛(wèi)星通信、深空通信、空間信息網(wǎng)絡(luò)。

    基金項目:國家自然科學(xué)基金(U21A20450,62171234,61971440);江蘇省前沿引領(lǐng)技術(shù)基礎(chǔ)研究專項(BK20192002);國家重點研發(fā)計劃(2022YFB2902600)

    久久这里有精品视频免费| 国产中年淑女户外野战色| 国产精品国产三级国产专区5o| 男人狂女人下面高潮的视频| 三级国产精品片| 乱码一卡2卡4卡精品| 国产伦理片在线播放av一区| www.av在线官网国产| 国产免费福利视频在线观看| 看非洲黑人一级黄片| 亚洲美女视频黄频| 久久影院123| 午夜福利视频精品| 黄色视频在线播放观看不卡| 丝瓜视频免费看黄片| 赤兔流量卡办理| 国产在线免费精品| 在线观看av片永久免费下载| 最新中文字幕久久久久| 亚洲国产精品专区欧美| 亚洲一级一片aⅴ在线观看| 国产精品欧美亚洲77777| 国产视频首页在线观看| 欧美另类一区| 午夜久久久在线观看| 亚洲欧洲日产国产| 男男h啪啪无遮挡| 嫩草影院入口| 久久韩国三级中文字幕| 国产亚洲精品久久久com| 亚洲欧美一区二区三区国产| 亚洲精品中文字幕在线视频 | 99精国产麻豆久久婷婷| 日韩不卡一区二区三区视频在线| 男男h啪啪无遮挡| 久久久久久久国产电影| 少妇高潮的动态图| 久久国产乱子免费精品| 免费观看无遮挡的男女| 三上悠亚av全集在线观看 | 亚洲,一卡二卡三卡| 欧美3d第一页| kizo精华| av免费观看日本| 人人妻人人添人人爽欧美一区卜| 建设人人有责人人尽责人人享有的| 久久久精品免费免费高清| 国产免费一级a男人的天堂| 国产精品人妻久久久影院| 免费观看的影片在线观看| 新久久久久国产一级毛片| h日本视频在线播放| 精品国产乱码久久久久久小说| a级毛片免费高清观看在线播放| 春色校园在线视频观看| 亚洲国产精品国产精品| 国产高清不卡午夜福利| 黑人高潮一二区| 在线观看免费高清a一片| 欧美日韩亚洲高清精品| 久久久久久久久久人人人人人人| 99久久精品热视频| 久久久久久久久久久久大奶| 又粗又硬又长又爽又黄的视频| 男男h啪啪无遮挡| 亚洲不卡免费看| 亚洲中文av在线| 如日韩欧美国产精品一区二区三区 | 男女啪啪激烈高潮av片| 国产精品秋霞免费鲁丝片| av国产精品久久久久影院| 99热6这里只有精品| 午夜视频国产福利| 黄片无遮挡物在线观看| av在线播放精品| 欧美成人精品欧美一级黄| 国产探花极品一区二区| 搡女人真爽免费视频火全软件| 久久久久久久精品精品| 亚洲国产精品国产精品| 久久久久国产网址| 久久久久久久久大av| 黄色毛片三级朝国网站 | 亚洲av福利一区| 热99国产精品久久久久久7| 成人美女网站在线观看视频| 亚洲av国产av综合av卡| 精品少妇久久久久久888优播| 婷婷色av中文字幕| h日本视频在线播放| 久久久久久久久大av| 免费观看的影片在线观看| 中文字幕久久专区| 亚洲一区二区三区欧美精品| 久久97久久精品| 国产91av在线免费观看| 国产高清不卡午夜福利| 国产亚洲午夜精品一区二区久久| 日本黄色片子视频| 国产黄片美女视频| 国产成人精品婷婷| 99热6这里只有精品| 伦精品一区二区三区| 中文欧美无线码| 天美传媒精品一区二区| 97在线人人人人妻| 肉色欧美久久久久久久蜜桃| 午夜激情久久久久久久| 国产精品久久久久久久久免| 久久久国产精品麻豆| 日本vs欧美在线观看视频 | 亚洲av国产av综合av卡| 国精品久久久久久国模美| 欧美日韩亚洲高清精品| 男女免费视频国产| 亚洲精品国产成人久久av| 国产精品一二三区在线看| 尾随美女入室| 国产一区二区在线观看av| 一区二区三区乱码不卡18| 亚洲国产毛片av蜜桃av| 亚洲精品亚洲一区二区| 日本爱情动作片www.在线观看| 久久久久久久国产电影| 全区人妻精品视频| 免费久久久久久久精品成人欧美视频 | 天堂中文最新版在线下载| 国产精品欧美亚洲77777| 国产中年淑女户外野战色| 久久精品国产亚洲av涩爱| 亚洲美女黄色视频免费看| 久久97久久精品| 久久久久久人妻| 久久女婷五月综合色啪小说| 爱豆传媒免费全集在线观看| 免费在线观看成人毛片| 欧美激情极品国产一区二区三区 | 女人久久www免费人成看片| www.色视频.com| 99久久精品国产国产毛片| 两个人免费观看高清视频 | 狠狠精品人妻久久久久久综合| 国产精品一区二区在线观看99| 久久久久久久久大av| 中文乱码字字幕精品一区二区三区| 亚洲精华国产精华液的使用体验| 亚洲av成人精品一二三区| 日韩免费高清中文字幕av| 国产成人精品一,二区| 精品人妻偷拍中文字幕| 91在线精品国自产拍蜜月| 久久久久久久久久久丰满| 精品国产露脸久久av麻豆| 自拍欧美九色日韩亚洲蝌蚪91 | 这个男人来自地球电影免费观看 | 在线 av 中文字幕| 国产在线视频一区二区| 777米奇影视久久| 最近的中文字幕免费完整| 亚洲久久久国产精品| 国内揄拍国产精品人妻在线| 成年人午夜在线观看视频| 国产在线男女| 精品久久久久久久久亚洲| 91精品国产国语对白视频| 久久久a久久爽久久v久久| 国产永久视频网站| 国精品久久久久久国模美| 新久久久久国产一级毛片| 日韩欧美 国产精品| 精品久久国产蜜桃| 自线自在国产av| 永久网站在线| 91在线精品国自产拍蜜月| 国产精品麻豆人妻色哟哟久久| 日韩欧美 国产精品| 9色porny在线观看| 精品久久久噜噜| 日产精品乱码卡一卡2卡三| 久久青草综合色| 丁香六月天网| 五月伊人婷婷丁香| 欧美日韩av久久| 亚洲国产日韩一区二区| 精品久久国产蜜桃| 中文字幕制服av| 欧美三级亚洲精品| 啦啦啦视频在线资源免费观看| 久久久a久久爽久久v久久| 汤姆久久久久久久影院中文字幕| av国产久精品久网站免费入址| 男女无遮挡免费网站观看| 欧美精品一区二区免费开放| 五月天丁香电影| 亚洲,欧美,日韩| 亚洲国产精品专区欧美| 国产男人的电影天堂91| 婷婷色av中文字幕| 午夜福利网站1000一区二区三区| 日本爱情动作片www.在线观看| 国产成人精品一,二区| 欧美精品一区二区免费开放| 能在线免费看毛片的网站| 成人午夜精彩视频在线观看| 亚洲国产精品专区欧美| 亚洲国产欧美日韩在线播放 | 97精品久久久久久久久久精品| 看免费成人av毛片| 97在线视频观看| 高清毛片免费看| 一级,二级,三级黄色视频| 久久97久久精品| 男女无遮挡免费网站观看| 国产熟女午夜一区二区三区 | 亚洲av成人精品一区久久| 久久99精品国语久久久| 女性被躁到高潮视频| 久久久久国产网址| 成人漫画全彩无遮挡| av黄色大香蕉| 蜜臀久久99精品久久宅男| 18禁在线无遮挡免费观看视频| 精品一品国产午夜福利视频| 一本色道久久久久久精品综合| 麻豆乱淫一区二区| 亚洲av欧美aⅴ国产| 男男h啪啪无遮挡| 久久久欧美国产精品| 日本wwww免费看| 国产黄色视频一区二区在线观看| 成人亚洲精品一区在线观看| 亚洲情色 制服丝袜| 亚洲国产毛片av蜜桃av| 国产一区二区在线观看av| 免费在线观看成人毛片| 欧美精品一区二区大全| 日韩欧美精品免费久久| 人妻制服诱惑在线中文字幕| 欧美精品国产亚洲| 亚洲国产精品国产精品| 午夜av观看不卡| 女性被躁到高潮视频| 女人久久www免费人成看片| av一本久久久久| 99热这里只有精品一区| 色94色欧美一区二区| 少妇精品久久久久久久| 水蜜桃什么品种好| 欧美变态另类bdsm刘玥| 日韩一本色道免费dvd| 菩萨蛮人人尽说江南好唐韦庄| 王馨瑶露胸无遮挡在线观看| 在线观看www视频免费| 免费看不卡的av| 蜜臀久久99精品久久宅男| 一级毛片久久久久久久久女| 在线观看www视频免费| 国产欧美亚洲国产| 亚洲欧美成人综合另类久久久| 亚洲国产最新在线播放| av视频免费观看在线观看| 只有这里有精品99| av国产精品久久久久影院| 又爽又黄a免费视频| 国产有黄有色有爽视频| 国产成人精品婷婷| 少妇人妻精品综合一区二区| 高清视频免费观看一区二区| 深夜a级毛片| 男男h啪啪无遮挡| 在线观看一区二区三区激情| av天堂中文字幕网| 一级毛片 在线播放| 中文字幕精品免费在线观看视频 | 亚洲一级一片aⅴ在线观看| 欧美成人精品欧美一级黄| 日韩一本色道免费dvd| 婷婷色综合www| 少妇高潮的动态图| av在线播放精品| 免费看不卡的av| 欧美精品国产亚洲| 亚洲第一av免费看| 色吧在线观看| 久久精品国产a三级三级三级| 欧美日韩在线观看h| 亚洲国产精品成人久久小说| 夜夜看夜夜爽夜夜摸| 少妇猛男粗大的猛烈进出视频| 97超碰精品成人国产| 久久人人爽人人爽人人片va| 久久国产精品男人的天堂亚洲 | 亚洲欧美精品专区久久| 51国产日韩欧美| 全区人妻精品视频| 国语对白做爰xxxⅹ性视频网站| 久久亚洲国产成人精品v| 欧美精品国产亚洲| 国产色婷婷99| 亚州av有码| 内地一区二区视频在线| 深夜a级毛片| 国产中年淑女户外野战色| 我的老师免费观看完整版| 国产精品久久久久久av不卡| www.av在线官网国产| 偷拍熟女少妇极品色| 久热这里只有精品99| 日日撸夜夜添| 国产亚洲av片在线观看秒播厂| 最黄视频免费看| 日韩一区二区三区影片| 亚洲激情五月婷婷啪啪| 久久久久网色| 欧美精品一区二区大全| 国产 精品1| 国产熟女欧美一区二区| 午夜福利视频精品| 亚洲国产日韩一区二区| 久久99热这里只频精品6学生| 成人午夜精彩视频在线观看| 欧美少妇被猛烈插入视频| 又爽又黄a免费视频| 视频区图区小说| 成人二区视频| 2022亚洲国产成人精品| 国产成人aa在线观看| 国产成人精品无人区| 精品久久久精品久久久| 精品国产一区二区三区久久久樱花| 人人澡人人妻人| kizo精华| 天堂8中文在线网| 另类亚洲欧美激情| 亚洲精品国产成人久久av| 亚洲丝袜综合中文字幕| 国产男女内射视频| 国产色爽女视频免费观看| 日本wwww免费看| 亚洲精品国产av蜜桃| 少妇人妻一区二区三区视频| 丝袜在线中文字幕| 国产极品粉嫩免费观看在线 | 在线观看人妻少妇| 97超视频在线观看视频| 少妇精品久久久久久久| 精品午夜福利在线看| av在线老鸭窝| 亚洲精品乱久久久久久| 成人综合一区亚洲| 一边亲一边摸免费视频| 80岁老熟妇乱子伦牲交| 97在线视频观看| 日韩精品免费视频一区二区三区 | 国产成人精品福利久久| 日韩大片免费观看网站| 亚洲精品中文字幕在线视频 | 97在线视频观看| 久久久久久久久久久丰满| 午夜免费观看性视频| av线在线观看网站| 男的添女的下面高潮视频| 丰满乱子伦码专区| 一级爰片在线观看| 在线观看av片永久免费下载| 日韩免费高清中文字幕av| 人人妻人人爽人人添夜夜欢视频 | 最黄视频免费看| 午夜老司机福利剧场| 日韩电影二区| 日韩精品有码人妻一区| 99久久精品国产国产毛片| 久久亚洲国产成人精品v| 亚州av有码| 天堂俺去俺来也www色官网| 看十八女毛片水多多多| 九色成人免费人妻av| 国产精品久久久久久久久免| 制服丝袜香蕉在线| 免费大片18禁| 伊人久久国产一区二区| av播播在线观看一区| 久久国内精品自在自线图片| 看免费成人av毛片| av一本久久久久| 高清在线视频一区二区三区| 国产精品偷伦视频观看了| 亚洲四区av| 男女无遮挡免费网站观看| 国产成人精品婷婷| 午夜免费男女啪啪视频观看| 在线观看国产h片| 另类亚洲欧美激情| 午夜免费鲁丝| 久久久国产精品麻豆| 王馨瑶露胸无遮挡在线观看| 国产av码专区亚洲av| 一本大道久久a久久精品| 精品国产露脸久久av麻豆| 国产在视频线精品| 日本-黄色视频高清免费观看| 亚洲欧美清纯卡通| 亚洲av欧美aⅴ国产| 狂野欧美白嫩少妇大欣赏| 日韩av不卡免费在线播放| 少妇人妻精品综合一区二区| 欧美日韩在线观看h| 久久久久久久大尺度免费视频| 国产 一区精品| 国产精品久久久久久久久免| 亚洲国产av新网站| 夜夜爽夜夜爽视频| 午夜久久久在线观看| 久久久欧美国产精品| 丰满迷人的少妇在线观看| 高清视频免费观看一区二区| 91精品伊人久久大香线蕉| 亚洲经典国产精华液单| 精品国产国语对白av| 亚洲成人av在线免费| 欧美 亚洲 国产 日韩一| 人妻夜夜爽99麻豆av| 丝袜脚勾引网站| 日韩一本色道免费dvd| 国产精品国产三级国产专区5o| 午夜av观看不卡| 国产欧美日韩一区二区三区在线 | 久久青草综合色| 日韩熟女老妇一区二区性免费视频| 最后的刺客免费高清国语| 成年美女黄网站色视频大全免费 | 成人午夜精彩视频在线观看| 久久综合国产亚洲精品| 黄色配什么色好看| 99九九在线精品视频 | 麻豆乱淫一区二区| 精华霜和精华液先用哪个| 黑人巨大精品欧美一区二区蜜桃 | 精华霜和精华液先用哪个| 欧美成人午夜免费资源| 中文精品一卡2卡3卡4更新| 两个人免费观看高清视频 | 精品国产露脸久久av麻豆| 国产精品熟女久久久久浪| 久久热精品热| 香蕉精品网在线| 精品久久久久久电影网| 久久久国产欧美日韩av| 2018国产大陆天天弄谢| 国产成人一区二区在线| 一区二区三区免费毛片| 免费观看av网站的网址| 免费在线观看成人毛片| 一本色道久久久久久精品综合| 国产精品一区www在线观看| 国产高清不卡午夜福利| 久久韩国三级中文字幕| 午夜影院在线不卡| 成人国产麻豆网| 国产成人aa在线观看| 老司机影院毛片| 欧美国产精品一级二级三级 | 一级毛片 在线播放| 美女视频免费永久观看网站| 成年美女黄网站色视频大全免费 | 欧美 亚洲 国产 日韩一| 国产成人精品福利久久| 亚洲av电影在线观看一区二区三区| 国产精品麻豆人妻色哟哟久久| 精品久久久久久电影网| 中文字幕亚洲精品专区| 亚洲精品,欧美精品| 午夜激情福利司机影院| 永久免费av网站大全| 大香蕉97超碰在线| av福利片在线观看| 国产av国产精品国产| 99精国产麻豆久久婷婷| 亚洲欧美一区二区三区国产| 五月伊人婷婷丁香| 丁香六月天网| 狂野欧美激情性xxxx在线观看| 亚洲国产精品999| 国产精品久久久久久av不卡| 国产女主播在线喷水免费视频网站| 丰满饥渴人妻一区二区三| 美女主播在线视频| 国产亚洲av片在线观看秒播厂| 男女边吃奶边做爰视频| 国产欧美亚洲国产| 色婷婷久久久亚洲欧美| 国产极品天堂在线| 成年av动漫网址| 一本色道久久久久久精品综合| 欧美变态另类bdsm刘玥| 久久韩国三级中文字幕| 永久网站在线| 亚洲av综合色区一区| 一区二区av电影网| 夜夜爽夜夜爽视频| 嫩草影院新地址| 成人亚洲精品一区在线观看| 五月天丁香电影| 久久久午夜欧美精品| 男人狂女人下面高潮的视频| 国产成人免费观看mmmm| av视频免费观看在线观看| 一区二区三区四区激情视频| 日本-黄色视频高清免费观看| 亚洲人与动物交配视频| 亚洲欧美成人综合另类久久久| 久久久国产欧美日韩av| 一级二级三级毛片免费看| 国产白丝娇喘喷水9色精品| 亚洲四区av| 只有这里有精品99| 国产日韩欧美亚洲二区| 亚洲国产最新在线播放| 天天躁夜夜躁狠狠久久av| 如何舔出高潮| 国产成人免费观看mmmm| 中文资源天堂在线| 一区二区三区免费毛片| 亚洲精品亚洲一区二区| 国产在线一区二区三区精| 午夜激情久久久久久久| 免费大片黄手机在线观看| 久久久久网色| 少妇被粗大猛烈的视频| h视频一区二区三区| 最近中文字幕高清免费大全6| 肉色欧美久久久久久久蜜桃| 亚洲精品成人av观看孕妇| 日本wwww免费看| 成年女人在线观看亚洲视频| 91在线精品国自产拍蜜月| 三级国产精品欧美在线观看| 女的被弄到高潮叫床怎么办| 国产精品免费大片| h日本视频在线播放| 国产黄色视频一区二区在线观看| 国内揄拍国产精品人妻在线| 国产av一区二区精品久久| 91精品伊人久久大香线蕉| 一级,二级,三级黄色视频| 中文字幕免费在线视频6| 亚洲精品乱久久久久久| 在线精品无人区一区二区三| 夜夜骑夜夜射夜夜干| 久久免费观看电影| 欧美日韩亚洲高清精品| 一级a做视频免费观看| 成年人免费黄色播放视频 | 国产精品.久久久| 中文乱码字字幕精品一区二区三区| 婷婷色综合大香蕉| 一个人免费看片子| 亚洲情色 制服丝袜| 永久网站在线| 欧美日韩亚洲高清精品| 多毛熟女@视频| 色94色欧美一区二区| 秋霞在线观看毛片| 一级毛片 在线播放| 中文资源天堂在线| 热99国产精品久久久久久7| 自拍偷自拍亚洲精品老妇| 老女人水多毛片| 成人影院久久| 国产av码专区亚洲av| 在线观看www视频免费| 搡女人真爽免费视频火全软件| av黄色大香蕉| 大码成人一级视频| 中文字幕人妻熟人妻熟丝袜美| 国产精品秋霞免费鲁丝片| 国产乱人偷精品视频| 少妇人妻 视频| 高清毛片免费看| 国产精品国产三级国产av玫瑰| 精品久久久噜噜| 一边亲一边摸免费视频| 国产精品国产av在线观看| 综合色丁香网| 国产精品.久久久| 国产精品久久久久成人av| 久久热精品热| 男人狂女人下面高潮的视频| 一级毛片 在线播放| 少妇被粗大的猛进出69影院 | h视频一区二区三区| 丰满迷人的少妇在线观看| 久久久久国产网址| 夜夜骑夜夜射夜夜干| 精品人妻偷拍中文字幕| 欧美丝袜亚洲另类| 人人妻人人澡人人看| 一个人免费看片子| 一级毛片电影观看| 国产一区二区三区综合在线观看 | 国产精品秋霞免费鲁丝片| 亚洲欧美成人精品一区二区| 亚洲四区av| 一区二区av电影网| 亚洲第一av免费看| 国产综合精华液| 精品酒店卫生间| 国产女主播在线喷水免费视频网站| 国产精品熟女久久久久浪| 国产成人精品久久久久久| 一本一本综合久久| 久久99精品国语久久久| 熟妇人妻不卡中文字幕| 黑人高潮一二区| 国产伦理片在线播放av一区|