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

    基于隱馬爾可夫模型的線譜跟蹤技術

    2014-03-08 06:43:40管景崇胡金華
    艦船科學技術 2014年5期
    關鍵詞:線譜馬爾可夫時刻

    管景崇,胡金華

    (海軍工程大學,湖北 武漢 430033)

    基于隱馬爾可夫模型的線譜跟蹤技術

    管景崇,胡金華

    (海軍工程大學,湖北 武漢 430033)

    線譜檢測和跟蹤是被動聲吶信號處理中的重要內容,本文給出一種基于隱馬爾可夫模型的線譜跟蹤方法。它采用前向后向算法對LOFAR譜圖上線譜進行狀態(tài)估計,然后根據連續(xù)檢驗對每根線譜的起始和終止時間進行檢測,實現對單根線譜和多根線譜的檢測與跟蹤。通過計算機仿真和海試數據處理,驗證了基于隱馬爾可夫模型的線譜跟蹤技術的有效性和穩(wěn)定性。

    隱馬爾可夫;被動聲吶;線譜檢測;線譜跟蹤

    0 引言

    隱馬爾可夫模型 (Hidden Markov Models,HMM)[1-2]是一種用參數表示,用于描述隨機過程統計特性的概率分布。它是一個雙重隨機過程,由馬爾可夫鏈和一般隨機過程組成。隱馬爾可夫模型適用于時間序列建模,并具有處理非線性時變信號的優(yōu)勢,在語音識別、圖像處理、機器視覺及生物醫(yī)學等領域得到廣泛應用。

    線譜跟蹤是被動聲吶信號處理的一個重要內容,對水下運動目標的檢測、跟蹤、識別及分類都具有重要作用,不僅有助于目標運動分析,而且還可用于對潛艇、魚雷等水下運動目標進行分類識別[3-5]。被動聲吶 (Low Frequency Analysis and Recording,LOFAR)譜圖包含豐富的目標信息,為提高目標檢測和識別性能,需要對LOFAR譜圖中的線譜進行自動檢測和跟蹤。

    線譜跟蹤可看作是對LOFAR圖上線譜的估計,并抵消線譜周圍的噪聲。當進行線譜跟蹤時,通常需要估計線譜數目、線譜的起止時間,并同時處理相互交叉的線譜。由于線譜具有時變特性,其狀態(tài)可用1階馬爾可夫鏈來描述,因此可采用HMM實現線譜自動檢測和跟蹤[6-15]。

    本文給出一種基于隱馬爾可夫模型的線譜跟蹤方法。它采用前向后向 (Forward Backward,FB)算法對LOFAR譜圖上線譜進行狀態(tài)估計,然后根據連續(xù)檢驗對每根線譜的起始和終止時間進行檢測。該方法屬于一種跟蹤置前檢測,能有效對多根線譜進行跟蹤,并能估計線譜的起止時間。通過計算機仿真和試驗數據處理,驗證了HMM線譜跟蹤的有效性和穩(wěn)定性。

    1 HMM 基本理論[2]

    HMM模型由若干個狀態(tài)組成,隨著時間的變化,各個狀態(tài)之間可以發(fā)生轉移,也可以在一個狀態(tài)內駐留,每個觀測序列對不同狀態(tài)都有相應的輸出的概率。由于被動聲吶LOFAR譜圖中每根線譜隨著時間變化,因此線譜對應的狀態(tài)可用馬爾可夫鏈來描述。

    1.1 HMM描述

    HMM模型通常采用三元組λ =(π,A,B)來描述。假設HMM為離散平穩(wěn)馬爾可夫鏈,其觀測序列為 Zk=(z1,…,zk),狀態(tài)序列為 Xk=(x1,…,xk),狀態(tài)集為 I。被動聲吶中觀測序列 Zk是LOFAR譜圖中第k時刻功率譜,而狀態(tài)序列Xk是頻率或者頻率的斜率。假定狀態(tài)序列Xk可以看作1階馬爾可夫鏈,則根據馬爾可夫鏈性質有:

    由于LOFAR譜圖上線譜可能存在交叉或彎曲,因此狀態(tài)除了考慮頻率外還需頻率的斜率,則線譜的狀態(tài)矢量可表示為:

    于是,線譜對應的各個狀態(tài)可由狀態(tài)方程描述為:

    式中:H為狀態(tài)矩陣;ηk為零均值加性噪聲;R為協方差矩陣。則狀態(tài)矩陣和協方差矩陣分別可表示為:

    由于初始狀態(tài)的先驗信息未知,因此常需要假定初始狀態(tài)分布為均勻分布,則初始概率分布為:

    根據式(3)可知,線譜在k時刻的狀態(tài)由頻率和頻率的斜率來描述。假定狀態(tài)轉移概率矩陣表示為A=(aj,s(i,r))M×M,則狀態(tài)轉移概率可表示為:

    式中:aj,s(i,r)為二維狀態(tài)轉移概率,表示k-1時刻狀態(tài)(j,s)轉移到k時刻狀態(tài)(i,r)的概率;det為行列式;c為歸一化因子。滿足下式:

    若似然概率矩陣表示為B=(bki)K×M,且第k時刻觀測矢量為zk,則似然概率bi(zk)=Pr(zk|xk=i)表示為k時刻處于狀態(tài)i的觀測矢量zk概率。當信噪比已知,根據貝葉斯準則可計算似然函數概率。然而實際情況信噪比是未知的,這時似然概率可利用功率譜概率積分非參數方法[16]來計算,用公式可表示為:

    式中Pk,i為到k時刻狀態(tài)i的功率譜。該式并不是似然概率的近似,而是概率統計特性。由于假定的模型不是確定的,而是1階馬爾可夫鏈,因此線譜的狀態(tài)序列對應著長時間功率譜積分的最大值。此外該式還可以減小計算量,節(jié)約計算時間。

    1.2 Forward Backward算法

    Forward Backward(FB)算法是一種局部最優(yōu)法,它是狀態(tài)的最大后驗概率估計。根據觀測數據Zk來估計線譜的可能狀態(tài)k,用公式表示為:

    由于FB算法采用αk(i)和βk(i)的遞歸計算,并且αk(i)和βk(i)的值都小于1,因此,αk(i)(隨著k的增加)和βk(i)(隨著k的減少)都迅速的趨向0,在計算過程中常常遇到溢出問題。為了解決溢出問題,通常采用增加比例因子的方法來對上述公式進行修正。根據HMM的前向和后向算法,前向遞推公式可表示為:

    2 HMM檢測和跟蹤

    傳統線譜跟蹤思想是先對線譜進行檢測即跟蹤初始化,然后再沿著時間對線譜進行跟蹤,可以概括為先檢測后跟蹤。HMM自動檢測跟蹤首先根據FB算法對線譜進行跟蹤,估計線譜所對應的狀態(tài),然后利用連續(xù)檢驗來判斷線譜是否存在。因此,在某種意義,HMM自動檢測和跟蹤可以看作是一種“跟蹤置前檢測”。

    2.1 線譜跟蹤

    一般,被動聲吶的LOFAR譜圖上同時存在多根線譜,而且目標線譜間可能相互交叉或彎曲。因此,實際所遇到的問題需要對多根線譜進行同時跟蹤處理。假設LOFAR譜圖上存在L根線譜,則L根線譜的最大后驗概率估計可表示為:

    當LOFAR譜圖上同時存在多根線譜時,多根線譜同時跟蹤處理一般較復雜,且處理過程中還必須避免產生“組合爆炸”問題。因此,本文采用互斥事件和馬爾可夫近似來解決上述的問題。為了避免FB算法返回相同L根線譜 (通常為最強的線譜),需將互斥事件作為前后向概率的條件。于是定義k時刻l事件為:

    多根線譜跟蹤處理是基于互斥事件和馬爾可夫近似假設。在這個假設前提下對FB算法進行修改,采用一種并行結構的FB算法對多根線譜進行同時跟蹤處理。假設L根線譜所對應的L個馬爾可夫鏈過程相互獨立。由于每根線譜對應狀態(tài)由頻率和頻率的斜率兩維矢量組成,因此HMM模型中相應的前向概率和后向概率需用二維來表示。根據互斥事件可知,FB算法中概率是以互斥事件為條件的概率,則多根線譜跟蹤的前向遞推公式為:

    多根線譜跟蹤處理時,為了消除多根線譜的前向概率或后向概率相互耦合,還需要對上述公式進行加權處理。定義2個變量如下:

    2.2 線譜檢測

    在LOFAR譜圖中,線譜并不是在整個觀測時間內存在。因此,線譜的跟蹤必須估計線譜在LOFAR譜圖上的起始和結束時刻,這就需要檢測所跟蹤的線譜在觀測的時間內是否存在。根據統計理論可知,可采用統計檢驗來判決線譜是否存在。若第k時刻所觀測數據為zk,假設H0表示為LOFAR譜圖不存在線譜,zk似然概率為P0(zk);假設H1表示為LOFAR譜圖存在線譜,zk似然概率為P1(zk),則LOFAR譜圖上線譜的有無可用二元假設來表示為:

    在上述2種假設條件下,若觀測數據是“無記憶”的,即只與當前觀測數據有關。根據貝葉斯準則,將所得的似然比與設定的門限比較,則可以判斷線譜是否存在。

    然而,實際中觀測數據在一段時間內持續(xù)。由于統計檢驗是一種瞬時檢驗,它只利用當前時刻的信息,而與當前時刻以前的信息無關。因此,本文采用連續(xù)檢驗來檢測目標的線譜在觀測時間內是否存在。若所跟蹤的線譜在第k-1時刻和第k+1時刻都存在,則根據連續(xù)檢驗,可判斷目標的線譜在第k時刻也存在。由于連續(xù)檢驗和FB算法機理相似,因此可將一種類似FB算法應用連續(xù)檢驗來判斷LOFAR譜圖上線譜是否存在。定義新的狀態(tài)dk,其中dk=0表示線譜不存在,dk=1表示線譜存在。若把線譜存在或不存在當作2種狀態(tài),則可將所處的狀態(tài)看作一階馬爾可夫鏈,則根據FB算法可估計線譜在觀測時間有無。由于目標線譜的狀態(tài)只有存在或不存在2種狀態(tài),因此狀態(tài)轉移概率矩陣可定義為:

    式中:DT為檢測門限;Λl(zk)為第L根線譜的似然比,與信噪比有關。然而,實際情況中信噪比是未知的,由式(31)很難得到似然比。因此,根據H0假設線譜統計獨立和H1假設線譜統計相關特性,將似然比定義為:

    根據HMM前向概率式(19)和后向概率式(20),經過數學推導,可將式(32)的似然比寫為:

    3 仿真和試驗數據處理

    3.1 仿真數據分析

    假定信號由單頻窄帶信號和零均值得高斯白噪聲組成,其中單頻信號的頻率為800 kHz,采樣頻率10 kHz,FFT的點數為8 192,信噪比為-20 dB。

    圖1 原始LOFAR圖Fig.1 Original LOFAR

    圖2 HMM線譜跟蹤Fig.2 HMM frequency line tracking

    圖1和圖2給出了單根線譜仿真數據的處理結果,其中圖1是原始的LOFAR譜圖。從圖中可以看出,線譜只在觀測時間50~150 s內存在。圖2是HMM線譜跟蹤進行估計處理結果,它利用FB算法對線譜進行跟蹤,然后采用連續(xù)檢驗對線譜起止時間進行檢測。從圖中可以看出,線譜僅在50~150 s觀測時間內存在,而其它觀測時間為0,這與原始LOFAR圖的線譜一致。

    假定信號由4個單頻窄帶信號和零均值得高斯白噪聲組成,且4個單頻窄帶信號的頻率分別為f1=600 Hz,f2=800 Hz,f3=820 Hz,f4=1 000 Hz,其中頻率為800 Hz的線譜僅在觀測時間50~150 s內存在,信噪比均為-20 dB,采樣頻率為10 kHz,FFT的點數為8 192,線譜跟蹤的參數為:L=10,vα=7,vβ=7,ω =0.001,DT=1。

    圖3 原始LOFAR圖Fig.3 Original LOFAR

    圖4 HMM線譜跟蹤Fig.4 HMM frequency line tracking

    圖3和圖4給出了多根線譜跟蹤仿真數據的處理結果。從圖3中可以看出,原始LOFAR圖中含有4根線譜,其中800 Hz的線譜僅在一段觀測時間內存在,而在其他觀測時間不存在。圖4是對圖3的LOFAR譜圖進行線譜跟蹤處理。其結果表明,HMM線譜跟蹤能有效地對圖3中的4根線譜進行跟蹤,并且能估計每根線譜的起止時間。

    3.2 試驗數據處理

    為了進一步驗證所提算法的有效性,本節(jié)給出1組試驗數據來進行分析。數據的采樣頻率為10 kHz,數據更新0.1 s,積分時間3 s,FFT的點數為8 192,滑動時間窗的窗長為10,線譜跟蹤參數為:L=10,vα=7,vβ=14,ω =0.005,DT=1.5。

    圖5是某試驗數據LOFAR分析所處理結果,從圖中可以看出,在觀測時間內同時存在多根線譜,其中高頻線譜成分較弱。圖6對圖5LOFAR譜圖上的線譜進行HMM自動檢測和跟蹤處理,其結果表明,HMM線譜跟蹤處理能有效地進行線譜檢測和跟蹤,而且更好地反映出原始LOFAR圖中線譜的信息,提高了線譜檢測能力和線譜質量,從而有利于目標檢測和分類識別。

    圖5 原始LOFAR圖Fig.5 Original LOFAR

    圖6 HMM線譜跟蹤Fig.6 HMM frequency line tracking

    5 結語

    本文通過將HMM應用到被動聲吶LOFAR譜圖中線譜的檢測和跟蹤,給出了一種基于HMM線譜自動檢測和跟蹤的方法。仿真和海試數據處理結果表明,該方法能有效地對多根線譜進行自動檢測和跟蹤,并能估計每根線譜的起止時間。為了避免多根線譜同時存在時產生的組合爆炸,本文對FB算法進行約束來同時處理多根線譜,消除線譜間相互影響,但其計算量隨著線譜數目增加而增加。

    [1]RABINER L R,JUANG B H.An introduction to hidden markov models[J].IEEE Acoustics,Speech,and Signal Processing Magazine,1986(3):4 -16.

    [2]RABINER L R.A tutorial on hidden markov models and selected application in speech recognition[J].Proceedings of the IEEE,1989,77(2):257 -286.

    [3]李啟虎,李敏,楊秀庭.水下目標輻射噪聲中單頻信號分量的檢測:理論分析[J].聲學學報,2008,33(3):193-196.

    [4]李啟虎,李敏,楊秀庭.水下目標輻射噪聲中單頻信號分量的檢測:數值仿真[J].聲學學報,2008,33(4):289 -293.

    [5]方世良,劉萍,賀今.水聲信號LOFAR譜圖的線譜跡圖象處理[J].聲學技術,2004,23(S):191 -193.

    [6]STREITR L,BARRETTR F.Frequency line tracking using hidden Markov models[J].IEEE Transactions on Acoustics,Speech,and Signal Processing,1990,38(4):586-598.

    [7]BARRETT R F,HOLDSWORTH D A.Frequency tracking using hidden Markov models with amplitude and phase information[J].IEEE Transactions on Signal Processing,1993,41(10):2965 -2976.

    [8]JAUFFRET C,BOUCHET D.Frequency line tracking on a lofargram:an efficient wedding between probabilistic data association modeling and dynamic programming technique[J].IEEE Conference Record of the Thirtieth Asilomar Conference on Signals,Systems and Computers,1996(1):486-490.

    [9]CAPPEL D V,ALINAT P.Frequency line extractor using multiple hidden Markov models[J].OCEANS’98.MTS/IEEE Conference Proceedings,1998(3):1481 -1485.

    [10]PAIRS S,JAUFFRET C.A new tracker for multiple frequency line[J].IEEE Proceedings of Aerospace Conference,2001(4):1771 -1782.

    [11]PAIRS S,JAUFFRET C.Frequency line tracking using HMM based schemes[J].IEEE Transaction on Aerospace and Electronic Systems,2003,39(2):439 -449.

    [12]GUNES T,ERDOL N.HMM based spectral frequency line tracking:improvements and new results[J].IEEE International Conference on Acoustics,Speech,and Signal Processing,2006:673 -676.

    [13]SANDBERG F,STRIDH M.Frequency tracking of atrial fibrillation using hidden Markov models[J].IEEE Transaction on Biomedical Engineering,2008,55(2):502-511.

    [14]KERELIUK C,DEPALLE P.Improved hidden Markov model partial tracking through time frequency analysis[J].Proceedings of the 11th International Conference on Digital Audio Effects,2008:1 -6.

    [15]黃海寧,尹力,張春華,等.基于HMM的微弱信號特征增強研究[C].中國聲學學會2002年全國聲學學術會議論文集,2002:475-476.

    [16]MARANDA BH,FAWCETT JS.Detection and localization of weak targets by space time integration[J].IEEE Journal of Oceanic Engineering,1991,16(2):189 -194.

    Research on freguency tine tracking based on HMM

    GUAN Jing-chong1,HU Jin-hua2
    (Naval University of Engineering,Wuhan 430033,China)

    Automatic detection and tracking of narrow band signals,which makes sonar operators easily detect,track and classify moving targets of interest,and improves their surveillance efficiency,is a crucial function of passive sonar.However,confronted with the estimation of the number of lines and the dates of birth and death in practical application,sonar system is difficult to achieve automatic frequency line tracking efficiently.In view of complex situations,a novelmethod of frequency line tracking based on HMM is presented.Each frequency line is modeled as a first order Markov chain,estimated by forward and backward algorithm,and detected using sequential test.The proposed method can efficiently achieve automatic multiple frequency line tracking,and estimate dates of birth and death of target tracks.The processing results from simulation and sea trials have validated the proposed method.

    HMM;passive sonar;sequential test;frequency line tracking

    TB556

    A

    1672-7649(2014)05-0110-06

    10.3404/j.issn.1672-7649.2014.05.023

    2014-01-24;

    2014-04-15

    管景崇(1981-),男,碩士研究生,從事裝備管理研究工作。

    猜你喜歡
    線譜馬爾可夫時刻
    UUV慣導系統多線譜振動抑制研究
    冬“傲”時刻
    捕獵時刻
    基于矢量線陣的目標低頻線譜提取方法
    應用聲學(2020年2期)2020-06-08 02:29:34
    保費隨機且?guī)в屑t利支付的復合馬爾可夫二項模型
    基于SOP的核電廠操縱員監(jiān)視過程馬爾可夫模型
    街拍的歡樂時刻到來了
    應用馬爾可夫鏈對品牌手機市場占有率進行預測
    “”維譜在艦船輻射噪聲線譜提取中的應用
    一天的時刻
    性色avwww在线观看| 国产精品av视频在线免费观看| 久久国产精品影院| 精品久久久久久久久久免费视频| 国产97色在线日韩免费| 最近在线观看免费完整版| 国产精品女同一区二区软件 | 狠狠狠狠99中文字幕| 在线免费观看的www视频| 久久精品人妻少妇| 日韩中文字幕欧美一区二区| 色老头精品视频在线观看| 亚洲国产精品合色在线| 真人做人爱边吃奶动态| 国产单亲对白刺激| 欧美色欧美亚洲另类二区| 国产高清有码在线观看视频| 欧美性猛交╳xxx乱大交人| 国产精品免费一区二区三区在线| 国产精品久久久久久亚洲av鲁大| 亚洲一区二区三区色噜噜| 网址你懂的国产日韩在线| 小蜜桃在线观看免费完整版高清| 99久国产av精品| 夜夜夜夜夜久久久久| 精品一区二区三区视频在线 | 级片在线观看| 午夜激情欧美在线| 中亚洲国语对白在线视频| 99国产精品一区二区蜜桃av| 90打野战视频偷拍视频| 亚洲精品粉嫩美女一区| 91麻豆精品激情在线观看国产| 精品人妻1区二区| av中文乱码字幕在线| 久久久久久九九精品二区国产| 国产精品98久久久久久宅男小说| 一区二区三区国产精品乱码| 免费人成在线观看视频色| 看免费av毛片| 欧美一级a爱片免费观看看| 日本三级黄在线观看| 亚洲av不卡在线观看| 国产69精品久久久久777片| 两人在一起打扑克的视频| 18禁黄网站禁片午夜丰满| 国产精品嫩草影院av在线观看 | 日韩高清综合在线| 日韩av在线大香蕉| 在线观看美女被高潮喷水网站 | 久久精品国产亚洲av涩爱 | 久久国产精品影院| 国产乱人伦免费视频| av福利片在线观看| 搞女人的毛片| 国产午夜精品论理片| www日本在线高清视频| 成人国产综合亚洲| 成人高潮视频无遮挡免费网站| 欧美+日韩+精品| 欧美精品啪啪一区二区三区| 白带黄色成豆腐渣| 亚洲精品一卡2卡三卡4卡5卡| 波野结衣二区三区在线 | 午夜福利在线在线| 此物有八面人人有两片| 午夜免费成人在线视频| 亚洲欧美日韩无卡精品| 免费看a级黄色片| 亚洲精品日韩av片在线观看 | 亚洲五月天丁香| 午夜福利在线观看吧| 午夜福利视频1000在线观看| 十八禁人妻一区二区| 嫁个100分男人电影在线观看| 亚洲国产色片| 国产中年淑女户外野战色| 51国产日韩欧美| 在线看三级毛片| 性欧美人与动物交配| 精品日产1卡2卡| 国产高清三级在线| 成人一区二区视频在线观看| 一级毛片高清免费大全| 小蜜桃在线观看免费完整版高清| 欧美日本视频| 日韩欧美在线乱码| 久久精品亚洲精品国产色婷小说| www日本在线高清视频| 免费av不卡在线播放| 国产精品久久久久久久电影 | 国产亚洲精品久久久com| 男人和女人高潮做爰伦理| 嫁个100分男人电影在线观看| 女警被强在线播放| 亚洲人成伊人成综合网2020| 成年免费大片在线观看| 一卡2卡三卡四卡精品乱码亚洲| 丰满人妻熟妇乱又伦精品不卡| 天堂网av新在线| 午夜日韩欧美国产| 看黄色毛片网站| 亚洲最大成人中文| 精品不卡国产一区二区三区| 日韩欧美国产在线观看| 99久久成人亚洲精品观看| 黄色片一级片一级黄色片| 成年免费大片在线观看| 亚洲av电影不卡..在线观看| 国产aⅴ精品一区二区三区波| 99热这里只有是精品50| 国产精品久久电影中文字幕| xxx96com| 窝窝影院91人妻| 伊人久久精品亚洲午夜| 国产精品1区2区在线观看.| 国产极品精品免费视频能看的| 91久久精品电影网| 美女免费视频网站| 岛国视频午夜一区免费看| 舔av片在线| av视频在线观看入口| 亚洲最大成人手机在线| 成年免费大片在线观看| 国产精品永久免费网站| 精品久久久久久久久久免费视频| 成人欧美大片| 成年版毛片免费区| 亚洲国产高清在线一区二区三| 婷婷精品国产亚洲av在线| 搡女人真爽免费视频火全软件 | 搞女人的毛片| 在线观看66精品国产| av在线蜜桃| 精华霜和精华液先用哪个| 一区福利在线观看| 亚洲内射少妇av| 免费看美女性在线毛片视频| 老鸭窝网址在线观看| 麻豆成人午夜福利视频| 久久久久九九精品影院| 国产97色在线日韩免费| 国产 一区 欧美 日韩| 色综合婷婷激情| 精品一区二区三区av网在线观看| 国产精品精品国产色婷婷| 熟妇人妻久久中文字幕3abv| 亚洲中文日韩欧美视频| 手机成人av网站| 国产亚洲精品一区二区www| 成年版毛片免费区| 日本五十路高清| 一进一出好大好爽视频| 欧美日韩中文字幕国产精品一区二区三区| а√天堂www在线а√下载| 一个人观看的视频www高清免费观看| 免费观看人在逋| 天天躁日日操中文字幕| 免费在线观看成人毛片| 中国美女看黄片| 法律面前人人平等表现在哪些方面| 午夜精品一区二区三区免费看| 天天一区二区日本电影三级| 一进一出抽搐gif免费好疼| 国产一区在线观看成人免费| 高清日韩中文字幕在线| 女人十人毛片免费观看3o分钟| 男人和女人高潮做爰伦理| www.熟女人妻精品国产| 少妇熟女aⅴ在线视频| 黄色日韩在线| 午夜影院日韩av| 五月玫瑰六月丁香| 精品免费久久久久久久清纯| 最后的刺客免费高清国语| 成人亚洲精品av一区二区| 91九色精品人成在线观看| 久久久久免费精品人妻一区二区| 免费在线观看日本一区| 国产精品国产高清国产av| 久久久久久九九精品二区国产| 一a级毛片在线观看| 国产精品一及| 天堂√8在线中文| 精品久久久久久久人妻蜜臀av| 亚洲成人久久性| 青草久久国产| 国产午夜精品久久久久久一区二区三区 | 国产色婷婷99| 日日夜夜操网爽| 在线观看66精品国产| 在线视频色国产色| 久久久久久久久大av| 女人十人毛片免费观看3o分钟| 精品熟女少妇八av免费久了| 99久久综合精品五月天人人| 国产激情偷乱视频一区二区| 国产精品一区二区三区四区免费观看 | 免费看十八禁软件| 很黄的视频免费| 国产综合懂色| 欧美日韩国产亚洲二区| 亚洲中文字幕一区二区三区有码在线看| 日本五十路高清| 欧美午夜高清在线| 搡老妇女老女人老熟妇| 毛片女人毛片| 99久久精品国产亚洲精品| 天堂影院成人在线观看| 国产精品一区二区三区四区免费观看 | 免费av观看视频| 国产探花极品一区二区| 国产一级毛片七仙女欲春2| 日韩 欧美 亚洲 中文字幕| 嫩草影院精品99| 欧美成狂野欧美在线观看| 亚洲一区高清亚洲精品| 人妻久久中文字幕网| 成人永久免费在线观看视频| 嫩草影视91久久| www.999成人在线观看| 母亲3免费完整高清在线观看| 久9热在线精品视频| 九色国产91popny在线| 日韩欧美在线乱码| 久久精品91蜜桃| 97超视频在线观看视频| 我要搜黄色片| 女警被强在线播放| 成人一区二区视频在线观看| av天堂中文字幕网| 少妇人妻一区二区三区视频| 久久久久久国产a免费观看| 亚洲av免费高清在线观看| 亚洲国产中文字幕在线视频| 少妇的逼水好多| 久久精品91蜜桃| eeuss影院久久| 亚洲,欧美精品.| 波多野结衣高清作品| 日韩av在线大香蕉| 99久久综合精品五月天人人| 国产精品永久免费网站| 国产伦在线观看视频一区| 大型黄色视频在线免费观看| 日韩欧美三级三区| 午夜福利在线观看吧| 亚洲第一欧美日韩一区二区三区| 成人精品一区二区免费| 日本黄色视频三级网站网址| 最新中文字幕久久久久| 90打野战视频偷拍视频| av女优亚洲男人天堂| 国产成人av激情在线播放| 一级黄色大片毛片| 久久精品国产亚洲av香蕉五月| 两个人视频免费观看高清| e午夜精品久久久久久久| 99久久精品一区二区三区| 亚洲一区二区三区色噜噜| 日韩免费av在线播放| 3wmmmm亚洲av在线观看| 一个人免费在线观看的高清视频| 给我免费播放毛片高清在线观看| 亚洲不卡免费看| 国产亚洲精品一区二区www| 嫩草影视91久久| 天堂√8在线中文| 国产精华一区二区三区| 亚洲国产精品久久男人天堂| 欧美乱色亚洲激情| 2021天堂中文幕一二区在线观| 亚洲精品国产精品久久久不卡| 亚洲熟妇中文字幕五十中出| 在线观看午夜福利视频| 蜜桃亚洲精品一区二区三区| 亚洲国产欧美网| 久久精品人妻少妇| 哪里可以看免费的av片| 青草久久国产| 久久久国产成人精品二区| 成人精品一区二区免费| 亚洲 欧美 日韩 在线 免费| www.色视频.com| 一个人观看的视频www高清免费观看| 老司机福利观看| 免费看光身美女| 中文字幕人成人乱码亚洲影| 免费高清视频大片| 男插女下体视频免费在线播放| 蜜桃久久精品国产亚洲av| 欧美日韩精品网址| 在线天堂最新版资源| 国产一区在线观看成人免费| 精品久久久久久久毛片微露脸| 狂野欧美白嫩少妇大欣赏| 韩国av一区二区三区四区| 波多野结衣高清无吗| 麻豆国产97在线/欧美| 99国产精品一区二区蜜桃av| 久久精品国产99精品国产亚洲性色| 亚洲国产日韩欧美精品在线观看 | 欧美一区二区亚洲| 欧美色欧美亚洲另类二区| 欧美+日韩+精品| 亚洲欧美日韩无卡精品| 51午夜福利影视在线观看| 亚洲成a人片在线一区二区| 欧美日本视频| av中文乱码字幕在线| 欧美性猛交黑人性爽| 99在线视频只有这里精品首页| 成人欧美大片| 伊人久久大香线蕉亚洲五| 欧美日韩中文字幕国产精品一区二区三区| 欧美色视频一区免费| 99热这里只有是精品50| 全区人妻精品视频| 88av欧美| 好男人在线观看高清免费视频| 最新美女视频免费是黄的| 在线观看66精品国产| 久久精品国产清高在天天线| 久久午夜亚洲精品久久| 久久伊人香网站| 一个人看的www免费观看视频| 亚洲乱码一区二区免费版| 亚洲av第一区精品v没综合| x7x7x7水蜜桃| 又黄又粗又硬又大视频| 免费av毛片视频| 1024手机看黄色片| 午夜免费男女啪啪视频观看 | 99久久综合精品五月天人人| 亚洲精品一卡2卡三卡4卡5卡| 国内揄拍国产精品人妻在线| 午夜福利在线观看吧| 九色国产91popny在线| 国产高清视频在线播放一区| 一二三四社区在线视频社区8| 狂野欧美白嫩少妇大欣赏| 欧美激情久久久久久爽电影| 日本 av在线| av视频在线观看入口| 久久久久精品国产欧美久久久| 白带黄色成豆腐渣| 色吧在线观看| 国产精华一区二区三区| 免费一级毛片在线播放高清视频| 成人国产一区最新在线观看| 亚洲熟妇中文字幕五十中出| 欧美色视频一区免费| 怎么达到女性高潮| 亚洲一区二区三区色噜噜| 长腿黑丝高跟| 亚洲人与动物交配视频| www日本在线高清视频| 搡女人真爽免费视频火全软件 | 美女黄网站色视频| 一个人免费在线观看电影| 亚洲成人久久爱视频| 色综合婷婷激情| 日韩成人在线观看一区二区三区| 深爱激情五月婷婷| aaaaa片日本免费| 韩国av一区二区三区四区| 亚洲精品日韩av片在线观看 | 99热6这里只有精品| 久久久久免费精品人妻一区二区| 色噜噜av男人的天堂激情| 美女cb高潮喷水在线观看| 人妻丰满熟妇av一区二区三区| 亚洲成人久久性| 成人特级av手机在线观看| 在线看三级毛片| www日本在线高清视频| 18禁国产床啪视频网站| 毛片女人毛片| 欧美高清成人免费视频www| 日韩欧美精品v在线| 国产高清有码在线观看视频| 色尼玛亚洲综合影院| 亚洲国产高清在线一区二区三| 国产黄a三级三级三级人| 最新在线观看一区二区三区| www日本在线高清视频| 97超视频在线观看视频| 三级毛片av免费| 精品福利观看| 国产欧美日韩精品亚洲av| 色视频www国产| 白带黄色成豆腐渣| 亚洲国产欧美人成| 99久久综合精品五月天人人| 亚洲中文字幕日韩| 亚洲国产欧美人成| 91麻豆av在线| 精品一区二区三区视频在线观看免费| 麻豆国产97在线/欧美| 亚洲欧美激情综合另类| 婷婷精品国产亚洲av| 国产欧美日韩精品亚洲av| 桃色一区二区三区在线观看| 亚洲av二区三区四区| 国产黄a三级三级三级人| 国产v大片淫在线免费观看| 日本黄色视频三级网站网址| 成人特级av手机在线观看| 色av中文字幕| 国内精品美女久久久久久| 最近视频中文字幕2019在线8| 国产成人啪精品午夜网站| 每晚都被弄得嗷嗷叫到高潮| 69人妻影院| 精品一区二区三区av网在线观看| 亚洲真实伦在线观看| 亚洲精品影视一区二区三区av| 麻豆一二三区av精品| 午夜福利欧美成人| a在线观看视频网站| 成人精品一区二区免费| 搡女人真爽免费视频火全软件 | av在线天堂中文字幕| 国产精品 欧美亚洲| 中文字幕av成人在线电影| 国产精品久久久久久人妻精品电影| 精品久久久久久久久久免费视频| 亚洲真实伦在线观看| 黄片大片在线免费观看| 在线播放无遮挡| 在线免费观看的www视频| 无人区码免费观看不卡| 在线观看66精品国产| 国产一区二区激情短视频| 人妻久久中文字幕网| 欧美乱码精品一区二区三区| 好男人电影高清在线观看| 亚洲aⅴ乱码一区二区在线播放| 精品一区二区三区av网在线观看| 美女大奶头视频| 一级作爱视频免费观看| 国产黄片美女视频| 国产一区二区亚洲精品在线观看| 精品久久久久久成人av| 国产三级中文精品| a在线观看视频网站| 国产伦一二天堂av在线观看| 亚洲国产欧洲综合997久久,| 亚洲人成网站在线播| 人妻夜夜爽99麻豆av| 女人十人毛片免费观看3o分钟| 成人精品一区二区免费| 久久久成人免费电影| 少妇人妻精品综合一区二区 | 亚洲精品色激情综合| 亚洲第一欧美日韩一区二区三区| 国产一区二区三区视频了| 欧美在线一区亚洲| a级一级毛片免费在线观看| 少妇丰满av| 亚洲人与动物交配视频| 国产精品98久久久久久宅男小说| 亚洲av日韩精品久久久久久密| 黄色丝袜av网址大全| 国产精品影院久久| 一边摸一边抽搐一进一小说| 国产伦精品一区二区三区四那| 久久这里只有精品中国| 99热6这里只有精品| 亚洲中文字幕日韩| 最近最新中文字幕大全免费视频| 波多野结衣巨乳人妻| 国产精品免费一区二区三区在线| 国产成人欧美在线观看| 日本免费一区二区三区高清不卡| 中亚洲国语对白在线视频| 91在线观看av| 亚洲av成人av| 国产精品永久免费网站| 最好的美女福利视频网| 午夜福利视频1000在线观看| 欧美+日韩+精品| 国产一区二区在线av高清观看| 国内久久婷婷六月综合欲色啪| 波多野结衣高清无吗| 成人无遮挡网站| 无限看片的www在线观看| 99在线人妻在线中文字幕| 日本在线视频免费播放| 黄色视频,在线免费观看| 制服人妻中文乱码| 人人妻人人澡欧美一区二区| 一本一本综合久久| 久久精品91蜜桃| 欧美性猛交黑人性爽| 老司机福利观看| 日韩精品青青久久久久久| 好男人在线观看高清免费视频| 在线十欧美十亚洲十日本专区| 久久久久久久久久黄片| 亚洲国产精品久久男人天堂| 午夜福利18| 亚洲人与动物交配视频| 日韩中文字幕欧美一区二区| 免费看美女性在线毛片视频| 天堂av国产一区二区熟女人妻| 国产精品一及| 亚洲欧美精品综合久久99| 此物有八面人人有两片| 亚洲国产精品成人综合色| 看黄色毛片网站| 男女床上黄色一级片免费看| 欧美另类亚洲清纯唯美| 亚洲精品日韩av片在线观看 | 欧美极品一区二区三区四区| 精品人妻1区二区| 欧美日韩一级在线毛片| 亚洲精品在线美女| 国产成人影院久久av| 日韩国内少妇激情av| a级毛片a级免费在线| 国产色爽女视频免费观看| 久久99热这里只有精品18| 亚洲av不卡在线观看| 亚洲成人精品中文字幕电影| 18禁黄网站禁片午夜丰满| 色哟哟哟哟哟哟| 欧美成人一区二区免费高清观看| 成人高潮视频无遮挡免费网站| 国产一区二区三区在线臀色熟女| xxx96com| 久久中文看片网| 法律面前人人平等表现在哪些方面| 亚洲精品一区av在线观看| 久9热在线精品视频| 国产极品精品免费视频能看的| 免费看十八禁软件| 欧美日韩精品网址| 成人午夜高清在线视频| 一本综合久久免费| 极品教师在线免费播放| 天堂av国产一区二区熟女人妻| 美女黄网站色视频| 国产av一区在线观看免费| 亚洲片人在线观看| 97超级碰碰碰精品色视频在线观看| 免费在线观看影片大全网站| 亚洲黑人精品在线| 手机成人av网站| av视频在线观看入口| 国产精品亚洲美女久久久| 嫩草影院精品99| 色综合婷婷激情| 久久久国产精品麻豆| 黑人欧美特级aaaaaa片| 国产主播在线观看一区二区| 色av中文字幕| 3wmmmm亚洲av在线观看| 国产精品久久久久久久久免 | 脱女人内裤的视频| 日本精品一区二区三区蜜桃| 亚洲欧美日韩高清在线视频| 亚洲,欧美精品.| 国产综合懂色| 97超视频在线观看视频| 欧美大码av| 日韩欧美一区二区三区在线观看| 精品99又大又爽又粗少妇毛片 | 美女高潮的动态| 国产精品久久久久久亚洲av鲁大| 给我免费播放毛片高清在线观看| 亚洲精品一卡2卡三卡4卡5卡| 色精品久久人妻99蜜桃| av福利片在线观看| 亚洲精华国产精华精| 老熟妇乱子伦视频在线观看| 女人十人毛片免费观看3o分钟| 老司机深夜福利视频在线观看| 超碰av人人做人人爽久久 | 日韩免费av在线播放| 精品熟女少妇八av免费久了| 欧美一区二区精品小视频在线| bbb黄色大片| 亚洲精品亚洲一区二区| 精品无人区乱码1区二区| 97碰自拍视频| 夜夜夜夜夜久久久久| 少妇人妻一区二区三区视频| 欧美zozozo另类| 看片在线看免费视频| 噜噜噜噜噜久久久久久91| 好看av亚洲va欧美ⅴa在| 免费观看人在逋| 久久精品国产亚洲av涩爱 | 成人欧美大片| 亚洲av电影在线进入| 欧美+日韩+精品| 99国产综合亚洲精品| 国产探花极品一区二区| 欧美zozozo另类| 变态另类成人亚洲欧美熟女| 成人无遮挡网站| 在线观看66精品国产| 桃红色精品国产亚洲av| 国产免费一级a男人的天堂| 久99久视频精品免费| 久久精品人妻少妇| 国产精品美女特级片免费视频播放器| 亚洲 国产 在线| 国产单亲对白刺激| 欧美bdsm另类| 中文字幕人妻熟人妻熟丝袜美 | 中出人妻视频一区二区| 国产蜜桃级精品一区二区三区| 黄色丝袜av网址大全|