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

    Costas序列短波MFSK通信時(shí)頻聯(lián)合同步方法

    2022-08-29 09:55:56李景慧李國軍葉昌榮王尊立
    關(guān)鍵詞:碼元示意圖頻譜

    李景慧,李國軍,葉昌榮,王尊立

    (1.重慶郵電大學(xué) 通信與信息工程學(xué)院,重慶 400065;2.重慶郵電大學(xué) 超視距可信信息傳輸研究所,重慶 400065;3.重慶郵電大學(xué) 光電工程學(xué)院,重慶 400065)

    0 引 言

    在現(xiàn)代數(shù)字通信系統(tǒng)中,同步性能是評(píng)估系統(tǒng)性能的關(guān)鍵指標(biāo),是可靠數(shù)字通信的重要保證。傳統(tǒng)同步方法一般通過檢測(cè)位于數(shù)據(jù)幀頭部的同步碼實(shí)現(xiàn),其性能與同步碼的長(zhǎng)度相關(guān)。同步碼太短,在低信噪比條件下將無法被識(shí)別;而同步碼太長(zhǎng),又將占用過多的信道資源,影響通信效率。因此,為惡劣環(huán)境下的短波通信系統(tǒng)設(shè)計(jì)更好的同步碼和更穩(wěn)定的檢測(cè)算法已經(jīng)成為研究同步的核心問題。實(shí)現(xiàn)同步的主要方法有相關(guān)法[1]、最大似然法(maximum likelihood, ML)[2-3]和似然比檢驗(yàn)(likelihood ratio test, LRT)[4]等3種方法。其中,相關(guān)法具有最低的復(fù)雜度,且最容易實(shí)現(xiàn)[5],但是在低信噪比條件下的性能最差;最大似然法和似然比檢驗(yàn)效果好于相關(guān)法,但計(jì)算比相關(guān)法復(fù)雜[6]。提高相關(guān)法在惡劣環(huán)境下的性能,降低算法實(shí)現(xiàn)難度,對(duì)于提高惡劣環(huán)境下的通信質(zhì)量具有重要意義。

    同步碼中使用較多的是偽隨機(jī)碼(PN碼)[7-8],延遲相關(guān)法和本地相關(guān)法是兩種常用的同步算法。延遲相關(guān)法具有較大的頻率偏移容限,但是在低信噪比和多徑信道下,其自相關(guān)函數(shù)曲線在峰值處變化平緩,無法精確指示同步起始點(diǎn)。本地相關(guān)法具有尖銳的自相關(guān)峰值,準(zhǔn)確度高于延遲相關(guān)法,能精確地指示同步起始點(diǎn),抗噪聲和多徑能力強(qiáng),但是對(duì)頻偏敏感。

    本文提出了一種基于Costas序列的時(shí)頻聯(lián)合同步方法,使用Costas序列作為同步碼,將傳統(tǒng)的滑動(dòng)窗口法擴(kuò)展到時(shí)頻域中,可以同時(shí)估計(jì)時(shí)延和頻偏,具有更好的抗噪聲和抗頻偏性能。

    1 Costas序列特性分析

    1.1 Costas序列定義及構(gòu)造方法

    任意每行每列有且僅有1個(gè)元素等于1、其余元素都為0的n×n階矩陣稱為置換矩陣。Costas序列是一類特殊的置換矩陣,它與自身任意方向的平移副本之間都至多有1個(gè)元素“1”重合,如矩陣A所示(其序列表示為[4,1,6,7,5,8,3,2])。

    (1)

    (2)

    由于置換矩陣A是有限維的,元素ai,j構(gòu)成的序列是有限長(zhǎng)的,因此在非循環(huán)差異函數(shù)C(r,s)中,當(dāng)i+r或j+s超出區(qū)間[1,N]時(shí),取ai+r,j+s為0。非循環(huán)差異函數(shù)滿足(3)式的n×n階置換矩陣稱為n階Costas序列[9]。

    (3)

    Costas序列構(gòu)造方法[9-11]分為計(jì)算機(jī)窮舉搜索法和基于有限域理論的構(gòu)造法兩種方法。計(jì)算機(jī)窮舉搜索法會(huì)窮舉出所有的置換矩陣,并檢測(cè)矩陣是否滿足Costas序列的定義,滿足則輸出該Costas序列。計(jì)算機(jī)窮舉搜索法可以輸出滿足所需階數(shù)的所有Costas序列,但由于窮舉法的復(fù)雜度較高,計(jì)算機(jī)窮舉法所能搜索的階數(shù)是有限制的[9],在構(gòu)造階數(shù)較大的Costas序列時(shí)需要使用基于有限域理論的構(gòu)造方法?;谟邢抻蚶碚摰臉?gòu)造法有Welch構(gòu)造法、Golomb構(gòu)造法和Lempel構(gòu)造法[9-10]等多種方法。

    本文使用8階Costas序列,計(jì)算機(jī)窮舉搜索法可以滿足使用。如無特別說明,本文中出現(xiàn)的Costas序列均由計(jì)算機(jī)窮舉搜索法生成。

    1.2 Costas序列相關(guān)特性分析

    對(duì)于一個(gè)Costas序列([4,1,6,7,5,8,3,2]),稱[4,1,6,7,5,8,3,2]為其序列表示,(1)式為其矩陣表示。Costas序列的自相關(guān)性體現(xiàn)在其矩陣表示上。為了分析矩陣的相關(guān)特性,定義函數(shù)RAB(a,b)為n×n階矩陣A和B的相關(guān)函數(shù)。計(jì)算相關(guān)函數(shù)的方法如圖1所示,矩陣A不動(dòng),矩陣B沿x軸和y軸進(jìn)行移動(dòng),其中a為矩陣B沿x軸移動(dòng)的單位數(shù)(-n+1≤a≤n-1),b為矩陣B沿y軸移動(dòng)的單位數(shù)(-n+1≤b≤n-1),則函數(shù)RAB(a,b)的值為矩陣B移動(dòng)后,兩矩陣中“1”單元格重合的個(gè)數(shù)。當(dāng)矩陣A和矩陣B相等時(shí),稱RAB(a,b)為自相關(guān)函數(shù)[12]。

    圖1 矩陣相關(guān)計(jì)算示意圖

    圖2是一個(gè)8階Costas序列的自相關(guān)函數(shù)計(jì)算結(jié)果,峰值幅度為其階數(shù)n,非峰值處的最大幅值為1。不同Costas序列做同步碼會(huì)表現(xiàn)出不同的性能。圖3和圖4分別為Costas序列1和序列2自相關(guān)函數(shù)示意圖。圖3中,序列1存在“相鄰1”,導(dǎo)致自相關(guān)函數(shù)存在臨峰(序列中相鄰兩個(gè)元素的差值為1的情況稱為“相鄰1”,序列自相關(guān)函數(shù)中與峰值相鄰的不為0的值稱為臨峰),在同步過程中,臨峰會(huì)對(duì)同步造成干擾。而圖4 中Costas序列2由于沒有“相鄰1”的存在,自相關(guān)函數(shù)中不存在臨峰,因此其在同步時(shí)受到的干擾會(huì)小于Costas序列1。

    圖2 8階Costas序列自相關(guān)函數(shù)

    圖3 Costas序列1自相關(guān)函數(shù)示意圖

    圖4 Costas序列2自相關(guān)函數(shù)示意圖

    2 基于Costas序列的MFSK信號(hào)時(shí)頻聯(lián)合同步方法

    2.1 粗同步方法

    Costas序列經(jīng)過FSK調(diào)制后轉(zhuǎn)換到時(shí)頻域,信號(hào)在時(shí)頻域的能量分布符合Costas序列自相關(guān)特性。本文通過離散短時(shí)傅里葉變換(discrete short time Fourier transform, DSTFT)將信號(hào)映射到時(shí)頻域中,如圖5所示。

    圖5 時(shí)頻譜轉(zhuǎn)換框圖

    在沒有信道干擾的情況下,x(t)的表達(dá)式為

    (4)

    (4)式中:A為幅度;g(t)為窗函數(shù);M為FSK調(diào)制階數(shù);M與Costas序列的階數(shù)n相同;Ts為碼元時(shí)間;fi為第i個(gè)碼元符號(hào)的調(diào)制頻率;θi為第i個(gè)碼元的初始相位。對(duì)x(t)采樣后得到x(n),采樣率為fs。將時(shí)域信號(hào)按照碼元時(shí)間Ts進(jìn)行分割,分割后的信號(hào)集合為{x1(k),x2(k),…,xM(k)},其中k為一個(gè)碼元的采樣點(diǎn)數(shù),k=Ts·fs。分別對(duì){x1(k),x2(k),…,xM(k)}進(jìn)行傅里葉變換可得

    Xi(W)=F[xi(k)]

    (5)

    (5)式中,F(xiàn)[·]表示傅里葉變換。將|Xi(W)|進(jìn)行組合得到信號(hào)x(t)的時(shí)頻譜,表示為

    [|X1(W)|,|X2(W)|,…,|XM(W)|]=

    TF[x(t)]

    (6)

    (6)式中,TF[·]表示時(shí)頻譜變換,其結(jié)果為矩陣形式,示意圖如圖6所示。

    圖6 時(shí)頻譜示意圖

    假設(shè)FTi,j為時(shí)頻譜在碼元時(shí)間i、頻率j處的能量值,定義在碼元時(shí)間t、頻率f處的判決量為

    (7)

    (7)式中:n為Costas序列階數(shù);Ci為Costas序列第i位的值。如果在時(shí)頻譜范圍內(nèi)判決量大于判決門限,則認(rèn)為信號(hào)存在(在一般情況下取最大相關(guān)檢測(cè)值的60%—90%作為判決門限)。圖7所示為iturHFMD信道下的判決量仿真結(jié)果,iturHFMD信道參數(shù)如表1所示。在信號(hào)嚴(yán)重衰落的情況下,通過(7)式定義的判決量依然可以準(zhǔn)確指示信號(hào)起止位置。

    圖7 判決量示意圖

    表1 信道參數(shù)

    定義第一個(gè)大于判決門限的自相關(guān)值坐標(biāo)為[tes,fes],則tes·Ts為時(shí)間粗估計(jì)值,fes·fs/N為頻率粗估計(jì)值,其中N為傅里葉變換點(diǎn)數(shù)。

    滑動(dòng)相關(guān)后的時(shí)間估計(jì)精度為Ts,頻率估計(jì)精度為fs/N=1/NΔT,其中ΔT為采樣間隔,NΔT為碼元時(shí)間。提高時(shí)延估計(jì)精度,需要減小碼元時(shí)間Ts,而Ts的減小會(huì)使得頻率估計(jì)精度減小。為了同時(shí)提高時(shí)間估計(jì)精度和頻偏估計(jì)精度,使短時(shí)傅里葉變換中窗函數(shù)的步進(jìn)長(zhǎng)度ΔL為Ts/K,其中K為正整數(shù)。同時(shí)使用過采樣(補(bǔ)零)提高頻譜分辨率,過采樣系數(shù)定義為E。K=4時(shí)的窗函數(shù)步進(jìn)示意圖如圖8所示。

    圖8 DSTFT步進(jìn)示意圖

    圖9為E=2、K=4的參數(shù)條件下,信號(hào)x(t)的時(shí)頻譜。此時(shí)的時(shí)間估計(jì)精度為Ts/K,頻率估計(jì)精度為1/ENΔT。通過這種方式提高估計(jì)精度會(huì)導(dǎo)致計(jì)算量成倍增加。因此,在適當(dāng)調(diào)整參數(shù)K和E后,使用離散頻譜分析方法進(jìn)行精估計(jì)。

    圖9 K=4, E=2的時(shí)頻譜示意圖

    2.2 精同步方法

    由2.1節(jié)的tes可以獲得一個(gè)碼元時(shí)間的頻譜|Xi(W)|,fes為|Xi(W)|的頻譜峰值坐標(biāo),示意圖如圖10所示。

    圖10 基于二倍過采樣的非線性頻偏估計(jì)示意圖

    圖10中,x為經(jīng)過傅里葉變換的頻譜峰值位置,x0為頻偏Δf影響后實(shí)際的頻譜峰值位置。假設(shè)窗函數(shù)為矩形窗,定義其功率譜函數(shù)表達(dá)式為[13]

    (8)

    根據(jù)窗函數(shù)的重心法則,在n趨近于無窮時(shí)有

    (9)

    根據(jù)(8)—(9)式可以計(jì)算出x0的位置坐標(biāo)為[13]

    (10)

    定義k為峰值譜線坐標(biāo),由(10)式可以推導(dǎo)出頻偏估計(jì)值為

    (11)

    最終的頻率估計(jì)為

    (12)

    在估計(jì)時(shí)延之前,使用Δf對(duì)信號(hào)頻譜進(jìn)行糾正。取一個(gè)碼元時(shí)間的信息進(jìn)行分析,對(duì)信號(hào)重新進(jìn)行時(shí)頻譜變換。由實(shí)驗(yàn)數(shù)據(jù)可得,在時(shí)頻譜的時(shí)域取8個(gè)數(shù)據(jù)點(diǎn),數(shù)據(jù)點(diǎn)的值大致服從線性分布,如圖11所示。

    圖11 基于1/4步進(jìn)滑動(dòng)窗口的線性時(shí)延估計(jì)示意圖

    假設(shè)峰值坐標(biāo)為x,次高值坐標(biāo)為x-1。定義x到x+3表示的直線函數(shù)為Y1(x),x-4到x-1表示的直線函數(shù)為Y2(x)。則Y1(x)與Y2(x)定義如下。

    (13)

    數(shù)據(jù)幀實(shí)際開始位置x0可以通過計(jì)算Y1(x)與Y2(x)的交點(diǎn)所得。

    即求解

    Y1(x)=Y2(x)

    (14)

    f(x)=2kx+b1-b2

    (15)

    所得唯一可行解即為x0,因此,時(shí)延精估計(jì)值為

    Δt=|x-x0|Ts/K

    (16)

    (16)式中:Ts為碼元時(shí)間;K為窗函數(shù)步進(jìn)系數(shù)。

    2.3 窗函數(shù)影響分析

    DSTFT中窗函數(shù)的不同會(huì)對(duì)同步性能產(chǎn)生影響。令窗函數(shù)步進(jìn)長(zhǎng)度ΔL為Ts,假設(shè)f1為碼元1所代表的頻率,f2為碼元2所代表的頻率。當(dāng)存在時(shí)延時(shí),第一個(gè)碼元符號(hào)的部分信息會(huì)延后到x2(k)中。

    x2(n)=A1sin[2πf1(n-τ)]+A2sin[2πf2(n-τ)]

    (17)

    (17)式中:A1表示第一個(gè)碼元時(shí)間的幅度;A2表示第二個(gè)碼元時(shí)間的幅度;τ為時(shí)間延遲。此時(shí)對(duì)x2(k)進(jìn)行傅里葉變換后,頻率f1所造成的頻譜泄漏現(xiàn)象會(huì)極大地干擾頻率f2的主瓣譜線,如圖12所示。這種會(huì)嚴(yán)重影響頻偏估計(jì)的干擾,可以通過選擇不同的窗函數(shù)或者全相位傅里葉變換來減小。

    圖12 時(shí)間延遲分析

    3 仿真驗(yàn)證

    仿真過程使用8階Costas序列[1,5,7,4,2,8,3,6]作為同步碼,窗函數(shù)為矩形窗函數(shù),使用MFSK(M=8)方式對(duì)傳輸序列進(jìn)行調(diào)制,采樣率設(shè)置為12 kHz,為2倍過采樣,窗函數(shù)步進(jìn)長(zhǎng)度為1/4碼元時(shí)間。根據(jù)國際電信聯(lián)盟給出的Watterson模型建議參數(shù)[14],選擇中緯度地區(qū)的3種短波信道仿真。3種信道參數(shù)如表1所示。

    圖13—圖15給出了本文方法在不同判決門限下虛警率和檢測(cè)率的ROC曲線。在AWGN信道條件下,隨著判決門限的逐漸增大,虛警率不斷降低,但檢測(cè)率變化很小。在短波信道中,不添加噪聲時(shí)的檢測(cè)性能與添加-20 dB噪聲干擾的檢測(cè)性能相近。綜上所述,本文同步方法對(duì)閾值不敏感,且具有優(yōu)秀的抗噪聲性能,雖在短波信道中的性能有所下降,但依然可以保證其在短波信道的可用性。

    圖13 ROC曲線(AWGN信道)

    圖14 ROC曲線(iturHF信道)

    圖15 ROC曲線(iturHF信道, -20 dB)

    以最大相關(guān)檢測(cè)值的70%作為判決門限,本文同步方法與傳統(tǒng)PN碼同步方法在iturHFMQ信道下的檢測(cè)性能如圖16所示。傳統(tǒng)的PN碼同步方法,在信噪比逐漸降低時(shí)性能出現(xiàn)陡降;本文同步方法在信噪比為-20 dB時(shí),依然具有可觀的檢測(cè)率,可以有效提高惡劣環(huán)境下短波通信的質(zhì)量。

    圖16 基于Costas序列同步方法的漏檢率

    圖17—圖18為本文同步方法在頻偏影響下的自相關(guān)峰值變化圖,仿真信道為iturHFMQ。由于自相關(guān)函數(shù)是關(guān)于時(shí)間和頻率的參數(shù),需要從兩個(gè)角度來分析其峰值變化,圖17為頻域角度的自相關(guān)峰值變化圖,圖18為時(shí)域角度的自相關(guān)峰值變化圖。本文同步方法的自相關(guān)峰值幅度主要受頻譜泄漏的影響,當(dāng)頻偏越接近1/2頻譜分辨率的奇數(shù)倍時(shí),影響越大。從時(shí)域分析,自相關(guān)峰值附近的值相對(duì)較低,出現(xiàn)了一個(gè)凹陷,驗(yàn)證了本文1.2節(jié)中所描述的“相鄰1”的影響。在小頻偏影響下,相關(guān)結(jié)果峰值幾乎不受頻偏的影響,在頻偏值造成的影響較大時(shí)(接近1/2頻譜分辨率的奇數(shù)倍),其相關(guān)結(jié)果峰值依然可以近似達(dá)到無頻偏時(shí)的70%左右。對(duì)比文獻(xiàn)[15]中的方法,本文同步方法可以有效減小大頻偏對(duì)同步造成的影響。

    圖17 同步相關(guān)峰示意圖(頻域)

    圖18 同步相關(guān)峰示意圖(時(shí)域)

    圖19為本文同步方法與PN碼相關(guān)法的定時(shí)精度對(duì)比。圖19中,PN碼長(zhǎng)度與本文使用的同步碼長(zhǎng)度相同(1.28 s)。由圖19可以看出,本文所提出的同步方法在信噪比較低時(shí)具有明顯的優(yōu)勢(shì)。

    圖19 定時(shí)誤差分析

    由于使用了離散頻譜分析方法進(jìn)行精估計(jì),本文同步方法的精度在低信噪比情況下呈現(xiàn)出明顯的周期性。設(shè)置信噪比為-17 dB,圖20和圖21為本文同步方法在iturHFMQ信道中的時(shí)間估計(jì)誤差和頻率估計(jì)誤差仿真圖。在低信噪比短波信道條件下,本文同步方法依然具有很高的估計(jì)精度,但估計(jì)誤差出現(xiàn)周期性變化的現(xiàn)象。這是由于在精同步過程中,依賴頻譜或自相關(guān)函數(shù)中最高值坐標(biāo)和次高值坐標(biāo)的判定,當(dāng)時(shí)間延遲或者頻率偏移接近其歸一化值-0.5、0、+0.5時(shí)會(huì)出現(xiàn)最高值坐標(biāo)或次高值坐標(biāo)判定錯(cuò)誤的情況。以歸一化頻率偏差接近±0.5為例,其峰值譜線的兩側(cè)譜線在實(shí)際中趨近于相等,由于受到噪聲的干擾,其值相對(duì)于真實(shí)值來說出現(xiàn)波動(dòng),會(huì)造成次高值判定錯(cuò)誤的情況,因此,估計(jì)誤差會(huì)出現(xiàn)一個(gè)周期性的變化(在本文的仿真條件中,頻譜分辨率為3.125 Hz,時(shí)間分辨率為0.04 s,歸一化是指一個(gè)分辨率)。

    圖20 時(shí)間估計(jì)誤差分析

    圖21 頻率估計(jì)誤差分析

    4 結(jié)束語

    本文提出了一種基于Costas序列的時(shí)頻聯(lián)合同步方法,經(jīng)仿真驗(yàn)證,該方法與基于PN碼的同步方法相比,具有更好的抗噪聲和頻偏能力,可以有效提高短波通信中的捕獲率和同步精度,且相同長(zhǎng)度下同步性能更好,更加適用于短波通信,可以有效提高惡劣信道條件下的同步性能,為提高短波通信的可靠性奠定基礎(chǔ)。

    猜你喜歡
    碼元示意圖頻譜
    一種用于深空探測(cè)的Chirp變換頻譜分析儀設(shè)計(jì)與實(shí)現(xiàn)
    LFM-BPSK復(fù)合調(diào)制參數(shù)快速估計(jì)及碼元恢復(fù)
    先畫示意圖再解答問題
    黔西南州旅游示意圖
    一種基于稀疏度估計(jì)的自適應(yīng)壓縮頻譜感知算法
    基于極大似然準(zhǔn)則的短猝發(fā)信號(hào)盲解調(diào)
    認(rèn)知無線電頻譜感知技術(shù)綜述
    兩張圖讀懂“青年之聲”
    一種基于功率限制下的認(rèn)知無線電的頻譜感知模型
    一種碼元同步時(shí)鐘信號(hào)的提取方法及單片機(jī)實(shí)現(xiàn)
    亚洲熟女精品中文字幕| 免费在线观看黄色视频的| 色94色欧美一区二区| 热99久久久久精品小说推荐| 欧美老熟妇乱子伦牲交| 黄色视频不卡| 叶爱在线成人免费视频播放| 色94色欧美一区二区| 大话2 男鬼变身卡| 亚洲欧美清纯卡通| 啦啦啦在线观看免费高清www| 狠狠婷婷综合久久久久久88av| 男女床上黄色一级片免费看| 国产精品.久久久| 晚上一个人看的免费电影| 在线免费观看不下载黄p国产| 一二三四在线观看免费中文在| 国产亚洲av高清不卡| 亚洲,一卡二卡三卡| 精品国产国语对白av| av不卡在线播放| 曰老女人黄片| 建设人人有责人人尽责人人享有的| 久久久精品区二区三区| 欧美日韩亚洲高清精品| 中文乱码字字幕精品一区二区三区| 人妻人人澡人人爽人人| 极品人妻少妇av视频| 久久午夜综合久久蜜桃| 国语对白做爰xxxⅹ性视频网站| 无限看片的www在线观看| 免费黄色在线免费观看| 一边摸一边做爽爽视频免费| www.自偷自拍.com| 考比视频在线观看| 黄色视频不卡| 亚洲一码二码三码区别大吗| 国产黄色免费在线视频| 男人操女人黄网站| 国产精品亚洲av一区麻豆 | 色精品久久人妻99蜜桃| xxxhd国产人妻xxx| 亚洲,一卡二卡三卡| 欧美乱码精品一区二区三区| 五月天丁香电影| 免费看av在线观看网站| 十八禁高潮呻吟视频| 亚洲伊人色综图| 精品一品国产午夜福利视频| 性高湖久久久久久久久免费观看| 女性生殖器流出的白浆| 国产精品 国内视频| 国产精品秋霞免费鲁丝片| 亚洲av成人精品一二三区| 久久久久国产一级毛片高清牌| 国产精品欧美亚洲77777| 看免费av毛片| 夜夜骑夜夜射夜夜干| 男人爽女人下面视频在线观看| 黄片播放在线免费| 嫩草影院入口| 最近手机中文字幕大全| 亚洲欧美一区二区三区国产| 国产无遮挡羞羞视频在线观看| 国产精品久久久久久精品古装| 老汉色∧v一级毛片| 国产成人系列免费观看| 蜜桃在线观看..| 久久女婷五月综合色啪小说| 男女边摸边吃奶| 久久久精品国产亚洲av高清涩受| 亚洲综合色网址| 精品第一国产精品| 美女福利国产在线| 五月天丁香电影| 韩国av在线不卡| 少妇人妻精品综合一区二区| 中文字幕高清在线视频| 亚洲国产最新在线播放| 欧美国产精品va在线观看不卡| 亚洲精品成人av观看孕妇| 黄色 视频免费看| 丝袜喷水一区| 欧美老熟妇乱子伦牲交| 久久国产精品大桥未久av| 男的添女的下面高潮视频| 美女脱内裤让男人舔精品视频| 男的添女的下面高潮视频| 91aial.com中文字幕在线观看| 成人手机av| 欧美日韩亚洲国产一区二区在线观看 | 成人黄色视频免费在线看| 人成视频在线观看免费观看| 国产一区二区三区av在线| 别揉我奶头~嗯~啊~动态视频 | 久久国产精品男人的天堂亚洲| 美女午夜性视频免费| 人体艺术视频欧美日本| 亚洲一码二码三码区别大吗| 天天操日日干夜夜撸| 91国产中文字幕| www.精华液| 中文字幕精品免费在线观看视频| 日韩一区二区三区影片| 免费看av在线观看网站| 国产精品久久久av美女十八| 在线精品无人区一区二区三| 久久99一区二区三区| 日韩制服骚丝袜av| 久久久久人妻精品一区果冻| 天天躁日日躁夜夜躁夜夜| 不卡av一区二区三区| 黄色 视频免费看| 日韩欧美一区视频在线观看| 90打野战视频偷拍视频| 亚洲av中文av极速乱| 又大又黄又爽视频免费| 日韩av免费高清视频| 一二三四中文在线观看免费高清| 九色亚洲精品在线播放| 免费黄网站久久成人精品| 日韩欧美一区视频在线观看| 免费女性裸体啪啪无遮挡网站| 黄频高清免费视频| 一区在线观看完整版| 久久久国产一区二区| 毛片一级片免费看久久久久| 国产精品99久久99久久久不卡 | 香蕉国产在线看| 自拍欧美九色日韩亚洲蝌蚪91| 成人漫画全彩无遮挡| 丝袜美足系列| 建设人人有责人人尽责人人享有的| 如何舔出高潮| 久久久欧美国产精品| 90打野战视频偷拍视频| 国产黄频视频在线观看| 大片电影免费在线观看免费| av一本久久久久| 丝袜美足系列| 999久久久国产精品视频| 国产一卡二卡三卡精品 | 日韩 欧美 亚洲 中文字幕| 久久狼人影院| 黄色视频不卡| 如日韩欧美国产精品一区二区三区| 2018国产大陆天天弄谢| 亚洲综合精品二区| 亚洲欧美日韩另类电影网站| 日韩 欧美 亚洲 中文字幕| 最近中文字幕2019免费版| 国产熟女午夜一区二区三区| 久久热在线av| 一级毛片我不卡| 国产一级毛片在线| 国产一卡二卡三卡精品 | 亚洲国产中文字幕在线视频| 啦啦啦视频在线资源免费观看| 五月天丁香电影| 欧美精品高潮呻吟av久久| 国产免费现黄频在线看| 国产 一区精品| 欧美日韩福利视频一区二区| 夫妻午夜视频| 亚洲精品国产av成人精品| 悠悠久久av| 精品人妻在线不人妻| 纵有疾风起免费观看全集完整版| 国产国语露脸激情在线看| 日本av手机在线免费观看| 亚洲自偷自拍图片 自拍| 免费黄频网站在线观看国产| 日韩,欧美,国产一区二区三区| av在线观看视频网站免费| av天堂久久9| 丁香六月欧美| 欧美亚洲日本最大视频资源| 涩涩av久久男人的天堂| 亚洲国产精品999| 久久国产精品大桥未久av| 亚洲精品自拍成人| 精品一区二区三区四区五区乱码 | av免费观看日本| 国产亚洲最大av| 老司机亚洲免费影院| 国产精品麻豆人妻色哟哟久久| 亚洲第一青青草原| 国产精品久久久久久人妻精品电影 | 天堂中文最新版在线下载| 操出白浆在线播放| 9色porny在线观看| 超碰97精品在线观看| 国产色婷婷99| 51午夜福利影视在线观看| 哪个播放器可以免费观看大片| 男女高潮啪啪啪动态图| 精品久久久精品久久久| 午夜福利视频在线观看免费| 亚洲成色77777| 国产精品国产三级专区第一集| 色综合欧美亚洲国产小说| 王馨瑶露胸无遮挡在线观看| 亚洲精品成人av观看孕妇| 免费看不卡的av| 欧美在线一区亚洲| 国产亚洲欧美精品永久| 国产精品嫩草影院av在线观看| av在线app专区| 一级爰片在线观看| 亚洲成色77777| 日韩不卡一区二区三区视频在线| 日日爽夜夜爽网站| 久久人人97超碰香蕉20202| 在线观看一区二区三区激情| 久久女婷五月综合色啪小说| 亚洲国产av新网站| 亚洲国产最新在线播放| 自线自在国产av| 午夜福利在线免费观看网站| 男人操女人黄网站| 日本av手机在线免费观看| 国产精品免费视频内射| 国产又色又爽无遮挡免| 亚洲一码二码三码区别大吗| 日本猛色少妇xxxxx猛交久久| 水蜜桃什么品种好| 少妇精品久久久久久久| 亚洲欧美精品自产自拍| 五月开心婷婷网| 性少妇av在线| 亚洲欧美清纯卡通| av在线app专区| 亚洲七黄色美女视频| 亚洲精品乱久久久久久| 日本午夜av视频| 久久影院123| 老司机在亚洲福利影院| 国产av一区二区精品久久| 国产精品国产三级专区第一集| 伦理电影免费视频| 建设人人有责人人尽责人人享有的| 午夜福利乱码中文字幕| 肉色欧美久久久久久久蜜桃| 国产一区亚洲一区在线观看| 久久国产亚洲av麻豆专区| 精品亚洲乱码少妇综合久久| 男女之事视频高清在线观看 | 色精品久久人妻99蜜桃| 制服诱惑二区| 午夜日韩欧美国产| 婷婷成人精品国产| 日韩熟女老妇一区二区性免费视频| 纵有疾风起免费观看全集完整版| 中文精品一卡2卡3卡4更新| 五月天丁香电影| 日韩精品免费视频一区二区三区| 两性夫妻黄色片| 国产一区二区三区综合在线观看| 少妇被粗大猛烈的视频| 亚洲 欧美一区二区三区| 国产日韩欧美亚洲二区| 国产成人精品福利久久| 国产精品偷伦视频观看了| 丁香六月欧美| 亚洲图色成人| 亚洲第一青青草原| av免费观看日本| 两个人看的免费小视频| 亚洲av综合色区一区| 啦啦啦 在线观看视频| 成人国语在线视频| 色播在线永久视频| 欧美av亚洲av综合av国产av | 午夜91福利影院| 男人添女人高潮全过程视频| 欧美激情高清一区二区三区 | 国产精品麻豆人妻色哟哟久久| 黄频高清免费视频| av一本久久久久| 国产精品久久久久久久久免| 十八禁网站网址无遮挡| 日本vs欧美在线观看视频| 午夜激情av网站| 十八禁网站网址无遮挡| 一本—道久久a久久精品蜜桃钙片| 美女国产高潮福利片在线看| 国产av一区二区精品久久| 国产精品无大码| 午夜免费鲁丝| 黑人巨大精品欧美一区二区蜜桃| 免费高清在线观看日韩| 大片电影免费在线观看免费| 国产 精品1| 午夜免费男女啪啪视频观看| 一区二区三区精品91| 欧美久久黑人一区二区| 中文乱码字字幕精品一区二区三区| 成年动漫av网址| 亚洲精品视频女| 久久人人爽av亚洲精品天堂| 久久久久久久国产电影| 男人添女人高潮全过程视频| 欧美亚洲 丝袜 人妻 在线| 男女边吃奶边做爰视频| 国产亚洲一区二区精品| 考比视频在线观看| 久久人妻熟女aⅴ| 日韩大片免费观看网站| 久久久久精品性色| 两个人免费观看高清视频| 18在线观看网站| 国产成人av激情在线播放| 晚上一个人看的免费电影| 亚洲在久久综合| 69精品国产乱码久久久| 免费观看性生交大片5| 91成人精品电影| 一级,二级,三级黄色视频| 久久国产亚洲av麻豆专区| 精品一区二区三卡| 午夜日本视频在线| 看免费av毛片| 国产爽快片一区二区三区| 在线看a的网站| 青青草视频在线视频观看| 午夜精品国产一区二区电影| 大话2 男鬼变身卡| 美女扒开内裤让男人捅视频| 妹子高潮喷水视频| 大香蕉久久网| 精品午夜福利在线看| 激情视频va一区二区三区| videosex国产| 国产97色在线日韩免费| 国产黄色视频一区二区在线观看| 国产伦理片在线播放av一区| av免费观看日本| 最近手机中文字幕大全| 男人添女人高潮全过程视频| 国产一区二区三区av在线| 91精品国产国语对白视频| 国产成人系列免费观看| 成人手机av| 2021少妇久久久久久久久久久| 一个人免费看片子| 别揉我奶头~嗯~啊~动态视频 | 精品一区二区三卡| 看非洲黑人一级黄片| 日本欧美国产在线视频| 宅男免费午夜| 看非洲黑人一级黄片| 亚洲四区av| 亚洲国产精品999| 桃花免费在线播放| 超碰97精品在线观看| 亚洲成人av在线免费| 色吧在线观看| 欧美成人午夜精品| 精品人妻在线不人妻| 性少妇av在线| 在线观看www视频免费| 国产伦人伦偷精品视频| 中文字幕最新亚洲高清| www.自偷自拍.com| 欧美日韩成人在线一区二区| 亚洲综合精品二区| 欧美老熟妇乱子伦牲交| 狂野欧美激情性xxxx| 黄网站色视频无遮挡免费观看| 90打野战视频偷拍视频| 久久久精品94久久精品| 三上悠亚av全集在线观看| 日韩一本色道免费dvd| 国产av国产精品国产| 美女中出高潮动态图| 亚洲美女视频黄频| 日韩 亚洲 欧美在线| 亚洲男人天堂网一区| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久久成人av| 19禁男女啪啪无遮挡网站| 2018国产大陆天天弄谢| 人人妻,人人澡人人爽秒播 | 精品酒店卫生间| 七月丁香在线播放| 一区二区三区激情视频| 欧美精品一区二区大全| 制服人妻中文乱码| 午夜影院在线不卡| 中文字幕人妻熟女乱码| 国产精品99久久99久久久不卡 | netflix在线观看网站| 宅男免费午夜| 一本色道久久久久久精品综合| 一级a爱视频在线免费观看| 丝袜美腿诱惑在线| 性色av一级| 免费观看人在逋| 丁香六月欧美| 人人妻人人爽人人添夜夜欢视频| 国产成人精品在线电影| 男女高潮啪啪啪动态图| 在线观看一区二区三区激情| 成人黄色视频免费在线看| 亚洲欧美清纯卡通| 精品人妻熟女毛片av久久网站| 婷婷色麻豆天堂久久| 色播在线永久视频| 欧美老熟妇乱子伦牲交| xxx大片免费视频| 亚洲欧美日韩另类电影网站| 免费在线观看完整版高清| 纵有疾风起免费观看全集完整版| 91精品伊人久久大香线蕉| 男人添女人高潮全过程视频| 国产成人91sexporn| 啦啦啦啦在线视频资源| 永久免费av网站大全| 卡戴珊不雅视频在线播放| 亚洲av综合色区一区| 日韩av不卡免费在线播放| 亚洲国产欧美网| 不卡视频在线观看欧美| 国产亚洲欧美精品永久| 午夜福利视频在线观看免费| 亚洲欧美成人综合另类久久久| 美女国产高潮福利片在线看| 丝袜美腿诱惑在线| 男男h啪啪无遮挡| 国产成人一区二区在线| 国产成人欧美在线观看 | 99久久精品国产亚洲精品| 婷婷色av中文字幕| 中文天堂在线官网| 欧美乱码精品一区二区三区| 日本黄色日本黄色录像| av一本久久久久| 男女无遮挡免费网站观看| 最新在线观看一区二区三区 | 国产精品欧美亚洲77777| 国产亚洲精品第一综合不卡| 国产一区二区三区av在线| 在线 av 中文字幕| 国产男女超爽视频在线观看| 看免费av毛片| 人妻 亚洲 视频| 18禁裸乳无遮挡动漫免费视频| 女的被弄到高潮叫床怎么办| 视频在线观看一区二区三区| 一本大道久久a久久精品| 成人免费观看视频高清| 91老司机精品| 亚洲,欧美精品.| 国产高清国产精品国产三级| 亚洲精品久久成人aⅴ小说| 色婷婷av一区二区三区视频| 不卡av一区二区三区| 国产精品免费视频内射| 在线看a的网站| 亚洲欧美清纯卡通| 欧美日韩综合久久久久久| 亚洲国产欧美在线一区| 久久精品久久久久久久性| av.在线天堂| 黄色毛片三级朝国网站| 欧美日本中文国产一区发布| 成人漫画全彩无遮挡| av在线app专区| 国产麻豆69| 99久久人妻综合| 亚洲天堂av无毛| 另类精品久久| 汤姆久久久久久久影院中文字幕| 免费日韩欧美在线观看| 久久人妻熟女aⅴ| 亚洲四区av| 黑丝袜美女国产一区| 亚洲国产最新在线播放| 19禁男女啪啪无遮挡网站| 久久青草综合色| 中文字幕人妻丝袜一区二区 | 国产xxxxx性猛交| 中国三级夫妇交换| 涩涩av久久男人的天堂| 99热网站在线观看| 韩国精品一区二区三区| 亚洲色图 男人天堂 中文字幕| 侵犯人妻中文字幕一二三四区| 国产av精品麻豆| 亚洲欧美一区二区三区久久| 18禁国产床啪视频网站| 精品一区在线观看国产| 国产午夜精品一二区理论片| 可以免费在线观看a视频的电影网站 | 欧美在线黄色| 国产精品人妻久久久影院| 日日啪夜夜爽| 男人操女人黄网站| 亚洲精品国产区一区二| 成人手机av| 亚洲成人一二三区av| 制服诱惑二区| 亚洲国产欧美日韩在线播放| 午夜av观看不卡| 国产精品久久久久久久久免| 狠狠婷婷综合久久久久久88av| 午夜免费鲁丝| 国产精品二区激情视频| 久久青草综合色| 国产成人精品福利久久| 久久 成人 亚洲| 欧美激情极品国产一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 美女午夜性视频免费| 美女扒开内裤让男人捅视频| 午夜免费观看性视频| 欧美少妇被猛烈插入视频| 人妻 亚洲 视频| 午夜久久久在线观看| h视频一区二区三区| 视频区图区小说| 黑人巨大精品欧美一区二区蜜桃| 免费观看a级毛片全部| 秋霞在线观看毛片| 国产成人精品无人区| av女优亚洲男人天堂| 电影成人av| 日韩一区二区三区影片| 亚洲精品久久久久久婷婷小说| 少妇人妻精品综合一区二区| 老鸭窝网址在线观看| 99热全是精品| 麻豆av在线久日| 天堂俺去俺来也www色官网| 亚洲精品aⅴ在线观看| 亚洲自偷自拍图片 自拍| 亚洲精品aⅴ在线观看| 高清欧美精品videossex| 欧美xxⅹ黑人| 欧美日韩av久久| 免费高清在线观看视频在线观看| 亚洲熟女毛片儿| 少妇被粗大猛烈的视频| 亚洲国产精品国产精品| 日韩一卡2卡3卡4卡2021年| 丝袜喷水一区| 免费不卡黄色视频| av在线app专区| 亚洲综合色网址| 亚洲第一av免费看| 亚洲精品乱久久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产精品国产精品| 国产精品秋霞免费鲁丝片| 国产精品麻豆人妻色哟哟久久| 亚洲国产欧美在线一区| 极品少妇高潮喷水抽搐| 国产亚洲欧美精品永久| 蜜桃在线观看..| 黑人欧美特级aaaaaa片| 久久鲁丝午夜福利片| 精品一区在线观看国产| 欧美激情 高清一区二区三区| 亚洲国产精品一区三区| 成年人免费黄色播放视频| 国产av国产精品国产| 久久精品人人爽人人爽视色| 精品国产国语对白av| 国产精品久久久久成人av| 久久久亚洲精品成人影院| 欧美久久黑人一区二区| 亚洲精品成人av观看孕妇| 男女边摸边吃奶| 最近的中文字幕免费完整| 欧美最新免费一区二区三区| 日本av手机在线免费观看| 一区福利在线观看| 国产激情久久老熟女| 精品亚洲成a人片在线观看| 色婷婷av一区二区三区视频| 中文字幕高清在线视频| 九草在线视频观看| 丰满迷人的少妇在线观看| 99国产精品免费福利视频| 久久ye,这里只有精品| 亚洲精华国产精华液的使用体验| 久久久国产一区二区| 亚洲精品久久成人aⅴ小说| 成年女人毛片免费观看观看9 | 国产乱人偷精品视频| 在线观看免费午夜福利视频| 老司机靠b影院| 欧美精品一区二区大全| 亚洲自偷自拍图片 自拍| av网站免费在线观看视频| 老司机在亚洲福利影院| 满18在线观看网站| 国产精品嫩草影院av在线观看| 肉色欧美久久久久久久蜜桃| 一级毛片我不卡| 日韩熟女老妇一区二区性免费视频| 国产免费视频播放在线视频| 午夜福利免费观看在线| 亚洲欧美一区二区三区久久| 色网站视频免费| 日日爽夜夜爽网站| 精品国产露脸久久av麻豆| 国产精品熟女久久久久浪| 你懂的网址亚洲精品在线观看| 欧美日韩国产mv在线观看视频| 久久女婷五月综合色啪小说| 丰满少妇做爰视频| 欧美激情高清一区二区三区 | 黄片小视频在线播放|