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

        ?

        電力線OFDM/OQAM通信系統(tǒng)信道相位預(yù)處理均衡算法*

        2018-05-28 09:25:07
        電訊技術(shù) 2018年5期
        關(guān)鍵詞:復(fù)數(shù)實(shí)數(shù)載波

        (重慶郵電大學(xué) 通信核心芯片、協(xié)議及系統(tǒng)應(yīng)用創(chuàng)新團(tuán)隊(duì),重慶400065)

        1 引 言

        目前,各種電力線通信(Power Line Communication,PLC)協(xié)議標(biāo)準(zhǔn)中規(guī)定使用的調(diào)制技術(shù)是正交頻分復(fù)用(Orthogonal Frequency Division Multiplexing,OFDM),主要是因?yàn)檫@種多載波調(diào)制技術(shù)可以較好地解決多徑信道所引起的頻率選擇性衰落,同時(shí)信道均衡也相對(duì)簡單。然而,CP-OFDM技術(shù)的缺點(diǎn)也很明顯,其中一個(gè)弊端是較差的頻率選擇性;另外,循環(huán)前綴(Cyclic Prefix,CP)的插入也降低了頻譜效率。本文提出將5G中的一種新波形即正交頻分復(fù)用/偏移正交幅度調(diào)制(Orthogonal Frequency Division Multiplexing/Offset Quadrature Amplitude Modulation,OFDM/OQAM)技術(shù)[1]引入到PLC系統(tǒng)中,代替了傳統(tǒng)的CP-OFDM技術(shù),能夠獲得更高的頻帶利用率,通過選用時(shí)頻聚焦性良好的濾波器,如升余弦濾波器、擴(kuò)展高斯濾波器(Extended Gaussian Function,EGF)和各向同性正交變換算法(Isotropic Orthogonal Transform Algorithm,IOTA)[2],使其時(shí)域和頻域都具有較快的帶外衰減速度,能夠較好地適應(yīng)PLC的頻率掩碼[3]。

        然而,OFDM/OQAM技術(shù)的引入也面臨著一些難題,其中最棘手的是OFDM/OQAM系統(tǒng)存在著固有干擾,這種固有干擾在實(shí)數(shù)域信道下可以通過簡單的迫零(Forced Zero,ZF)均衡后取實(shí)部操作來消除,但是PLC復(fù)數(shù)信道所引入的時(shí)域符號(hào)間干擾(Inter-symbol Interference,ISI)和頻域載波間干擾(Inter-carrier Interference,ICI)并不能通過簡單的均衡后取實(shí)部操作來消除,這將嚴(yán)重影響電力線通信系統(tǒng)的性能,因此采用行之有效的均衡算法來抑制這些干擾是很有必要的。文獻(xiàn)[4]提出了一種兩步最小均方誤差(Minimum Mean Square Error,MMSE)均衡算法,首先運(yùn)用MMSE均衡并通過初始判決來消除某個(gè)子載波相鄰位置間的干擾,然后再用MMSE均衡來消除每個(gè)子載波上殘余的符號(hào)間干擾。文獻(xiàn)[5]提出了針對(duì)多路傳輸?shù)淖赃m應(yīng)正弦/余弦調(diào)制濾波器組均衡器(Adaptive Sine/Cosine-modulated Filter Bank Equalizer for Transmultiplexer,ASCET)來實(shí)現(xiàn)OFDM/OQAM系統(tǒng)的均衡,該方法使用了3個(gè)抽頭的ZF均衡器,具有比單抽頭的ZF均衡更好的誤碼性能。文獻(xiàn)[6]在分析了OFDM/OQAM系統(tǒng)固有干擾成因的基礎(chǔ)上,提出了一種尋找常數(shù)值γ的盲均衡算法,該算法γ值的選擇成為最后成功的關(guān)鍵,但其選擇十分困難。文獻(xiàn)[7]提出了一種迭代干擾消除均衡算法,首先通過ZF均衡對(duì)數(shù)據(jù)進(jìn)行初始判決,然后依據(jù)判決值重構(gòu)干擾項(xiàng),最后將干擾消除。研究結(jié)果表明,上述的均衡算法能夠取得比ZF均衡更好的誤碼性能,但是都存在復(fù)雜度過高的問題。

        本文在ZF均衡的基礎(chǔ)上提出一種基于信道相位預(yù)處理的均衡算法,性能較ZF均衡有一定提升,復(fù)雜度略高于ZF均衡,但卻遠(yuǎn)低于上述其他均衡算法。該均衡算法的基本思想是在均衡前將接收信號(hào)乘上一個(gè)相位因子,盡可能使等效信道從復(fù)數(shù)域向?qū)崝?shù)域靠近,然后通過取實(shí)操作就可以消除一部分純虛數(shù)干擾,在均衡時(shí)再除以信道函數(shù)與該相位因子乘積的實(shí)部,這樣就可以有效減小摻雜在有用信號(hào)中的干擾。本文的主要貢獻(xiàn)包括:將OFDM/OQAM技術(shù)應(yīng)用到PLC系統(tǒng)中,并和傳統(tǒng)的CP-OFDM技術(shù)進(jìn)行了性能對(duì)比;提出了一種新穎的信道均衡方法,能獲得復(fù)雜度與性能的折中,并與傳統(tǒng)的ZF均衡進(jìn)行了性能對(duì)比。

        2 電力線多徑信道OFDM/OQAM解調(diào)信號(hào)模型

        2.1 OFDM/OQAM系統(tǒng)模型

        OFDM/OQAM技術(shù)沿用傳統(tǒng)的CP-OFDM技術(shù),不同之處在于OQAM調(diào)制和濾波器的選擇。圖1是基于快速傅里葉逆變換(Inverse Fast Fourier Transform,IFFT)和快速傅里葉變換(Fast Fourier Transform,FFT)運(yùn)算的OFDM/OQAM的系統(tǒng)框圖。

        圖1 OFDM/OQAM系統(tǒng)結(jié)構(gòu)框圖Fig.1 Block diagram of OFDM/OQAM system

        2.2 OFDM/OQAM發(fā)送信號(hào)及應(yīng)用在PLC系統(tǒng)中滿足的條件

        OFDM/OQAM發(fā)送信號(hào)[8]可以表示為

        (1)

        式中:N表示子載波個(gè)數(shù);am,n表示第n個(gè)符號(hào)、第m個(gè)子載波傳輸?shù)姆?hào)數(shù)據(jù),它來自于QAM映射后符號(hào)的實(shí)部和虛部;υ0表示子載波間隔;τ0表示符號(hào)實(shí)部和虛部之間的偏移間隔,τ0=T/2,且滿足τ0υ0=1/2;gm,n(t)表示時(shí)頻格點(diǎn)坐標(biāo)為(m,n)的基函數(shù)。與傳統(tǒng)CP-OFDM系統(tǒng)在復(fù)數(shù)域不同的是,OFDM/OQAM系統(tǒng)僅僅滿足在實(shí)數(shù)域嚴(yán)格正交,實(shí)數(shù)域正交性可表示如下:

        R{j(m+n-m0-n0)+(m-m0)(n+n0)·

        Ag((n0-n)τ0,(m-m0)υ0)}=

        δm,m0δn,n0。

        (2)

        式中:當(dāng)m=m0時(shí),δm,m0=1;否則δm,m0=0。實(shí)數(shù)域正交性條件可進(jìn)一步解釋為當(dāng)m=m0,n=n0時(shí),〈gm,n,gm0,n0〉為實(shí)數(shù)符號(hào);不滿足這一條件時(shí),〈gm,n,gm0,n0〉為純虛數(shù)符號(hào)。Ag(τ,υ)為模糊函數(shù),表示為

        (3)

        然而,PLC是基帶傳輸系統(tǒng),若要將OFDM/OQAM應(yīng)用在PLC系統(tǒng)中,發(fā)送信號(hào)必須為實(shí)數(shù),文獻(xiàn)[9]提出了OFDM/OQAM的共軛對(duì)稱(Hermitian Symmetry,HS)形式,即HS-OQAM,并論證了電力線系統(tǒng)中實(shí)現(xiàn)HS-OQAM的條件為

        a0,n=aM,n=0,

        (4)

        am,n=aN-m,n(-1)D-M-ne-jφ0。

        (5)

        式中:M代表離散時(shí)間偏移,M=N/2;D=L-1,L為原型濾波器長度;φ0為附加相位值,通常為了方便,取值為0。

        2.3 OFDM/OQAM解調(diào)信號(hào)

        PLC信道通常是多徑衰落信道,典型的有4徑和15徑,根據(jù)著名的Zimmermann信道模型[10]可建模如下:

        (6)

        式中:gi為第i條路徑的加權(quán)系數(shù),表示沿著這個(gè)路徑的反射和傳輸因子,一條路徑上的反射越多,加權(quán)因子gi就越小;參數(shù)a0、a1和k都是由測(cè)量得到;τi表示第i條路徑的延時(shí)。該模型代表了N條不同路徑信號(hào)的疊加。

        經(jīng)過PLC多徑信道加噪聲之后,OFDM/OQAM接收信號(hào)可表示為

        r(t)=h(t)?s(t)+η(t)=

        6 Effect of continuous quality improvement on renal function in patients with chronic kidney disease of stage 3-4

        (7)

        式中:h(t,τ)表示多徑信道沖激響應(yīng),Δ為信道的最大時(shí)延擴(kuò)展。為了便于分析,此處暫不考慮噪聲的影響,則式(7)可進(jìn)一步化簡為

        g(t-τ-nτ0)e-j2πmυ0τdτ。

        (8)

        假設(shè)多徑信道的時(shí)延擴(kuò)展Δ遠(yuǎn)小于一個(gè)符號(hào)周期,則在0≤τ≤Δ時(shí)間內(nèi),可近似認(rèn)為g(t-τ-nτ0)≈g(t-nτ0),這樣式(8)可簡化為

        (9)

        (10)

        接收信號(hào)在經(jīng)過了接收濾波器之后,時(shí)頻格點(diǎn)(n0,m0)上的解調(diào)信號(hào)可表示為

        (11)

        3 傳統(tǒng)ZF均衡算法

        假設(shè)接收端已獲得理想的信道頻域響應(yīng)矩陣,接收信號(hào)在經(jīng)過ZF均衡后,可得到

        (12)

        上式等號(hào)右邊第二項(xiàng)為待解調(diào)數(shù)據(jù)受到的周圍符號(hào)的干擾。對(duì)上式進(jìn)行取實(shí)部操作,結(jié)果為

        (13)

        上文已分析過濾波器基函數(shù)的特性,當(dāng)m≠m0,n≠n0時(shí),〈gm,n(t),gm0,n0(t)〉為純虛數(shù),進(jìn)一步分析可發(fā)現(xiàn):

        (1)若為實(shí)數(shù)域信道,則干擾項(xiàng)為純虛數(shù),通過取實(shí)部操作可完全消除周圍數(shù)據(jù)的干擾,恢復(fù)出原發(fā)送信號(hào),即

        (14)

        (2)若為復(fù)數(shù)域信道,則干擾項(xiàng)為復(fù)數(shù),通過取實(shí)部操作也不能消除周圍數(shù)據(jù)的干擾,不可準(zhǔn)確恢復(fù)出原發(fā)送信號(hào),即

        (15)

        電力線OFDM/OQAM通信系統(tǒng)中不同子載波的頻率響應(yīng)各不相同且都為復(fù)數(shù),故僅僅通過ZF均衡并不能消除周圍數(shù)據(jù)符號(hào)的干擾。

        4 基于信道相位預(yù)處理的均衡算法

        前面已分析過在頻域?qū)崝?shù)信道下,OFDM/OQAM系統(tǒng)能夠通過取實(shí)部操作消除干擾,但是在頻域復(fù)數(shù)信道下并不能。由此帶來的啟示是:只要在系統(tǒng)接收濾波和信道均衡之前采取一定的方法對(duì)接收信號(hào)進(jìn)行處理,使接收信號(hào)等效信道盡量逼近實(shí)數(shù),即讓等效信道的虛部盡量小,就能夠通過取實(shí)部操作消除掉大部分干擾。通常在計(jì)算干擾分量時(shí)只考慮待解調(diào)數(shù)據(jù)周圍一階鄰域格點(diǎn)上的干擾,如圖2所示。

        圖2 一階鄰域格點(diǎn)干擾分布圖Fig.2 The disturbance distribution of first order neighborhood lattice sites

        這是因?yàn)镺FDM/OQAM系統(tǒng)選用的濾波器一般都具有良好的時(shí)頻聚焦性,距離較遠(yuǎn)的符號(hào)造成的干擾已經(jīng)非常小了,可忽略不計(jì)。此時(shí),式(11)可改寫為

        (16)

        為描述方便,這里將復(fù)數(shù)信道頻率響應(yīng)表示為幅度與相位的形式:Hm=|Hm|ejφm,φm表示復(fù)數(shù)信道的相位。

        對(duì)任意的子載波而言,在接收信號(hào)上乘以一個(gè)適當(dāng)?shù)南辔徽{(diào)整因子后再進(jìn)行接收濾波,此時(shí)信號(hào)變?yōu)?/p>

        rm0,n0e-jφ=〈r(t)e-jφ,gm0,n0(t)〉=

        am0,n0Hm0e-jφ+

        (17)

        分析式(17)可以發(fā)現(xiàn),應(yīng)盡可能選擇相位調(diào)整因子e-jφ,使得頻點(diǎn)為m0-1、m0和m0+1的等效信道為實(shí)數(shù)。實(shí)際上不可能使上述3個(gè)頻點(diǎn)的等效信道同時(shí)為實(shí)數(shù),但可以使某一頻點(diǎn)的等效信道為實(shí)數(shù),例如選擇相鄰前一子載波的相位作為相位調(diào)整因子e-jφm0-1,則頻點(diǎn)為m0-1的等效信道就為實(shí)數(shù)。將e-jφm0-1代入式(17),有

        rm0,n0e-jφm0-1=am0,n0Hm0e-jφm0-1+

        am0,n0Hm0e-jφm0-1+

        (18)

        此時(shí)再進(jìn)行取實(shí)部操作,可得

        R{rm0,n0e-jφm0-1}=am0,n0R{Hm0e-jφm0-1}+

        (19)

        式(19)等號(hào)右邊第一項(xiàng)為有用信號(hào),第二項(xiàng)為等效干擾。觀察發(fā)現(xiàn)接收信號(hào)乘以相位調(diào)整因子e-jφm0-1并取實(shí)部之后,等效干擾分量中來自第m0-1個(gè)子載波的ICI就被消除了,此時(shí)等效干擾分量中僅剩下ISI和來自第m0+1個(gè)子載波的ICI,如圖3所示。

        圖3 相位預(yù)處理后一階鄰域格點(diǎn)干擾分布圖Fig.3 The disturbance distribution of first order neighborhood lattice after phase pretreatment

        在進(jìn)行信道相位預(yù)處理之后,再采用ZF均衡,此時(shí)應(yīng)除以信道頻率響應(yīng)與相位調(diào)整因子乘積的實(shí)部,即

        am0,n0+

        (20)

        對(duì)比式(12)和式(20)可以看出,信道相位預(yù)處理操作減小了解調(diào)數(shù)據(jù)符號(hào)周圍的干擾分量,而計(jì)算量方面,僅僅是在接收信號(hào)的基礎(chǔ)上乘以一個(gè)相位調(diào)整因子,若僅統(tǒng)計(jì)一個(gè)符號(hào)上所有子載波信道均衡的復(fù)雜度,則相比于ZF均衡,增加的復(fù)雜度主要包括計(jì)算式(18)的8N次復(fù)數(shù)乘法和計(jì)算式(20)中均衡系數(shù)時(shí)的N次復(fù)數(shù)乘法,但是由于最后在均衡時(shí)去除了第m0-1個(gè)子載波上的干擾,此時(shí)其實(shí)又減少了3N次復(fù)數(shù)乘法運(yùn)算,所以信道相位預(yù)處理均衡算法相較于ZF均衡算法,總共增加了8N+N-3N=5N次復(fù)數(shù)乘法運(yùn)算。由于并沒有出現(xiàn)指數(shù)級(jí)或冪次級(jí)運(yùn)算,因此所提算法增加的復(fù)雜度相對(duì)于性能的提升還是可以接受的?;谛诺老辔活A(yù)處理的均衡算法原理框圖如圖4所示。

        圖4 基于信道相位預(yù)處理的均衡算法原理框圖Fig.4 Principle block diagram of equalization algorithm based on channel phase preprocessing

        5 性能仿真

        5.1 參數(shù)設(shè)置

        本文利用Matlab工具來進(jìn)行仿真驗(yàn)證,在對(duì)比了OFDM/OQAM技術(shù)和傳統(tǒng)的CP-OFDM技術(shù)應(yīng)用在PLC系統(tǒng)中的性能之后,接著驗(yàn)證了所提的基于信道相位預(yù)處理均衡算法抵抗干擾的有效性,并與傳統(tǒng)的ZF均衡進(jìn)行了性能對(duì)比。仿真是以國家電網(wǎng)公司發(fā)布的《電力線通信解決方案-低壓電力線寬帶載波通信技術(shù)規(guī)范第4-1部分:物理層通信協(xié)議(報(bào)批稿)》[11]為依據(jù),并在文獻(xiàn)[12]中Zimmermann提出的15徑慢時(shí)變信道下進(jìn)行的。仿真幀數(shù)為1 000幀,其他參數(shù)如表1所示。

        表1 仿真參數(shù)Tab.1 Simulation parameter

        5.2 仿真結(jié)果及分析

        圖5所示的是PLC系統(tǒng)中OFDM/OQAM技術(shù)和CP-OFDM技術(shù)的性能對(duì)比曲線,兩種技術(shù)是在同等的仿真條件下進(jìn)行,采用的都是ZF均衡。觀察圖5可以發(fā)現(xiàn),OFDM/OQAM系統(tǒng)的性能始終要好于CP-OFDM系統(tǒng),這是因?yàn)镺FDM/OQAM系統(tǒng)選用了時(shí)頻聚焦性良好的IOTA濾波器,使得系統(tǒng)有較強(qiáng)的抗ICI和ISI的能力,而CP-OFDM系統(tǒng)盡管有循環(huán)前綴的保護(hù),但是PLC時(shí)域信道很長的拖尾仍然會(huì)導(dǎo)致ISI的存在。以誤碼率1.0×10-6為例,OFDM/OQAM系統(tǒng)與CP-OFDM系統(tǒng)相比,大約有1 dB的性能提升。

        圖5 OFDM/OQAM和CP-OFDM性能對(duì)比曲線圖Fig.5 Performance comparison between OFDM/OQAM and CP-OFDM

        圖6所示的是本文提出的基于信道相位預(yù)處理的均衡算法與傳統(tǒng)的ZF均衡算法性能對(duì)比曲線圖,從圖中可以看出,信道相位預(yù)處理的均衡算法性能要稍好于ZF均衡,這是因?yàn)樗峋馑惴ㄏ啾扔赯F均衡算法,將前一子載波對(duì)當(dāng)前子載波產(chǎn)生的載波間干擾進(jìn)行了消除,還有不同符號(hào)同一頻點(diǎn)上的子載波對(duì)當(dāng)前子載波產(chǎn)生的符號(hào)間干擾以及后一子載波對(duì)當(dāng)前子載波產(chǎn)生的載波間干擾都沒有進(jìn)行處理,所以性能沒有提升太多。同樣以誤碼率1.0×10-6為例,基于信道相位預(yù)處理的均衡算法與ZF均衡算法相比,大約有0.5 dB的性能提升。

        圖6 所提均衡算法與傳統(tǒng)ZF均衡算法性能對(duì)比曲線圖Fig.6 Performance comparison between proposed equalization algorithm and traditional ZF equalization algorithm

        6 結(jié)束語

        本文將OFDM/OQAM技術(shù)成功應(yīng)用到PLC系統(tǒng)中,取得了比傳統(tǒng)CP-OFDM技術(shù)更好的性能優(yōu)勢(shì)。但是OFDM/OQAM技術(shù)引入的濾波器僅在實(shí)數(shù)域正交,經(jīng)過PLC復(fù)數(shù)域信道會(huì)產(chǎn)生嚴(yán)重的自干擾。本文在詳細(xì)分析了OFDM/OQAM技術(shù)采用ZF均衡會(huì)存在殘余干擾的原理后,提出了一種基于信道相位預(yù)處理的均衡算法。通過PLC15徑信道的仿真,驗(yàn)證了所提算法比傳統(tǒng)均衡算法具有更好的誤比特率性能。需要指出,本文所提算法只消除了部分子載波的干擾,未來的研究工作可圍繞進(jìn)一步消除剩余子載波的干擾來展開。

        [1] ACHAICHIA P,BOT M L,SIOHAN P.OFDM/OQAM:a solution to efficiently increase the capacity of future PLC networks[J].IEEE Transactions on Power Delivery,2011,26(4):2443-2455.

        [2] RAZAVI R,XIAO P,TAFAZOLLI R.Information theoretic analysis of OFDM/OQAM with utilized intrinsic interference[J].IEEE Signal Processing Letters,2014,22(5):618-622.

        [3] ALESSANDRO D S,TONELLO A M,LAMPE L.Adaptive pulse-shaped OFDM with application to in-home power line communications[J].Telecommunication Systems,2012,51(1):3-13.

        [4] IKHLEF A,LOUVEAUX J.An enhanced MMSE per subchannel equalizer for highly frequency selective channels for FBMC/OQAM systems[C]//Proceedings of 2009 Workshop on Signal Processing Advances in Wireless Communications.Perugia,Italy:IEEE,2009:186-190.

        [5] LIN H,SIOHAN P,TELECOM F.HS-OQAM PLC:long filter or equalizer?[C]//Proceedings of the Third Workshop on Power Line Communications(WSPLC).Udine, Italy:IEEE,2009:1-4.

        [6] SAVAUX V,BADER F,NAOUES M.Blind equalization using constant modulus algorithm adapted to OFDM/OQAM modulation[C]//Proceedings of 2016 International Symposium on Wireless Communication Systems.Poznan,Poland:IEEE,2016:28-33.

        [7] 孫希東.OFDM/OQAM系統(tǒng)在瑞利信道下的信道估計(jì)與均衡方法[J].自動(dòng)化技術(shù)與應(yīng)用,2016,35(3):36-41.

        SUN Xidong.Channel estimation and equalization method for OFDM/OQAM systems over Rayleigh channels[J].Automation Technology and Applications,2016,35(3):36-41.(in Chinese)

        [8] LIN B,TANG X,FANG X,et al.Efficient frequency-domain channel equalisation methods for OFDM/OQAM-PON with intensity modulation and direct detection[J].IET Communications,2017,11(6):872-877.

        [9] LIN H,SIOHAN P.Capacity analysis for indoor PLC using different multi-carrier modulation schemes[J].IEEE Transactions on Power Delivery,2010,25(1):113-124.

        [10] ABENOV R,POKAMESTOV D,GELTSER A.Multipath powerline communications channel(PLC) modelling[C]//Proceedings of 2015 IEEE International Conference on Microwaves,Communications,Antennas and Electronic Systems.Tel Aviv,Israel:IEEE,2015:1-4.

        [11] 國家電網(wǎng)公司.Q/GDW 374.3-2016低壓電力線寬帶載波通信技術(shù)規(guī)范 第4-1部分:物理層通信協(xié)議(報(bào)批稿)[S].北京:國家電網(wǎng)公司,2016.

        [12] ZIMMERMANN M,DOSTERT K.A multipath model for the powerline channel[J].IEEE Transactions on Communications,2002,50(4):553-559.

        猜你喜歡
        復(fù)數(shù)實(shí)數(shù)載波
        “實(shí)數(shù)”實(shí)戰(zhàn)操練
        評(píng)析復(fù)數(shù)創(chuàng)新題
        求解復(fù)數(shù)模及最值的多種方法
        數(shù)系的擴(kuò)充和復(fù)數(shù)的引入
        復(fù)數(shù)
        認(rèn)識(shí)實(shí)數(shù)
        1.1 實(shí)數(shù)
        比較實(shí)數(shù)的大小
        應(yīng)急廣播系統(tǒng)中副載波的構(gòu)建與應(yīng)用
        低壓載波通訊測(cè)試儀的開發(fā)與應(yīng)用
        国产经典免费视频在线观看 | 亚洲综合色区另类av| 一本色综合亚洲精品蜜桃冫 | 国产伦一区二区三区久久| 青青草成人在线免费视频| 亚洲中文字幕乱码一二三| 免费乱理伦片在线观看| 国产zzjjzzjj视频全免费| 免青青草免费观看视频在线| 青青草手机在线免费视频| 欧洲熟妇色xxxx欧美老妇软件 | 亚洲欧洲AV综合色无码| 99久久精品人妻一区| 日本一道综合久久aⅴ免费| 久久国产精久久精产国| 国产丰满乱子伦无码专| 日本按摩偷拍在线观看| 国产md视频一区二区三区| 欧洲色综合| 亚洲av乱码一区二区三区观影| 大香焦av一区二区三区| 中文字幕v亚洲日本| 99久热re在线精品99 6热视频| 青青草99久久精品国产综合| 蜜桃视频一区二区三区四| 久久久久九九精品影院| 最新亚洲精品国偷自产在线| 日韩中文字幕不卡网站| 精品av一区二区在线| 一区二区视频中文字幕| 精品人妻无码视频中文字幕一区二区三区| 国产啪精品视频网站丝袜| 99国产精品欲av麻豆在线观看| 少妇高潮久久蜜柚av| 亚洲色成人www永久在线观看| 亚洲国产美女精品久久| 国产免费的视频一区二区| 户外精品一区二区三区| 亚洲成av人片在线观看ww| 国内自拍偷拍亚洲天堂| 亚洲国产日韩一区二区三区四区 |