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

    高質(zhì)量LMSCT時頻分析算法及其在雷達(dá)信號目標(biāo)檢測中的應(yīng)用

    2022-03-08 02:38:06郝國成張必超張雅冰石光耀王盼盼
    上海交通大學(xué)學(xué)報 2022年2期
    關(guān)鍵詞:時頻雜波分析方法

    郝國成, 張必超, 鍋 娟, 張雅冰, 石光耀, 王盼盼, 張 薇

    (1. 中國地質(zhì)大學(xué)(武漢) 機(jī)械與電子信息學(xué)院,武漢 430074; 2. 杜克大學(xué) 數(shù)學(xué)系,美國Durham 27708; 3. 中國地質(zhì)大學(xué)(武漢) 復(fù)雜系統(tǒng)先進(jìn)控制與智能自動化湖北重點實驗室, 武漢 430074)

    時頻分析是處理非平穩(wěn)信號的重要工具之一,能夠?qū)⒁痪S時序信號變換到時間頻率坐標(biāo)軸并提供信號頻率隨時間變化的聯(lián)合分布關(guān)系[1].目前,時頻分析方法已經(jīng)廣泛應(yīng)用于地震信號[2-4]、旋轉(zhuǎn)機(jī)械故障信號[5]、金屬破裂信號[6]等.海雜波信號也具有非平穩(wěn)隨機(jī)的特征,如何提高海面小目標(biāo)(如浮冰、小船、蛙人以及飛機(jī)殘骸)的檢測能力是一項長期且艱巨的任務(wù)[7].傳統(tǒng)時頻分析方法包括短時Fourier變換(STFT)[8]、小波變換(WT)[9]及Wigner-Ville分布(WVD)[10]等.其中,STFT和WT屬于線性時頻分析方法,根據(jù)海森堡不確定準(zhǔn)則,線性時頻方法不能同時實現(xiàn)時間和頻率的高分辨率.二次型時頻分析方法WVD在處理單分量信號時能夠得到高分辨率的時頻分布,但是在處理多分量信號會出現(xiàn)交叉項干擾,這限制了WVD在現(xiàn)實生活中的應(yīng)用.

    近些年,針對傳統(tǒng)時頻分析方法的不足,許多優(yōu)秀的時頻分析方法被提出以解決上述的問題.根據(jù)選擇的窗函數(shù)如果和待分析信號的內(nèi)部特征相近,就會得到聚集度較高的時頻表示這種思想.參數(shù)化時頻分析方法朝著識別信號的內(nèi)部特征和構(gòu)造窗函數(shù)的方向發(fā)展.線調(diào)頻小波變換(CT)[11]通過引入額外的調(diào)頻率(CR)參數(shù),在處理線性調(diào)制信號時能夠得到聚集度較高的時頻分布.但是,當(dāng)CR與信號的瞬時頻率差別較大時,通過該方法獲得的時頻表示存在時頻能量模糊現(xiàn)象,不利于信號時頻特征的提取[12].后處理時頻分析方法通過對時頻譜的進(jìn)一步處理,使得傳統(tǒng)時頻分析方法的能量譜聚集度得到極大的改善.再賦值方法(RM)[13]通過在時間軸和頻率軸方向重排時頻點,實現(xiàn)了時頻分布的高分辨率.然而,RM是完全依賴時頻譜的,導(dǎo)致該方法失去了信號重構(gòu)能力.文獻(xiàn)[14]首先提出同步壓縮理論并應(yīng)用到連續(xù)小波變換(CWT).文獻(xiàn)[15]將同步壓縮理論應(yīng)用到STFT中,基于STFT的同步壓縮變換(SST)極大地改善了STFT時頻模糊現(xiàn)象并且保留了信號的重構(gòu)能力.SST在處理強(qiáng)時變信號時,其時頻分布結(jié)果會出現(xiàn)分辨率降低的現(xiàn)象.為了克服這一問題,涌現(xiàn)了許多新的時頻分析方法,比如2階SST[16-17]、3階SST[18-19]、高階SST[20]等.隨著階數(shù)的提高,這些方法的計算量也在不斷地增加[21].文獻(xiàn)[21]提出一種新的時頻分析方法,局部最大值同步變換(LMSST).該算法通過在頻率方向檢測時頻譜的局部最大值,構(gòu)造一個全新的頻率再賦值規(guī)則.通過這種再賦值操作,時頻分布的分辨率與SST和RM相比進(jìn)一步得到提高.針對CT存在的問題,本文將LMSST的頻率重排規(guī)則推廣到CT中,提出局部最大值同步壓縮線調(diào)頻小波變換(LMSCT)算法,解決了調(diào)頻率參數(shù)選取對CT時頻表示分辨率的影響.同時,仿真實驗表明,LMSCT比CT具有更好的抗噪性能.最后,將LMSCT應(yīng)用到海雜波信號中,為海雜波背景下的小目標(biāo)檢測提供一個可參考的技術(shù)手段.

    1 算法介紹

    1.1 同步壓縮理論

    一個具有調(diào)幅和調(diào)頻規(guī)律的多分量信號f(t)可以表示為

    (1)

    式中:Ak(t)和φk(t)分別為第k個分量的瞬時幅度和瞬時相位;t為時間變量.

    對于一個待分析信號f滿足f∈L2(R),L為勒貝格積分,對應(yīng)的窗函數(shù)g為實函數(shù)且滿足g∈L2(R),那么此信號的STFT可以定義為

    (2)

    式中:ω為信號的瞬時頻率;u為時間變量.

    SST通過獲得STFT的瞬時頻率信息,然后進(jìn)行頻率點壓縮來提高時頻譜的聚集度.要想獲得STFT的瞬時頻率,首先求取式(2)對時間t的1階導(dǎo)數(shù),則有:

    iωS(t,ω)-Sg′(t,ω)

    (3)

    STFT的瞬時頻率信息ω0(t,ω)可由如下公式求得:

    (4)

    (5)

    式中:δ(·)為狄拉克函數(shù).通過這一后處理操作,可以得到比原STFT聚集度更高的時頻分布.

    1.2 LMSST算法

    作為SST算法的改進(jìn)方法,LMSST算法定義了一個新的頻率再賦值規(guī)則,規(guī)則如下式所示:

    ωm(t,ω)=

    (6)

    式中:Δ為離散頻率間隔;|S(t,ω)|為STFT的時頻譜.假如任意兩個分量相隔足夠的距離,由于窗函數(shù)的Fourier變換在0處達(dá)到最大值,式(6)可以簡化為

    ωm(t,ω)=

    (7)

    為了獲得理想時頻分布,所有的模糊時頻系數(shù)都應(yīng)該沿頻率方向分配到時頻軌跡.因此,LMSST可以實現(xiàn)更高的時頻能量聚集度,可以表示為

    (8)

    1.3 LMSCT算法

    CT同樣是一個非常有效的時頻分析方法,通過引入CR參數(shù)β,其本身就成為了STFT的推廣形式,公式定義為[22]

    (9)

    如果參數(shù)β能夠很好地匹配待分析信號的調(diào)頻率,就能夠獲得能量聚集度較高的時頻分布結(jié)果.但是對于非線性調(diào)頻信號而言,此時參數(shù)β不能夠很好地匹配變化的調(diào)頻率,因此CT算法在處理這類信號時就會受到限制.

    根據(jù)LMSST方法的思想,通過對CT的結(jié)果進(jìn)一步處理,在頻率方向重新分配時頻系數(shù),再賦值規(guī)則可以表示為

    ωL(t,ω)=

    (10)

    因此,提出一個新的時頻分析方法并命名為局部最大值同步壓縮線調(diào)頻小波變換,算法的定義如以下式所示:

    TLMSCT(t,η)=

    (11)

    傳統(tǒng)高階SST隨著階數(shù)的增加,計算量也在相應(yīng)增加.原始SST需要執(zhí)行1次STFT步驟,當(dāng)階數(shù)增加到4時,執(zhí)行STFT的次數(shù)就變?yōu)?1次.而LMSCT方法只需要執(zhí)行1次CT,計算量也相對較小.

    2 數(shù)值實驗

    為了驗證本文提出算法的有效性,需要進(jìn)行2組仿真信號的實驗 ,將本文提出的算法與前文引用的STFT、CT和SST方法進(jìn)行對比.首先將文獻(xiàn)[23]中的仿真信號進(jìn)行驗證,信號1的函數(shù)表達(dá)式為

    f1(t)=sin[2π(10t+5t2/4+t3/9-t4/160)]

    (12)

    信號1的瞬時頻率是時變的并且是單分量的,實驗中該信號的抽樣頻率為100 Hz.根據(jù)式(12),在圖1(a)中繪制了信號的真實頻率分布,圖1(b)展示的是該信號STFT結(jié)果,STFT存在嚴(yán)重的能量擴(kuò)散,不適合分析信號頻率隨時間變化的關(guān)系.在圖2中分別使用CT、SST和LMSCT對信號1做進(jìn)一步處理.

    圖1 仿真信號1的時頻分布Fig.1 Time-frequency distribution of simulation Signal 1

    圖2(a)的時頻分布圖由CT得到,且β=4π.同STFT的結(jié)果相比,CT在一定程度上提高了該信號時頻分布的聚集度.但是從圖2(a)中的兩處標(biāo)記可以看出,在瞬時頻率較小標(biāo)記處,CR的值與信號的瞬時頻率相近,其時頻分布的聚集度相對較好,在瞬時頻率較大標(biāo)記處,CR的值與信號瞬時頻率的差別較大,時頻能量擴(kuò)散幅度較大.這兩處的對比驗證了前文提到的CR值的選取對信號時頻分布的影響.圖2(b) 和2(c)分別由SST與LMSCT得到.與SST的結(jié)果相比,LMSCT的聚集度更高.為了更加直觀地比較,分別截取了SST和LMSCT的局部效果圖進(jìn)行展示,截取部分為圖2(b)和2(c)中的標(biāo)記,局部效果如圖3所示.圖3(a)和3(b)分別為SST與LMSCT的局部展示.通過局部結(jié)果的展示可知,LMSCT的聚集度更高.

    圖2 不同算法時頻分布Fig.2 Time frequency distributions of different algorithms

    圖3 時頻分布局部效果Fig.3 Local effects of time-frequency distribution

    為了系統(tǒng)地比較幾種常見的時頻分析方法聚集度的高低,在這里引入Rényi熵的概念[24].文獻(xiàn)[25]提出將3階Rényi熵用于提供時頻分布的信息度量,文獻(xiàn)[24]詳細(xì)研究了Rényi熵的性質(zhì)和潛在的應(yīng)用價值.Rényi熵是評價時頻分布能量聚集度的指標(biāo)之一.Rényi熵越小代表時頻聚集度越高,相反Rényi熵越大,代表時頻聚集度越低.α階Rényi熵的計算公式如下式所示:

    (13)

    式中:T(t,ω)為信號時頻分布.一般情況下,α默認(rèn)為3.其中,信號1在不同算法下的Rényi熵值如表1所示.

    表1 不同算法處理信號1得到的Rényi熵

    接下來選取一個二分量信號進(jìn)行仿真驗證,信號2的函數(shù)表達(dá)式為

    f2(t)=sin{2π[40t+sin(1.5t)]}+

    sin{2π[17t+6sin(1.5t)]}

    (14)

    仿真實驗中,設(shè)置信號抽樣頻率為100 Hz,抽樣時間為4 s.信號實際頻率分布如圖4(a)所示,圖4(b)為WVD的時頻分布結(jié)果.從圖4中可以觀察到,對于多分量信號,由于交叉項干擾的存在,不能有效提取頻率隨時間變化的信息,所以WVD不適合處理這類型的信號.將主流的幾種時頻分析方法與本文提出的方法進(jìn)行比較,對比結(jié)果如圖5所示.由圖5可知,由LMSCT算法得到的時頻分布的能量聚集度比STFT、CT和SST更高.將圖5的時頻分布轉(zhuǎn)換到三維中,如圖6所示,其中:A為幅值;E為能量譜.由圖6可以清晰地觀察到,STFT和CT的時頻分布存在能量擴(kuò)散的情況;SST與LMSCT能量聚集度較高,然而相較于SST的結(jié)果,LMSCT的幅度比較一致,結(jié)果更加準(zhǔn)確.

    圖4 仿真信號2時頻結(jié)果Fig.4 Time-frequency distribution of simulation Signal 2

    圖5 信號2時頻分布Fig.5 Time-frequency distribution of Signal 2

    圖6 時頻分布的三維展示Fig.6 Three-dimensional representation of time-frequency distribution

    信號2在不同算法下的Rényi熵值如表2所示.為了更直觀地比較,在圖7中繪制了Rényi熵值的散點圖,其中:R為Rényi熵.

    表2 不同算法處理信號2得到的Rényi熵

    圖7 不同算法的Rényi熵值Fig.7 Rényi entropies of different algorithms

    通過表2和圖7列出的幾種算法Rényi熵值大小的比較,可以確定本文提出的算法具有更高的時頻聚集度.通常情況下,現(xiàn)實生活中的真實信號都會含有噪聲,將會影響信號時頻分布的可讀性.接下來對信號2添加高斯白噪聲,然后再進(jìn)行時頻分析算法處理,評估不同算法的抗噪聲性能.

    在信號2中加入信噪比為7 dB的加性高斯白噪聲,不同算法的處理結(jié)果如圖8所示.從圖8(a)和8(b)可以看出,受噪聲影響,STFT與CT的時頻分布不是連續(xù)的,嚴(yán)重影響信號時頻變化特征的提取.圖8(c)和8(d)分別為SST與LMSCT的結(jié)果,同STFT、CT結(jié)果相比,時頻譜能量聚集度都有較大的改善.

    圖8 不同算法抗噪效果Fig.8 Anti-noise effects of different algorithms

    為了進(jìn)一步比較SST與LMSCT的抗噪性能,分別截取SST和LMSCT的部分時頻表示進(jìn)行對比,截取部分為圖8(c)和8(d)標(biāo)記部分,局部結(jié)果如圖9所示.圖9(a)和9(b)分別為SST與LMSCT的局部效果.通過對比可以看出,在相同的頻率點范圍,LMSCT的頻率變化曲線比SST更加平滑,而且沒有出現(xiàn)能量消失的情況,受噪聲影響較小.綜上所述,LMSCT算法的時頻聚集度更高,且抗噪性能更好.

    圖9 抗噪性能局部圖Fig.9 Local diagram of anti-noise performance

    3 IPIX雷達(dá)數(shù)據(jù)分析

    IPIX雷達(dá)由加拿大McMaster大學(xué)Adaptive Systems實驗室設(shè)計,Haykin教授帶領(lǐng)其團(tuán)隊分別于1993年和1998年利用IPIX采集并公開了大量高分辨海雜波數(shù)據(jù),該數(shù)據(jù)已經(jīng)成為測試?yán)走_(dá)檢測算法的重要基準(zhǔn)數(shù)據(jù),在海雜波特性研究方面也做出了重要的貢獻(xiàn)[26].IPIX雷達(dá)是一個全相參的X波段雷達(dá),具有I通道和Q通道兩路收發(fā)信號.

    實驗中,雷達(dá)架設(shè)在加拿大東海岸一個高出海平面25 ~ 30 m的固定位置,雷達(dá)朝大西洋海面照射,待檢測目標(biāo)是被鋁絲包裹直徑1 m的漂浮圓球.雷達(dá)工作頻率為9.3 GHz,波束寬度為0.9°,距離分辨率為30 m.雷達(dá)工作在駐留模式,連續(xù)接收來自某一確定方向的海面回波,脈沖重復(fù)頻率為 1 kHz,駐留時間約為131 s,每組數(shù)據(jù)包含14個距離單元的回波信號.由于雷達(dá)以低掠射角照射目標(biāo),目標(biāo)物體隨海面起伏和擺動導(dǎo)致目標(biāo)能量擴(kuò)散,并且在進(jìn)行數(shù)據(jù)采集時采取了距離過采樣,所以目標(biāo)所在單元周圍的臨近單元會受到目標(biāo)能量的影響,記為受影響單元.每個數(shù)據(jù)文件都由純海雜波、目標(biāo)所在單元回波和受影響目標(biāo)單元回波組成.

    測量該數(shù)據(jù)時的雷達(dá)環(huán)境參數(shù)及數(shù)據(jù)組成如表3所示.回波信號幅度隨距離和時間的變化關(guān)系如圖10所示,其中:U為距離單元.每個距離單元由 131 072 個時域數(shù)據(jù)點組成,目標(biāo)位于第8個單元,距離為2.66 km,由于目標(biāo)的起伏和漂移,7、9、10這3個單元成為受影響單元.為了進(jìn)一步分析,本文選取水平發(fā)射水平接收(HH)極化下的I通道數(shù)據(jù)進(jìn)行處理與分析.

    表3 IPIX雷達(dá)54#文件環(huán)境參數(shù)及數(shù)據(jù)組成

    圖10 海雜波時間-距離-能量圖Fig.10 Didgram of sea clutter time-distance-energy

    單元1不是單元8的鄰近單元,在數(shù)據(jù)采集過程中不會受到目標(biāo)起伏和漂移的影響.首先處理該單元的數(shù)據(jù),分析純凈海雜波信號的時頻分布結(jié)果.單元1純凈海雜波使用3種算法的處理結(jié)果如圖11所示,其中:Sp為抽樣點.從圖11(a)和11(b)中可以看出,時頻分布結(jié)果包含了大量的背景噪聲,并且從圖11(c)的標(biāo)記部分展示的結(jié)果可知,LMSCT算法較好地抑制了噪聲的影響,可以為后續(xù)的含噪聲信號的處理提供參考依據(jù).接下來將對目標(biāo)單元及其鄰近單元的數(shù)據(jù)進(jìn)行處理,對目標(biāo)單元和其他單元數(shù)據(jù)的時頻分布結(jié)果進(jìn)行分析,提取有效的時變特征信息.

    圖11 純凈海雜波(單元1)Fig.11 Pure sea clutter (Unit 1)

    圖12、13、15和16為單元8的鄰近單元處理結(jié)果,分別為單元6、單元7、單元9和單元10.圖14為含有目標(biāo)的單元8處理結(jié)果.由圖14可知,對含目標(biāo)的單元8進(jìn)行處理后,與單元1的處理結(jié)果相比,可以觀察到在0頻附近出現(xiàn)清晰的頻率曲線,并且通過對比圖14(a)~14(c)的時頻結(jié)果發(fā)現(xiàn),LMSCT的頻率曲線更加精細(xì),更有利于時頻變化特征的提取.從圖13、15和16的處理結(jié)果中可以看出,由于距離目標(biāo)單元較近,所以也會出現(xiàn)相似的頻率曲線.由圖12可知,同樣距離目標(biāo)單元較近的單元6卻沒有受到影響.這個現(xiàn)象說明,此影響并不是均勻向兩邊擴(kuò)散.若不能提前確定單元8是含目標(biāo)的數(shù)據(jù),很難從單元7、8、9、10中確定出目標(biāo)單元.為了能夠進(jìn)一步確定目標(biāo)所在單元的位置,采用前文提到的Rényi熵.熵代表信號的無序性,也就是信息量越大,不確定性越小,對應(yīng)熵值越小.時頻分布Rényi熵越小的單元,其含有的信息量越大,即出現(xiàn)目標(biāo)的可能性越大.單元1、6、7、8、9和10分別使用LMSCT算法得到的時頻分布后,計算取得的Rényi熵值如圖17所示.從每個單元Rényi熵的散點圖可以得出,在單元8的Rényi熵值最小,此單元含有的信息量最大,因此可以判定單元8就是含目標(biāo)的單元.

    圖12 純凈海雜波(單元6)Fig.12 Pure sea clutter (Unit 6)

    圖13 受影響海雜波(單元7)Fig.13 Affected sea clutter (Unit 7)

    圖14 含目標(biāo)海雜波(單元8)Fig.14 Target sea clutter (Unit 8)

    圖15 受影響海雜波(單元9)Fig.15 Affected sea clutter (Unit 9)

    圖16 受影響海雜波(單元10)Fig.16 Affected sea clutter (Unit 10)

    圖17 由LMSCT得到的Rényi熵Fig.17 Rényi entropies obtained by LMSCT

    4 結(jié)語

    為提高信號時頻分布的聚集度及抑制噪聲干擾的性能,本文提出一種新的高質(zhì)量的時頻分析算法,即LMSCT.該算法主要思想是將LMSST與CT算法結(jié)合,首先對待分析信號進(jìn)行CT處理,采用局部最大值同步壓縮規(guī)則對CT結(jié)果進(jìn)行頻率點的重新分配,進(jìn)而可得到目標(biāo)信號較為精確的時頻分布圖.數(shù)值實驗結(jié)果表明,LMSCT與其他時頻分析方法相比,該算法時頻分布的能量聚集度更高且抗噪聲性能更好.實際IPIX雷達(dá)信號分析應(yīng)用中,與目前主流的時頻分析算法相比,LMSCT有較小的Rényi熵與瞬時頻率聚集度更高的優(yōu)點以及大幅抑制噪聲干擾,通過比較目標(biāo)信號的Rényi熵,可以有效地確定目標(biāo)所在的距離單元.

    猜你喜歡
    時頻雜波分析方法
    STAR2000型空管一次雷達(dá)雜波抑制淺析
    基于EMD的MEMS陀螺儀隨機(jī)漂移分析方法
    一種角接觸球軸承靜特性分析方法
    中國設(shè)立PSSA的可行性及其分析方法
    中國航海(2019年2期)2019-07-24 08:26:40
    密集雜波環(huán)境下確定性退火DA-HPMHT跟蹤算法
    相關(guān)廣義復(fù)合分布雷達(dá)海雜波仿真
    遙測遙控(2015年2期)2015-04-23 08:15:22
    基于時頻分析的逆合成孔徑雷達(dá)成像技術(shù)
    對采樣數(shù)據(jù)序列進(jìn)行時頻分解法的改進(jìn)
    雙線性時頻分布交叉項提取及損傷識別應(yīng)用
    核安全設(shè)備疲勞分析方法與步驟
    91大片在线观看| 精品久久久久久,| 在线天堂中文资源库| 亚洲国产欧美一区二区综合| 欧美中文日本在线观看视频| 国产精品野战在线观看 | 久久人妻福利社区极品人妻图片| 久久这里只有精品19| 高潮久久久久久久久久久不卡| 水蜜桃什么品种好| 在线视频色国产色| 天堂中文最新版在线下载| 午夜精品久久久久久毛片777| 中文字幕最新亚洲高清| 少妇 在线观看| 亚洲国产看品久久| 男女之事视频高清在线观看| 久久草成人影院| 麻豆一二三区av精品| 999精品在线视频| 中文亚洲av片在线观看爽| 美女午夜性视频免费| 91av网站免费观看| 国产精品亚洲av一区麻豆| 亚洲av成人一区二区三| 999久久久国产精品视频| 亚洲av成人av| 一进一出抽搐动态| 韩国av一区二区三区四区| 999久久久国产精品视频| 成年人免费黄色播放视频| 国产精品偷伦视频观看了| 国产av一区在线观看免费| 色老头精品视频在线观看| 久久国产乱子伦精品免费另类| 精品久久久久久久毛片微露脸| 精品熟女少妇八av免费久了| 真人做人爱边吃奶动态| 少妇粗大呻吟视频| 欧美日韩中文字幕国产精品一区二区三区 | 色播在线永久视频| 国产精品98久久久久久宅男小说| 男人舔女人下体高潮全视频| 一个人观看的视频www高清免费观看 | 日本撒尿小便嘘嘘汇集6| 国产成人精品久久二区二区免费| 精品人妻1区二区| videosex国产| 亚洲一区高清亚洲精品| 手机成人av网站| 精品一区二区三卡| 1024视频免费在线观看| 国产精品永久免费网站| 久久久久久免费高清国产稀缺| 宅男免费午夜| 男人操女人黄网站| 精品国产乱子伦一区二区三区| 久久精品91蜜桃| 免费在线观看完整版高清| 久久久久九九精品影院| 国产精品电影一区二区三区| 亚洲 国产 在线| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲五月天丁香| 夫妻午夜视频| 亚洲精品国产区一区二| 视频区欧美日本亚洲| 99久久国产精品久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 真人一进一出gif抽搐免费| 精品一区二区三区av网在线观看| 老司机靠b影院| 一二三四社区在线视频社区8| 久久九九热精品免费| 精品国产乱子伦一区二区三区| 黄色 视频免费看| 老汉色av国产亚洲站长工具| 三级毛片av免费| 亚洲av片天天在线观看| 男女午夜视频在线观看| 色精品久久人妻99蜜桃| 50天的宝宝边吃奶边哭怎么回事| 国产精品综合久久久久久久免费 | 大型av网站在线播放| 人妻久久中文字幕网| 午夜激情av网站| 国产蜜桃级精品一区二区三区| 1024视频免费在线观看| 国产1区2区3区精品| www.自偷自拍.com| 交换朋友夫妻互换小说| www.熟女人妻精品国产| 日韩中文字幕欧美一区二区| 国产黄a三级三级三级人| 热99re8久久精品国产| 亚洲av成人一区二区三| 国产高清激情床上av| 亚洲黑人精品在线| 国产亚洲欧美98| 久久 成人 亚洲| av免费在线观看网站| 日韩欧美国产一区二区入口| 精品国产美女av久久久久小说| 国产午夜精品久久久久久| 中文字幕人妻丝袜制服| 一级,二级,三级黄色视频| 高清毛片免费观看视频网站 | 亚洲一区二区三区不卡视频| 亚洲五月色婷婷综合| 亚洲专区中文字幕在线| 国产真人三级小视频在线观看| 免费高清视频大片| 99精国产麻豆久久婷婷| 一级片免费观看大全| 日本vs欧美在线观看视频| 国产精品一区二区免费欧美| 亚洲一区中文字幕在线| 女同久久另类99精品国产91| 久久午夜亚洲精品久久| 日韩有码中文字幕| 午夜成年电影在线免费观看| 国产免费男女视频| 日本五十路高清| 色综合站精品国产| 欧美日韩一级在线毛片| 久久人人爽av亚洲精品天堂| 国产亚洲精品综合一区在线观看 | 国产一区在线观看成人免费| 国产免费现黄频在线看| 香蕉国产在线看| 黑人猛操日本美女一级片| 亚洲专区字幕在线| 好男人电影高清在线观看| 深夜精品福利| 国产精品久久视频播放| 精品日产1卡2卡| 亚洲一区二区三区不卡视频| 中文字幕人妻丝袜制服| 十八禁人妻一区二区| 人妻久久中文字幕网| 搡老熟女国产l中国老女人| 国产成人欧美| 1024香蕉在线观看| 亚洲国产欧美一区二区综合| 久久草成人影院| 91九色精品人成在线观看| 精品电影一区二区在线| www.999成人在线观看| 中亚洲国语对白在线视频| 人妻久久中文字幕网| 男男h啪啪无遮挡| 亚洲欧美激情在线| 少妇粗大呻吟视频| 可以在线观看毛片的网站| 国产精品久久久久成人av| 国产精品电影一区二区三区| 老熟妇乱子伦视频在线观看| 欧美国产精品va在线观看不卡| 国产亚洲av高清不卡| 我的亚洲天堂| 亚洲午夜精品一区,二区,三区| 999久久久精品免费观看国产| 精品久久久久久,| 美女高潮喷水抽搐中文字幕| 高清黄色对白视频在线免费看| 视频区图区小说| 亚洲情色 制服丝袜| 不卡av一区二区三区| 三上悠亚av全集在线观看| 久久精品国产亚洲av香蕉五月| 国产精品自产拍在线观看55亚洲| 久久九九热精品免费| 欧美另类亚洲清纯唯美| 一本综合久久免费| 麻豆久久精品国产亚洲av | 成人手机av| 男女高潮啪啪啪动态图| 亚洲国产精品合色在线| 亚洲美女黄片视频| 99精国产麻豆久久婷婷| a在线观看视频网站| 国产熟女午夜一区二区三区| 琪琪午夜伦伦电影理论片6080| 久久久久九九精品影院| 久久久久精品国产欧美久久久| 水蜜桃什么品种好| 色在线成人网| 一夜夜www| 老司机午夜福利在线观看视频| 国产亚洲精品一区二区www| 女人高潮潮喷娇喘18禁视频| 亚洲人成77777在线视频| 欧美日韩视频精品一区| 婷婷丁香在线五月| 在线观看免费高清a一片| 不卡一级毛片| 欧美老熟妇乱子伦牲交| 啪啪无遮挡十八禁网站| 免费在线观看亚洲国产| 免费一级毛片在线播放高清视频 | 国产亚洲精品久久久久久毛片| 久久人人精品亚洲av| 亚洲精品粉嫩美女一区| 又大又爽又粗| 精品午夜福利视频在线观看一区| 黄色毛片三级朝国网站| 动漫黄色视频在线观看| 国产高清激情床上av| 香蕉丝袜av| 精品无人区乱码1区二区| 国产激情欧美一区二区| 热re99久久精品国产66热6| 久久精品国产亚洲av香蕉五月| 欧美成人性av电影在线观看| 夜夜夜夜夜久久久久| 午夜福利影视在线免费观看| 成人特级黄色片久久久久久久| 黑人操中国人逼视频| 男女午夜视频在线观看| 久久精品国产亚洲av高清一级| 午夜成年电影在线免费观看| 中文亚洲av片在线观看爽| 欧美激情极品国产一区二区三区| 国产av又大| 亚洲色图综合在线观看| 亚洲欧美日韩无卡精品| 欧美日韩亚洲综合一区二区三区_| 在线免费观看的www视频| 国产精品久久久人人做人人爽| 国内毛片毛片毛片毛片毛片| 日本撒尿小便嘘嘘汇集6| 精品久久久精品久久久| 午夜免费观看网址| 自拍欧美九色日韩亚洲蝌蚪91| 老汉色∧v一级毛片| 91av网站免费观看| 在线国产一区二区在线| 欧美激情 高清一区二区三区| 十分钟在线观看高清视频www| 久久久久久免费高清国产稀缺| 久久精品国产亚洲av高清一级| 可以在线观看毛片的网站| 国产高清激情床上av| 无人区码免费观看不卡| 变态另类成人亚洲欧美熟女 | 国产黄色免费在线视频| 一边摸一边抽搐一进一出视频| 国产精品九九99| 伦理电影免费视频| av在线播放免费不卡| bbb黄色大片| 侵犯人妻中文字幕一二三四区| 成人亚洲精品一区在线观看| 久久青草综合色| 国产精品香港三级国产av潘金莲| 十分钟在线观看高清视频www| 老司机靠b影院| a级毛片在线看网站| 欧美日韩亚洲高清精品| 一区二区三区精品91| 99国产精品一区二区三区| 国产三级黄色录像| 国产高清国产精品国产三级| 精品一品国产午夜福利视频| 亚洲第一欧美日韩一区二区三区| 最好的美女福利视频网| 99久久国产精品久久久| 国产亚洲精品久久久久久毛片| 18禁美女被吸乳视频| 真人做人爱边吃奶动态| 国产黄a三级三级三级人| 天堂影院成人在线观看| 999久久久国产精品视频| 国产黄色免费在线视频| av免费在线观看网站| 91av网站免费观看| 91九色精品人成在线观看| 亚洲精品美女久久av网站| √禁漫天堂资源中文www| 女人高潮潮喷娇喘18禁视频| 亚洲一卡2卡3卡4卡5卡精品中文| 精品国产亚洲在线| 视频在线观看一区二区三区| 欧美日韩av久久| 久久中文字幕人妻熟女| 99久久99久久久精品蜜桃| 亚洲熟妇中文字幕五十中出 | 国产午夜精品久久久久久| 欧美大码av| 啦啦啦免费观看视频1| 很黄的视频免费| 女人被躁到高潮嗷嗷叫费观| 男人操女人黄网站| 999精品在线视频| 久久久久国内视频| 国产精品1区2区在线观看.| 中文字幕另类日韩欧美亚洲嫩草| 99久久99久久久精品蜜桃| 日韩三级视频一区二区三区| 怎么达到女性高潮| 长腿黑丝高跟| 精品国内亚洲2022精品成人| 国产精品一区二区三区四区久久 | 国产精品免费视频内射| 久久青草综合色| 三上悠亚av全集在线观看| 亚洲人成77777在线视频| 校园春色视频在线观看| 亚洲五月色婷婷综合| 欧美日本中文国产一区发布| 午夜老司机福利片| 久久久久国产精品人妻aⅴ院| 久久中文看片网| 女警被强在线播放| 国产精品一区二区精品视频观看| 51午夜福利影视在线观看| 精品国产国语对白av| 女生性感内裤真人,穿戴方法视频| 成人影院久久| 91麻豆精品激情在线观看国产 | 在线视频色国产色| 国产单亲对白刺激| 亚洲国产精品999在线| 日本wwww免费看| 老汉色∧v一级毛片| 大香蕉久久成人网| 一区二区三区激情视频| 超色免费av| 国产精品免费视频内射| 人成视频在线观看免费观看| 水蜜桃什么品种好| 母亲3免费完整高清在线观看| www日本在线高清视频| 免费观看人在逋| 热re99久久国产66热| 欧美激情极品国产一区二区三区| 久久 成人 亚洲| 天堂√8在线中文| 9色porny在线观看| 一二三四社区在线视频社区8| 久久热在线av| 99精品欧美一区二区三区四区| 亚洲精品国产色婷婷电影| 日韩精品免费视频一区二区三区| 亚洲男人的天堂狠狠| 在线观看一区二区三区激情| 欧美黄色淫秽网站| 成人三级黄色视频| 高清av免费在线| 国产在线观看jvid| 成人三级做爰电影| 国产av一区二区精品久久| 高清av免费在线| 久久狼人影院| 国产精品永久免费网站| a级片在线免费高清观看视频| 国产精品永久免费网站| 欧美老熟妇乱子伦牲交| 亚洲欧美日韩高清在线视频| 亚洲熟妇熟女久久| 热re99久久国产66热| 成人亚洲精品一区在线观看| 好看av亚洲va欧美ⅴa在| 黄色视频,在线免费观看| 国产成人精品久久二区二区免费| 极品教师在线免费播放| 日韩国内少妇激情av| 日韩欧美免费精品| 欧美日韩瑟瑟在线播放| 亚洲av熟女| 老熟妇乱子伦视频在线观看| 久久精品成人免费网站| 日本黄色视频三级网站网址| 久久午夜亚洲精品久久| 大型黄色视频在线免费观看| 女人精品久久久久毛片| 无遮挡黄片免费观看| 首页视频小说图片口味搜索| 国产不卡一卡二| 亚洲精品国产一区二区精华液| 午夜a级毛片| 国产精品 欧美亚洲| 久热这里只有精品99| 亚洲成国产人片在线观看| 精品电影一区二区在线| 丝袜人妻中文字幕| 亚洲国产欧美网| 美女午夜性视频免费| 欧美日韩视频精品一区| 激情在线观看视频在线高清| 女性生殖器流出的白浆| 国产av一区在线观看免费| 国产亚洲精品第一综合不卡| 18禁裸乳无遮挡免费网站照片 | 国产一区二区在线av高清观看| 校园春色视频在线观看| 91字幕亚洲| 免费看十八禁软件| 国产亚洲精品一区二区www| 中文字幕人妻丝袜一区二区| 成熟少妇高潮喷水视频| 999久久久国产精品视频| 18美女黄网站色大片免费观看| 老司机亚洲免费影院| 日韩视频一区二区在线观看| 欧美日韩国产mv在线观看视频| 黄色视频,在线免费观看| 又黄又粗又硬又大视频| 国产精品 国内视频| 国产视频一区二区在线看| 国产色视频综合| 亚洲一区二区三区色噜噜 | 丰满饥渴人妻一区二区三| 在线观看免费日韩欧美大片| 天天躁夜夜躁狠狠躁躁| 在线av久久热| 日本 av在线| 搡老乐熟女国产| 男女床上黄色一级片免费看| 中文字幕高清在线视频| 免费观看人在逋| 亚洲午夜精品一区,二区,三区| 国产精品久久视频播放| 欧美亚洲日本最大视频资源| 精品国产一区二区久久| 啦啦啦免费观看视频1| 多毛熟女@视频| 亚洲精品一二三| 亚洲欧美日韩高清在线视频| 国产一区二区三区在线臀色熟女 | 高清毛片免费观看视频网站 | 午夜免费激情av| 国产亚洲av高清不卡| 久久久久久久久久久久大奶| 欧美 亚洲 国产 日韩一| 电影成人av| 免费久久久久久久精品成人欧美视频| 亚洲,欧美精品.| 两个人看的免费小视频| 国产精品九九99| 国产熟女xx| 亚洲精品中文字幕一二三四区| 国产精品一区二区精品视频观看| tocl精华| 欧美日本亚洲视频在线播放| 国产精品久久电影中文字幕| 欧美黑人欧美精品刺激| 色综合婷婷激情| 日本a在线网址| 国产乱人伦免费视频| 精品福利永久在线观看| 女性被躁到高潮视频| 国产伦人伦偷精品视频| 精品少妇一区二区三区视频日本电影| 电影成人av| 欧美精品啪啪一区二区三区| 亚洲精品中文字幕一二三四区| 少妇 在线观看| 欧美人与性动交α欧美精品济南到| 男人舔女人的私密视频| 精品午夜福利视频在线观看一区| a级片在线免费高清观看视频| 亚洲精品av麻豆狂野| 黑丝袜美女国产一区| 激情在线观看视频在线高清| 久久午夜亚洲精品久久| 国产高清视频在线播放一区| 在线观看免费日韩欧美大片| 亚洲中文字幕日韩| 又大又爽又粗| 又紧又爽又黄一区二区| 亚洲中文av在线| 搡老乐熟女国产| 国产欧美日韩一区二区三| 亚洲男人的天堂狠狠| 国产精品久久久av美女十八| 免费在线观看亚洲国产| ponron亚洲| 丝袜人妻中文字幕| 午夜激情av网站| 天堂中文最新版在线下载| 丝袜美足系列| 亚洲中文av在线| 窝窝影院91人妻| 久久久水蜜桃国产精品网| 欧美久久黑人一区二区| a级片在线免费高清观看视频| 最近最新中文字幕大全电影3 | 亚洲欧美日韩另类电影网站| 久久精品91无色码中文字幕| av免费在线观看网站| 亚洲精华国产精华精| 激情在线观看视频在线高清| 狠狠狠狠99中文字幕| 国产成+人综合+亚洲专区| 亚洲人成77777在线视频| 中文字幕高清在线视频| 琪琪午夜伦伦电影理论片6080| 国产激情欧美一区二区| 亚洲熟妇熟女久久| 窝窝影院91人妻| 久久久国产精品麻豆| netflix在线观看网站| av网站免费在线观看视频| 在线观看免费视频日本深夜| 麻豆av在线久日| 国产亚洲av高清不卡| 国产精品偷伦视频观看了| 色哟哟哟哟哟哟| 99久久国产精品久久久| 久久久国产欧美日韩av| 人妻久久中文字幕网| 久久精品国产亚洲av香蕉五月| 国产精品一区二区免费欧美| 女生性感内裤真人,穿戴方法视频| 丝袜美足系列| 日日干狠狠操夜夜爽| 午夜激情av网站| 丝袜人妻中文字幕| 999精品在线视频| 精品久久久久久电影网| 亚洲精品一区av在线观看| 老司机午夜福利在线观看视频| 精品一品国产午夜福利视频| 国产人伦9x9x在线观看| 黄色女人牲交| 精品一区二区三区av网在线观看| 国产欧美日韩精品亚洲av| 国产亚洲精品综合一区在线观看 | 日本三级黄在线观看| 国产精品国产av在线观看| 国产aⅴ精品一区二区三区波| 岛国视频午夜一区免费看| 精品一区二区三区av网在线观看| 啦啦啦 在线观看视频| 亚洲av成人一区二区三| 亚洲精品中文字幕一二三四区| 超碰成人久久| 久久亚洲真实| 欧美日韩精品网址| av超薄肉色丝袜交足视频| 亚洲精品成人av观看孕妇| 一区二区三区精品91| 欧美亚洲日本最大视频资源| 99香蕉大伊视频| 好男人电影高清在线观看| 女人精品久久久久毛片| 国产主播在线观看一区二区| 中文字幕高清在线视频| 国产黄色免费在线视频| 久久人人精品亚洲av| 91老司机精品| 久久精品国产亚洲av高清一级| 丝袜在线中文字幕| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品粉嫩美女一区| 亚洲欧美精品综合一区二区三区| 国产成人一区二区三区免费视频网站| 国产蜜桃级精品一区二区三区| netflix在线观看网站| 黄片大片在线免费观看| 一级a爱片免费观看的视频| 亚洲色图综合在线观看| 一级片免费观看大全| 国产精品影院久久| 久久热在线av| 狠狠狠狠99中文字幕| 国产欧美日韩一区二区三区在线| 免费av中文字幕在线| av在线天堂中文字幕 | av超薄肉色丝袜交足视频| 日韩欧美国产一区二区入口| 国产人伦9x9x在线观看| 黄频高清免费视频| 三上悠亚av全集在线观看| 高清在线国产一区| 国内毛片毛片毛片毛片毛片| 黄色 视频免费看| 午夜成年电影在线免费观看| 成在线人永久免费视频| 天堂俺去俺来也www色官网| 性欧美人与动物交配| 免费高清视频大片| 国产精品一区二区精品视频观看| 97超级碰碰碰精品色视频在线观看| 亚洲av美国av| 亚洲国产精品合色在线| 精品国产一区二区久久| 亚洲av美国av| 9热在线视频观看99| 美女 人体艺术 gogo| 免费av毛片视频| 无限看片的www在线观看| 99在线视频只有这里精品首页| 久久精品aⅴ一区二区三区四区| 黑人操中国人逼视频| 如日韩欧美国产精品一区二区三区| 男女下面进入的视频免费午夜 | 亚洲成国产人片在线观看| 19禁男女啪啪无遮挡网站| 国产精品免费视频内射| 国产单亲对白刺激| 女性生殖器流出的白浆| 天堂俺去俺来也www色官网| 一级片'在线观看视频| 日韩视频一区二区在线观看| 亚洲一区二区三区色噜噜 | 天天影视国产精品| 亚洲精品av麻豆狂野| 高清毛片免费观看视频网站 | 女人精品久久久久毛片| 亚洲熟妇熟女久久| 一边摸一边做爽爽视频免费| 色综合欧美亚洲国产小说| videosex国产|