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

    基于高階累積量和譜分析的數(shù)字調(diào)制信號(hào)識(shí)別

    2017-01-10 03:53:12孫汝峰劉順蘭
    關(guān)鍵詞:類間譜線特征參數(shù)

    孫汝峰,劉順蘭

    (杭州電子科技大學(xué) 通信工程學(xué)院,浙江 杭州 310018)

    基于高階累積量和譜分析的數(shù)字調(diào)制信號(hào)識(shí)別

    孫汝峰,劉順蘭

    (杭州電子科技大學(xué) 通信工程學(xué)院,浙江 杭州 310018)

    提出了一種基于高階累積量和譜分析識(shí)別多種數(shù)字調(diào)制信號(hào)的算法。首先根據(jù)各調(diào)制信號(hào)四階和八階累積量的不同,定義一個(gè)特征參數(shù)實(shí)現(xiàn)信號(hào)的類間識(shí)別;其次根據(jù)不同調(diào)制信號(hào)二次方譜與四次方譜的不同,提取出相應(yīng)的特征參數(shù),從而實(shí)現(xiàn)信號(hào)的類內(nèi)識(shí)別。仿真實(shí)驗(yàn)結(jié)果表明,該方法在較低信噪比條件下可以對(duì)2/4/8PSK、2/4/8FSK信號(hào)實(shí)現(xiàn)識(shí)別,且識(shí)別率較高,具有很強(qiáng)的實(shí)用性。

    調(diào)制識(shí)別;高階累積量;二次方譜;四次方譜

    0 引言

    通信信號(hào)的調(diào)制識(shí)別在電子偵察和無(wú)線電監(jiān)控等領(lǐng)域占據(jù)著十分重要的地位,主要任務(wù)是在調(diào)制信息未知的情況下確定調(diào)制信號(hào)的調(diào)制方式以及估計(jì)信號(hào)的一些參數(shù)(如載波頻率、波特率等),為之后的信號(hào)分析處理提供依據(jù)。

    近年來(lái)國(guó)內(nèi)外提出了很多比較有效的調(diào)制信號(hào)識(shí)別方法,主要可以分為決策理論方法和特征參數(shù)模式識(shí)別方法[1-4],其中決策論方法需要的先驗(yàn)知識(shí)較多[1],相比較下特征參數(shù)模式識(shí)別的方法更為實(shí)用,常用的特征參數(shù)模式分類特征有信號(hào)的瞬時(shí)特征、高階累積量特征、小波變換特征、譜相關(guān)特征等。參考文獻(xiàn)[2]基于信號(hào)瞬時(shí)特征識(shí)別對(duì)噪聲比較敏感,在低信噪比的環(huán)境下識(shí)別率比較低。參考文獻(xiàn)[3]利用四階累積量實(shí)現(xiàn)2PSK、4PSK和8PSK信號(hào)的分類,并證明高階累積量對(duì)信號(hào)星座圖的平移,尺度和相位旋轉(zhuǎn)等變換具有不變性。這種方法對(duì)于識(shí)別少數(shù)的信號(hào)具有不錯(cuò)的性能,但是對(duì)于識(shí)別較多的調(diào)制信號(hào)僅僅基于高階累積量特征的識(shí)別需要信號(hào)序列的統(tǒng)計(jì)信息,需要較大的樣本空間才能得到良好的識(shí)別效果,無(wú)疑增加了算法復(fù)雜度,并且MFSK信號(hào)的高階累積量值都相同,直接運(yùn)用高階累積量無(wú)法識(shí)別出MFSK信號(hào)。參考文獻(xiàn)[4]基于小波變換,則需要信號(hào)精確的相位信息。參考文獻(xiàn)[5]通過(guò)分析信號(hào)的譜相關(guān)平面圖,提取出一組譜相關(guān)特征參數(shù)來(lái)實(shí)現(xiàn)對(duì)信號(hào)的識(shí)別。

    本文提出一種基于高階累積量與譜分析的調(diào)制識(shí)別方法。利用高階累積量提取更少的特征參數(shù)對(duì)信號(hào)進(jìn)行類間識(shí)別。利用不同信號(hào)經(jīng)過(guò)非線性變換后譜線特征的不同提取信號(hào)的特征參數(shù)對(duì)信號(hào)進(jìn)行類內(nèi)識(shí)別,在低信噪比環(huán)境下實(shí)現(xiàn)更高的識(shí)別率。計(jì)算機(jī)仿真結(jié)果驗(yàn)證了算法的實(shí)用性和有效性。

    1 信號(hào)模型

    一般情況下接收到的受噪聲污染過(guò)的數(shù)字調(diào)制信號(hào)的模型可表示為[5]:

    (1)

    其中,k=1,2,3,...,N,N為發(fā)送碼元序列的長(zhǎng)度,Ik表示碼元序列,p(t)為基帶碼元波形,T為碼元寬度,fc為載波頻率,θc為載波相位,Es為發(fā)送碼元波形的能量,n(t)為零均值的高斯白噪聲。

    接收端對(duì)接收信號(hào)進(jìn)行預(yù)處理,假設(shè)載波頻率、相位、定時(shí)同步,下變頻后的復(fù)基帶信號(hào)可以表示為[6]:

    (2)

    其中,Δθc為載波相位差。

    調(diào)制信號(hào)的類型不同,Ik的表達(dá)形式也不一樣,具體如下:

    對(duì)MPSK信號(hào):

    Ik=ak+jbk∈{ej2π(m-1)/M,m=1,2,...,M}

    (3)

    對(duì)MFSK信號(hào):

    Ik=ej2πfkt

    (4)

    其中,fk∈{(2m-1-M)Δf,m=1,2,...,M} ,Δf為信號(hào)的頻率間隔。

    2 高階累積量和高階矩

    對(duì)于一個(gè)具有零均值的復(fù)隨機(jī)過(guò)程X(t),其高階矩定義為[7]:

    Mpq=E[X(t)(p-q)X*(t)q]

    (5)

    其中,E[·]表示求期望運(yùn)算。

    累積量定義為[8]:

    Cpq=Cum{X(t),…,X(t),X*(t),…,X*(t)}

    (6)

    其中,X(t)為p-q項(xiàng),X*(t)為q項(xiàng)。Cum為累積矩,* 表示函數(shù)的共軛。

    常用的累積量與矩的關(guān)系如下[9]:

    C20=M20

    (7)

    C21=M21

    (8)

    (9)

    (10)

    (11)

    (12)

    (13)

    因發(fā)送信號(hào)s(t)與高斯白噪聲n(t)兩者獨(dú)立,根據(jù)累積量的性質(zhì)由式(1)可得:

    Cum(x(t))=Cum(s(t))+Cum(n(t))

    (14)

    由于零均值高斯白噪聲大于二階的累積量值為零,則式(14)也可表示為:Cum(x(t))=Cum(s(t))。即接收信號(hào)的高階累積量值等于發(fā)送信號(hào)的高階累積量值,據(jù)此可以消除高斯白噪聲的影響。

    假設(shè)符號(hào)發(fā)送概率相等,信號(hào)能量為Es,且無(wú)高斯噪聲影響,采用總體平均代替統(tǒng)計(jì)平均的方法[10],可以得到各種調(diào)制信號(hào)的高階累積量理論值,具體如表1所示。

    表1 數(shù)字調(diào)制信號(hào)的理論高階累積量值

    調(diào)制方式C20C21C40C42C60C63C802FSK0Es0-E2s04E3s04FSK0Es0-E2S04E3s08FSK0Es0-E2s04E3s02PSKEsEs-2E2s-2E2s16E3s16E3s-272E4s4PSK0EsE2s-E2s04E3s-34E4s8PSK0Es0-E2s04E3sE4s

    3 基于信號(hào)高階累積量實(shí)現(xiàn)調(diào)制信號(hào)的類間識(shí)別

    分析表1可以看出,MFSK與MPSK的C80值存在著較大的差別,可以利用C80這個(gè)值實(shí)現(xiàn)MFSK與MPSK的類間識(shí)別。本文定義特征參數(shù)F如下:

    (15)

    特征參數(shù)選取信號(hào)的八階和四階累積量,可以減小噪聲的影響,采用絕對(duì)值形式可以減小相位抖動(dòng)對(duì)特征參數(shù)的影響,采用比值的形式可以消除幅度對(duì)特征參數(shù)的影響。根據(jù)表1可以得到各調(diào)制信號(hào)的F值,具體如表2所示。由表2可以看出,MFSK信號(hào)的F值均為0,MPSK信號(hào)的F值均為非0值,利用特征參數(shù)F可以實(shí)現(xiàn)MFSK和MPSK的類間識(shí)別。

    表2 各調(diào)制信號(hào)的特征參數(shù)F的理論

    4 調(diào)制信號(hào)的譜線特征

    由于不同的調(diào)制信號(hào)經(jīng)過(guò)非線性變換后,其頻譜會(huì)呈現(xiàn)不一樣的譜線特征,因此利用譜線特征可以實(shí)現(xiàn)信號(hào)的類內(nèi)識(shí)別。

    4.1 MPSK的二次方譜與四次方譜的譜線特征

    *p(t-kT)p(t-nT)]

    (16)

    [11]分析了信號(hào)的譜線生成特性,當(dāng)E[s2(t)]具有周期時(shí)變性時(shí),其頻域(二次方譜)會(huì)產(chǎn)生出離散的譜線。而對(duì)4PSK信號(hào),由式(16)可知E[s2(t)]=0,因此不具有周期時(shí)變性,所以其二次方譜不存在離散譜線。

    同理,4PSK四次方形式的統(tǒng)計(jì)期望值為:

    (17)

    式(17)是一個(gè)以T為周期的函數(shù),其頻域會(huì)產(chǎn)生離散譜線,相應(yīng)的傅里葉變換形式(信號(hào)的四次方譜)為:

    *A(f)δ(f-k/T)

    (18)

    同理可以分析得到,2PSK信號(hào)的二次方譜和四次方譜均存在譜線,8PSK信號(hào)的二次與四次方譜均不存在譜線。圖1給出了無(wú)噪聲干擾下4PSK的譜線特征圖。

    圖1 4PSK的譜線特征圖

    4.2 MFSK的二次方譜的譜線特征

    根據(jù)之前分析,MFSK的等效復(fù)低通信號(hào)形式可用式(2)表示,其中Ik用式(4)表示。不失一般性,p(t)為升余弦型脈沖[12],滾降系數(shù)α取0<α<0.5。以2FSK信號(hào)為例,2FSK信號(hào)做二次方譜為:

    (19)

    其中,A(f)=F{p2(t)},由于α<0.5的限制,由式(19)可知2FSK信號(hào)的二次方譜在±2Δf處有兩條明顯譜線。同理分析可得,4FSK信號(hào)的二次方譜在±2Δf、±6Δf處有4條譜線,8FSK信號(hào)的二次方譜在±2Δf,±6Δf,±10Δf,±14Δf處有8條譜線。

    5 基于信號(hào)的譜線特征實(shí)現(xiàn)調(diào)制信號(hào)的類內(nèi)識(shí)別

    對(duì)于MPSK信號(hào),由于其不同階數(shù)的PSK信號(hào)特征參數(shù)F的值不同,可以通過(guò)本文定義的F值實(shí)現(xiàn)類內(nèi)識(shí)別。但是本文建議采用MPSK信號(hào)的二次方譜和四次方譜零頻處是否存在譜線作為識(shí)別特征值。如果信號(hào)的二次方譜在零頻位置處存在譜線,令其特征參數(shù)f1=1,否則f1=0。如果信號(hào)的四次方譜在零頻位置處存在譜線,令其特征參數(shù)f2=1,否則f2=0。定義N_p=[f1,f2],據(jù)之前的分析,MPSK信號(hào)的N_p值表示如下:

    (20)

    由式(20)可以實(shí)現(xiàn)MPSK信號(hào)的類內(nèi)識(shí)別。

    對(duì)于MFSK信號(hào),由于其二次方譜的譜線個(gè)數(shù)對(duì)應(yīng)信號(hào)的M值,因此這里提取特征參數(shù)N_f,如果其所測(cè)信號(hào)二次方譜中譜線的個(gè)數(shù)N_f=2,則為2FSK信號(hào);如果N_f=4,則為4FSK信號(hào);如果N_f=8,則為8FSK信號(hào)。

    本文建議的基于高階累積量和譜分析的數(shù)字調(diào)制信號(hào)識(shí)別流程圖如圖2所示。

    圖2 調(diào)制信號(hào)類間與類內(nèi)識(shí)別流程圖

    6 性能仿真與分析

    設(shè)MPSK、MFSK信號(hào)的載波頻率為5 kHz,采樣頻率為40 kHz,碼元速率為1 000 b/s,且MFSK的頻率間隔為5 kHz,數(shù)據(jù)長(zhǎng)度N=2 000,加性噪聲為高斯白噪聲。在相同的信噪比環(huán)境下,對(duì)信號(hào)進(jìn)行100次獨(dú)立試驗(yàn)[13-14]。

    由于特征參數(shù)F的幅度范圍比較大,為觀察方便,將其劃分為兩段([0-5],[20-80])。分別繪制出調(diào)制信號(hào)的特征參數(shù)F隨信噪比變化曲線,如圖3所示。

    圖3 不同信噪比情況下調(diào)制信號(hào)的特征參數(shù)F值

    從圖3可以看出,MFSK與MPSK信號(hào)的F值都比較接近于前面計(jì)算的理論值,兩者的F值存在著較為明顯的差別,這就說(shuō)明了可以利用F值來(lái)實(shí)現(xiàn)兩者的類間識(shí)別。選取F=1作為類間識(shí)別的閾值,當(dāng)信號(hào)的F值大于1的時(shí)候,信號(hào)判定為MPSK信號(hào),否則判定為MFSK信號(hào)。

    在實(shí)現(xiàn)MFSK與MPSK信號(hào)的類間識(shí)別仿真之后,利用上述所提的特征參數(shù)根據(jù)算法流程對(duì)兩者進(jìn)行類內(nèi)識(shí)別的仿真。同一信噪比,對(duì)所有信號(hào)進(jìn)行500次獨(dú)立試驗(yàn),取識(shí)別正確的次數(shù)和識(shí)別總次數(shù)的比值為信號(hào)的識(shí)別率,取信號(hào)的類間識(shí)別率與類內(nèi)識(shí)別率的乘積作為信號(hào)的總識(shí)別率,仿真結(jié)果如圖4所示。

    圖4 MPSK和MFSK信號(hào)的識(shí)別率

    由圖4可知,采用本文建議的方法,MPSK信號(hào)在信噪比為1 dB情況下,2PSK、4PSK、8PSK信號(hào)的識(shí)別率分別為100%、85%、100%;當(dāng)信噪比大于2 dB時(shí),3種PSK信號(hào)的識(shí)別率都達(dá)到了100%。MFSK信號(hào)在信噪比為1 dB的情況下,2FSK、4FSK、8FSK信號(hào)的識(shí)別率分別為54%、67%、45%;當(dāng)信噪比等于3 dB時(shí),3種FSK信號(hào)信號(hào)的識(shí)別率達(dá)到了85%以上;當(dāng)信噪比大于等于4 dB時(shí),3種FSK信號(hào)信號(hào)的識(shí)別率達(dá)到了100%。由此可見(jiàn),本文建議的方法在低信噪比下也取得了較高的識(shí)別率。MFSK信號(hào)的識(shí)別率在信噪比較低的情況下比MPSK信號(hào)的識(shí)別率要低,分析原因主要是因?yàn)槭茉肼曈绊戭愰g識(shí)別的特征值F在信噪比低的情況下要比理論值偏大,在對(duì)MFSK信號(hào)識(shí)別時(shí),會(huì)出現(xiàn)大于閾值的情況,導(dǎo)致判斷出錯(cuò)。

    參考文獻(xiàn)[9]中分別采用高階累積量的特征值和微分后的累積量構(gòu)造特征值的方法來(lái)識(shí)別MPSK和MFSK信號(hào),在相同的仿真實(shí)驗(yàn)環(huán)境和信號(hào)參數(shù)設(shè)置,本文方法與參考文獻(xiàn)[9]進(jìn)行性能對(duì)比試驗(yàn),其對(duì)比結(jié)果如表3所示。從表3可以看出,相比參考文獻(xiàn)[9],在相同條件下本文所提方法的識(shí)別性能顯著提高。

    表3 不同方法的識(shí)別率對(duì)比(%)

    7 結(jié)論

    本文基于高階累積量和譜分析的理論知識(shí),定義了一個(gè)基于信號(hào)八階累積量和四階累積量的特征參數(shù)F,用來(lái)實(shí)現(xiàn)MFSK和MPSK信號(hào)的類間識(shí)別。由于高斯白噪聲高于二階的累積量值為零,因此此方法具有很好的抗噪聲性能。根據(jù)不同信號(hào)的二次方譜與四次方譜的譜線特性不同,分別提取出特征參數(shù)N_p和N_f對(duì)MPSK和MFSK信號(hào)實(shí)現(xiàn)類內(nèi)識(shí)別。計(jì)算機(jī)仿真結(jié)果表明,此方法在低信噪比的情況下可以取得理想的識(shí)別率,證明了此方法的有效性,抗噪聲性能較強(qiáng)。

    參考文獻(xiàn)

    [1] NANDI A K,AZZOUZ E E.Algorithms for automatic modulation recognition of communication signals[J].IEEE Transactions on Communication, 1998, 46(4):431-436.

    [2] 位小記,謝紅,郭慧.基于瞬時(shí)特征參數(shù)的數(shù)字調(diào)制識(shí)別算法[J]. 傳感器與微系統(tǒng),2011,30(1):127-130.

    [3] 陳衛(wèi)東,楊紹全.利用累量不變量對(duì)MPSK信號(hào)分類[J].西安電子科技大學(xué)學(xué)報(bào),2002,29(2):229-232.

    [4] 胡建偉,湯建龍,楊紹全.使用小波變換的MPSK信號(hào)調(diào)制類型識(shí)別[J].電路與系統(tǒng)學(xué)報(bào),2006,11(3),130-134.

    [5] 韓國(guó)棟,蔡斌,鄔江興.調(diào)制分析與識(shí)別的譜相關(guān)方法[J].系統(tǒng)工程與電子技術(shù),2001,23(3):34-36.

    [6] 徐曉瑤,劉娟,楊東.多徑信道下MPSK 信號(hào)調(diào)制識(shí)別算法的研究[J].電子技術(shù)應(yīng)用,2010,36(2):103-106.

    [7] PROAKIS J G. Digital communicatio (4th edition)[M]. NewYork: McGraw-HillBook Co.,2001.

    [8] 劉獻(xiàn)玲.基于累量的OFDM信號(hào)調(diào)制識(shí)別[D].西安:西安科技大學(xué),2007.

    [9] 包錫銳,吳瑛,周欣. 基于高階累積量的數(shù)字調(diào)制信號(hào)識(shí)別算法[J]. 信息工程大學(xué)學(xué)報(bào), 2007,80(4):463-467.

    [10] SHERMEH A E,GHAZALIAN R. Recognition of communication signal types using genetic algorithm and support vector machines based on the higher-order statistics[J].Digital Signal Processing, 2010, 20(6):1748-1757.

    [11] 楊琳,徐小東,路友榮,等.常見(jiàn)數(shù)字通信信號(hào)的譜線特征分析[J]. 電子與信息學(xué)報(bào),2009,31(5):1067-1071.

    [12] 丁革媛,高寶芹,孫強(qiáng),等.無(wú)線通信技術(shù)的發(fā)展研究[J].微型機(jī)與應(yīng)用,2014,33(10):1-3,6.

    [13] 吳小紅.LabVIEW中調(diào)用MATLAB進(jìn)行信號(hào)處理方法研究[J].微型機(jī)與應(yīng)用,2014,33(10):66-68.

    [14] 秦明偉,李湘魯,梁亞均,等.一種PSK符號(hào)率盲估計(jì)新方法[J].電子技術(shù)應(yīng)用,2014,40(8):105-108.

    Digital modulation signals recognition based on high order cumulant and spectral analysis

    Sun Rufeng,Liu Shunlan

    (School of Communication Engineering,Hangzhou Dianzi University,Hangzhou 310018,China)

    This paper proposes a digital modulation signal recognition technique based on higher order cumulant and spectrum analysis. Firstly, a feature based on the difference of signal’s four order and eight order cumulant is defined to realize the inter class recognition .Then another characteristic parameters based on the difference of signals’ square spectrum and quartic spectrum are extracted to realize the within class identification. The simulation results show that the proposed technique can recognize 2/4/8PSK, 2/4/8FSK at low SNR, achieves high recognition rate, and is very practical.

    modulation recognition; high order cumulant; square spectrum; quartic spectrum

    TN911.72

    A

    1674-7720(2016)04-0076-04

    孫汝峰,劉順蘭.基于高階累積量和譜分析的數(shù)字調(diào)制信號(hào)識(shí)別[J] .微型機(jī)與應(yīng)用,2016,35(4):76-79,84.

    2015-08-20)

    孫汝峰(1990-),男,碩士研究生,主要研究方向:數(shù)字調(diào)制信號(hào)的調(diào)制識(shí)別。

    劉順蘭(1965-),女,碩士,教授,主要研究方向:信號(hào)處理,軟件無(wú)線電。

    猜你喜歡
    類間譜線特征參數(shù)
    故障診斷中信號(hào)特征參數(shù)擇取方法
    基于特征參數(shù)化的木工CAD/CAM系統(tǒng)
    基于HITRAN光譜數(shù)據(jù)庫(kù)的合并譜線測(cè)溫仿真研究
    基于OTSU改進(jìn)的布匹檢測(cè)算法研究
    基于貝葉斯估計(jì)的多類間方差目標(biāo)提取*
    基于類間相對(duì)均勻性的紙張表面缺陷檢測(cè)
    鐵合金光譜譜線分離實(shí)驗(yàn)研究
    基于PSO-VMD的齒輪特征參數(shù)提取方法研究
    基于改進(jìn)最大類間方差法的手勢(shì)分割方法研究
    鍶原子光鐘鐘躍遷譜線探測(cè)中的程序控制
    国产淫语在线视频| 久久久色成人| 亚洲av成人精品一二三区| 亚洲精品aⅴ在线观看| 亚洲精华国产精华液的使用体验| 99热6这里只有精品| 青青草视频在线视频观看| 久久国产乱子免费精品| 精品熟女少妇av免费看| 欧美 日韩 精品 国产| 中文字幕免费在线视频6| 欧美3d第一页| 国产成人aa在线观看| 狂野欧美激情性bbbbbb| av国产精品久久久久影院| 国产精品精品国产色婷婷| 国产亚洲午夜精品一区二区久久| 大码成人一级视频| 亚州av有码| 亚洲av男天堂| www.色视频.com| 国产成人91sexporn| 中文天堂在线官网| 欧美区成人在线视频| 在线观看www视频免费| 男女免费视频国产| 欧美97在线视频| 天天影视国产精品| 人妻 亚洲 视频| 国产av一区二区精品久久| 美女高潮到喷水免费观看| videosex国产| 亚洲av电影在线进入| 岛国毛片在线播放| 大陆偷拍与自拍| 青草久久国产| 精品欧美一区二区三区在线| 成人18禁高潮啪啪吃奶动态图| 精品一区在线观看国产| 国语对白做爰xxxⅹ性视频网站| 99热全是精品| 免费黄频网站在线观看国产| 国产成人一区二区在线| 99国产综合亚洲精品| 成人亚洲欧美一区二区av| 成人手机av| 一本色道久久久久久精品综合| 成人影院久久| 国产成人a∨麻豆精品| 欧美人与性动交α欧美软件| 国产亚洲一区二区精品| 国产爽快片一区二区三区| 国产精品麻豆人妻色哟哟久久| 午夜免费男女啪啪视频观看| 国产精品一区二区在线观看99| 久久精品国产亚洲av高清一级| 最新的欧美精品一区二区| 亚洲欧美成人综合另类久久久| 婷婷丁香在线五月| 我的亚洲天堂| 中文字幕精品免费在线观看视频| 侵犯人妻中文字幕一二三四区| 韩国精品一区二区三区| 人妻 亚洲 视频| 精品免费久久久久久久清纯 | 免费看不卡的av| 精品少妇黑人巨大在线播放| 大陆偷拍与自拍| 男女免费视频国产| 亚洲av成人精品一二三区| 中国美女看黄片| 波多野结衣一区麻豆| 久久久久国产一级毛片高清牌| 黑人欧美特级aaaaaa片| 亚洲欧美一区二区三区久久| 丰满饥渴人妻一区二区三| 国产欧美日韩一区二区三区在线| 久久久精品免费免费高清| 成年人免费黄色播放视频| 成年美女黄网站色视频大全免费| a级毛片在线看网站| 久久这里只有精品19| 午夜久久久在线观看| 国产日韩一区二区三区精品不卡| 日韩免费高清中文字幕av| 日韩一区二区三区影片| 狠狠精品人妻久久久久久综合| 精品视频人人做人人爽| www.999成人在线观看| 性色av乱码一区二区三区2| 国产精品成人在线| 操美女的视频在线观看| 国产亚洲午夜精品一区二区久久| 久久 成人 亚洲| 中文乱码字字幕精品一区二区三区| 另类精品久久| 波多野结衣av一区二区av| 国产欧美日韩一区二区三 | 男人添女人高潮全过程视频| 精品少妇黑人巨大在线播放| 纵有疾风起免费观看全集完整版| 久9热在线精品视频| 高清av免费在线| 亚洲精品av麻豆狂野| 脱女人内裤的视频| 久久久久网色| 亚洲国产精品国产精品| 久久精品国产亚洲av涩爱| 久久人妻福利社区极品人妻图片 | 国产亚洲av高清不卡| 在线观看人妻少妇| 捣出白浆h1v1| 精品福利永久在线观看| 99热网站在线观看| 久久影院123| 国产精品99久久99久久久不卡| 欧美日韩黄片免| 少妇精品久久久久久久| www.999成人在线观看| 久久99一区二区三区| 久久精品亚洲av国产电影网| 日韩制服骚丝袜av| 欧美少妇被猛烈插入视频| 啦啦啦中文免费视频观看日本| 日韩 亚洲 欧美在线| 三上悠亚av全集在线观看| 日本91视频免费播放| 黄片小视频在线播放| 成年人免费黄色播放视频| 老鸭窝网址在线观看| 看十八女毛片水多多多| 深夜精品福利| 国产成人精品久久二区二区91| 操出白浆在线播放| 亚洲欧洲日产国产| 热re99久久国产66热| 侵犯人妻中文字幕一二三四区| 亚洲伊人色综图| 国产成人一区二区在线| 国产精品 国内视频| 国产成人91sexporn| 脱女人内裤的视频| 亚洲欧美日韩另类电影网站| 晚上一个人看的免费电影| 另类亚洲欧美激情| 九草在线视频观看| 超色免费av| 日本vs欧美在线观看视频| 青春草亚洲视频在线观看| www日本在线高清视频| 一级黄片播放器| 免费黄频网站在线观看国产| 蜜桃在线观看..| 一级片免费观看大全| 爱豆传媒免费全集在线观看| 91精品国产国语对白视频| 日韩中文字幕欧美一区二区 | av又黄又爽大尺度在线免费看| 99国产精品99久久久久| 亚洲欧美激情在线| 国产淫语在线视频| 免费观看av网站的网址| 久9热在线精品视频| 成人影院久久| 女性被躁到高潮视频| 亚洲少妇的诱惑av| 久久久久精品人妻al黑| 1024视频免费在线观看| 久久久久精品国产欧美久久久 | 欧美国产精品一级二级三级| 精品卡一卡二卡四卡免费| 一二三四在线观看免费中文在| 成年av动漫网址| 欧美变态另类bdsm刘玥| 男女免费视频国产| 国产成人免费观看mmmm| 亚洲伊人久久精品综合| 亚洲情色 制服丝袜| 一本大道久久a久久精品| 亚洲av综合色区一区| 久久狼人影院| 黑丝袜美女国产一区| 欧美人与性动交α欧美精品济南到| 一边亲一边摸免费视频| 亚洲中文日韩欧美视频| 曰老女人黄片| 国产深夜福利视频在线观看| 国产伦理片在线播放av一区| 欧美成狂野欧美在线观看| 久久热在线av| 精品国产超薄肉色丝袜足j| 国产欧美日韩一区二区三 | 深夜精品福利| 一边亲一边摸免费视频| 午夜免费成人在线视频| 日韩av在线免费看完整版不卡| 久久 成人 亚洲| avwww免费| 男人操女人黄网站| 狂野欧美激情性bbbbbb| 成人亚洲欧美一区二区av| 国产视频首页在线观看| 无遮挡黄片免费观看| 国产黄色免费在线视频| 国产1区2区3区精品| 亚洲国产看品久久| 久久久精品国产亚洲av高清涩受| 国产精品免费视频内射| 久久人妻福利社区极品人妻图片 | 精品久久久久久久毛片微露脸 | 午夜日韩欧美国产| 高清视频免费观看一区二区| 成年女人毛片免费观看观看9 | 免费看十八禁软件| 国产精品熟女久久久久浪| 欧美激情 高清一区二区三区| 亚洲国产av影院在线观看| 国产精品熟女久久久久浪| 国产在线视频一区二区| 亚洲欧洲精品一区二区精品久久久| 成年av动漫网址| 久久天躁狠狠躁夜夜2o2o | 美女大奶头黄色视频| 男男h啪啪无遮挡| 一区二区日韩欧美中文字幕| 日韩大片免费观看网站| 婷婷色av中文字幕| 在线看a的网站| 久久久久久久久免费视频了| 欧美变态另类bdsm刘玥| 91麻豆av在线| 亚洲国产精品国产精品| 亚洲美女黄色视频免费看| 午夜福利影视在线免费观看| 亚洲男人天堂网一区| 成人国产av品久久久| 欧美精品一区二区免费开放| 高清欧美精品videossex| 亚洲av美国av| 国产精品 欧美亚洲| 国产爽快片一区二区三区| 中国美女看黄片| 一本综合久久免费| 大香蕉久久网| 精品少妇久久久久久888优播| 精品视频人人做人人爽| 成年女人毛片免费观看观看9 | 成人手机av| 在现免费观看毛片| 男女床上黄色一级片免费看| 亚洲国产av影院在线观看| 亚洲成国产人片在线观看| 久久久久久久久久久久大奶| 19禁男女啪啪无遮挡网站| 成人黄色视频免费在线看| 日本av手机在线免费观看| 另类精品久久| 亚洲,欧美,日韩| 精品国产一区二区久久| 国产在视频线精品| 久久久国产欧美日韩av| 国产成人精品久久久久久| √禁漫天堂资源中文www| 亚洲欧美日韩另类电影网站| 日本黄色日本黄色录像| 国产成人一区二区三区免费视频网站 | 男女下面插进去视频免费观看| 国产精品 欧美亚洲| 天堂俺去俺来也www色官网| 99国产精品免费福利视频| 少妇人妻 视频| av电影中文网址| 欧美黑人精品巨大| 婷婷成人精品国产| 天天躁狠狠躁夜夜躁狠狠躁| 国产女主播在线喷水免费视频网站| 亚洲综合色网址| 手机成人av网站| 一区二区日韩欧美中文字幕| 汤姆久久久久久久影院中文字幕| 国产精品一区二区在线不卡| 在线观看www视频免费| 成年人免费黄色播放视频| 国产一区二区在线观看av| 午夜福利影视在线免费观看| 久久精品熟女亚洲av麻豆精品| 久久天躁狠狠躁夜夜2o2o | 午夜影院在线不卡| 婷婷色麻豆天堂久久| 精品一区二区三区av网在线观看 | 国产精品人妻久久久影院| 好男人电影高清在线观看| 欧美精品一区二区免费开放| 久久精品久久久久久久性| 日韩人妻精品一区2区三区| 中文字幕精品免费在线观看视频| 香蕉丝袜av| 国产成人a∨麻豆精品| 99国产精品免费福利视频| 91麻豆av在线| 下体分泌物呈黄色| 99九九在线精品视频| 999久久久国产精品视频| 亚洲九九香蕉| 国产日韩一区二区三区精品不卡| 飞空精品影院首页| 成年人免费黄色播放视频| 亚洲av成人精品一二三区| 美国免费a级毛片| 国产日韩欧美在线精品| 日韩av在线免费看完整版不卡| 免费在线观看日本一区| 十八禁高潮呻吟视频| 51午夜福利影视在线观看| 男女国产视频网站| 久久久久久免费高清国产稀缺| 午夜福利免费观看在线| av网站在线播放免费| 亚洲熟女毛片儿| 国产免费现黄频在线看| 18禁裸乳无遮挡动漫免费视频| 亚洲精品第二区| 亚洲中文av在线| 少妇的丰满在线观看| 啦啦啦啦在线视频资源| 最近中文字幕2019免费版| 精品一区在线观看国产| 乱人伦中国视频| 欧美精品人与动牲交sv欧美| 下体分泌物呈黄色| 色精品久久人妻99蜜桃| 久久人人爽人人片av| 人成视频在线观看免费观看| 亚洲午夜精品一区,二区,三区| 成年美女黄网站色视频大全免费| 五月天丁香电影| 亚洲欧美精品自产自拍| 亚洲国产精品成人久久小说| 一级片免费观看大全| 亚洲一区中文字幕在线| 国产成人精品久久二区二区91| 欧美性长视频在线观看| 亚洲国产精品一区二区三区在线| 亚洲国产精品成人久久小说| 久久免费观看电影| 又紧又爽又黄一区二区| 亚洲激情五月婷婷啪啪| 亚洲五月色婷婷综合| 日韩av不卡免费在线播放| 欧美精品亚洲一区二区| 久久天躁狠狠躁夜夜2o2o | 香蕉丝袜av| 咕卡用的链子| 天堂中文最新版在线下载| 高清欧美精品videossex| 欧美日韩亚洲综合一区二区三区_| 午夜91福利影院| 夫妻性生交免费视频一级片| 91九色精品人成在线观看| 黑丝袜美女国产一区| av国产久精品久网站免费入址| 成人免费观看视频高清| 男男h啪啪无遮挡| videosex国产| 午夜福利影视在线免费观看| 在线看a的网站| 中文乱码字字幕精品一区二区三区| 亚洲av日韩在线播放| 女性被躁到高潮视频| 高清黄色对白视频在线免费看| 夜夜骑夜夜射夜夜干| 波野结衣二区三区在线| 精品久久久精品久久久| 亚洲av日韩在线播放| 亚洲一区中文字幕在线| 久久精品成人免费网站| 激情视频va一区二区三区| 超色免费av| 热re99久久精品国产66热6| 日本av免费视频播放| 国产成人欧美| 在线观看免费午夜福利视频| av在线app专区| 一本—道久久a久久精品蜜桃钙片| 精品亚洲成a人片在线观看| 两人在一起打扑克的视频| 国产在线一区二区三区精| 操美女的视频在线观看| 久久国产精品男人的天堂亚洲| av网站在线播放免费| 国产主播在线观看一区二区 | 黄色视频不卡| 国产免费一区二区三区四区乱码| 黄频高清免费视频| 国产精品免费视频内射| 国产一区二区三区av在线| 欧美日韩黄片免| 99精品久久久久人妻精品| 中文字幕制服av| 18禁黄网站禁片午夜丰满| 91字幕亚洲| 亚洲欧美成人综合另类久久久| 精品人妻1区二区| 国产精品亚洲av一区麻豆| 欧美在线黄色| av天堂久久9| 人妻一区二区av| 亚洲色图综合在线观看| 国产精品一区二区在线观看99| 人妻人人澡人人爽人人| 男的添女的下面高潮视频| 狂野欧美激情性bbbbbb| 男人操女人黄网站| 精品亚洲乱码少妇综合久久| 大片电影免费在线观看免费| 亚洲第一青青草原| 美女福利国产在线| 亚洲欧美一区二区三区黑人| 一级黄片播放器| av在线播放精品| 免费在线观看完整版高清| 久久精品久久精品一区二区三区| 亚洲成色77777| 丰满少妇做爰视频| 亚洲七黄色美女视频| 欧美国产精品va在线观看不卡| 亚洲av日韩在线播放| 国产精品免费视频内射| 亚洲一卡2卡3卡4卡5卡精品中文| 涩涩av久久男人的天堂| e午夜精品久久久久久久| 久久久久久久国产电影| 国产亚洲精品久久久久5区| 精品少妇黑人巨大在线播放| 十分钟在线观看高清视频www| 久久国产亚洲av麻豆专区| 久热这里只有精品99| 国产精品一区二区在线不卡| 国产亚洲欧美在线一区二区| 午夜福利影视在线免费观看| 久久天躁狠狠躁夜夜2o2o | 老鸭窝网址在线观看| 热re99久久精品国产66热6| 又大又黄又爽视频免费| 丝袜喷水一区| 欧美少妇被猛烈插入视频| 91精品三级在线观看| 男男h啪啪无遮挡| 欧美xxⅹ黑人| 一区二区三区激情视频| 久久久精品国产亚洲av高清涩受| 美女中出高潮动态图| 亚洲国产精品一区三区| 欧美人与性动交α欧美精品济南到| 午夜两性在线视频| 国产极品粉嫩免费观看在线| 自线自在国产av| 国产精品一国产av| 久久精品久久久久久久性| 在线亚洲精品国产二区图片欧美| 国产精品国产三级专区第一集| 秋霞在线观看毛片| 国产亚洲精品久久久久5区| 91精品国产国语对白视频| 91麻豆精品激情在线观看国产 | 在线观看免费视频网站a站| 日本猛色少妇xxxxx猛交久久| 美女中出高潮动态图| 日本vs欧美在线观看视频| 新久久久久国产一级毛片| 性少妇av在线| 午夜视频精品福利| 国产欧美日韩综合在线一区二区| 18禁观看日本| 日本wwww免费看| 黄色 视频免费看| 久久九九热精品免费| 性色av一级| 日韩av在线免费看完整版不卡| 曰老女人黄片| 满18在线观看网站| 咕卡用的链子| 免费高清在线观看日韩| 丝袜美足系列| 久久国产精品大桥未久av| 精品少妇内射三级| 老司机影院毛片| 日日摸夜夜添夜夜爱| 久热这里只有精品99| 日韩中文字幕视频在线看片| 韩国精品一区二区三区| 日本午夜av视频| 少妇粗大呻吟视频| 日韩人妻精品一区2区三区| 水蜜桃什么品种好| 亚洲国产精品一区二区三区在线| 欧美乱码精品一区二区三区| 国产高清不卡午夜福利| 久久 成人 亚洲| 日本黄色日本黄色录像| 50天的宝宝边吃奶边哭怎么回事| 日本黄色日本黄色录像| 亚洲情色 制服丝袜| 欧美黄色片欧美黄色片| 成人18禁高潮啪啪吃奶动态图| 少妇被粗大的猛进出69影院| 精品久久蜜臀av无| 国产精品久久久av美女十八| 777米奇影视久久| 亚洲国产看品久久| 美女扒开内裤让男人捅视频| 精品视频人人做人人爽| 国产精品欧美亚洲77777| 午夜福利一区二区在线看| 伦理电影免费视频| 亚洲第一青青草原| 91老司机精品| 欧美人与善性xxx| 久久99热这里只频精品6学生| 欧美激情极品国产一区二区三区| 免费一级毛片在线播放高清视频 | 婷婷色麻豆天堂久久| 久久人人爽人人片av| 在线观看免费日韩欧美大片| 亚洲熟女毛片儿| 后天国语完整版免费观看| 婷婷成人精品国产| 日本一区二区免费在线视频| 成年av动漫网址| 女人精品久久久久毛片| 日本欧美视频一区| 大片电影免费在线观看免费| 久久精品aⅴ一区二区三区四区| 免费av中文字幕在线| 亚洲国产毛片av蜜桃av| av又黄又爽大尺度在线免费看| 午夜影院在线不卡| 中文字幕人妻熟女乱码| 侵犯人妻中文字幕一二三四区| 亚洲精品自拍成人| 女性被躁到高潮视频| 国产一区二区三区综合在线观看| 欧美精品av麻豆av| 中文字幕人妻丝袜一区二区| 国产精品秋霞免费鲁丝片| 日韩制服丝袜自拍偷拍| 七月丁香在线播放| 午夜影院在线不卡| 91麻豆av在线| 精品人妻一区二区三区麻豆| 一本色道久久久久久精品综合| 欧美97在线视频| 黑人猛操日本美女一级片| 在线观看免费视频网站a站| 多毛熟女@视频| 亚洲精品美女久久久久99蜜臀 | 男女无遮挡免费网站观看| 黑人猛操日本美女一级片| 亚洲视频免费观看视频| 人人澡人人妻人| 国产一区二区三区综合在线观看| 中文字幕av电影在线播放| 色94色欧美一区二区| 丰满少妇做爰视频| 肉色欧美久久久久久久蜜桃| 欧美国产精品va在线观看不卡| cao死你这个sao货| 97人妻天天添夜夜摸| 免费看十八禁软件| 一区二区三区激情视频| 欧美日韩一级在线毛片| av在线老鸭窝| 搡老乐熟女国产| 黄色视频不卡| 女人爽到高潮嗷嗷叫在线视频| 成年美女黄网站色视频大全免费| 欧美人与性动交α欧美软件| 亚洲av日韩在线播放| 麻豆乱淫一区二区| 亚洲一码二码三码区别大吗| 亚洲欧美精品综合一区二区三区| 亚洲激情五月婷婷啪啪| 国产黄色视频一区二区在线观看| 肉色欧美久久久久久久蜜桃| a级毛片黄视频| 亚洲一区二区三区欧美精品| 亚洲精品国产区一区二| 性色av乱码一区二区三区2| av在线播放精品| 悠悠久久av| av福利片在线| 国产爽快片一区二区三区| 午夜老司机福利片| 黑人巨大精品欧美一区二区蜜桃| 18禁裸乳无遮挡动漫免费视频| 青春草视频在线免费观看| 91老司机精品| 日韩熟女老妇一区二区性免费视频| 午夜福利视频精品| 国产精品久久久人人做人人爽| 伦理电影免费视频| 一二三四社区在线视频社区8| 国产免费一区二区三区四区乱码| 男女免费视频国产| 久久人人爽人人片av| 波多野结衣一区麻豆| 乱人伦中国视频| 国产精品人妻久久久影院| 久久精品aⅴ一区二区三区四区| 激情五月婷婷亚洲| 欧美老熟妇乱子伦牲交| 50天的宝宝边吃奶边哭怎么回事| 日韩,欧美,国产一区二区三区| 99久久99久久久精品蜜桃| 男女免费视频国产|