亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        復(fù)雜環(huán)境下的UWB脈沖成形算法

        2016-12-24 07:16:16孫希延紀(jì)元法
        關(guān)鍵詞:超寬帶誤碼率時域

        孫希延, 劉 健, 紀(jì)元法

        (1. 桂林電子科技大學(xué)信息與通信學(xué)院, 廣西 桂林 541004;2. 廣西精密導(dǎo)航技術(shù)與應(yīng)用重點(diǎn)實驗室, 廣西 桂林 541004)

        ?

        復(fù)雜環(huán)境下的UWB脈沖成形算法

        孫希延1,2, 劉 健1, 紀(jì)元法1,2

        (1. 桂林電子科技大學(xué)信息與通信學(xué)院, 廣西 桂林 541004;2. 廣西精密導(dǎo)航技術(shù)與應(yīng)用重點(diǎn)實驗室, 廣西 桂林 541004)

        針對超寬帶定位系統(tǒng)在多用戶和弱信號環(huán)境下接收端易出現(xiàn)碼間串?dāng)_(inter symbol interference, ISI)和多用戶干擾(multiuser interference, MUI)等問題,提出一種基于修正最小均方誤差估計(minimum mean square error estimation, MMSE)的自適應(yīng)迭代算法進(jìn)行超寬帶(ultra-wideband,UWB)脈沖設(shè)計。該算法選取修正Hermite多項式(modefied Hermite polynomial, MHP)作為脈沖設(shè)計基函數(shù)。通過分析MHP的時頻特性對最小均方誤差估計準(zhǔn)則下組合脈沖的各階MHP系數(shù)進(jìn)行修正,在此基礎(chǔ)上根據(jù)組合脈沖功率譜密度與美國聯(lián)邦通訊委員會(federal communications commission, FCC)輻射掩蔽之間的擬合程度設(shè)計自適應(yīng)算法,對組合脈沖進(jìn)行自適應(yīng)調(diào)整實現(xiàn)脈沖設(shè)計。仿真結(jié)果表明該方案獲得的脈沖有較高功率利用率,同時在多用戶接入和弱信號檢測方面也具備很強(qiáng)的適用性。

        修正Hermite多項式; 自適應(yīng)迭代; 多用戶干擾; 弱信號檢測

        0 引 言

        超寬帶通信的信號載體是時域很窄且占空比很低的脈沖。已有多種脈沖波形被提出用來設(shè)計超寬帶(ultra-wideband,UWB)脈沖,如Laplace脈沖、高斯脈沖、Rayleigh脈沖、Hermite脈沖、扁長橢球波脈沖等。但在美國聯(lián)邦通訊委員會(federal communications commission, FCC)規(guī)定的UWB頻段3.1~10.6 GHz[1]內(nèi)單一脈沖頻譜利用率低,限制了設(shè)計脈沖的功率。通過多種脈沖組合設(shè)計UWB脈沖是一種有效方案。

        文獻(xiàn)[2]通過組合多個Schlotz脈沖設(shè)計UWB脈沖,但組合脈沖時域?qū)挾容^大,通信過程中易發(fā)生碼間串?dāng)_(inter symbol interference,ISI),系統(tǒng)誤碼率偏高。文獻(xiàn)[3-4]將3.1~10.6 GHz頻段分塊,在各頻塊內(nèi)分別選取適合的高斯脈沖逼近FCC掩蔽,最后將各頻塊的脈沖組合成UWB脈沖。該方法設(shè)計參數(shù)易于獲取,但高頻段輻射功率偏低。文獻(xiàn)[5]分析了以扁長橢球波脈沖為基函數(shù)的UWB脈沖的多址性能,但該脈沖在低頻段的頻譜利用率不足。文獻(xiàn)[6]通過調(diào)整前三階Rayleigh脈沖的成形因子來設(shè)計脈沖,雖然該脈沖頻譜利用率優(yōu)于高斯脈沖,但多址性能并不突出。文獻(xiàn)[7-8]選用修正Hermite多項式(modified hermite polynominal, MHP)設(shè)計脈沖,因其基函數(shù)僅限于偶數(shù)階MHP,限制了脈沖功率的靈活性,且多址接入性能仍有待改進(jìn)。文獻(xiàn)[9-10]分別利用正余弦波窄帶信號設(shè)計UWB脈沖,所得脈沖在低信噪比環(huán)境中易受干擾。文獻(xiàn)[11]將粒子群優(yōu)化算法應(yīng)用于超寬帶脈沖設(shè)計,所得脈沖具有一定抗多址干擾(multiple access interference, MUI)性能,但用戶數(shù)大于50時,MUI性能嚴(yán)重下降。另外文獻(xiàn)[12]結(jié)合認(rèn)知無線電理論提出了基于特征值分解的脈沖整形方案來抑制窄帶通信系統(tǒng)對UWB通信的干擾,并以5.2 GHz 無線局域網(wǎng)絡(luò)(wireless local area networks,WLAN)和8.5 GHz無線電定位服務(wù)為干擾源進(jìn)行性能分析,但未對MUI性能進(jìn)行分析。

        為解決上述問題,本文選取各階相互正交的MHP作為UWB脈沖設(shè)計的基函數(shù),通過分析MHP時頻特性,提出一種基于修正最小均方誤差估計準(zhǔn)則(minimum mean square error estimation, MMSE)的自適應(yīng)迭代算法。該算法能兼容更多基函數(shù)并對脈沖功率譜密度(power spectral density, PSD)進(jìn)行自適應(yīng)調(diào)整,通過與廣泛使用的脈沖設(shè)計方案進(jìn)行仿真對比,驗證了該算法設(shè)計脈沖的良好性能。

        1 修正Hermite多項式

        集成電路或雙極性晶體管生成的超寬帶脈沖[13]可用Hermite多項式近似逼近。各階MHP相互正交且無直流分量,其表達(dá)式為

        h0(t)=e-t2/4τ2

        (1)

        (2)

        (3)

        式(2)、式(3)的頻域表達(dá)式分別為

        (4)

        (5)

        式中,H0(f)、Hn+1(f)分別為h0(t)、hn+1(t)的傅里葉變換(Fourier transform, FT)形式。

        前4階MHP時域波形如圖1所示。

        圖1 前4階Hermite脈沖時域波形(τ=0.8e-10)Fig.1 Time domain waveform of first 4 orders

        由于衰減因子e-t2/4的衰減速度大于MHP增長速度,故各階MHP時域脈寬趨于一致[14]。選取合理τ值,即可獲得合適的脈沖寬度降低ISI。

        由式(5)可推出

        Hn(f)∝(j)nfne-4π2f2τ2

        (6)

        式中,f為頻率。進(jìn)一步推導(dǎo)可得峰值頻率fpeak與τ的關(guān)系為

        (7)

        式(6)說明偶數(shù)階MHP的FT為實函數(shù),奇數(shù)階的FT為虛函數(shù)。式(7)表明隨著MHP階數(shù)的增高,其峰值頻率被搬移到更高頻段。從圖2的仿真結(jié)果可以看出:隨著MHP階數(shù)的增加,能量譜密度(energy spectral density,ESD)峰值將向更高頻率端移動。

        圖2 前4階Hermite脈沖的ESDFig.2 ESD of first 4 orders

        圖3和圖4分別表示τ對MHP時域脈寬和頻域ESD的影響:τ正比于MHP時域?qū)挾?反比于ESD帶寬。為保證設(shè)計脈沖的脈寬在1 ns以內(nèi),本文選取的各階MHP均滿足τ≤0.8e-10。

        圖3 τ對脈沖持續(xù)時間的影響Fig.3 Effect of τ on pulse duration

        圖4 τ對Hermite脈沖ESD的影響Fig.4 Effect of τ on ESD

        2 UWB脈沖自適應(yīng)迭代算法

        2.1 基函數(shù)選取

        記各階MHP的權(quán)系數(shù)組成的向量為B=[β1,β2,β3,…βL]T。各階MHP的FT組成的向量為G=[G1,τ1,G2,τ2,G3,τ3,…,GL,τL]T,其中L為設(shè)計脈沖中MHP的最大階數(shù),τ1,τ2,…,τL均小于0.8e-10。相應(yīng)設(shè)計脈沖的功率譜密度為

        (8)

        (9)

        式中,1≤i≤L。由于各階MHP具有正交性,式(9)可近似表示[8]為

        (10)

        2.2 MMSE準(zhǔn)則設(shè)計方案分析

        MMSE準(zhǔn)則是以誤差平方和最小為前提尋找最佳擬合曲線的方法。其表達(dá)式為

        (11)

        式中,f(t)為目標(biāo)函數(shù)。式(11)可改寫為

        (12)

        按MMSE準(zhǔn)則將設(shè)計脈沖的PSD向S(f)擬合,即求式(12)取得最小值時,設(shè)計脈沖中各階MHP的權(quán)重向量為

        Bmmse=[-1.726 4i,-0.578 1,1.053 3i,2.509 7,-2.785 1i,-1.712 4,1.101 9i,1.575 5,

        1.860 0i,-1.334 1,1.026 1i,1.502 2,-1.805 7i, 1.279 1,0.788 5i]

        圖5為MMSE準(zhǔn)則下設(shè)計脈沖PSD的擬合曲線。由于MMSE無法在具體頻點(diǎn)對功率值加以約束,故圖5中設(shè)計脈沖的低頻段PSD大于目標(biāo)函數(shù)。此外,由于奇數(shù)階MHP的FT是虛函數(shù),導(dǎo)致Bmmse中奇數(shù)項為虛數(shù),因此奇數(shù)階MHP實際上也無法用來進(jìn)行實際的脈沖設(shè)計[8]。

        圖5 MMSE組合脈沖的PSD擬合曲線Fig.5 PSD fitting curve of MMSE combined pulse

        2.3 基于修正MMSE的自適應(yīng)迭代算法

        為對圖5中低頻部分PSD進(jìn)行修正并克服Bmmse中奇數(shù)項不能參與脈沖設(shè)計的問題,本文設(shè)計一種新的迭代選擇算法。算法過程描述如下:

        步驟 1 參數(shù)初始化:設(shè)定BF空間φ(本文φ由前15階MHP組成)、迭代次數(shù)N、權(quán)重區(qū)間[α,β]、搜索步長Δ。

        步驟 2 PSD曲線擬合:用隨機(jī)數(shù)法和MMSE準(zhǔn)則分別向FCC輻射掩蔽擬合,生成系數(shù)向量Brandom和Bmmse。

        步驟 3 系數(shù)向量修正:生成新的系數(shù)向量E,其中E(i)=Bmmse(i),i=2,4,6,8,…;E(i)=Brandom(i),i=1,3,5,7,…。

        步驟 4 對E進(jìn)行自適應(yīng)迭代:以E為初始權(quán)重生成初始組合脈沖,在[α,β]區(qū)間內(nèi)對各BF的系數(shù)進(jìn)行迭代,即對各頻點(diǎn)做如下討論:

        首先定義迭代過程的約束條件為

        ε(f)=S(f)-Sp(f)>0

        (13)

        步驟 4.1 式(13)成立則記錄該組權(quán)值En,n為獲得En時的迭代次數(shù)。按步驟1設(shè)定的步長Δ繼續(xù)搜索,權(quán)值超出[α,β]或迭代次數(shù)達(dá)到N,則迭代終止。

        步驟 4.2 式(13)不成立,則按式(14)給出的自適應(yīng)公式對迭代步長進(jìn)行自適應(yīng)調(diào)整后繼續(xù)迭代

        (14)

        式中,自適應(yīng)因子μ的約束條件為μ≥1;ε(f)的約束條件由式(13)給出。

        圖6為步驟4所述初始脈沖的PSD,可以看到修正后的組合脈沖仍未滿足式(13),為此迭代過程應(yīng)側(cè)重降低低階MHP權(quán)重,并調(diào)節(jié)高階MHP在低頻段的分布,即通過式(14)增大低階MHP的搜索步長同時降低高階MHP在低頻段搜索步長。圖7為算法流程圖。

        圖6 系數(shù)替換后的初始脈沖的PSDFig.6 PSD of initial pulse

        圖7 基于MMSE修正的迭代算法流程圖Fig.7 Flow chart of iterative algorithm

        圖8 設(shè)計脈沖的PSDFig.8 PSD of designed pulse

        3 脈沖性能仿真

        超寬帶跳時(time hopping,TH)調(diào)制模型可表示為

        (15)

        r(t)=ru(t)+rmui(t)+n(t)

        (16)

        式中,ru(t)和rmui(t)分別代表接收機(jī)收到的有用信號和MUI噪聲。在不同信道環(huán)境中用TH模型對接收端的多用戶誤碼率(bit error rate, BER)和定時同步性能進(jìn)行分析[15-16]。對比脈沖為文獻(xiàn)[2]的正交Gaussian脈沖和文獻(xiàn)[6]的Reyleigh脈沖。仿真參數(shù)設(shè)置為:數(shù)據(jù)流速率Rb=20 Mb/s,δ=0.5 ns,Tf=6。

        3.1 接收端信號的誤碼率分析

        用戶數(shù)分別設(shè)置為10、20、100;信噪比仿真范圍為[-5 dB,30 dB]。仿真結(jié)果如圖9、圖10所示。圖9為不同信噪比下設(shè)計脈沖與Reyleigh脈沖的MUI性能比較??梢钥闯鲈谟脩魯?shù)達(dá)到100時,設(shè)計脈沖的誤碼率僅為Reyleigh脈沖的15%。圖10為設(shè)計脈沖與正交Gaussian脈沖的MUI性能對比。當(dāng)用戶數(shù)為100時,設(shè)計脈沖的誤碼率僅為正交Gaussian脈沖的60%;用戶數(shù)為50時,設(shè)計脈沖的BER僅為正交Gaussian脈沖的40%;用戶數(shù)為10時,設(shè)計脈沖的BER僅為正交Gaussian脈沖的2.7%。

        圖9 Reyleigh脈沖與設(shè)計脈沖誤碼率與信噪比關(guān)系對比圖Fig.9 Comparison of Reyleigh pulse and design pulse

        圖10 正交Gaussian脈沖與設(shè)計脈沖誤碼率與信噪比關(guān)系對比圖Fig.10 Comparison of orthogonal Gaussian pulse and design pulse

        3.2 弱信號環(huán)境下同步性能分析

        在信噪比為-20 dB的加性高斯白噪聲 (additive white gaussian noise, AWGN)信道中,將一導(dǎo)頻序列分別用本文設(shè)計脈沖、正交Gaussian脈沖和Reyleigh脈沖發(fā)射,在接收端檢測信號的到達(dá)時間(time of arrival, TOA)估計,以比較各脈沖的定時同步性能。其中TOA真實值為:21.344 ns。各脈沖經(jīng)信道傳播后的TOA測量值如表1所示。

        表1 各脈沖的TOA測量值

        如表1所示,設(shè)計脈沖的定時同步誤差低于另外兩種脈沖,這說明設(shè)計脈沖在復(fù)雜信道中具有更好的魯棒性。為對本文設(shè)計脈沖的相關(guān)性進(jìn)行分析,將上述3種脈沖輸入相關(guān)器,各脈沖的相關(guān)器輸出如圖11~圖13所示。

        圖11 Hermit脈沖的相關(guān)性能Fig.11 Correlator output of Hermit pulse

        圖12 正交高斯脈沖的相關(guān)性能Fig.12 Correlator output of quadrature Gauss pulse

        圖13 Reyleigh脈沖的相關(guān)性能Fig.13 Correlator output of Reyleigh pulse

        圖11~圖13分別為設(shè)計脈沖、正交Gaussian脈沖和Reyleigh脈沖的相關(guān)器輸出情況,通過分析可得:①在弱信號環(huán)境下設(shè)計脈沖的相關(guān)峰值高于另外兩種脈沖,易于檢測;②設(shè)計脈沖的噪底較薄,說明設(shè)計脈沖較另外兩種脈沖具有更好的自相關(guān)和互相關(guān)性能。

        4 結(jié)束語

        針對UWB信號多址干擾和碼間串?dāng)_問題,本文選用1到15階Hermite脈沖作為BF進(jìn)行UWB脈沖設(shè)計,通過分析MHP頻域特性,選取了合適的τ值來控制脈寬,減小ISI干擾。在設(shè)計脈沖PSD曲線向輻射掩蔽擬合的過程中提出用隨機(jī)系數(shù)來修正MMSE權(quán)重的方法,以此擴(kuò)充參與脈沖設(shè)計的BF數(shù)量。采用自適應(yīng)迭代算法來克服MMSE不能在具體頻點(diǎn)上對設(shè)計脈沖PSD進(jìn)行約束的問題。通過自適應(yīng)地調(diào)整迭代步長來保證各頻點(diǎn)上設(shè)計脈沖的PSD達(dá)到最佳擬合值。改變算法中自適應(yīng)系數(shù)即可在算法計算量和擬合曲線NESP間取得平衡。最后從通信鏈路誤碼率入手,在多用戶環(huán)境下對設(shè)計脈沖的性能進(jìn)行仿真并與常用脈沖設(shè)計方案進(jìn)行對比分析。仿真結(jié)果表明設(shè)計脈沖在獲得較高的NESP的同時能保證通信鏈路具有較強(qiáng)魯棒性。

        [1] Keshavarz S N, Hajizadeh S, Hamidi M, et al. A novel UWB pulse waveform design method[C]∥Proc.oftheInternationalConferenceonNextGenerationMobileapplications,ServicesandTechnologies, 2010:168-173.

        [2] Li L, Wang P, Wu X D, et al. Improved UWB pulse shaping method based on Gaussian derivatives[C]∥Proc.oftheIETInternationalCommunicationConferenceonWirelessMobileandComputing, 2011:438-442.

        [3] Menon M B, Gopakumar A, Iqbal N V. A hybrid approach for UWB pulse shaping[C]∥Proc.ofthe2ndInternationalConerenceonElectronics&CommunicationSystems, 2015:263-270.

        [4] Bowen B W, Nelson W S, Avise J C. Revision of part 15 of the commission’s rules regaring ultra-wideband transmission system[C]∥Proc.oftheNationalAcademyofSciencesoftheUnitedStatesofAmerica, 1993, 90(12):5574-5577.

        [5] Chen L L, Yu X, Dou Z. Research on the performance of the multiband impulse radio UWB communication system based on PSWF[J].JournalofHarbinEngineeringUniversity,2014, 35(4):499-503.(陳麗麗,于欣,竇崢.采用橢球波脈沖的多波段超寬帶性能研究[J].哈爾濱工程大學(xué)學(xué)報,2014,35(4):499-503.)

        [6] Chen R, Zeng X, Dong D. Algorithm for UWB pulse design and its performance analysis[C]∥Proc.oftheIEEEInternationalConferenceonIntegrationTechnology, 2007:74-77.

        [7] Wu X L, Sha X J, Zhang N T. Modified Hermite function based pulse shaping algorithm for ultra wideband communications[C]∥Proc.oftheRadioandWirelessSysmposium, 2007:395-398.

        [8] Lin D, Wu X L, Sha X J. Analysis on performance improvement in multiuser UWB system under imperfect power control[C]∥Proc.ofthe4thInternationalConferenceonWirelessCommunication,NetworkingandMobileComputing, 2008:1-4.

        [9] Shi X. UWB waveform design method for arrowband interference suppression based on raised cosine pulse[C]∥Proc.oftheIEEEInternationalConferenceonComputerScienceandAutomationEngineering, 2011:646-649.

        [10] Chen H, Jia Z, Cai X. Waveform design for UWB communication based on combined sinusoid Gaussian pulse[C]∥Proc.oftheInternationalConferenceonComputer,Mechatronics,andControntrolElectronicEngineering, 2010:267-271.

        [11] Jia Z B, Chen H. A new UWB pulse design method and performance analysis[J].JournalofYunnanUniversity(NaturalSciencesEdition), 2013,27(2):162-166. (賈占彪, 陳紅. 一種新的超寬帶脈沖設(shè)計方法與性能分析[J].云南大學(xué)學(xué)報(自然科學(xué)版),2013,27(2):162-166.)

        [12] El-Khamy R, Shaaban S, Ghaleb I, et al. An interference aware cognitive radio UWB(IA-CR-UWB) system using eigen pulse shaping[C]∥Proc.oftheInternationalConferenceonTechnologicalAdvancesinElectrical,ElectronicsandComputerEngineering, 2013:16-21.

        [13] Gu L A. Research and design of ultra wideband pulse generation circuit[J].ModernElectronicsTechnique, 2016, 39(3): 131-134.(顧麗愛. 超寬帶脈沖產(chǎn)生電路的研究與設(shè)[J].現(xiàn)代電子技術(shù), 2016, 39(3): 131-134.)

        [14] Wu X L, Qiu X, Ning X Y, et al. Orthogonal Hermite functions based orthogonal pulse shaping method[C]∥Proc.oftheElectricalandComputerEngineering, 2009:496-499.

        [15] Mishra S, Rajesh A, Bora P K. Performance of pulse shape modulation of UWB signals using composite Hermite pulses[C]∥Proc.oftheNationalConferenceonCommunications, 2012: 1-5.

        [16] Hidayat R, Miyanaga Y. IR-UWB pulse position modulation and pulse shape modulation through S-V channel model[C]∥Proc.oftheInternationalConferenceonCommunicationSoftwareandNetworks, 2010:214-217.

        UWB pulse shaping algorithm in complex environment

        SUN Xi-yan1,2, LIU Jian1, JI Yuan-fa1,2

        (1.SchoolofInformationandCommunication,GuilinUniversityofElectronicTechnology,Guilin541004,China;2.GuangxiKeyLaboratoryofPrecisionNavigationTechnologyandApplication,Guilin541004,China)

        In order to reduce the inter symbol interference (ISI) and multiuser interference (MUI) which exist in the ultra wideband positioning system under the multiuser and weak signal environment, an adaptive iterative algorithm based on minimum mean square error estimation (MMSE) for ultra-wideband (UWB) pulse design is proposed. Modified Hermite polynomial (MHP) is selected as the basis function of the designed pulse. By analyzing the time-frequency characteristics of MHP, coefficients of combined pulse are modified by the MMSE criterion. According to the fitting degree of power spectrum density between combined pulse and Federal communications commission (FCC) radiation masking, an adaptive algorithm is designed to realize pulse design. Simulation results show that the proposed scheme can obtain higher power utilization, and also obtain strong applicability in multiuser access and weak signal detection.

        modified Hermite polynomial(MHP); adaptive iteration; multiuser interference (MUI); weak signal detection

        2016-04-11;

        2016-09-13;網(wǎng)絡(luò)優(yōu)先出版日期:2016-10-08。

        國家自然科學(xué)基金(61162007,61271284,61362005,61561016);廣西自然科學(xué)基金(2013GXNSFA019004,2014GXNSFAA118352,2014GXNSFBA118280);廣西無線寬帶通信與信號處理重點(diǎn)實驗室主任基金;廣西信息科學(xué)實驗中心資助課題

        TN 914, TN 781

        A

        10.3969/j.issn.1001-506X.2016.12.05

        孫希延(1973-),女,研究員,博士,主要研究方向為衛(wèi)星導(dǎo)航和電子對抗。

        E-mail: sunxiyan1@163.com

        劉 健(1991-),男,碩士研究生,主要研究方向為超寬帶室內(nèi)定位。

        E-mail: sdjkjsdh@163.com

        紀(jì)元法(1975-),男,教授,博士,主要研究方向為信號處理、衛(wèi)星導(dǎo)航。

        E-mail: jiyuanfa@163.com

        網(wǎng)絡(luò)優(yōu)先出版地址:http://www.cnki.net/kcms/detail/11.2422.TN.20161008.1500.006.html

        猜你喜歡
        超寬帶誤碼率時域
        基于多模諧振器的超寬帶濾波器設(shè)計
        面向通信系統(tǒng)的誤碼率計算方法
        基于時域信號的三電平逆變器復(fù)合故障診斷
        超寬帶雷達(dá)建筑物結(jié)構(gòu)稀疏成像
        一種基于射頻開關(guān)切換的偽隨機(jī)編碼超寬帶多發(fā)多收雷達(dá)設(shè)計
        一種脊波導(dǎo)超寬帶濾波器
        電子制作(2016年11期)2016-11-07 08:43:34
        基于極大似然準(zhǔn)則與滾動時域估計的自適應(yīng)UKF算法
        基于時域逆濾波的寬帶脈沖聲生成技術(shù)
        基于時域波形特征的輸電線雷擊識別
        電測與儀表(2015年2期)2015-04-09 11:28:50
        泰克推出BERTScope誤碼率測試儀
        国产精品三级国产精品高| 成人午夜性a级毛片免费| 精品人妻无码一区二区三区蜜桃一 | 视频一区视频二区自拍偷拍| 蜜臀一区二区三区精品 | 久精品国产欧美亚洲色aⅴ大片| 亚洲综合av在线在线播放| 日本岛国精品中文字幕| 国产三级一区二区三区在线观看| 国产成人高清在线观看视频| 扒开腿狂躁女人爽出白浆| 天堂一区人妻无码| 国产精品久久久久免费a∨不卡| 人妻露脸国语对白字幕| 亚洲视频高清一区二区| 人妻丰满熟妇av无码区| 天天操夜夜操| 激情亚洲综合熟女婷婷| 亚洲中文字幕第15页| 欧美成人午夜免费影院手机在线看 | 99riav国产精品视频| 无码人妻av一区二区三区蜜臀| 国产精品视频一区日韩丝袜| 中文字幕视频二区三区| 亚洲一区二区国产激情| 久久99精品久久水蜜桃| 国产精品白浆一区二小说| 精品熟妇av一区二区三区四区| 亚洲熟妇av一区二区在线观看| 国内精品久久久人妻中文字幕| 亚洲精品久久久无码av片软件| 精品无码人妻久久久一区二区三区| 国产精品高清视亚洲一区二区 | 内射人妻少妇无码一本一道| 国产免费午夜a无码v视频| 激情综合五月天开心久久| 香港三级日本三韩级人妇久久| 后入到高潮免费观看| 抽搐一进一出试看60秒体验区| 国产美女黄性色av网站| 亚洲第一女人的天堂av|