• <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)視過程馬爾可夫模型
    街拍的歡樂時刻到來了
    應用馬爾可夫鏈對品牌手機市場占有率進行預測
    “”維譜在艦船輻射噪聲線譜提取中的應用
    一天的時刻
    久久久亚洲精品成人影院| 欧美日韩国产mv在线观看视频| 亚洲av国产av综合av卡| 久久亚洲国产成人精品v| 免费黄网站久久成人精品| 在线 av 中文字幕| 青春草国产在线视频| 日本爱情动作片www.在线观看| 日韩av在线免费看完整版不卡| av福利片在线观看| 麻豆成人av视频| 久久影院123| 深夜a级毛片| 九草在线视频观看| 久久精品国产鲁丝片午夜精品| 99久久综合免费| 日韩强制内射视频| 亚洲精品乱久久久久久| 国产一区二区三区av在线| 久久久久久久久久久久大奶| 国产有黄有色有爽视频| 日韩,欧美,国产一区二区三区| 大香蕉久久网| 精品少妇黑人巨大在线播放| 国产精品麻豆人妻色哟哟久久| 观看免费一级毛片| 成年av动漫网址| 国产69精品久久久久777片| av国产精品久久久久影院| 亚洲内射少妇av| 久久久久人妻精品一区果冻| 国产成人精品婷婷| 国产亚洲欧美精品永久| 精品午夜福利在线看| 精品一区二区免费观看| 视频区图区小说| 少妇被粗大猛烈的视频| 午夜av观看不卡| 一区在线观看完整版| 国产一区二区三区av在线| 深夜a级毛片| 人人妻人人澡人人看| 国产一区二区三区av在线| 欧美日韩视频精品一区| 国产无遮挡羞羞视频在线观看| 国产男女内射视频| 免费大片18禁| 久久97久久精品| 激情五月婷婷亚洲| 人人妻人人添人人爽欧美一区卜| 丰满乱子伦码专区| 搡老乐熟女国产| 视频中文字幕在线观看| 成人影院久久| 成人特级av手机在线观看| 久久综合国产亚洲精品| 综合色丁香网| 国产免费视频播放在线视频| 美女内射精品一级片tv| 久久久久网色| 亚洲四区av| 国内精品宾馆在线| 午夜91福利影院| 国产免费视频播放在线视频| 国产一区二区在线观看日韩| 精品久久久久久久久亚洲| 亚洲激情五月婷婷啪啪| 边亲边吃奶的免费视频| 男女啪啪激烈高潮av片| 男的添女的下面高潮视频| 久久久久久久久久久免费av| 国产精品一区二区在线观看99| av卡一久久| 能在线免费看毛片的网站| 亚洲av男天堂| 日韩av在线免费看完整版不卡| 国产一区二区三区av在线| 久久99精品国语久久久| 一区在线观看完整版| 免费播放大片免费观看视频在线观看| 少妇人妻一区二区三区视频| 我的老师免费观看完整版| 永久免费av网站大全| 亚洲精品视频女| 国产精品.久久久| 国产黄片美女视频| 欧美xxⅹ黑人| 人体艺术视频欧美日本| 欧美日本中文国产一区发布| 人人妻人人看人人澡| 亚洲,欧美,日韩| 精品久久久噜噜| 精品人妻一区二区三区麻豆| 妹子高潮喷水视频| 在线观看一区二区三区激情| 日韩av免费高清视频| www.色视频.com| 亚洲精品色激情综合| 人妻夜夜爽99麻豆av| 精品酒店卫生间| 另类亚洲欧美激情| 波野结衣二区三区在线| 欧美日韩综合久久久久久| 全区人妻精品视频| 美女大奶头黄色视频| 99国产精品免费福利视频| 欧美变态另类bdsm刘玥| av线在线观看网站| 久久精品国产鲁丝片午夜精品| 少妇被粗大的猛进出69影院 | 午夜福利,免费看| a级一级毛片免费在线观看| 亚洲欧美成人精品一区二区| 亚洲av二区三区四区| 寂寞人妻少妇视频99o| 22中文网久久字幕| 日韩一本色道免费dvd| 欧美亚洲 丝袜 人妻 在线| 69精品国产乱码久久久| 欧美日韩在线观看h| 欧美精品国产亚洲| av一本久久久久| 王馨瑶露胸无遮挡在线观看| 亚洲成人手机| 麻豆成人午夜福利视频| 免费看日本二区| 亚洲美女搞黄在线观看| 大又大粗又爽又黄少妇毛片口| a级片在线免费高清观看视频| 人妻制服诱惑在线中文字幕| 亚洲一区二区三区欧美精品| 亚洲精品aⅴ在线观看| 最黄视频免费看| 大码成人一级视频| 久久精品久久精品一区二区三区| 国产精品久久久久久精品电影小说| 日本欧美国产在线视频| 国产无遮挡羞羞视频在线观看| 精品久久久噜噜| 亚洲丝袜综合中文字幕| 丁香六月天网| 免费黄网站久久成人精品| 午夜免费观看性视频| 亚洲中文av在线| 我的老师免费观看完整版| av在线观看视频网站免费| 极品人妻少妇av视频| 亚洲av电影在线观看一区二区三区| 国模一区二区三区四区视频| 日韩 亚洲 欧美在线| 亚洲人成网站在线播| 亚洲三级黄色毛片| 国产男女超爽视频在线观看| 男的添女的下面高潮视频| 欧美bdsm另类| 最近最新中文字幕免费大全7| 美女xxoo啪啪120秒动态图| 午夜老司机福利剧场| 免费看av在线观看网站| 中文字幕久久专区| 精品亚洲成国产av| 在线播放无遮挡| 黑丝袜美女国产一区| 男的添女的下面高潮视频| 人人妻人人澡人人看| 多毛熟女@视频| 国产黄片视频在线免费观看| 51国产日韩欧美| 91成人精品电影| 午夜视频国产福利| 免费观看av网站的网址| 久久精品久久久久久噜噜老黄| 在线观看三级黄色| 日韩成人av中文字幕在线观看| 国产精品熟女久久久久浪| 亚洲欧美日韩另类电影网站| 中文字幕免费在线视频6| 久久婷婷青草| 性高湖久久久久久久久免费观看| 夜夜看夜夜爽夜夜摸| 国产精品国产三级专区第一集| 精品一区二区免费观看| 成人漫画全彩无遮挡| 大片免费播放器 马上看| 超碰97精品在线观看| 日韩视频在线欧美| 伦理电影大哥的女人| 我的老师免费观看完整版| 亚洲精品自拍成人| 观看av在线不卡| 又粗又硬又长又爽又黄的视频| 99九九线精品视频在线观看视频| 五月天丁香电影| 亚洲精品中文字幕在线视频 | 国产免费福利视频在线观看| 汤姆久久久久久久影院中文字幕| 少妇被粗大猛烈的视频| 少妇被粗大猛烈的视频| 高清av免费在线| 丰满少妇做爰视频| 在线 av 中文字幕| 少妇熟女欧美另类| 亚洲精品456在线播放app| 亚洲丝袜综合中文字幕| 女人久久www免费人成看片| 22中文网久久字幕| 国产成人91sexporn| av天堂久久9| 人妻人人澡人人爽人人| 免费黄色在线免费观看| 久久久久久人妻| av视频免费观看在线观看| 久久久欧美国产精品| 国产精品99久久久久久久久| 亚洲在久久综合| 亚洲欧美日韩东京热| 人妻一区二区av| 日本-黄色视频高清免费观看| 色吧在线观看| 丰满迷人的少妇在线观看| 国产在视频线精品| 香蕉精品网在线| 亚洲国产毛片av蜜桃av| 精品一品国产午夜福利视频| 久久女婷五月综合色啪小说| 欧美日韩综合久久久久久| 老司机影院毛片| 久久精品国产自在天天线| 日本猛色少妇xxxxx猛交久久| 99九九在线精品视频 | 日韩av不卡免费在线播放| 我的女老师完整版在线观看| 免费观看的影片在线观看| 日韩电影二区| 日本-黄色视频高清免费观看| 成人国产av品久久久| 99re6热这里在线精品视频| 免费观看a级毛片全部| 亚洲欧美日韩卡通动漫| 老熟女久久久| 欧美xxⅹ黑人| 国产男女内射视频| 中文字幕精品免费在线观看视频 | 精品卡一卡二卡四卡免费| 欧美成人精品欧美一级黄| 成年人免费黄色播放视频 | 午夜老司机福利剧场| 国产成人精品一,二区| 制服丝袜香蕉在线| 高清欧美精品videossex| 99久久精品一区二区三区| 交换朋友夫妻互换小说| 一级a做视频免费观看| 亚洲人与动物交配视频| 免费看光身美女| 成人国产av品久久久| www.av在线官网国产| 少妇人妻久久综合中文| 久久国产亚洲av麻豆专区| 国产成人aa在线观看| 亚洲精华国产精华液的使用体验| 亚洲av中文av极速乱| 多毛熟女@视频| 美女内射精品一级片tv| 精品亚洲成国产av| 寂寞人妻少妇视频99o| 男女边吃奶边做爰视频| 国产日韩欧美视频二区| 十八禁网站网址无遮挡 | 亚洲欧美一区二区三区国产| tube8黄色片| h日本视频在线播放| www.色视频.com| 日韩免费高清中文字幕av| 80岁老熟妇乱子伦牲交| 男人舔奶头视频| 伦理电影免费视频| 人妻系列 视频| 久久久久国产网址| 亚洲自偷自拍三级| 久久久久网色| 亚洲国产精品999| 3wmmmm亚洲av在线观看| 王馨瑶露胸无遮挡在线观看| 日本91视频免费播放| 男女啪啪激烈高潮av片| 91在线精品国自产拍蜜月| 亚洲成人一二三区av| 国产伦在线观看视频一区| 日本色播在线视频| 亚洲情色 制服丝袜| 能在线免费看毛片的网站| 精品视频人人做人人爽| 最新的欧美精品一区二区| 久久久国产一区二区| 亚洲av综合色区一区| 亚洲成人av在线免费| av免费在线看不卡| 大码成人一级视频| 久久精品久久久久久久性| 国产男女超爽视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 夜夜骑夜夜射夜夜干| 亚洲第一av免费看| 欧美xxⅹ黑人| 国产国拍精品亚洲av在线观看| 偷拍熟女少妇极品色| 亚洲国产欧美日韩在线播放 | 久久久久人妻精品一区果冻| 国产中年淑女户外野战色| 国产淫语在线视频| 久久99蜜桃精品久久| 色视频在线一区二区三区| 亚洲国产欧美日韩在线播放 | 一二三四中文在线观看免费高清| 免费看av在线观看网站| 少妇高潮的动态图| 久久久久精品久久久久真实原创| av在线播放精品| 成年av动漫网址| 免费看光身美女| 人体艺术视频欧美日本| 日韩三级伦理在线观看| 男女边吃奶边做爰视频| 99久久精品国产国产毛片| 久久国产精品男人的天堂亚洲 | 美女福利国产在线| 99热6这里只有精品| av有码第一页| 少妇被粗大猛烈的视频| 久久午夜福利片| 又黄又爽又刺激的免费视频.| 美女大奶头黄色视频| 成人18禁高潮啪啪吃奶动态图 | 久久综合国产亚洲精品| 香蕉精品网在线| 久久毛片免费看一区二区三区| 深夜a级毛片| 午夜91福利影院| 美女中出高潮动态图| 这个男人来自地球电影免费观看 | 日韩欧美 国产精品| 各种免费的搞黄视频| 下体分泌物呈黄色| 人人妻人人爽人人添夜夜欢视频 | 久久 成人 亚洲| 伊人久久国产一区二区| 永久免费av网站大全| 久久午夜综合久久蜜桃| 亚洲av不卡在线观看| 亚洲中文av在线| 十分钟在线观看高清视频www | 亚洲图色成人| 三级国产精品欧美在线观看| 国产精品欧美亚洲77777| 国产无遮挡羞羞视频在线观看| 黄色日韩在线| 18禁裸乳无遮挡动漫免费视频| 麻豆成人午夜福利视频| 精品一区二区三卡| 少妇人妻一区二区三区视频| 嫩草影院入口| 成人美女网站在线观看视频| 22中文网久久字幕| 欧美日韩视频高清一区二区三区二| 久热久热在线精品观看| 人人妻人人澡人人爽人人夜夜| 亚州av有码| 国产精品一区www在线观看| 亚洲第一区二区三区不卡| 亚洲精品亚洲一区二区| 91在线精品国自产拍蜜月| 涩涩av久久男人的天堂| 亚洲成色77777| 日本色播在线视频| 婷婷色综合www| 观看免费一级毛片| 欧美xxⅹ黑人| 又粗又硬又长又爽又黄的视频| 亚洲精品第二区| 午夜福利,免费看| av福利片在线| 国语对白做爰xxxⅹ性视频网站| 日韩精品免费视频一区二区三区 | 又大又黄又爽视频免费| 亚洲婷婷狠狠爱综合网| 日本欧美视频一区| 91久久精品电影网| 成人影院久久| 99热这里只有精品一区| 久久久精品免费免费高清| 国产精品国产三级专区第一集| 一区二区三区四区激情视频| 亚洲精品乱码久久久v下载方式| 欧美最新免费一区二区三区| 一级a做视频免费观看| 黄色视频在线播放观看不卡| 久久狼人影院| 国产黄频视频在线观看| 人人澡人人妻人| 成年av动漫网址| 寂寞人妻少妇视频99o| 精品一区在线观看国产| 婷婷色av中文字幕| 亚洲经典国产精华液单| 最近最新中文字幕免费大全7| 国产高清有码在线观看视频| 午夜精品国产一区二区电影| 一级毛片黄色毛片免费观看视频| 亚洲欧洲精品一区二区精品久久久 | 精品酒店卫生间| 亚洲av成人精品一二三区| 亚洲四区av| 欧美激情国产日韩精品一区| 国产伦精品一区二区三区视频9| 亚洲国产精品成人久久小说| 人妻一区二区av| 99热网站在线观看| av女优亚洲男人天堂| 十八禁网站网址无遮挡 | 欧美日韩亚洲高清精品| 熟女电影av网| 午夜免费鲁丝| 美女中出高潮动态图| 亚洲欧美精品专区久久| 日韩免费高清中文字幕av| 国产伦在线观看视频一区| 丰满人妻一区二区三区视频av| 少妇熟女欧美另类| av福利片在线观看| 欧美少妇被猛烈插入视频| 亚洲国产成人一精品久久久| 大话2 男鬼变身卡| 啦啦啦在线观看免费高清www| 免费看日本二区| 大话2 男鬼变身卡| 日本爱情动作片www.在线观看| 观看免费一级毛片| 一级,二级,三级黄色视频| 美女福利国产在线| 日韩欧美一区视频在线观看 | 最近的中文字幕免费完整| 免费大片黄手机在线观看| 精品久久国产蜜桃| 性色av一级| 久久久国产精品麻豆| 久久精品国产自在天天线| 美女cb高潮喷水在线观看| 日日啪夜夜爽| 插逼视频在线观看| 少妇裸体淫交视频免费看高清| 男女无遮挡免费网站观看| 国产真实伦视频高清在线观看| 乱系列少妇在线播放| 国产精品一区二区在线观看99| 人妻夜夜爽99麻豆av| 激情五月婷婷亚洲| 丝瓜视频免费看黄片| 日韩人妻高清精品专区| 午夜91福利影院| 国产欧美日韩精品一区二区| 国内揄拍国产精品人妻在线| 老司机影院成人| 国产成人freesex在线| 丰满人妻一区二区三区视频av| 亚洲av不卡在线观看| 国产精品国产三级专区第一集| 中文字幕精品免费在线观看视频 | 亚洲精品久久午夜乱码| 新久久久久国产一级毛片| 亚洲欧美日韩另类电影网站| 亚洲不卡免费看| freevideosex欧美| 久久久久网色| 免费av不卡在线播放| 97在线视频观看| videos熟女内射| 少妇被粗大猛烈的视频| 久久人人爽人人片av| 欧美区成人在线视频| 午夜免费鲁丝| 久久毛片免费看一区二区三区| 黑人高潮一二区| 国产深夜福利视频在线观看| 国产真实伦视频高清在线观看| 少妇 在线观看| 亚洲欧美精品自产自拍| 成人影院久久| 国产精品久久久久久久久免| 一级,二级,三级黄色视频| 国内精品宾馆在线| 国产av精品麻豆| 久久精品久久精品一区二区三区| 国产在线男女| 哪个播放器可以免费观看大片| 国产黄色视频一区二区在线观看| 九九在线视频观看精品| 午夜免费鲁丝| 欧美变态另类bdsm刘玥| 十分钟在线观看高清视频www | 在线观看免费视频网站a站| 久久精品国产自在天天线| 天堂俺去俺来也www色官网| 边亲边吃奶的免费视频| 美女cb高潮喷水在线观看| 国产黄片美女视频| 久久99热这里只频精品6学生| 少妇被粗大的猛进出69影院 | 一区二区三区免费毛片| www.av在线官网国产| 久久午夜综合久久蜜桃| 人人妻人人看人人澡| 在线观看人妻少妇| 国产极品粉嫩免费观看在线 | 国产精品人妻久久久影院| 高清毛片免费看| 国产极品粉嫩免费观看在线 | 久久韩国三级中文字幕| 精品人妻熟女毛片av久久网站| 欧美性感艳星| 亚洲精品一二三| av一本久久久久| 另类亚洲欧美激情| 啦啦啦在线观看免费高清www| 精品亚洲成国产av| 国产成人免费观看mmmm| 熟妇人妻不卡中文字幕| 亚洲天堂av无毛| 91精品国产九色| 国产一区二区在线观看日韩| 国产精品伦人一区二区| 最新的欧美精品一区二区| 99热6这里只有精品| freevideosex欧美| 国产日韩欧美视频二区| 国产精品国产三级专区第一集| 亚洲精品第二区| 国产永久视频网站| 色婷婷久久久亚洲欧美| 人妻夜夜爽99麻豆av| 精品一区二区三区视频在线| 欧美xxⅹ黑人| 女性被躁到高潮视频| 寂寞人妻少妇视频99o| 18禁在线无遮挡免费观看视频| 国产精品.久久久| 美女福利国产在线| 国产黄片美女视频| 日韩欧美精品免费久久| 高清午夜精品一区二区三区| 一本一本综合久久| 永久网站在线| 人妻少妇偷人精品九色| 久久人妻熟女aⅴ| 国模一区二区三区四区视频| 国产黄片美女视频| 亚洲婷婷狠狠爱综合网| 久久精品国产亚洲av涩爱| 9色porny在线观看| 欧美激情国产日韩精品一区| 精品少妇久久久久久888优播| 久久精品国产亚洲av涩爱| 国产成人a∨麻豆精品| 国产白丝娇喘喷水9色精品| 国产日韩一区二区三区精品不卡 | 国产成人免费无遮挡视频| 亚洲av二区三区四区| 一级二级三级毛片免费看| 亚洲精品一二三| 亚洲av福利一区| 欧美 日韩 精品 国产| 永久免费av网站大全| 热99国产精品久久久久久7| 高清视频免费观看一区二区| 噜噜噜噜噜久久久久久91| 久热这里只有精品99| 国产成人91sexporn| 丰满饥渴人妻一区二区三| 汤姆久久久久久久影院中文字幕| 久久久久精品性色| 欧美日本中文国产一区发布| 黑人猛操日本美女一级片| 中国美白少妇内射xxxbb| 色视频在线一区二区三区| 国产精品99久久久久久久久| 视频中文字幕在线观看| 欧美日本中文国产一区发布| 国产高清有码在线观看视频| 日韩熟女老妇一区二区性免费视频| 毛片一级片免费看久久久久| 久久久久久久久久人人人人人人| 插阴视频在线观看视频| 一级毛片黄色毛片免费观看视频| 日韩欧美精品免费久久| 午夜福利在线观看免费完整高清在| 美女cb高潮喷水在线观看| 成人毛片60女人毛片免费| 国产精品久久久久久精品电影小说| 日本黄大片高清| 一本一本综合久久| 午夜日本视频在线| 麻豆乱淫一区二区| 欧美日韩在线观看h| 国产91av在线免费观看| 日日撸夜夜添| av在线老鸭窝| 国产精品免费大片| a级片在线免费高清观看视频| 一区二区三区四区激情视频| 人人妻人人澡人人爽人人夜夜| videossex国产| 成年美女黄网站色视频大全免费 | av女优亚洲男人天堂| 亚洲精品视频女| 日本黄色片子视频|