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

    部分連接寬帶毫米波全雙工MIMO系統(tǒng)混合波束成形設(shè)計(jì)

    2024-04-26 03:15:40王庚善楊志家李志航
    關(guān)鍵詞:全雙工干擾信號(hào)波束

    王庚善,楊志家,李志航,4

    1(中國科學(xué)院 沈陽自動(dòng)化研究所 機(jī)器人學(xué)國家重點(diǎn)實(shí)驗(yàn)室,沈陽 110016) 2(中國科學(xué)院 網(wǎng)絡(luò)化控制系統(tǒng)重點(diǎn)實(shí)驗(yàn)室,沈陽 110016) 3(中國科學(xué)院 機(jī)器人與智能制造創(chuàng)新研究院,沈陽 110169) 4(中國科學(xué)院大學(xué),北京 100049)

    0 引 言

    近年來,無線通信技術(shù)迅猛發(fā)展,5G規(guī)?;瘧?yīng)用不斷推進(jìn),6G相關(guān)研究逐步展開,進(jìn)入萬物互聯(lián)時(shí)代,在傳輸速率、時(shí)延、功耗等方面提出了新的更高要求[1-3].為此,許多新興的物理層技術(shù)被廣泛提出并深入研究,其中包括毫米波、大規(guī)模多輸入多輸出(Multiple-Input Multiple-Output,MIMO)和全雙工無線通信等[3-5].毫米波是指頻率在26.5-300GHz的電磁波[6],該頻段可以提供豐富的頻譜資源用于無線通信,從根本上解決Sub-6 GHz頻段的頻譜資源緊張問題,從而大幅提升傳輸速率[7];大規(guī)模MIMO利用多天線波束成形實(shí)現(xiàn)空間復(fù)用和高度定向收發(fā)數(shù)據(jù),可以降低通信干擾,提高通信速率[8,9];全雙工無線通信是指同時(shí)同頻進(jìn)行數(shù)據(jù)發(fā)送與接收,理論上最大可以將頻譜效率提升為傳統(tǒng)時(shí)分或頻分雙工的兩倍,同時(shí)也可以降低傳輸時(shí)延[10].這3種技術(shù)分別從不同角度來提高傳輸速率.由于毫米波波長較短,雖然路徑損耗嚴(yán)重,但是相比于Sub-6 GHz頻段,在相同的物理尺寸內(nèi),可以配備更多天線,因此毫米波全雙工MIMO系統(tǒng)可以將上述3種新興的物理層技術(shù)有效地結(jié)合起來,是進(jìn)一步提升傳輸速率的可行方法[11,12].

    由于全雙工系統(tǒng)可以同時(shí)同頻發(fā)送和接收信號(hào),所以其自身的發(fā)送信號(hào)也會(huì)被接收端接收,稱之為“自干擾”信號(hào)[10].由于全雙工節(jié)點(diǎn)自身接收與發(fā)送天線之間的距離遠(yuǎn)小于接收天線與遠(yuǎn)端節(jié)點(diǎn)發(fā)送天線之間的距離,因此自干擾信號(hào)的強(qiáng)度遠(yuǎn)高于有效接收信號(hào).如果無法有效抑制自干擾信號(hào),系統(tǒng)的頻譜效率將會(huì)顯著降低.同時(shí),由于毫米波射頻鏈路的成本和功耗較高[7],因此實(shí)現(xiàn)毫米波全雙工MIMO通信的關(guān)鍵問題是在可接受的硬件功耗和復(fù)雜度前提下,如何抑制自干擾信號(hào)、提升頻譜效率.自干擾信號(hào)通常需要采用被動(dòng)抑制與主動(dòng)抑制方法相結(jié)合,在空域、模擬域和數(shù)字域進(jìn)行聯(lián)合抑制,使剩余自干擾信號(hào)的功率在合理范圍[10,13].文獻(xiàn)[14]搭建了第一套完整的毫米波全雙工系統(tǒng),通過放置天線信號(hào)自反射器、設(shè)計(jì)模擬域?qū)S眯酒约皵?shù)字域算法,實(shí)現(xiàn)了超過80 dB的自干擾抑制,頻譜效率是半雙工系統(tǒng)的1.7倍左右.該研究證明了毫米波全雙工通信的可行性和優(yōu)越性,但是僅考慮單輸入單輸出系統(tǒng).由于所提方法需要自反射器和專用芯片,其硬件和計(jì)算復(fù)雜度會(huì)隨著全雙工節(jié)點(diǎn)天線數(shù)目的二次方增長,所以,該方法很難擴(kuò)展到毫米波全雙工MIMO系統(tǒng).在Sub-6 GHz頻段的全雙工MIMO系統(tǒng)中,利用天線選擇、波束選取、零空間投影以及發(fā)送和接收波束成形等不需要任何輔助硬件電路的方法進(jìn)行自干擾抑制和頻譜效率優(yōu)化被廣泛研究[15-17].這些方法通常基于全數(shù)字波束成形來實(shí)現(xiàn),而全數(shù)字波束成形要求每根天線連接一條單獨(dú)的射頻鏈路,考慮到毫米波射頻鏈路較高的成本和功耗,故基于全數(shù)字波束成形的方法在毫米波全雙工MIMO系統(tǒng)中不可行.為了進(jìn)一步降低系統(tǒng)功耗和硬件復(fù)雜度,采用模擬波束成形的毫米波全雙工MIMO系統(tǒng)被提出并研究[11,18].模擬波束成形只需要一條射頻鏈路,所有天線通過移相器與射頻鏈路相連,是硬件復(fù)雜度最低、最簡單的毫米波波束成形方法[7].但是,模擬波束成形無法生成多個(gè)波束、無法支持多數(shù)據(jù)流、多用戶通信,嚴(yán)重限制了毫米波全雙工通信的性能和應(yīng)用.

    為了進(jìn)一步提升系統(tǒng)性能同時(shí)保持較低功耗和復(fù)雜度,基于混合波束成形的毫米波全雙工MIMO系統(tǒng)得到越來越多的關(guān)注與研究[19-24].混合波束成形包括模擬和數(shù)字波束成形兩部分,其中模擬波束成形通常由毫米波頻段的移相器網(wǎng)絡(luò)實(shí)現(xiàn),數(shù)字波束成形則在基帶完成[7].相比于全數(shù)字波束成形,混合波束成形需要更少的射頻鏈路,功耗和復(fù)雜度低;相比于模擬波束成形,混合波束成形支持多數(shù)據(jù)流、多用戶通信,性能更優(yōu).根據(jù)天線與射頻鏈路的連接方式,混合波束成形結(jié)構(gòu)又分為全連接結(jié)構(gòu)(每根天線通過移相器與每條射頻鏈路相連)和部分連接結(jié)構(gòu)(每根天線通過移相器僅與一條射頻鏈路相連).全連接結(jié)構(gòu)的性能更優(yōu),但是部分連接的功耗和硬件復(fù)雜度更低.當(dāng)前關(guān)于毫米波全雙工MIMO系統(tǒng)混合波束成形的研究多集中于全連接結(jié)構(gòu)[19-22].文獻(xiàn)[19]以最小化自干擾信號(hào)的功率作為優(yōu)化目標(biāo),利用拉格朗日乘子法和子空間分解方法求出混合波束成形矩陣,所提方法能夠?qū)⒆愿蓴_信號(hào)抑制超過30 dB,但是該研究僅考慮窄帶系統(tǒng)和無限精度移相器.文獻(xiàn)[20]提出了基于稀疏空間編碼和零空間投影的混合波束成形方法來抑制自干擾信號(hào),能夠有效提升頻譜效率.在后續(xù)研究中,該方法又被進(jìn)一步擴(kuò)展至寬帶毫米波全雙工系統(tǒng)[21].文獻(xiàn)[22]考慮采用有限精度移相器的毫米波全雙工MIMO系統(tǒng),以最大化系統(tǒng)頻譜效率作為優(yōu)化目標(biāo),采用懲罰對(duì)偶分解框架進(jìn)行混合波束成形設(shè)計(jì),得到近似最優(yōu)解,然而該研究僅考慮單天線用戶和窄帶系統(tǒng).針對(duì)部分連接毫米波全雙工MIMO系統(tǒng),當(dāng)前研究相對(duì)較少[23,24].文獻(xiàn)[23]提出了基于半正定松弛(Semidefinite Relaxation,SDR)、交替最小化和零空間投影的自干擾抑制混合波束成形方法,在特定條件下,該方法實(shí)現(xiàn)的頻譜效率要優(yōu)于半雙工通信.該研究證明了部分連接毫米波全雙工通信的可行性,但是計(jì)算復(fù)雜度較高、僅考慮窄帶系統(tǒng)和無限精度移相器.文獻(xiàn)[24]考慮寬帶毫米波全雙工系統(tǒng),提出了波束成形和自干擾抑制聯(lián)合設(shè)計(jì)方法.該方法能夠有效抑制自干擾信號(hào)、提高接收功率,但是在建模時(shí)并沒有考慮與節(jié)點(diǎn)之間距離相關(guān)的上下行信道的實(shí)際路徑損耗.

    針對(duì)上述相關(guān)研究中存在的不足,本文對(duì)部分連接寬帶毫米波全雙工MIMO系統(tǒng)的混合波束成形問題進(jìn)行研究,主要?jiǎng)?chuàng)新點(diǎn)和貢獻(xiàn)總結(jié)如下:

    1)考慮更加符合實(shí)際通信場景的部分連接毫米波全雙工MIMO系統(tǒng),即有限精度移相器、寬帶毫米波系統(tǒng)和上下行信道的實(shí)際路徑損耗.

    2)提出了一種可行的三階段混合波束成形方法.首先,不考慮自干擾信號(hào),利用上行與下行各子信道的奇異值分解和注水功率分配算法得到各子帶最優(yōu)全數(shù)字發(fā)送和接收波束成形矩陣;然后,根據(jù)得到的全數(shù)字波束成形矩陣,利用矩陣分解方法求出滿足部分連接結(jié)構(gòu)的模擬波束成形矩陣和各子帶數(shù)字波束成形矩陣;最后,提出利用最小均方誤差(Minimum Mean Square Error,MMSE)準(zhǔn)則更新各子帶數(shù)字接收波束成形矩陣來處理剩余自干擾信號(hào),降低其引起的頻譜效率損失.

    3)理論分析表明,所提方法的計(jì)算復(fù)雜度要低于相關(guān)算法;同時(shí),仿真結(jié)果表明,所提方法的頻譜效率要高于相關(guān)對(duì)比方法和理想半雙工毫米波通信.當(dāng)移相器精度為6位、信噪比為20 dB時(shí),采用所提方法的部分連接寬帶毫米波全雙工MIMO系統(tǒng)實(shí)現(xiàn)的頻譜效率大約是理想半雙工毫米波系統(tǒng)的1.41倍.

    本文所用數(shù)學(xué)符號(hào)介紹如下:粗體小寫和大寫字母分別表示向量和矩陣,AH,A-1,A?,Tr{A}和‖A‖F(xiàn)分別表示矩陣A的共軛轉(zhuǎn)置、逆、廣義逆、跡和Frobenius范數(shù),A(i,:),A(:,j)和A(i,j)分別表示矩陣A的第i行、第j列和位于第i行、第j列的元素.IN表示N×N維單位矩陣,0N(0M×N)表示N×N(M×N)維零矩陣,M×N表示M×N維復(fù)數(shù)矩陣構(gòu)成的線性空間.CN(μ,Σ)表示服從均值為μ、方差為Σ的循環(huán)對(duì)稱復(fù)高斯隨機(jī)向量的分布,N(μ,σ2)表示服從均值為μ、方差為σ2的實(shí)高斯隨機(jī)變量的分布.blkdiag{a1,a2,…,aN}表示由向量a1,a2,…,aN構(gòu)成的塊狀對(duì)角矩陣,|·|表示標(biāo)量的絕對(duì)值或矩陣的行列式,「·表示向上取整即取超過參數(shù)的最小整數(shù),{·}和Re{·}分別表示取參數(shù)的期望和實(shí)部.

    1 系統(tǒng)模型與問題描述

    1.1 系統(tǒng)模型

    本文考慮單蜂窩寬帶毫米波全雙工MIMO系統(tǒng),包括一個(gè)全雙工基站BS和兩個(gè)半雙工移動(dòng)端MS1和MS2,如圖1所示.其中,基站和移動(dòng)端都采用部分連接結(jié)構(gòu),寬帶信號(hào)采用正交頻分復(fù)用(Orthogonal Frequency Division Multiplexing,OFDM)實(shí)現(xiàn).全雙工基站BS配備NT根發(fā)送天線、NR根接收天線,發(fā)送和接收射頻鏈路的數(shù)目均為NRF;移動(dòng)端MS1配備MT根發(fā)送天線和MRF條發(fā)送射頻鏈路,移動(dòng)端MS2配備MR根接收天線和MRF條接收射頻鏈路.OFDM子載波的數(shù)目為K,每個(gè)子載波的平行數(shù)據(jù)流數(shù)目均為Ns.由毫米波混合波束成形結(jié)構(gòu)中數(shù)據(jù)流、射頻鏈路和天線數(shù)目之間的關(guān)系可知,Ns≤NRF≤min{NT,NR},Ns≤MRF≤min{MT,MR}.正常工作時(shí),全雙工基站BS在上行鏈路接收移動(dòng)端MS1的發(fā)送信號(hào),同時(shí)在下行鏈路向移動(dòng)端MS2發(fā)送信號(hào).由于毫米波信道的路徑損耗嚴(yán)重,為了處理簡單,本文假設(shè)兩移動(dòng)端MS1和MS2之間的距離足夠遠(yuǎn),因此,可以忽略二者之間的相互干擾.

    圖1 部分連接寬帶毫米波全雙工MIMO系統(tǒng)Fig.1 Wideband mmWave full-duplex MIMO system with partially-connected structures

    在上行鏈路中,移動(dòng)端MS1第k個(gè)子帶的發(fā)送符號(hào)向量為sU[k]∈Ns×1,滿足首先,發(fā)送符號(hào)向量sU[k]經(jīng)過該子帶的數(shù)字發(fā)送波束成形矩陣FBBU[k]∈MRF×Ns預(yù)編碼,再通過K點(diǎn)快速傅里葉逆變換(Inverse Fast Fourier Transform,IFFT)變?yōu)闀r(shí)域信號(hào),并增加循環(huán)前綴(Cyclic Prefix,CP),以保證子載波間的正交性同時(shí)消除符號(hào)間干擾;然后,該信號(hào)經(jīng)發(fā)送射頻鏈路上變頻至毫米波頻段,通過模擬發(fā)送波束成形矩陣FRFU∈MT×MRF處理之后,得到移動(dòng)端MS1的最終發(fā)送信號(hào);最后,該發(fā)送信號(hào)經(jīng)過上行信道發(fā)送至基站BS.同理,在下行鏈路中,基站第k個(gè)子帶的發(fā)送符號(hào)向量sD[k]∈Ns×1,滿足經(jīng)過該子帶的數(shù)字發(fā)送波束成形矩陣FBBD[k]∈NRF×Ns、K點(diǎn)IFFT增加CP以及模擬發(fā)送波束成矩陣FRFD∈NT×NRF處理之后,得到基站BS的最終發(fā)送信號(hào),該信號(hào)經(jīng)過下行信道發(fā)送至移動(dòng)端MS2.

    由于基站BS工作于全雙工模式,因此,不但能夠通過上行信道接收來自移動(dòng)端MS1的信號(hào),而且能夠通過自干擾信道接收自干擾信號(hào).該接收信號(hào)首先經(jīng)過模擬接收波束成形矩陣WRFU∈NR×NRF處理,再由接收射頻鏈路下變頻至基帶;然后,通過移除CP并進(jìn)行快速傅里葉變換(Fast Fourier Transform,FFT)變?yōu)轭l域信號(hào),最后,經(jīng)過數(shù)字接收波束成形矩陣WBBU[k]∈NRF×Ns,k=1,2,…,K,處理之后,得到每個(gè)子帶的最終基帶接收信號(hào).基站BS第k個(gè)子帶的基帶接收信號(hào)可以具體表示為:

    (1)

    其中,HU[k]∈NR×MT和HSI[k]∈NR×NT分別表示第k個(gè)子帶的上行信道和自干擾信道,表示基站BS接收端第k個(gè)子帶的加性高斯白噪聲.WU[k]=WRFUWBBU[k]表示基站接收端第k個(gè)子帶的整體混合接收波束成形矩陣,FU[k]=FRFUFBBU[k]和FD[k]=FRFDFBBD[k]則分別表示移動(dòng)端MS1和基站BS發(fā)送端第k個(gè)子帶的整體混合發(fā)送波束成形矩陣.在基站BS第k個(gè)子帶的接收信號(hào)即公式(1)中,第1項(xiàng)表示從移動(dòng)端MS1接收到的有效信號(hào),第2項(xiàng)表示基帶噪聲信號(hào),第3項(xiàng)表示自干擾信號(hào).同理,半雙工移動(dòng)端MS2的接收信號(hào)經(jīng)過模擬接收波束成形矩陣WRFD∈MR×MRF、移除CP和K點(diǎn)FFT以及各子帶數(shù)字接收波束成形矩陣WBBD[k]∈MRF×Ns,k=1,2,…,K,處理之后,得到第k個(gè)子帶的基帶接收信號(hào)為:

    (2)

    其中,HD[k]∈MR×NT表示第k個(gè)子帶的下行信道,表示移動(dòng)端MS2第k個(gè)子帶的加性高斯白噪聲,WD[k]=WRFDWBBD[k]則表示移動(dòng)端MS2第k個(gè)子帶的整體混合接收波束成形矩陣.

    1.2 信道模型

    由于基站BS與移動(dòng)端MS1(MS2)之間的距離足夠遠(yuǎn),所以,上行和下行信道充分滿足遠(yuǎn)場條件.為了準(zhǔn)確描述毫米波MIMO-OFDM信道特征,本文采用抽頭延遲線寬帶信道模型[25],其中每個(gè)抽頭信道又建模為Saleh-Valenzuela(SV)空間信道模型[26,27].具體地,第d個(gè)抽頭的信道模型如下:

    (3)

    (4)

    其中,D表示CP的長度.因此,可以根據(jù)上行與下行鏈路中發(fā)送和接收天線數(shù)目、散射簇和路徑數(shù)目以及到達(dá)角、離開角等來更改式(3)中的參數(shù),建模各抽頭上行和下行信道;然后,根據(jù)式(4)得到上行與下行各子帶信道.

    由于全雙工基站BS發(fā)送和接收天線陣列距離較近,無法充分滿足遠(yuǎn)場條件,所以自干擾信道模型必須將近場效應(yīng)考慮在內(nèi),通常建模為萊斯衰減信道[11,21,24],其第d個(gè)抽頭的表達(dá)式為:

    (5)

    (6)

    其中,ρ為歸一化因子,保證為載波波長,rmn表示基站BS第m根接收天線與第n根發(fā)送天線之間的距離,可以根據(jù)文獻(xiàn)[11]中的圖2,利用余弦定理求出.表示自干擾信道第d個(gè)抽頭的非視距傳播部分,滿足遠(yuǎn)場條件,所以,可用式(3)表示的SV空間信道模型進(jìn)行建模.對(duì)時(shí)域抽頭自干擾信道進(jìn)行DFT,如式(4)所示,可得頻域各子帶自干擾信道.

    圖2 第2階段混合波束成形的收斂性Fig.2 Convergence of the hybrid beamforming design in the second stage

    1.3 部分連接結(jié)構(gòu)

    1.4 問題描述

    假設(shè)每個(gè)子帶的發(fā)送信號(hào)均為高斯信號(hào),則第k個(gè)子帶的上行頻譜效率RU[k]和下行頻譜效率RD[k]分別定義為:

    本文假設(shè)上行、下行和自干擾信道的狀態(tài)信息在基站BS已知且接收射頻鏈路模數(shù)轉(zhuǎn)換器(Analog-to-Digital Converters,ADCs)的動(dòng)態(tài)范圍無限大[19-21],考慮發(fā)送功率約束、模擬波束成形矩陣非零元素的有限精度約束、單位模長約束以及由部分連接結(jié)構(gòu)引入的約束,旨在通過聯(lián)合設(shè)計(jì)全雙工基站BS和半雙工移動(dòng)端MS1、MS2的混合波束成形矩陣來最大化系統(tǒng)的整體頻譜效率.該問題的數(shù)學(xué)表述如下:

    (7a)

    (7b)

    (7c)

    FRFU(i,j)∈D,|FRFU(i,j)|=1,?(i,j)∈FU

    (7d)

    FRFD(i,j)∈D,|FRFD(i,j)|=1,?(i,j)∈FD

    (7e)

    WRFU(i,j)∈D,|WRFU(i,j)|=1,?(i,j)∈WU

    (7f)

    WRFD(i,j)∈D,|WRFD(i,j)|=1,?(i,j)∈WD

    (7g)

    2 三階段混合波束成形方法

    針對(duì)非凸優(yōu)化問題(7),本文提出了一種可行的混合波束成形方法.雖然所提方法并不能求出最優(yōu)解,但可以求出一組考慮系統(tǒng)性能的可行解.所提混合波束成形方法共包括3個(gè)階段,詳述如下.

    2.1 不考慮自干擾信號(hào)的最優(yōu)全數(shù)字波束成形

    (8a)

    (8b)

    其中,PU[k]∈Ns×Ns為對(duì)角矩陣,每個(gè)對(duì)角線元素表示在移動(dòng)端MS1通過注水算法分配給第k個(gè)子帶中對(duì)應(yīng)數(shù)據(jù)流的發(fā)送功率,滿足同理,若下行第k個(gè)子信道HD[k]的奇異值分解為其中DD[k]∈MR×NT表示對(duì)角矩陣,其對(duì)角線元素為遞減排列的奇異值,UD[k]∈MR×MR和VD[k]∈NT×NT分別表示由左奇異向量和右奇異向量構(gòu)成的酉陣,則基站BS發(fā)送端第k個(gè)子帶的最優(yōu)全數(shù)字發(fā)送波束成形矩陣和移動(dòng)端MS2第k個(gè)子帶的最優(yōu)全數(shù)字接收波束成形矩陣分別為:

    (9a)

    (9b)

    其中PD[k]∈Ns×Ns為對(duì)角矩陣,其對(duì)角線元素表示在基站BS發(fā)送端利用注水算法分配給第k個(gè)子帶中對(duì)應(yīng)數(shù)據(jù)流的發(fā)送功率,滿足

    2.2 不考慮自干擾信號(hào)的混合波束成形

    (10a)

    (10b)

    FRFD(i,j)∈D,|FRFD(i,j)|=1,?(i,j)∈FD

    (10c)

    2.2.1 更新模擬發(fā)送波束成形矩陣FRFD

    當(dāng)固定其他變量同時(shí)移除常數(shù)項(xiàng)時(shí),問題(10)關(guān)于模擬發(fā)送波束成形矩陣FRFD的子問題為:

    (11)

    s.t.(10c)

    (12)

    s.t.(10c)

    (13)

    由于各子帶的數(shù)字發(fā)送波束成形矩陣相互獨(dú)立,所以,當(dāng)其他變量固定時(shí),問題(10)關(guān)于第k個(gè)子帶的數(shù)字發(fā)送波束成形矩陣FBBD[k]的子問題為:

    (14)

    (15)

    同理,可求出其他各子帶數(shù)字發(fā)送波束成形矩陣.交替迭代更新模擬和各子帶數(shù)字發(fā)送波束成形矩陣,直至滿足某一終止條件.最后,對(duì)數(shù)字發(fā)送波束成形矩陣FBBD[k],k=1,2,…,K,進(jìn)行如下功率歸一化,以滿足發(fā)送功率約束:

    (16)

    2.3 MMSE準(zhǔn)則更新各子帶數(shù)字接收波束成形矩陣

    前兩階段求得的混合波束成形設(shè)計(jì)并沒有考慮自干擾信號(hào),所以,剩余自干擾信號(hào)存在,并且會(huì)造成較大的頻譜效率損失.在考慮復(fù)雜度前提下,線性MMSE接收機(jī)是實(shí)現(xiàn)系統(tǒng)最大傳輸速率的最優(yōu)數(shù)字接收機(jī)[28,31].因此,為了有效降低自干擾信號(hào)對(duì)頻譜效率的影響,本階段提出在已經(jīng)求得的混合波束成形設(shè)計(jì)的基礎(chǔ)上,利用MMSE準(zhǔn)則來更新基站BS接收端和移動(dòng)端MS2各子帶的數(shù)字接收波束成形矩陣.首先,基站BS接收端和移動(dòng)端MS2第k個(gè)子帶的均方差矩陣EU[k]和ED[k]分別定義如下:

    EU[k]=

    ED[k]=

    根據(jù)上述定義,當(dāng)其他變量保持不變同時(shí)移除常數(shù)項(xiàng)時(shí),關(guān)于數(shù)字接收波束成形矩陣WBBU[k]的最小均方誤差問題如下:

    (17)

    其中,

    (18)

    至此,針對(duì)部分連接寬帶毫米波全雙工MIMO系統(tǒng),本文所提出的三階段混合波束成形設(shè)計(jì)完成,其中,第2階段不考慮自干擾信號(hào)的混合波束成形設(shè)計(jì)總結(jié)于算法1.由于第3階段采用MMSE準(zhǔn)則更新各子帶數(shù)字接收波束成形矩陣,所以,將本文提出的混合波束成形方法命名為WBFD_PCS_MMSE,其具體細(xì)節(jié)總結(jié)于算法2.

    算法1.第2階段不考慮自干擾信號(hào)的混合波束成形設(shè)計(jì)

    1.初始化XRF∈nant×nrf為塊狀對(duì)角矩陣,非零元素滿足有限精度和單位模長約束,設(shè)置最大迭代次數(shù)L和誤差精度δ,l=0;

    2.nsub=nant/nrf

    3.repeat

    4.i=0;

    5.repeat

    7.until(i>nant)

    11.until(|Δ1-Δ2|<δ或l>L)

    算法2.部分連接寬帶毫米波全雙工MIMO系統(tǒng)三階段混合波束成形方法WBFD_PCS_MMSE

    3 算法復(fù)雜度分析

    本節(jié)以所需復(fù)數(shù)乘法的次數(shù)作為算法復(fù)雜度的量度,對(duì)本文提出的部分連接寬帶毫米波全雙工MIMO系統(tǒng)混合波束成形方法WBFD_PCS_MMSE的算法復(fù)雜度進(jìn)行分析.所提方法包括3階段,各階段算法復(fù)雜度具體分析如下:

    (19)

    其中,(a)成立的條件是假設(shè)基站BS發(fā)送和接收天線的數(shù)目要遠(yuǎn)大于移動(dòng)端MS1和MS2的天線數(shù)目,同時(shí)基站BS的射頻鏈路數(shù)目大于等于移動(dòng)端MS1和MS2的射頻鏈路數(shù)目.

    4 仿真結(jié)果與分析

    本節(jié)將利用蒙特卡洛仿真方法對(duì)提出的部分連接寬帶毫米波全雙工MIMO系統(tǒng)混合波束成形方法WBFD_PCS_MMSE的性能進(jìn)行評(píng)估.為了更好地展現(xiàn)所提方法的性能,考慮如下對(duì)比:

    1)Ideal FD:理想全雙工系統(tǒng),假設(shè)沒有自干擾信號(hào)同時(shí)采用全數(shù)字波束成形設(shè)計(jì);

    2)FD_BFC:文獻(xiàn)[21]中提出的基于自干擾抑制的全連接寬帶毫米波全雙工系統(tǒng)混合波束成形設(shè)計(jì);

    3)FD_MinSIP:將文獻(xiàn)[19]中提出的基于自干擾信號(hào)功率最小化的混合波束成形方法擴(kuò)展至全連接寬帶毫米波全雙工系統(tǒng);

    4)FD_PCS_JBIC:文獻(xiàn)[24]中提出的部分連接寬帶毫米波全雙工系統(tǒng)聯(lián)合優(yōu)化接收功率和自干擾抑制的混合波束成形方法;

    5)FD_PCS_SDRNSP:將文獻(xiàn)[23]中提出的部分連接毫米波全雙工系統(tǒng)的自干擾抑制混合波束成形方法擴(kuò)展到對(duì)應(yīng)的部分連接寬帶毫米波全雙工MIMO系統(tǒng);

    6)Ideal HD:理想半雙工系統(tǒng),基站BS所有天線均用于半雙工通信,同時(shí)采用全數(shù)字波束成形設(shè)計(jì).

    4.1 收斂性

    4.2 頻譜效率

    圖3比較了不同方法的頻譜效率隨信噪比的變化,其中,WBFD_PCS_NoMMSE表示只采用所提混合波束成形設(shè)計(jì)的前兩階段所實(shí)現(xiàn)的頻譜效率.從圖中可以看出,本文提出的三階段混合波束成形方法WBFD_PCS_MMSE的頻譜效率要明顯優(yōu)于WBFD_PCS_NoMMSE,說明所提方法的第3階段必不可少,該階段可以有效抑制剩余自干擾信號(hào)、提升頻譜效率.從圖中也可以觀察到,本文所提方法WBFD_PCS_MMSE的頻譜效率要比方法WBFD_BFC差,主要原因是后者為全連接毫米波全雙工系統(tǒng),能夠提供更好的波束成形增益;雖然方法WBFD_MinSIP也考慮全連接毫米波全雙工系統(tǒng),但是由于沒有對(duì)矩陣分解引起的剩余自干擾信號(hào)進(jìn)行處理,所以該方法的頻譜效率要比本文所提方法差.同時(shí),可以看到本文所提方法實(shí)現(xiàn)的頻譜效率要優(yōu)于WBFD_PCS_JBIC和WBFD_PCS_SDRNSP這兩種針對(duì)部分連接寬帶毫米波全雙工MIMO系統(tǒng)的混合波束成形方法.通過第4節(jié)分析可知,這兩種對(duì)比方法的實(shí)際計(jì)算復(fù)雜度要高于本文所提方法,因此本文所提方法WBFD_PCS_MMSE具有一定的優(yōu)越性.此外,當(dāng)信噪比高于0dB時(shí),本文所提方法的頻譜效率要高于理想半雙工系統(tǒng),證明了部分連接寬帶毫米波全雙工MIMO通信的可行性;當(dāng)信噪比為20dB時(shí),本文所提方法的頻譜效率大約是理想半雙工通信的1.41倍.

    圖3 頻譜效率隨信噪比的變化Fig.3 Spectral efficiency versussignal-to-noise ratio

    圖4對(duì)比了不同方法的頻譜效率隨基站BS射頻鏈路數(shù)目的變化,其中信噪比等于10dB.對(duì)于部分連接結(jié)構(gòu),為了保證天線數(shù)目可以被射頻鏈路數(shù)目整除,基站BS射頻鏈路數(shù)目NRF={4,8,16}.隨著基站BS射頻鏈路數(shù)目的增加,模擬和數(shù)字波束成形的自由度增加,可以更好地抑制自干擾信號(hào)同時(shí)提高波束成形增益,因此,圖4中采用本文所提方法WBFD_PCS_MMSE的部分連接寬帶毫米波全雙工MIMO系統(tǒng)的頻譜效率上升.同時(shí),在任何射頻鏈路數(shù)目條件下,本文所提方法WBFD_PCS_MMSE的頻譜效率都要優(yōu)于考慮相同系統(tǒng)的對(duì)比方法WBFD_PCS_JBIC和WBFD_PCS_SDRNSP.當(dāng)射頻鏈路數(shù)目較少時(shí),如NRF=4,由于仿真實(shí)驗(yàn)中的平行數(shù)據(jù)流數(shù)目Ns=3,所以,沒有充足的自由度同時(shí)進(jìn)行有效自干擾抑制和數(shù)據(jù)流傳輸.因此,本文所提方法WBFD_PCS_MMSE的頻譜效率要低于基于方法WBFD_MinSIP的全連接毫米波全雙工系統(tǒng),甚至要低于理想半雙工系統(tǒng).當(dāng)射頻鏈路數(shù)目大于等于6時(shí),所提方法的頻譜效率要高于方法WBFD_MinSIP和理想半雙工系統(tǒng),但是仍低于方法WBFD_BFC,具體原因與圖3類似,這里不再贅述.

    圖4 頻譜效率隨基站射頻鏈路數(shù)目的變化Fig.4 Spectral efficiency versus the number ofradio-frequency chains at the base station

    圖5描述了不同方法的頻譜效率隨基站BS天線數(shù)目的變化,其中信噪比等于10dB.圖中橫坐標(biāo)表示基站BS總的天線數(shù)目,其中一半是發(fā)送天線、另一半是接收天線.從圖中可以觀察到,隨著天線數(shù)目的增加,本文所提方法WBFD_PCS_MMSE的頻譜效率上升,這是因?yàn)楦嗟奶炀€能夠提供更高的天線多樣性和自由度,從而可以更好地抑制自干擾信號(hào),同時(shí)提高波束成形增益.相比于考慮相同系統(tǒng)的混合波束成形方法WBFD_PCS_JBIC和WBFD_PCS_SDRNSP,本文所提方法WBFD_PCS_MMSE能夠?qū)崿F(xiàn)更高的頻譜效率,但與方法WBFD_PCS_JBIC的頻譜效率差距逐漸減小.本文所提方法的頻譜效率比基于方法WBFD_BFC的全連接毫米波全雙工系統(tǒng)差,但是二者之間的頻譜效率差距先上升后下降,在天線數(shù)目為128時(shí)最大,大約為9.1 bits/s/Hz.同時(shí),本文所提方法的頻譜效率要明顯高于采用方法WBFD_MinSIP的全連接毫米波全雙工系統(tǒng)和理想半雙工系統(tǒng).

    圖5 頻譜效率隨基站天線數(shù)目的變化Fig.5 Spectral efficiency versus the number of antennas at the base station

    圖6提供了不同方法的頻譜效率隨子載波數(shù)目的變化,其中信噪比等于10dB.由圖可知,隨著子載波數(shù)目的增加,本文所提方法WBFD_PCS_MMSE的頻譜效率略微下降.這是因?yàn)樽虞d波數(shù)目的增加會(huì)引起信道的頻率選擇性上升,從而更難設(shè)計(jì)出滿足所有子載波信道的模擬發(fā)送和接收波束成形矩陣,在一定程度上降低了自干擾抑制能力和波束成形增益,造成頻譜效率下降.此外,從圖6中可以得到與圖3和圖5中類似的結(jié)論,即本文所提方法WBFD_PCS_MMSE的頻譜效率比方法WBFD_BFC差,但是要明顯優(yōu)于考慮相同系統(tǒng)的對(duì)比方法WBFD_PCS_JBIC和WBFD_PCS_SDRNSP,以及WBFD_MinSIP和理想半雙工系統(tǒng).

    圖6 頻譜效率隨子載波數(shù)目的變化Fig.6 Spectral efficiency versus the number of subcarriers

    圖7仿真了本文所提方法WBFD_PCS_MMSE采用不同精度移相器實(shí)現(xiàn)的頻譜效率隨信噪比的變化.從圖中可以看出,隨著移相器精度的提高,所提方法的頻譜效率先上升后逐漸飽和,原因是高精度的移相器可以提供更豐富的碼本,從而實(shí)現(xiàn)更優(yōu)的模擬波束成形.當(dāng)信噪比高于10 dB時(shí),所提方法即使采用1位精度移相器實(shí)現(xiàn)的頻譜效率也要高于理想半雙工系統(tǒng),這與文獻(xiàn)[22]中針對(duì)全連接毫米波全雙工系統(tǒng)的結(jié)論一致;當(dāng)移相器精度超過4位時(shí),所提方法能夠?qū)崿F(xiàn)接近于無窮精度移相器的頻譜效率,再提高移相器精度對(duì)頻譜效率的提升不明顯.因此,權(quán)衡成本、功耗和性能,對(duì)于采用本文所提方法的部分連接寬帶毫米波全雙工MIMO系統(tǒng)建議使用4位精度移相器.當(dāng)信噪比等于20 dB時(shí),所提方法采用4位精度移相器可以實(shí)現(xiàn)理想全雙工系統(tǒng)大約76.2%的頻譜效率,同時(shí)大約是理想半雙工系統(tǒng)頻譜效率的1.41倍.

    圖7 不同精度移相器的頻譜效率隨信噪比的變化Fig.7 Spectral efficiency using phase shifters with different resolutions versus SNR

    5 結(jié) 論

    本文對(duì)采用有限精度移相器的部分連接寬帶毫米波全雙工MIMO系統(tǒng)的混合波束成形問題進(jìn)行研究,提出了一種可行的三階段混合波束成形方法,該方法可以有效抑制自干擾信號(hào)同時(shí)實(shí)現(xiàn)較優(yōu)的頻譜效率.分析表明,所提方法的算法復(fù)雜度要低于相關(guān)對(duì)比方法;仿真結(jié)果表明,所提方法具有較好的收斂性,能夠?qū)崿F(xiàn)比相關(guān)研究和理想半雙工毫米波通信更高的頻譜效率.

    在后續(xù)工作中,將進(jìn)一步對(duì)考慮ADC有限動(dòng)態(tài)范圍和用戶間相互干擾的部分連接寬帶毫米波全雙工MIMO系統(tǒng)的混合波束成形問題進(jìn)行研究.

    猜你喜歡
    全雙工干擾信號(hào)波束
    對(duì)抗全雙工主動(dòng)竊聽的安全高效D2D通信策略
    無線電工程(2024年5期)2024-07-20 00:00:00
    正弦采樣信號(hào)中單一脈沖干擾信號(hào)的快速剔除實(shí)踐方法
    基于粒子群算法的光纖通信干擾信號(hào)定位方法
    毫米波大規(guī)模陣列天線波束掃描研究*
    圓陣多波束測角探究
    電子測試(2018年6期)2018-05-09 07:31:54
    RLS算法在同時(shí)同頻全雙工系統(tǒng)中的自干擾抑制性能分析
    Helix陣匹配場三維波束形成
    淺析監(jiān)控干擾信號(hào)的優(yōu)化處置措施
    同頻同時(shí)全雙工數(shù)字自干擾抵消算法
    5G全雙工技術(shù)淺析
    亚洲婷婷狠狠爱综合网| 成年av动漫网址| 春色校园在线视频观看| 自拍欧美九色日韩亚洲蝌蚪91 | 国产男人的电影天堂91| 99热国产这里只有精品6| 国产成人精品无人区| 久久久久国产网址| 看免费成人av毛片| 国产成人精品一,二区| 国产欧美日韩一区二区三区在线 | 国产精品久久久久久av不卡| 九色成人免费人妻av| 亚洲怡红院男人天堂| 日韩中字成人| 成人毛片a级毛片在线播放| .国产精品久久| 久久久久久久久久成人| 伦精品一区二区三区| 建设人人有责人人尽责人人享有的| 亚洲欧美中文字幕日韩二区| 精品国产国语对白av| 伦理电影免费视频| 久久人人爽人人爽人人片va| 草草在线视频免费看| 久久精品国产自在天天线| 纵有疾风起免费观看全集完整版| 亚洲国产日韩一区二区| 18+在线观看网站| 美女中出高潮动态图| 久久精品国产鲁丝片午夜精品| 亚洲av综合色区一区| 水蜜桃什么品种好| 六月丁香七月| 国产亚洲最大av| 国产成人精品无人区| 欧美区成人在线视频| 久久国产精品男人的天堂亚洲 | 国产视频内射| av女优亚洲男人天堂| 日韩免费高清中文字幕av| 亚洲,一卡二卡三卡| 国产精品女同一区二区软件| 日韩中文字幕视频在线看片| 99热网站在线观看| 成人国产麻豆网| 欧美bdsm另类| 观看av在线不卡| 亚洲精品亚洲一区二区| 美女大奶头黄色视频| 成年美女黄网站色视频大全免费 | 国产精品99久久久久久久久| 两个人免费观看高清视频 | 久久国产乱子免费精品| 一级,二级,三级黄色视频| 永久网站在线| 永久网站在线| 精品国产一区二区三区久久久樱花| 一个人看视频在线观看www免费| 新久久久久国产一级毛片| 亚洲图色成人| 国产欧美日韩一区二区三区在线 | 国产精品一区www在线观看| 日韩亚洲欧美综合| 一本久久精品| 免费大片黄手机在线观看| 亚洲精品中文字幕在线视频 | 丁香六月天网| 中文在线观看免费www的网站| 自线自在国产av| 18禁动态无遮挡网站| 久久久久久久大尺度免费视频| 啦啦啦啦在线视频资源| 晚上一个人看的免费电影| 啦啦啦啦在线视频资源| 国产精品伦人一区二区| 亚洲精品中文字幕在线视频 | 丝袜喷水一区| 极品少妇高潮喷水抽搐| 新久久久久国产一级毛片| 亚洲美女搞黄在线观看| 黑人巨大精品欧美一区二区蜜桃 | 欧美亚洲 丝袜 人妻 在线| 亚洲精品视频女| 久久免费观看电影| 91午夜精品亚洲一区二区三区| 亚洲av欧美aⅴ国产| 91午夜精品亚洲一区二区三区| 国产亚洲5aaaaa淫片| 日本与韩国留学比较| 又大又黄又爽视频免费| 亚洲国产日韩一区二区| 免费av不卡在线播放| 黄色一级大片看看| 黄色日韩在线| 99精国产麻豆久久婷婷| 插逼视频在线观看| 色94色欧美一区二区| 国产高清不卡午夜福利| 精品国产露脸久久av麻豆| 成年美女黄网站色视频大全免费 | 麻豆成人av视频| 亚洲在久久综合| 女性生殖器流出的白浆| 伦精品一区二区三区| 女人久久www免费人成看片| 国产高清不卡午夜福利| 制服丝袜香蕉在线| 伊人亚洲综合成人网| 日产精品乱码卡一卡2卡三| 如日韩欧美国产精品一区二区三区 | 国产黄色视频一区二区在线观看| 2021少妇久久久久久久久久久| 一本色道久久久久久精品综合| 日日撸夜夜添| 五月玫瑰六月丁香| 国产精品免费大片| 久久人人爽av亚洲精品天堂| 精品一区在线观看国产| 亚洲精品国产成人久久av| 久久久久人妻精品一区果冻| 国产精品女同一区二区软件| 99九九线精品视频在线观看视频| 美女xxoo啪啪120秒动态图| 成人影院久久| 9色porny在线观看| 成人漫画全彩无遮挡| 观看美女的网站| 精品少妇内射三级| 一本大道久久a久久精品| 久久国内精品自在自线图片| 在线 av 中文字幕| 亚洲人与动物交配视频| 另类精品久久| 成年av动漫网址| 97在线人人人人妻| 欧美97在线视频| h视频一区二区三区| 久久久欧美国产精品| 欧美高清成人免费视频www| 两个人免费观看高清视频 | 我的老师免费观看完整版| 国产精品一区二区性色av| 久久99蜜桃精品久久| 久久精品久久久久久噜噜老黄| 日韩制服骚丝袜av| 久久久国产一区二区| 国产精品一区二区在线不卡| 嫩草影院新地址| 一级毛片我不卡| 少妇 在线观看| 纯流量卡能插随身wifi吗| 成年av动漫网址| 老女人水多毛片| 欧美一级a爱片免费观看看| 日韩视频在线欧美| 18+在线观看网站| 国产探花极品一区二区| 五月开心婷婷网| 久久久久久久久久久免费av| 最近中文字幕2019免费版| 日本黄大片高清| 亚洲精品视频女| 日韩中字成人| 久久久久久久久大av| xxx大片免费视频| 亚洲国产精品一区三区| 欧美日韩视频高清一区二区三区二| 人人妻人人添人人爽欧美一区卜| 一本大道久久a久久精品| 建设人人有责人人尽责人人享有的| 欧美成人精品欧美一级黄| 中文字幕人妻熟人妻熟丝袜美| 精品一区二区三卡| 18+在线观看网站| 人妻制服诱惑在线中文字幕| 欧美三级亚洲精品| 看非洲黑人一级黄片| 18禁裸乳无遮挡动漫免费视频| 亚洲欧洲日产国产| 国产av国产精品国产| 一本大道久久a久久精品| 国产中年淑女户外野战色| 国产精品伦人一区二区| 免费黄网站久久成人精品| h视频一区二区三区| 日本黄色日本黄色录像| 日本免费在线观看一区| 精品人妻偷拍中文字幕| 一区在线观看完整版| 大码成人一级视频| 国产综合精华液| 啦啦啦中文免费视频观看日本| 亚洲婷婷狠狠爱综合网| 黄色视频在线播放观看不卡| 香蕉精品网在线| 国产精品麻豆人妻色哟哟久久| 男女无遮挡免费网站观看| 日本黄色日本黄色录像| 久久人妻熟女aⅴ| 91aial.com中文字幕在线观看| 一级毛片 在线播放| 亚洲欧洲精品一区二区精品久久久 | av不卡在线播放| a级毛色黄片| 欧美+日韩+精品| 国内少妇人妻偷人精品xxx网站| 少妇丰满av| 国产高清国产精品国产三级| 国语对白做爰xxxⅹ性视频网站| 三级国产精品片| 国产成人精品一,二区| 免费观看在线日韩| 久久久久精品久久久久真实原创| 黑人猛操日本美女一级片| 大香蕉97超碰在线| 国产91av在线免费观看| 国产在线男女| 一级毛片黄色毛片免费观看视频| 久久99蜜桃精品久久| 十分钟在线观看高清视频www | 狂野欧美激情性bbbbbb| 国产一区有黄有色的免费视频| www.av在线官网国产| 最新的欧美精品一区二区| 我的女老师完整版在线观看| 国产日韩一区二区三区精品不卡 | 亚洲精品乱码久久久久久按摩| 九九久久精品国产亚洲av麻豆| 黑人猛操日本美女一级片| 免费观看无遮挡的男女| 久久久国产一区二区| 国产高清不卡午夜福利| 成人午夜精彩视频在线观看| 蜜桃在线观看..| 亚洲无线观看免费| 久久av网站| av有码第一页| av.在线天堂| 亚洲av日韩在线播放| 国产黄频视频在线观看| 久久精品国产亚洲av涩爱| 国产有黄有色有爽视频| 亚洲精品aⅴ在线观看| 男人舔奶头视频| 日本午夜av视频| 69精品国产乱码久久久| 成人黄色视频免费在线看| 国产 一区精品| 人人妻人人澡人人看| 各种免费的搞黄视频| 中文精品一卡2卡3卡4更新| 成人影院久久| 精品卡一卡二卡四卡免费| 99久久精品国产国产毛片| 久久人人爽人人片av| 最近的中文字幕免费完整| 免费黄频网站在线观看国产| 秋霞伦理黄片| 欧美人与善性xxx| 色网站视频免费| 亚洲激情五月婷婷啪啪| 亚洲美女黄色视频免费看| 亚洲欧洲国产日韩| 亚洲怡红院男人天堂| 久久99热这里只频精品6学生| 美女内射精品一级片tv| 曰老女人黄片| 亚洲欧美一区二区三区国产| 亚洲欧美日韩另类电影网站| 极品教师在线视频| 亚洲精品一区蜜桃| kizo精华| 国产精品一区二区性色av| 久久久午夜欧美精品| 亚洲伊人久久精品综合| 成人特级av手机在线观看| 在线观看免费高清a一片| 亚洲av综合色区一区| av在线app专区| 赤兔流量卡办理| 女性生殖器流出的白浆| 欧美精品人与动牲交sv欧美| 街头女战士在线观看网站| 亚洲av福利一区| 亚洲国产精品一区三区| 美女cb高潮喷水在线观看| 欧美丝袜亚洲另类| 久久精品国产a三级三级三级| 中文字幕亚洲精品专区| 免费大片18禁| 亚洲va在线va天堂va国产| 国产成人a∨麻豆精品| 少妇被粗大猛烈的视频| 国产精品久久久久久久电影| 天堂中文最新版在线下载| 中文天堂在线官网| 国产av码专区亚洲av| 22中文网久久字幕| 欧美日韩精品成人综合77777| 18禁裸乳无遮挡动漫免费视频| 精品少妇久久久久久888优播| 亚洲av综合色区一区| 91成人精品电影| 久久av网站| 精品酒店卫生间| 亚洲欧美日韩东京热| 国产精品不卡视频一区二区| 亚洲av中文av极速乱| 黑人猛操日本美女一级片| 一区二区三区乱码不卡18| 国产亚洲一区二区精品| av在线老鸭窝| 色94色欧美一区二区| 亚洲av成人精品一二三区| 2018国产大陆天天弄谢| 色视频在线一区二区三区| 国产乱来视频区| 伊人久久国产一区二区| 久久亚洲国产成人精品v| 五月伊人婷婷丁香| 热re99久久国产66热| 狠狠精品人妻久久久久久综合| 少妇被粗大的猛进出69影院 | 亚洲国产成人一精品久久久| 在线播放无遮挡| 国产av精品麻豆| 日本色播在线视频| 亚洲欧洲精品一区二区精品久久久 | 制服丝袜香蕉在线| 国产真实伦视频高清在线观看| a级毛片免费高清观看在线播放| 王馨瑶露胸无遮挡在线观看| 性色av一级| 午夜免费观看性视频| 久久久国产精品麻豆| 2022亚洲国产成人精品| 亚洲国产欧美在线一区| 免费观看性生交大片5| 国产老妇伦熟女老妇高清| 亚洲av在线观看美女高潮| 免费大片黄手机在线观看| 亚洲精品亚洲一区二区| 亚洲伊人久久精品综合| 一级a做视频免费观看| 亚洲欧美日韩另类电影网站| 中文在线观看免费www的网站| 亚洲无线观看免费| 久久99一区二区三区| 免费久久久久久久精品成人欧美视频 | 男的添女的下面高潮视频| 成人国产麻豆网| 噜噜噜噜噜久久久久久91| 午夜91福利影院| av在线app专区| 日韩,欧美,国产一区二区三区| 国产女主播在线喷水免费视频网站| 中文欧美无线码| 日本猛色少妇xxxxx猛交久久| 亚洲va在线va天堂va国产| 极品少妇高潮喷水抽搐| 久久久久国产网址| 国产精品成人在线| 国产日韩一区二区三区精品不卡 | 国产黄色免费在线视频| 亚洲av.av天堂| 黄色毛片三级朝国网站 | 肉色欧美久久久久久久蜜桃| 91aial.com中文字幕在线观看| 视频中文字幕在线观看| 午夜激情福利司机影院| 亚洲国产毛片av蜜桃av| 成年美女黄网站色视频大全免费 | 一级毛片黄色毛片免费观看视频| 国产一区二区在线观看日韩| 狂野欧美激情性xxxx在线观看| 建设人人有责人人尽责人人享有的| 青春草视频在线免费观看| 天堂中文最新版在线下载| 免费久久久久久久精品成人欧美视频 | 亚洲成色77777| 亚洲精品自拍成人| 一级a做视频免费观看| av不卡在线播放| 国产欧美日韩精品一区二区| 午夜福利网站1000一区二区三区| 日韩精品有码人妻一区| 黄色怎么调成土黄色| 国产美女午夜福利| 97超碰精品成人国产| 大话2 男鬼变身卡| 少妇人妻久久综合中文| 亚洲,一卡二卡三卡| 色视频www国产| 国产精品不卡视频一区二区| 最新的欧美精品一区二区| 激情五月婷婷亚洲| 亚洲精品国产av成人精品| 久久精品国产鲁丝片午夜精品| 免费看光身美女| 亚洲精品,欧美精品| 在线观看国产h片| 亚洲精品第二区| 日韩一区二区视频免费看| 免费大片黄手机在线观看| 男的添女的下面高潮视频| 国产男人的电影天堂91| 一本久久精品| 午夜精品国产一区二区电影| 欧美激情国产日韩精品一区| 肉色欧美久久久久久久蜜桃| 精品一区二区免费观看| 青青草视频在线视频观看| 青春草国产在线视频| 免费少妇av软件| 亚洲av国产av综合av卡| 一级片'在线观看视频| 欧美3d第一页| 人人妻人人爽人人添夜夜欢视频 | 日日撸夜夜添| 性色avwww在线观看| 成人二区视频| 内射极品少妇av片p| 久久热精品热| 在线观看三级黄色| 亚洲av免费高清在线观看| 男人和女人高潮做爰伦理| 欧美激情国产日韩精品一区| 99久久中文字幕三级久久日本| 亚洲高清免费不卡视频| 十八禁高潮呻吟视频 | 丝袜喷水一区| 全区人妻精品视频| 久久精品国产自在天天线| 91久久精品国产一区二区成人| 亚洲图色成人| 全区人妻精品视频| 狂野欧美激情性xxxx在线观看| 久久6这里有精品| 成人综合一区亚洲| 最近中文字幕2019免费版| 国产高清国产精品国产三级| 韩国高清视频一区二区三区| 97超碰精品成人国产| 黄片无遮挡物在线观看| 在线观看国产h片| 十八禁高潮呻吟视频 | 亚洲怡红院男人天堂| 男男h啪啪无遮挡| 伊人久久国产一区二区| av福利片在线观看| 99久久精品热视频| 亚洲国产最新在线播放| videos熟女内射| 欧美三级亚洲精品| 少妇的逼好多水| 国产美女午夜福利| 80岁老熟妇乱子伦牲交| 亚洲av欧美aⅴ国产| 久久人人爽人人爽人人片va| 一本久久精品| 狠狠精品人妻久久久久久综合| 国产精品国产av在线观看| 一级,二级,三级黄色视频| 日韩人妻高清精品专区| 国产精品秋霞免费鲁丝片| 午夜老司机福利剧场| 国产乱来视频区| 中文资源天堂在线| 久久久a久久爽久久v久久| 久久免费观看电影| 久久 成人 亚洲| 亚洲美女搞黄在线观看| 久久久国产精品麻豆| 视频区图区小说| 大片免费播放器 马上看| 中文字幕人妻熟人妻熟丝袜美| 全区人妻精品视频| 国产亚洲午夜精品一区二区久久| 国产精品久久久久久av不卡| 久热久热在线精品观看| a级片在线免费高清观看视频| 免费人妻精品一区二区三区视频| 久久久久人妻精品一区果冻| 国产亚洲一区二区精品| 伊人久久精品亚洲午夜| 亚洲欧美成人综合另类久久久| 热re99久久国产66热| 国产黄色视频一区二区在线观看| 99久久精品一区二区三区| 日本免费在线观看一区| 日韩精品有码人妻一区| 久久97久久精品| 国产精品一区www在线观看| 久久久欧美国产精品| 日韩av不卡免费在线播放| 少妇裸体淫交视频免费看高清| 韩国高清视频一区二区三区| 麻豆成人av视频| 久久免费观看电影| 精品久久国产蜜桃| 免费在线观看成人毛片| 97超碰精品成人国产| 免费播放大片免费观看视频在线观看| 欧美日韩国产mv在线观看视频| 日本欧美视频一区| 这个男人来自地球电影免费观看 | 大香蕉久久网| 亚洲av.av天堂| 国产精品一区www在线观看| 亚洲欧美成人精品一区二区| 亚洲天堂av无毛| 毛片一级片免费看久久久久| 久久久久久久久久久丰满| 久久人妻熟女aⅴ| 大码成人一级视频| 插阴视频在线观看视频| 美女中出高潮动态图| 免费人妻精品一区二区三区视频| 黄色毛片三级朝国网站 | 精品人妻一区二区三区麻豆| 日本黄大片高清| 高清av免费在线| 亚洲欧美一区二区三区国产| av专区在线播放| 自拍偷自拍亚洲精品老妇| 国产成人a∨麻豆精品| 亚洲精品国产av蜜桃| av.在线天堂| 最近中文字幕高清免费大全6| 国产日韩欧美亚洲二区| 伊人亚洲综合成人网| 18禁在线播放成人免费| 亚洲精品中文字幕在线视频 | 80岁老熟妇乱子伦牲交| 汤姆久久久久久久影院中文字幕| 国产视频首页在线观看| 极品教师在线视频| 国产黄片美女视频| 精品午夜福利在线看| www.色视频.com| 嫩草影院入口| 乱码一卡2卡4卡精品| 亚洲无线观看免费| av福利片在线| av不卡在线播放| 久久亚洲国产成人精品v| 少妇丰满av| 最近的中文字幕免费完整| 三级国产精品欧美在线观看| 在线观看国产h片| 免费看日本二区| av在线观看视频网站免费| 好男人视频免费观看在线| 欧美老熟妇乱子伦牲交| 免费久久久久久久精品成人欧美视频 | 亚洲国产精品一区二区三区在线| 一级二级三级毛片免费看| 99热这里只有是精品50| 国产精品成人在线| 三级国产精品片| 亚洲成人av在线免费| 日韩强制内射视频| av黄色大香蕉| 男女边摸边吃奶| 国产成人aa在线观看| 成年人免费黄色播放视频 | 婷婷色综合大香蕉| 精品国产国语对白av| av福利片在线观看| 国产日韩一区二区三区精品不卡 | 免费人妻精品一区二区三区视频| videossex国产| 日韩电影二区| 69精品国产乱码久久久| 中文乱码字字幕精品一区二区三区| 18+在线观看网站| 日日撸夜夜添| 妹子高潮喷水视频| 久久久久久久久久成人| 久久综合国产亚洲精品| 亚洲国产毛片av蜜桃av| 最黄视频免费看| 色婷婷久久久亚洲欧美| 少妇高潮的动态图| 有码 亚洲区| 国产成人精品一,二区| 久久久久精品久久久久真实原创| av女优亚洲男人天堂| 男女无遮挡免费网站观看| 亚洲内射少妇av| 久久久久久久久久久丰满| 日本av免费视频播放| 性色av一级| 欧美精品一区二区大全| 大片免费播放器 马上看| 中文资源天堂在线| 18禁在线播放成人免费| 国产色爽女视频免费观看| 亚洲精品亚洲一区二区| 少妇熟女欧美另类| 亚洲欧美一区二区三区黑人 | 另类精品久久| 成人亚洲欧美一区二区av| 国语对白做爰xxxⅹ性视频网站| 亚洲,欧美,日韩| 免费观看a级毛片全部| 久热久热在线精品观看| 纵有疾风起免费观看全集完整版| 久久精品久久精品一区二区三区| 国产欧美日韩精品一区二区| 欧美精品国产亚洲| 日日啪夜夜爽| 欧美日韩在线观看h| 亚洲高清免费不卡视频|