• <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)
    女人爽到高潮嗷嗷叫在线视频| 丝袜美腿诱惑在线| 青青草视频在线视频观看| 成人三级做爰电影| 日本一区二区免费在线视频| a级片在线免费高清观看视频| 丝袜在线中文字幕| 日韩免费高清中文字幕av| 丁香六月天网| 精品一区二区三区av网在线观看 | 一个人免费在线观看的高清视频| 亚洲色图av天堂| www.精华液| 久久精品亚洲av国产电影网| 亚洲成人国产一区在线观看| 黄色 视频免费看| 久久青草综合色| 日本一区二区免费在线视频| 日韩熟女老妇一区二区性免费视频| 老司机靠b影院| 18禁裸乳无遮挡动漫免费视频| 日韩欧美一区视频在线观看| 女性被躁到高潮视频| 午夜老司机福利片| 99re6热这里在线精品视频| 啦啦啦视频在线资源免费观看| 高清黄色对白视频在线免费看| 丁香欧美五月| 久久久久网色| 欧美久久黑人一区二区| 欧美+亚洲+日韩+国产| 亚洲avbb在线观看| 亚洲精品中文字幕一二三四区 | 亚洲精华国产精华精| 午夜福利欧美成人| 激情在线观看视频在线高清 | 日本vs欧美在线观看视频| 国产精品一区二区在线观看99| 韩国精品一区二区三区| 免费日韩欧美在线观看| 狂野欧美激情性xxxx| 一级,二级,三级黄色视频| 一二三四社区在线视频社区8| 日韩一区二区三区影片| 天堂8中文在线网| 国产麻豆69| 男女下面插进去视频免费观看| 成年动漫av网址| 日韩中文字幕欧美一区二区| 免费一级毛片在线播放高清视频 | 国产国语露脸激情在线看| 俄罗斯特黄特色一大片| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品一区二区精品视频观看| 我的亚洲天堂| 精品久久久久久电影网| 国产av精品麻豆| 男女无遮挡免费网站观看| 亚洲三区欧美一区| 美女福利国产在线| 中文字幕色久视频| 国产精品 国内视频| 欧美黑人精品巨大| 不卡一级毛片| 黄色怎么调成土黄色| 99精品久久久久人妻精品| 国产成人精品在线电影| 亚洲成人免费av在线播放| 搡老熟女国产l中国老女人| 99热网站在线观看| 亚洲精品中文字幕一二三四区 | 黄色片一级片一级黄色片| 亚洲人成77777在线视频| 久久精品国产a三级三级三级| 一区二区三区乱码不卡18| 国产97色在线日韩免费| 亚洲精华国产精华精| 午夜福利乱码中文字幕| 亚洲五月色婷婷综合| 久久久久精品国产欧美久久久| 大香蕉久久成人网| 久久毛片免费看一区二区三区| 婷婷成人精品国产| 精品久久久精品久久久| 久久久久精品国产欧美久久久| 精品国产乱码久久久久久小说| 免费观看人在逋| 亚洲精品久久成人aⅴ小说| 久久精品亚洲精品国产色婷小说| 精品久久蜜臀av无| 后天国语完整版免费观看| 欧美日韩视频精品一区| 一进一出抽搐动态| 亚洲第一欧美日韩一区二区三区 | 亚洲成人免费电影在线观看| 日韩制服丝袜自拍偷拍| 欧美精品高潮呻吟av久久| 国产精品久久久人人做人人爽| 女人久久www免费人成看片| 亚洲熟妇熟女久久| 天堂动漫精品| 欧美日本中文国产一区发布| 人人妻人人爽人人添夜夜欢视频| 另类亚洲欧美激情| 欧美老熟妇乱子伦牲交| 精品乱码久久久久久99久播| 国产一区二区三区视频了| xxxhd国产人妻xxx| 精品欧美一区二区三区在线| 国产福利在线免费观看视频| 免费在线观看完整版高清| 捣出白浆h1v1| 男人舔女人的私密视频| 国产精品1区2区在线观看. | 成人手机av| 国产精品 国内视频| 99riav亚洲国产免费| 欧美乱妇无乱码| 成人黄色视频免费在线看| 精品国产一区二区久久| 久9热在线精品视频| 欧美精品一区二区免费开放| 免费看十八禁软件| 老鸭窝网址在线观看| 人成视频在线观看免费观看| 精品免费久久久久久久清纯 | 亚洲久久久国产精品| 国产免费av片在线观看野外av| 巨乳人妻的诱惑在线观看| 一个人免费看片子| 免费看a级黄色片| 国产精品影院久久| 自线自在国产av| 国产欧美日韩一区二区三区在线| 国产在线一区二区三区精| 亚洲专区中文字幕在线| 我要看黄色一级片免费的| 啦啦啦中文免费视频观看日本| 免费久久久久久久精品成人欧美视频| 久久久国产成人免费| 亚洲精华国产精华精| 亚洲av第一区精品v没综合| www.999成人在线观看| 另类亚洲欧美激情| 亚洲第一av免费看| 国产在线一区二区三区精| 在线观看www视频免费| 国产精品国产高清国产av | 亚洲 欧美一区二区三区| 一级片'在线观看视频| 精品久久久久久电影网| 久9热在线精品视频| 十八禁人妻一区二区| 电影成人av| 欧美激情高清一区二区三区| 欧美激情极品国产一区二区三区| 国产91精品成人一区二区三区 | 日本vs欧美在线观看视频| 高清毛片免费观看视频网站 | 欧美 亚洲 国产 日韩一| 女人被躁到高潮嗷嗷叫费观| 12—13女人毛片做爰片一| 成人特级黄色片久久久久久久 | 五月开心婷婷网| 成人手机av| 欧美精品av麻豆av| 在线 av 中文字幕| 飞空精品影院首页| 少妇精品久久久久久久| 男女无遮挡免费网站观看| 午夜福利一区二区在线看| 满18在线观看网站| 在线av久久热| 精品少妇内射三级| 精品国产超薄肉色丝袜足j| 午夜福利乱码中文字幕| 亚洲黑人精品在线| 国产精品秋霞免费鲁丝片| 欧美久久黑人一区二区| 国产视频一区二区在线看| 久久精品aⅴ一区二区三区四区| 国产三级黄色录像| 欧美人与性动交α欧美精品济南到| 中文字幕最新亚洲高清| 欧美 亚洲 国产 日韩一| 日韩免费av在线播放| 亚洲第一青青草原| 久久人妻熟女aⅴ| 国产一区二区在线观看av| 最新美女视频免费是黄的| 色婷婷久久久亚洲欧美| 久久人人97超碰香蕉20202| 欧美黄色片欧美黄色片| 国产av精品麻豆| 欧美日韩一级在线毛片| 人人澡人人妻人| 欧美日韩国产mv在线观看视频| 不卡一级毛片| 国产精品麻豆人妻色哟哟久久| 高清视频免费观看一区二区| 国产成人精品无人区| 日本wwww免费看| 亚洲 欧美一区二区三区| 老熟妇乱子伦视频在线观看| 亚洲美女黄片视频| 亚洲熟妇熟女久久| 国产精品成人在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲专区字幕在线| 十八禁高潮呻吟视频| 人妻一区二区av| a在线观看视频网站| 十八禁网站网址无遮挡| 免费观看av网站的网址| 99香蕉大伊视频| 国产在线精品亚洲第一网站| 欧美日韩精品网址| 啦啦啦免费观看视频1| 欧美日韩亚洲综合一区二区三区_| 麻豆国产av国片精品| 丰满人妻熟妇乱又伦精品不卡| 欧美 亚洲 国产 日韩一| 亚洲全国av大片| 高清视频免费观看一区二区| cao死你这个sao货| 首页视频小说图片口味搜索| 精品国产乱子伦一区二区三区| 国产精品久久久久久精品电影小说| 欧美+亚洲+日韩+国产| www.999成人在线观看| 一边摸一边抽搐一进一小说 | 如日韩欧美国产精品一区二区三区| 制服诱惑二区| 免费av中文字幕在线| 欧美久久黑人一区二区| 美女国产高潮福利片在线看| 999精品在线视频| h视频一区二区三区| 国产午夜精品久久久久久| 亚洲五月婷婷丁香| 欧美日韩亚洲国产一区二区在线观看 | 超碰97精品在线观看| 久久精品国产a三级三级三级| 久久久国产精品麻豆| 免费久久久久久久精品成人欧美视频| 老司机靠b影院| 性色av乱码一区二区三区2| 两个人看的免费小视频| 男女午夜视频在线观看| 国产黄色免费在线视频| 午夜日韩欧美国产| 在线永久观看黄色视频| 久久久国产成人免费| 欧美国产精品va在线观看不卡| 美女高潮喷水抽搐中文字幕| 国产精品亚洲av一区麻豆| 日韩中文字幕欧美一区二区| 国产av精品麻豆| 亚洲欧洲精品一区二区精品久久久| 韩国精品一区二区三区| 中亚洲国语对白在线视频| 在线观看免费高清a一片| 国产又色又爽无遮挡免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 色婷婷久久久亚洲欧美| 亚洲精品美女久久av网站| 精品午夜福利视频在线观看一区 | 最近最新免费中文字幕在线| 亚洲精品av麻豆狂野| 亚洲国产看品久久| 久9热在线精品视频| 欧美激情久久久久久爽电影 | 色老头精品视频在线观看| 十八禁网站免费在线| 777米奇影视久久| 老司机靠b影院| 纵有疾风起免费观看全集完整版| 一级毛片精品| 精品视频人人做人人爽| 欧美日韩一级在线毛片| av免费在线观看网站| 老司机影院毛片| 在线播放国产精品三级| 嫩草影视91久久| 色精品久久人妻99蜜桃| 在线观看舔阴道视频| 国产不卡av网站在线观看| 欧美日韩福利视频一区二区| 看免费av毛片| 成人永久免费在线观看视频 | 国产高清国产精品国产三级| 欧美黄色淫秽网站| 成人国语在线视频| 色婷婷久久久亚洲欧美| 老汉色av国产亚洲站长工具| 亚洲精品乱久久久久久| 精品第一国产精品| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品av麻豆狂野| 欧美精品啪啪一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 国产一区二区激情短视频| 国产单亲对白刺激| 两个人看的免费小视频| 一本综合久久免费| av欧美777| 国产黄频视频在线观看| 亚洲av成人一区二区三| 男人舔女人的私密视频| 岛国毛片在线播放| 成人国语在线视频| 久久精品亚洲精品国产色婷小说| 啦啦啦免费观看视频1| 成人国产一区最新在线观看| 精品久久久久久久毛片微露脸| 大香蕉久久成人网| √禁漫天堂资源中文www| 黄网站色视频无遮挡免费观看| 精品人妻熟女毛片av久久网站| 国产欧美日韩一区二区精品| 精品高清国产在线一区| 丝袜喷水一区| 黑丝袜美女国产一区| 可以免费在线观看a视频的电影网站| 国产精品99久久99久久久不卡| 两个人免费观看高清视频| 色精品久久人妻99蜜桃| 国产精品一区二区免费欧美| 变态另类成人亚洲欧美熟女 | 国产欧美日韩一区二区三| 国产男女超爽视频在线观看| 99热国产这里只有精品6| 一级,二级,三级黄色视频| 两人在一起打扑克的视频| 国产亚洲精品一区二区www | 亚洲av成人不卡在线观看播放网| 老汉色av国产亚洲站长工具| 亚洲一码二码三码区别大吗| 久久精品国产99精品国产亚洲性色 | 狠狠婷婷综合久久久久久88av| 中文字幕人妻熟女乱码| 国精品久久久久久国模美| 亚洲欧美激情在线| 久久中文看片网| 国产成人影院久久av| 久久久久视频综合| 手机成人av网站| 国产精品久久久人人做人人爽| 国产精品 国内视频| 熟女少妇亚洲综合色aaa.| 黄色怎么调成土黄色| 妹子高潮喷水视频| 亚洲精华国产精华精| 久久精品亚洲精品国产色婷小说| 成人免费观看视频高清| 午夜免费鲁丝| 午夜福利一区二区在线看| 日韩大片免费观看网站| 久久久久久久大尺度免费视频| 老鸭窝网址在线观看| 成人永久免费在线观看视频 | 真人做人爱边吃奶动态| 精品少妇黑人巨大在线播放| 午夜成年电影在线免费观看| 亚洲欧美色中文字幕在线| 国产欧美日韩综合在线一区二区| 建设人人有责人人尽责人人享有的| 交换朋友夫妻互换小说| av电影中文网址| 美女高潮到喷水免费观看| 淫妇啪啪啪对白视频| 成年动漫av网址| 成年女人毛片免费观看观看9 | 日韩成人在线观看一区二区三区| 最黄视频免费看| 亚洲精品在线美女| a在线观看视频网站| 欧美日韩av久久| 亚洲精品中文字幕一二三四区 | 久久久久久久久久久久大奶| 国产精品久久久久成人av| 免费高清在线观看日韩| av国产精品久久久久影院| 亚洲男人天堂网一区| 男女高潮啪啪啪动态图| 免费黄频网站在线观看国产| 日本一区二区免费在线视频| 亚洲国产欧美日韩在线播放| 免费看十八禁软件| 国产精品1区2区在线观看. | 后天国语完整版免费观看| 在线观看免费高清a一片| 亚洲 国产 在线| 在线观看免费视频网站a站| 欧美日韩亚洲综合一区二区三区_| 国产成人av教育| 不卡av一区二区三区| 亚洲,欧美精品.| 黄片大片在线免费观看| 无限看片的www在线观看| 国产一卡二卡三卡精品| 99久久精品国产亚洲精品| 757午夜福利合集在线观看| 国产精品二区激情视频| 夜夜夜夜夜久久久久| 国产精品98久久久久久宅男小说| 天堂动漫精品| 大码成人一级视频| 久久久久久久国产电影| 99国产精品一区二区蜜桃av | 国产激情久久老熟女| 曰老女人黄片| 两个人免费观看高清视频| 欧美日韩一级在线毛片| 欧美 亚洲 国产 日韩一| 日本欧美视频一区| 免费人妻精品一区二区三区视频| 757午夜福利合集在线观看| 欧美日韩亚洲综合一区二区三区_| 王馨瑶露胸无遮挡在线观看| 女人久久www免费人成看片| 中文字幕av电影在线播放| 成年人免费黄色播放视频| 50天的宝宝边吃奶边哭怎么回事| 国产伦理片在线播放av一区| 亚洲精品国产色婷婷电影| 久久亚洲精品不卡| 国产亚洲欧美精品永久| 国产日韩欧美在线精品| 777米奇影视久久| 国产成人系列免费观看| 十分钟在线观看高清视频www| 亚洲少妇的诱惑av| 亚洲国产av新网站| 国产亚洲一区二区精品| 亚洲av第一区精品v没综合| 亚洲欧洲精品一区二区精品久久久| 国产高清视频在线播放一区| 亚洲精品自拍成人| 亚洲av片天天在线观看| 国产欧美日韩一区二区三| 免费不卡黄色视频| 不卡av一区二区三区| 欧美变态另类bdsm刘玥| 波多野结衣一区麻豆| 久久 成人 亚洲| 久久狼人影院| 三级毛片av免费| 黑人巨大精品欧美一区二区mp4| 天天添夜夜摸| 亚洲欧洲日产国产| 亚洲专区中文字幕在线| 国产欧美日韩精品亚洲av| 国产福利在线免费观看视频| 国产在线精品亚洲第一网站| 极品少妇高潮喷水抽搐| 啪啪无遮挡十八禁网站| 91大片在线观看| 国产在线一区二区三区精| 久久天躁狠狠躁夜夜2o2o| 精品少妇久久久久久888优播| 国产亚洲精品第一综合不卡| 日韩大码丰满熟妇| 精品少妇内射三级| 久久精品国产亚洲av香蕉五月 | 又紧又爽又黄一区二区| www.999成人在线观看| 亚洲综合色网址| www.自偷自拍.com| 麻豆av在线久日| 91九色精品人成在线观看| 人人妻人人爽人人添夜夜欢视频| 亚洲免费av在线视频| 侵犯人妻中文字幕一二三四区| 国产精品.久久久| 国产高清videossex| 18禁观看日本| 高清毛片免费观看视频网站 | 黑人欧美特级aaaaaa片| 亚洲精品美女久久久久99蜜臀| 中文字幕av电影在线播放| 精品少妇黑人巨大在线播放| 久久久精品94久久精品| 老熟女久久久| 男女边摸边吃奶| 精品久久久精品久久久| 十八禁网站免费在线| 国产亚洲欧美精品永久| 欧美+亚洲+日韩+国产| tocl精华| 欧美精品亚洲一区二区| 久久国产精品大桥未久av| 精品国内亚洲2022精品成人 | 国产亚洲精品一区二区www | 久久精品91无色码中文字幕| 色婷婷久久久亚洲欧美| 精品国产乱码久久久久久男人| 亚洲精品久久午夜乱码| 国产淫语在线视频| 老司机午夜福利在线观看视频 | 老司机影院毛片| xxxhd国产人妻xxx| 蜜桃在线观看..| 国产精品免费大片| 涩涩av久久男人的天堂| 久久九九热精品免费| 水蜜桃什么品种好| 男女之事视频高清在线观看| 亚洲中文日韩欧美视频| 久久av网站| 一级片'在线观看视频| 三上悠亚av全集在线观看| 久久久久久久久久久久大奶| 啦啦啦视频在线资源免费观看| 久久亚洲精品不卡| 国产成人一区二区三区免费视频网站| 国产精品98久久久久久宅男小说| 精品国产一区二区三区久久久樱花| 美女福利国产在线| 欧美在线黄色| 欧美人与性动交α欧美精品济南到| 国产成+人综合+亚洲专区| 1024香蕉在线观看| 国产成人欧美在线观看 | 国产日韩欧美视频二区| 精品久久久久久久毛片微露脸| 成人影院久久| 天天影视国产精品| 99热国产这里只有精品6| 国产精品免费视频内射| 国产精品自产拍在线观看55亚洲 | svipshipincom国产片| 老熟女久久久| 国产色视频综合| 激情在线观看视频在线高清 | 啦啦啦 在线观看视频| 欧美另类亚洲清纯唯美| 成人国产一区最新在线观看| 最新的欧美精品一区二区| 欧美人与性动交α欧美软件| 国产一区二区三区综合在线观看| 国产精品免费大片| 丝瓜视频免费看黄片| 亚洲人成电影观看| 又黄又粗又硬又大视频| 亚洲男人天堂网一区| 亚洲国产欧美在线一区| 自线自在国产av| 视频在线观看一区二区三区| 黑人操中国人逼视频| 99国产精品免费福利视频| 欧美一级毛片孕妇| 久久久久久人人人人人| 十八禁人妻一区二区| 女性被躁到高潮视频| 大香蕉久久网| 欧美久久黑人一区二区| 日韩精品免费视频一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久久精品古装| 亚洲黑人精品在线| netflix在线观看网站| 深夜精品福利| 一级毛片电影观看| 久久久国产欧美日韩av| 如日韩欧美国产精品一区二区三区| 午夜激情av网站| 一夜夜www| 91九色精品人成在线观看| 日韩三级视频一区二区三区| 欧美国产精品一级二级三级| 少妇精品久久久久久久| 美女主播在线视频| 中文亚洲av片在线观看爽 | av有码第一页| 精品第一国产精品| 国产又色又爽无遮挡免费看| 黄色片一级片一级黄色片| 国产无遮挡羞羞视频在线观看| 69av精品久久久久久 | av有码第一页| 91老司机精品| 久久毛片免费看一区二区三区| 亚洲综合色网址| 男女午夜视频在线观看| 欧美日韩一级在线毛片| 久久久久久免费高清国产稀缺| 欧美精品人与动牲交sv欧美| 我要看黄色一级片免费的| 午夜福利在线免费观看网站| 我的亚洲天堂| 久久毛片免费看一区二区三区| 国产亚洲欧美精品永久| 一个人免费看片子| 日本撒尿小便嘘嘘汇集6| 91字幕亚洲| 国产欧美日韩一区二区精品| 中文亚洲av片在线观看爽 | 久久久久久久久久久久大奶| 亚洲va日本ⅴa欧美va伊人久久| 99国产综合亚洲精品| 色综合婷婷激情| 欧美激情极品国产一区二区三区| 高潮久久久久久久久久久不卡| 一进一出好大好爽视频| 国产精品久久电影中文字幕 | 757午夜福利合集在线观看| 亚洲av日韩在线播放| 老司机在亚洲福利影院| 国产黄色免费在线视频| 日韩中文字幕欧美一区二区| 制服诱惑二区| 亚洲国产欧美日韩在线播放| 国产欧美日韩综合在线一区二区| 黑人巨大精品欧美一区二区蜜桃|