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

    基于雙源IR-UWB生物雷達的強反射雜波抑制

    2017-01-19 07:47:20張芫蓓王帥杰王昊祁富貴王健琪李釗
    中國醫(yī)療設(shè)備 2016年11期
    關(guān)鍵詞:雙源雜波天線

    張芫蓓,王帥杰,王昊,祁富貴,王健琪,李釗

    第四軍醫(yī)大學 生物醫(yī)學工程系,陜西西安 710032

    基于雙源IR-UWB生物雷達的強反射雜波抑制

    張芫蓓,王帥杰,王昊,祁富貴,王健琪,李釗

    第四軍醫(yī)大學 生物醫(yī)學工程系,陜西西安 710032

    在基于生物雷達的生命探測技術(shù)研究中,強反射雜波干擾對雷達回波的信噪比影響十分重要,普通的信號處理方法很難將其去除,從而使雷達系統(tǒng)發(fā)生誤判或者漏判。本文在現(xiàn)有超寬帶生物雷達技術(shù)的基礎(chǔ)上,提出了一種基于雙源IR-UWB生物雷達的強反射雜波抑制方法。利用兩種不同中心頻率的天線同時進行探測,首先對不同頻率雷達回波信號分別進行通道內(nèi)處理,實現(xiàn)人體微弱生命信號的增強;然后進行通道間處理,并采用自適應(yīng)雜波消除技術(shù)抑制雷達回波中的強反射雜波干擾。在穿透磚墻和廢墟兩種不同場景下進行實驗,結(jié)果表明本文所提出的方法可以明顯提高雷達探測的準確率。

    超寬帶;生物雷達;生命探測;自適應(yīng)雜波消除技術(shù);雜波抑制

    引言

    超寬帶(Ultra-Wide Band,UWB)生物雷達融合了雷達技術(shù)和生物醫(yī)學工程技術(shù),通過發(fā)射出的電磁波穿透非金屬介質(zhì)(木板、墻壁、廢墟等),檢測到生命體目標的體動,或者呼吸、心跳等生命體征引起的體表微動,并以此為依據(jù)探測和識別生命體目標。與基于紅外、光學和超聲的探測技術(shù)相比,UWB生物雷達因其具有非接觸、穿透性強、能夠獲得目標距離信息、抗干擾能力強等特點,可廣泛應(yīng)用于醫(yī)學、軍事、反恐等領(lǐng)域,特別是在地震、塌方等災(zāi)害發(fā)生后的應(yīng)急救援中,具有不可替代的優(yōu)勢[1-5]。

    目前的UWB生物雷達技術(shù)多采用沖激脈沖(Impulse Radio,IR)雷達,它是以固定的脈沖重復(fù)頻率發(fā)射超短脈沖來實現(xiàn)超寬帶的[6-7]。IR-UWB生物雷達在地震、塌方等災(zāi)后應(yīng)急救援場景中,通常采用小時窗分段探測技術(shù)來提高雷達的探測精度和回波信噪比[8]。但當目標壓埋較深,距離雷達較遠時,使用這種技術(shù)就會增加探測時間,影響探測效率。因此,為了能夠更加快速準確地發(fā)現(xiàn)幸存人員,應(yīng)當使用大時窗探測技術(shù)。

    但是當生物雷達進行大時窗探測時,由于探測環(huán)境十分復(fù)雜,雷達發(fā)射的電磁波通常需要穿透多種不同介質(zhì)才能照射到被壓埋的人體目標。當IR-UWB生物雷達發(fā)射的電磁波穿過兩種介電常數(shù)不同的介質(zhì)(如磚和空氣)時,由于雷達硬件電路不可能工作在理想條件下,發(fā)射脈沖必然會產(chǎn)生細微抖動,原本應(yīng)該靜止的不同介質(zhì)交界面處的雷達回波,將會隨時間不斷變化,產(chǎn)生強反射雜波[9-10]。有時這種強反射雜波和人體的呼吸信號十分相近,普通的雷達回波信號處理方法很難將其去除,從而導(dǎo)致雷達系統(tǒng)發(fā)生誤判或者漏判,浪費寶貴的救援時間和資源,影響對幸存人員的搜救[11]。

    在以往的IR-UWB生物雷達探測研究中,為了消除這種由強反射障礙物引起的非靜態(tài)雜波,有研究通過在人體目標身后放置吸波材料來消除磚墻、地面等帶來的強雜波干擾[12-13],但是這種方法只能在實驗研究中運用,無法應(yīng)用在實際的災(zāi)后搜救場景中;還有研究通過改變雷達回波的起始位置將這種非靜態(tài)雜波干擾調(diào)整到雷達接收時窗之外[14],但這種方法會導(dǎo)致探測時出現(xiàn)盲區(qū),影響雷達系統(tǒng)的實用性;還有些研究利用特殊的雷達回波信號處理算法,對這種非靜態(tài)雜波的消除具有一定的效果,但是當障礙物距離人體目標非常近時(如人體靠墻站立或平躺在地面上),由障礙物發(fā)射產(chǎn)生的非靜態(tài)雜波和人體目標的呼吸信號就變得更加難以區(qū)分[15-16]。

    本文提出了一種基于雙源IR-UWB生物雷達的強反射雜波抑制方法,利用兩種不同中心頻率的天線,通過等效時間采樣和分時復(fù)用技術(shù)同時對目標進行探測,再配合自適應(yīng)雜波消除算法,能夠有效地抑制這種強反射雜波干擾,提高雷達判斷的準確率。

    1 基于雙源IR-UWB生物雷達的通道內(nèi)處理算法

    雙源IR-UWB生物雷達的發(fā)射天線發(fā)射雷達脈沖,雷達脈沖被人體目標反射,再通過雙源IR-UWB生物雷達的接收天線接收到被人體目標反射后的雷達脈沖,即雷達原始回波信號。雙源IR-UWB生物雷達包含兩對具有不同中心頻率的天線,每對天線均包括一個發(fā)射天線和一個接收天線,兩個接收天線接收到兩個具有不同中心頻率的雷達原始回波信號R1和R2。分別對這兩個通道的回波信號做通道內(nèi)處理,具體步驟如下:

    1.1 距離累積

    用R表示任意一個通道內(nèi)的雷達原始回波信號,R可以表示為二維矩陣R(m,n)。雷達原始回波信號二維矩陣示意圖,見圖1。圖 中橫軸表示探測時間,稱為“慢時間”,沿著探測時間的方向為慢時間方向,單位是s;縱軸表示探測距離,稱為“快時間”,沿著探測距離的方向為快時間方向,單位是ns。在某一時刻,沿著快時間方向的信號,即二維矩陣的列向量,叫做“道信號”;在某一距離點上,沿著慢時間方向的信號,即二維矩陣的行向量,叫做“點信號”。由于快時間方向上鄰近距離點處的雷達回波的調(diào)制方式大致相同,且具有一定的相關(guān)性,因此,在不影響有用信息的前提下,首先對二維矩陣R(m,n)在空間域沿著快時間方向進行距離累積,即:

    式中,R1(l,n) (l=1,2,...,L)為距離累積后的雷達回波數(shù)據(jù),Q為沿著快時間方向累積的窗寬,L為累積后在快時間方向的距離點數(shù),且其中表示向下取整,設(shè)定窗寬Q=40,采樣點數(shù)M為8192,則L為200,從而大大降低了雷達數(shù)據(jù)處理的運算量,減少了探測所需的運算時間,提高了廢墟搜救幸存人員的工作效率。

    圖1 雷達原始回波信號二維矩陣示意圖

    1.2 歸一化

    距離累積后的雷達回波信號R1在時間域沿著慢時間方向進行歸一化處理,得到歸一化后的雷達回波信號R2。為了增強距離雷達較遠處目標信號的幅度,補償信號在穿透和傳播過程中的衰減,提高雷達回波信號的信噪比,對距離累積后的信號R1,在時間域沿著慢時間方向進行點信號的歸一化處理,即:

    式中,R2(l,n)為歸一化后的雷達回波數(shù)據(jù),N為R1中包含的道信號的個數(shù)。雷達回波數(shù)據(jù)R1(l,n)經(jīng)過歸一化處理后,每個距離上的點信號的幅值范圍都在-1到1之間。

    1.3 線性趨勢消除

    對雷達回波信號R2進行線性趨勢消除,得到雷達回波信號R3。由于IR-UWB生物雷達的硬件不可能在理想條件下工作,采集數(shù)據(jù)過程中往往伴隨著回波基線的漂移。因此,本文采用線性趨勢消除(Linear Trend Subtraction,LTS)來去除雷達回波信號R2中的靜態(tài)雜波以及線性漂移,LTS通過線性最小二乘擬合估計出雷達回波信號R2在慢時間方向上的直流分量和線性趨勢后,再從回波數(shù)據(jù)中減去,即:

    式中,R3表示LTS后的雷達回波信號,

    y = [n/N,IN],n = [0,1,2,...,N-1]T,IN為一個長度是N且元素都是1的列向量,N為R2中包含的道信號的個數(shù)。

    1.4 低通濾波

    對雷達回波信號R3進行低通濾波,得到濾波后的雷達回波信號R4。由于IR-UWB生物雷達的硬件在工作過程中不可避免的會產(chǎn)生高頻噪聲,而且人體目標的呼吸信號又是一個窄帶的準周期信號,因此,為了進一步提高雷達回波的信噪比,就需要在慢時間方向上對雷達回波信號R3的點信號進行低通濾波:

    式中,R4(l,n)為濾波后的雷達回波數(shù)據(jù),“*”表示卷積運算,h(t)為有限沖擊響應(yīng)(Finite Impulse Response, FIR)濾波器的沖擊函數(shù)。根據(jù)人體的呼吸頻率,該算法中低通濾波器的截止頻率為0.5 Hz。

    按照上述4個步驟,分別對兩個通道內(nèi)的雷達原始回波信號R1和R2進行通道內(nèi)處理,得到處理后的雷達回波信號R41和R42。

    2 基于雙源IR-UWB生物雷達的通道間處理算法

    本文中雙源IR-UWB生物雷達包含中心頻率分別為270 MHz和400 MHz的兩對蝶型天線,采用等效時間采樣、多通道分時復(fù)用技術(shù),實現(xiàn)不同中心頻率天線每發(fā)射出一道信號后交替工作,從而使兩種中心頻率的天線在工作時相互不產(chǎn)生干擾。由于兩對不同中心頻率的天線分時復(fù)用交替探測,且交替時間間隔非常短,可近似認為兩對天線同時對目標進行探測,兩對不同頻率的雷達回波中人體呼吸信號具有一定的相關(guān)性;但是兩對天線的脈沖發(fā)射源是相互獨立的,因此產(chǎn)生的脈沖抖動是隨機的,不具有相關(guān)性。

    R41和R42分別是中心頻率為270 MHz和400 MHz通道內(nèi)處理后的雷達回波信號。對R41和R42進行自適應(yīng)雜波消除,得到人體呼吸信號的估計值S(k)。基于雙源IRUWB生物雷達的通道間處理具體算法的原理,見圖2。主輸入信號y(k)為提取的雷達回波信號R41的點信號,由人體目標的呼吸信號S1(k)和非靜態(tài)雜波n1(k)組成,而參考信號x(k)為提取的雷達回波信號R42的點信號,由人體目標的呼吸信號S2(k)和非靜態(tài)雜波n2(k)組成。由于兩對不同中心頻率的天線同時對人體目標進行探測,因此呼吸信號S1(k)和S2(k)具有一定的相關(guān)性,又由于兩對天線的脈沖發(fā)射源是相互獨立的,因此產(chǎn)生的脈沖抖動是隨機的,非靜態(tài)雜波n1(k)和n2(k)不具有相關(guān)性。利用誤差信號e(k),通過最小均方算法調(diào)整FIR濾波器的參數(shù),使濾波器的輸出信號S(k)為相關(guān)的人體呼吸信號的估計值,將雷達回波信號中的強反射雜波等非靜態(tài)雜波去除,計算步驟如下:

    其中,F(xiàn)IR濾波器的階數(shù)L=50,步長因子μ=10-5,權(quán)系數(shù)向量ωi(k)的初始值全是0,k表示輸出信號S(k)中的數(shù)據(jù)序號。

    圖2 自適應(yīng)雜波消除算法的原理框圖

    3 雙源IR-UWB生物雷達探測實驗

    3.1 穿墻實驗

    穿墻探測場景,見圖3。房間內(nèi)兩堵平行的磚墻厚度均為24 cm,將雙源IR-UWB生物雷達系統(tǒng)擺放在其中一堵墻后,對無人的房間進行探測。探測完成后,一名健康成年男子緊靠另一堵墻且正對雷達站立,平靜呼吸,身體保持靜止,雷達與人體目標之間的距離約為3 m,對其進行探測。

    圖3 穿墻探測場景示意圖

    為了能夠更加直觀地識別人體目標的呼吸信號及其所在位置,對處理后的雷達回波數(shù)據(jù)沿慢時間方向計算每個距離上的點信號的方差,得到雷達回波信號的能量譜。在有人體目標和無人體目標兩種情況時,采用普通IR-UWB生物雷達處理方法(基于單源IR-UWB生物雷達的常規(guī)處理算法[16])得到的雷達回波數(shù)據(jù)能量譜,見圖4。從圖4(a)中可以看出,3 m附近的能量值明顯大于其余位置的能量值,這正是由人體目標的呼吸引起的。而圖4(b)中,3 m附近也出現(xiàn)了很明顯的能量峰值,這是因為IR-UWB雷達在探測過程中,人體身后的墻產(chǎn)生了強反射雜波干擾,它與人體目標的呼吸信號十分相似,無法通過普通的IRUWB生物雷達處理方法濾除。因此,通過普通的生物雷達探測技術(shù),有人體目標和無人體目標兩種情況的雷達回波數(shù)據(jù)都被判斷為探測區(qū)域內(nèi)有人體目標,且目標距離都約為3 m,圖4(b)的回波數(shù)據(jù)產(chǎn)生誤判。

    圖4 穿墻場景采用普通方法處理的雷達回波數(shù)據(jù)能量譜

    在有人體目標和無人體目標兩種情況時,采用基于雙源IR-UWB生物雷達的強反射雜波抑制方法得到的雷達回波數(shù)據(jù)能量譜,見圖5。從圖5(a)中可以看出,只有3 m附近的能量峰值較強,其最大值在2.88 m處,正好是人體目標所站的位置。而圖5(b)中有多處較強的能量峰值,這是由于能量譜進行歸一化后,原本很小的峰值都變得比較明顯。實際上,出現(xiàn)在圖4(b)中近似呼吸的強反射雜波已經(jīng)被較好的抑制。因此,圖5(b)中雷達回波數(shù)據(jù)的判斷結(jié)果為探測區(qū)域內(nèi)無人體目標,與實際情況相符。

    圖5 穿墻場景采用基于雙源IR-UWB生物雷達的強反射雜波抑制方法處理的雷達回波數(shù)據(jù)能量譜

    3.2 穿廢墟實驗

    穿廢墟探測場景,見圖6。模擬廢墟高度約為3 m,包括兩個10 cm厚的混凝土預(yù)制板,兩段厚度約為1 m的磚塊堆,還有一個高度約為80 cm的廢墟空洞。IR-UWB雷達系統(tǒng)放置在模擬廢墟頂端,方向朝下進行探測。探測完成后,一名健康成年男子平躺在廢墟空洞中,胸壁正對雷達,平靜呼吸,身體保持靜止,雷達與人體目標之間的距離接近3 m,對其進行探測。

    在有人體目標和無人體目標兩種情況時,采用普通IRUWB生物雷達處理方法得到的雷達回波數(shù)據(jù)能量譜,見圖7。兩幅圖的能量波形很相似,都存在多處較強的能量峰值,這主要是因為IR-UWB雷達系統(tǒng)在穿透模擬廢墟探測過程中,由于雷達脈沖的抖動,在兩種介質(zhì)(預(yù)制板和空氣、磚和空氣)的交界面處原本應(yīng)該靜止的雷達回波卻隨慢時間不停變化,產(chǎn)生近似呼吸的強反射雜波干擾。因此,兩幅圖的雷達回波數(shù)據(jù)判斷結(jié)果都是探測區(qū)域內(nèi)無人體目標,圖7(b)的探測結(jié)果發(fā)生了漏判。

    圖6 穿廢墟探測場景示意圖

    圖7 穿廢墟場景采用普通方法處理的雷達回波能量譜

    在有人體目標和無人體目標兩種情況時,采用基于雙源IR-UWB生物雷達的強反射雜波抑制方法得到的雷達回波數(shù)據(jù)能量譜,見圖8。雖然圖8(a)中仍存在多處較強的能量峰值,但實際上由廢墟交界面產(chǎn)生的近似呼吸的強反射雜波已經(jīng)被大大的削弱了,只是能量歸一化使得原本幅度很小的能量峰值又變得比較明顯。因此,這幅圖的回波數(shù)據(jù)判別結(jié)果為探測區(qū)域內(nèi)無人體目標,與實際情況相符。圖8(b)中只剩下一處較強的能量峰值,這正是人體目標的呼吸能量,其能量最大值在4.05 m處。由于該坐標距離計算的是電磁波在空氣中的傳播距離,而電磁波在該模擬廢墟中的衰減很大,經(jīng)過距離校準后的實際位置約為2.83 m,與人體目標所躺的位置相符。

    4 結(jié)語

    在基于IR-UWB生物雷達的生命探測技術(shù)研究中,強反射雜波干擾對雷達回波的信噪比影響十分重要。本文采用一種雙源IR-UWB生物雷達配合自適應(yīng)雜波消除方法能夠較好地抑制不同場景中的強反射雜波干擾,提高了生物雷達系統(tǒng)判斷的準確率,為該領(lǐng)域的進一步發(fā)展提供了新的研究思路和方法。

    圖8 穿廢墟場景采用基于雙源IR-UWB生物雷達的強反射雜波抑制方法處理的雷達回波能量譜

    [1] Kumar A,Li Z,Liang Q,et al.Experimental study of throughwall human detection using ultra wideband radar sensors[J]. Measurement,2014,47(1):869-879.

    [2] Lou J,Wang Y.Implement of UWB Radar for Human Life Detection[J].Appl Mech Mater,2012,239-240:348-352.

    [3] Liu L,Liu Z,Xie H,et al.Numerical simulation of UWB impulse radar vital sign detection at an earthquake disaster site[J].Ad Hoc Networks,2014,13(1):34-41.

    [4] Wang Y,Yu X,Zhang Y,et al.Detecting and monitoring the micro-motions of trapped people hidden by obstacles based on wavelet entropy with low centre-frequency UWB radar[J].Int J Remote Sens,2015,36(5):1349-1366.

    [5] Ren L,Yun SK,Wang H,et al.Noncontact Multiple Heartbeats Detection and Subject Localization Using UWB Impulse Doppler Radar[J].IEEE Microw Wireless Compon Lett,2015, 25(10):690-692.

    [6] Gezici S.Theoretical Limits for Estimation of Periodic Movements in Pulse-Based UWB Systems[J].IEEE J Sel Top Signal Process,2007,1(3):405-417.

    [7] Lazaro A,Girbau D,Villarino R.Analysis of Vital Signs Monitoring Using an Ir-UWB Radar[J].Prog Electromagn Res, 2010,100(1):265-284.

    [8] 安強,李釗,呂昊,等.UWB搜救生物雷達分段時窗探測系統(tǒng)的實現(xiàn)[J].中國醫(yī)療設(shè)備,2015,30(11):18-21.

    [9] Nezirovic A.Stationary clutter-and linear-trend suppression in impulse-radar-based respiratory motion detection[A].Ultra-Wideband (ICUWB),2011 IEEE International Conference on[C].New York:IEEE,2011:331-335.

    [10] Crabbe S.Detection and localization of persons behind obstacles using M-sequence through-the-wall radar[J].Proc Spie,2006(6201):62010I-62010I-12.

    [11] 梁福來,李釗,安強,等.脈沖式超寬帶生命探測雷達時間抖動抑制[J].雷達學報,2015,4(4):439-444.

    [12] Nezirovic A,Yarovoy AG,Ligthart LP.Experimental study on human being detection using UWB radar[A].2006 International Radar Symposium[C].New York:IEEE,2006:1-4.

    [13] Kiriazi JE,Boric-Lubecke O,Lubecke VM.Dual-Frequency Technique for Assessment of Cardiopulmonary Effective RCS and Displacement[J].IEEE Sensors J,2012,12(3):574-582.

    [14] Tenenbaum B,Karshan M,Ziskin D,et al.Adaptive Ground Clutter Removal Algorithm for Ground Penetrating Radar Applications in Harsh Environments[J].Sens Imag Int J,2006, 7(3):71-89.

    [15] Singh S,Liang Q,Chen D,et al.Sense through wall human detection using UWB radar[J].EURASIP JWCN,2011,(1):1-11.

    [16] Lv H,Lu GH,Jing XJ,et al.A new ultra-wideband radar for detecting survivors buried under earthquake rubbles[J].Microw Opt Technol Lett,2010,52(11):2621-2624.

    Suppression of Strong Clutter Based on the Dual-Frequency IR-UWB Bioradar

    ZHANG Yuan-bei, WANG Shuai-jie, WANG Hao, QI Fu-gui, WANG Jian-qi, LI Zhao
    Faculty of Biomedical Engineering, the Fourth Military Medical University, Xi'an Shaanxi 710032, China

    In the study of life detection based on bioradar, strong clutter has a very important influence on the signal-to-noise ratio of radar echo. And it is difficult to be removed by using traditional signal processing methods, so that false alarm or missing result usually occurs in the complex environment. In the light of existing ultra-wideband (UWB) radar technology, this paper put forward a suppression method of strong clutter based on the dual-frequency IR-UWB radar, which could detect target at the same time by using two pairs of antennas with different center frequencies. At first, the radar echo signal was respectively processed in different channels to enhance the weak vital signal of human. Then in order to suppress the interference of strong clutters, the data of two channels were processed together with adaptive clutter elimination technology. Experiments were conducted in two different scenarios of brick wall and simulation ruins, and the results showed that the accuracy of radar detection could be significantly improved by using the method proposed in this paper.

    ultra-wide band; bioradar; life detection; adaptive clutter elimination technology; clutter suppression

    R318;TN957.51

    A

    10.3969/j.issn.1674-1633.2016.11.006

    1674-1633(2016)11-0028-05

    2016-09-18

    2016-09-23

    陜西省自然科學基礎(chǔ)研究計劃青年人才項目(2014JQ2-6034)。

    作者郵箱:lizhaofmmu@fmmu.edu.cn

    猜你喜歡
    雙源雜波天線
    STAR2000型空管一次雷達雜波抑制淺析
    雙源采購跳躍-擴散庫存控制模型
    自動化學報(2018年2期)2018-04-12 05:46:16
    ETC相控陣天線與普通天線應(yīng)用對比分析
    ALLESS轉(zhuǎn)動天線射頻旋轉(zhuǎn)維護與改造
    電子制作(2016年1期)2016-11-07 08:42:54
    理論宣講要上接天線、下接地氣
    學習月刊(2016年19期)2016-07-11 01:59:46
    密集雜波環(huán)境下確定性退火DA-HPMHT跟蹤算法
    雙源信號深組合導(dǎo)航在新一代運載火箭制導(dǎo)系統(tǒng)中應(yīng)用初探
    相關(guān)廣義復(fù)合分布雷達海雜波仿真
    遙測遙控(2015年2期)2015-04-23 08:15:22
    微波雷達海雜波混沌特性分析
    彈載北斗抗干擾天線系統(tǒng)分析與設(shè)計
    精品久久久久久久久亚洲 | 此物有八面人人有两片| 成年女人永久免费观看视频| 99久久精品国产国产毛片| 国产精品伦人一区二区| 搡老熟女国产l中国老女人| 国产高潮美女av| 99热网站在线观看| 欧美成人一区二区免费高清观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美日韩东京热| 久久久久久九九精品二区国产| www.色视频.com| av视频在线观看入口| 国内精品久久久久久久电影| 亚洲国产日韩欧美精品在线观看| av天堂在线播放| 色哟哟·www| 亚洲精品影视一区二区三区av| 极品教师在线视频| 欧美不卡视频在线免费观看| 99热网站在线观看| 欧美成人一区二区免费高清观看| 国产在线精品亚洲第一网站| 我的老师免费观看完整版| 久久久久久久久大av| 啪啪无遮挡十八禁网站| 国内久久婷婷六月综合欲色啪| 人妻制服诱惑在线中文字幕| 亚洲成人中文字幕在线播放| 男女之事视频高清在线观看| 欧美极品一区二区三区四区| 国产精品1区2区在线观看.| 一级a爱片免费观看的视频| 国产精品自产拍在线观看55亚洲| 久久精品国产99精品国产亚洲性色| 亚洲精华国产精华液的使用体验 | 在线观看66精品国产| 精品久久久久久久久亚洲 | 国产极品精品免费视频能看的| 国产免费av片在线观看野外av| 男女视频在线观看网站免费| 日韩欧美国产在线观看| av天堂中文字幕网| 88av欧美| 国产美女午夜福利| 国模一区二区三区四区视频| 亚洲av成人精品一区久久| 日韩中文字幕欧美一区二区| 免费电影在线观看免费观看| 嫁个100分男人电影在线观看| 高清在线国产一区| 又紧又爽又黄一区二区| 精品乱码久久久久久99久播| 一边摸一边抽搐一进一小说| 欧美国产日韩亚洲一区| 日本黄色视频三级网站网址| 亚洲熟妇中文字幕五十中出| 在线免费观看不下载黄p国产 | 久久久久久久精品吃奶| 又爽又黄a免费视频| 国产精品98久久久久久宅男小说| 午夜a级毛片| 一夜夜www| 在线看三级毛片| 18+在线观看网站| 此物有八面人人有两片| 91在线精品国自产拍蜜月| 蜜桃亚洲精品一区二区三区| 亚洲欧美日韩东京热| 3wmmmm亚洲av在线观看| 日日干狠狠操夜夜爽| 久久草成人影院| 一个人看视频在线观看www免费| 亚洲国产精品sss在线观看| 99九九线精品视频在线观看视频| 久久精品国产亚洲网站| 一个人看视频在线观看www免费| 18禁黄网站禁片免费观看直播| 观看美女的网站| 成人二区视频| 在线a可以看的网站| 丰满的人妻完整版| 在现免费观看毛片| 男女那种视频在线观看| 黄色一级大片看看| 99riav亚洲国产免费| 3wmmmm亚洲av在线观看| 亚洲无线在线观看| 在现免费观看毛片| 成人亚洲精品av一区二区| 欧美+日韩+精品| 国产精品久久久久久亚洲av鲁大| 亚洲精品在线观看二区| 黄色女人牲交| av视频在线观看入口| 日韩欧美在线二视频| 小说图片视频综合网站| 久久精品综合一区二区三区| 老司机午夜福利在线观看视频| 看免费成人av毛片| 久久久久久久久久黄片| 美女 人体艺术 gogo| 久久久久免费精品人妻一区二区| 久久精品国产亚洲av香蕉五月| 亚洲精品久久国产高清桃花| 亚洲av成人精品一区久久| 我的老师免费观看完整版| 麻豆国产97在线/欧美| 亚洲成人久久爱视频| 少妇的逼水好多| 日韩精品青青久久久久久| 我要搜黄色片| 国产 一区 欧美 日韩| 欧美日韩黄片免| 哪里可以看免费的av片| 久久久国产成人精品二区| 婷婷亚洲欧美| 亚洲熟妇熟女久久| 国产成人aa在线观看| 在线播放无遮挡| 亚洲精品日韩av片在线观看| 亚洲精华国产精华精| 十八禁国产超污无遮挡网站| 国产精品综合久久久久久久免费| 国内毛片毛片毛片毛片毛片| 老师上课跳d突然被开到最大视频| av.在线天堂| 国产精品日韩av在线免费观看| 亚洲欧美激情综合另类| 波野结衣二区三区在线| 成年女人永久免费观看视频| 亚洲久久久久久中文字幕| 国产高清不卡午夜福利| 精品久久久久久久久av| 欧美极品一区二区三区四区| 永久网站在线| 最近视频中文字幕2019在线8| 91精品国产九色| 黄色日韩在线| 欧美一区二区国产精品久久精品| 老司机福利观看| 久久香蕉精品热| АⅤ资源中文在线天堂| 女的被弄到高潮叫床怎么办 | 一区二区三区四区激情视频 | 91麻豆av在线| 国产亚洲精品久久久久久毛片| 亚洲av成人av| 老司机福利观看| 日本黄色视频三级网站网址| 欧美日韩乱码在线| 午夜精品一区二区三区免费看| 波多野结衣高清无吗| 制服丝袜大香蕉在线| 国产精品永久免费网站| 波多野结衣高清无吗| 亚洲国产精品成人综合色| 99热这里只有是精品50| 婷婷精品国产亚洲av在线| 色播亚洲综合网| 亚洲av成人精品一区久久| 日日啪夜夜撸| 嫁个100分男人电影在线观看| 小蜜桃在线观看免费完整版高清| 精品久久久久久久久亚洲 | 欧美日本亚洲视频在线播放| 啪啪无遮挡十八禁网站| 国产探花在线观看一区二区| 久久精品夜夜夜夜夜久久蜜豆| 午夜免费激情av| 精品久久久久久久人妻蜜臀av| 久久国产乱子免费精品| 简卡轻食公司| 在线观看免费视频日本深夜| 18禁在线播放成人免费| 国产成人a区在线观看| 国产av在哪里看| 亚洲精华国产精华精| 国产av麻豆久久久久久久| 日韩欧美精品v在线| a在线观看视频网站| 国产精品人妻久久久久久| 乱码一卡2卡4卡精品| av天堂中文字幕网| 国国产精品蜜臀av免费| 女的被弄到高潮叫床怎么办 | 亚洲一区高清亚洲精品| 99久久久亚洲精品蜜臀av| 免费大片18禁| 免费人成在线观看视频色| 亚洲欧美日韩无卡精品| 可以在线观看的亚洲视频| 一个人看视频在线观看www免费| 日韩中字成人| 三级毛片av免费| 国产精品国产高清国产av| 国产精品人妻久久久影院| 最新在线观看一区二区三区| 女同久久另类99精品国产91| 国产伦人伦偷精品视频| 小蜜桃在线观看免费完整版高清| 国产高清激情床上av| 99精品久久久久人妻精品| 少妇高潮的动态图| 欧美性猛交╳xxx乱大交人| 欧美国产日韩亚洲一区| 美女大奶头视频| 日本欧美国产在线视频| 国产成人一区二区在线| 哪里可以看免费的av片| 乱码一卡2卡4卡精品| 人人妻人人澡欧美一区二区| 1024手机看黄色片| 三级男女做爰猛烈吃奶摸视频| 久久久久久国产a免费观看| 欧美高清性xxxxhd video| 成人av一区二区三区在线看| av视频在线观看入口| 欧美潮喷喷水| 我的女老师完整版在线观看| 亚洲男人的天堂狠狠| 大型黄色视频在线免费观看| 国内精品宾馆在线| 熟妇人妻久久中文字幕3abv| 国产精品三级大全| 国产av一区在线观看免费| 国产av在哪里看| 免费观看的影片在线观看| 亚洲无线观看免费| 久久久久久久久久黄片| 亚洲男人的天堂狠狠| 男人舔奶头视频| 18禁裸乳无遮挡免费网站照片| 欧美日韩国产亚洲二区| 亚洲国产欧洲综合997久久,| 欧美国产日韩亚洲一区| 欧美日本亚洲视频在线播放| 亚洲最大成人手机在线| h日本视频在线播放| 亚洲狠狠婷婷综合久久图片| 国产欧美日韩精品一区二区| avwww免费| 51国产日韩欧美| av福利片在线观看| 丝袜美腿在线中文| 日韩国内少妇激情av| 有码 亚洲区| 亚洲图色成人| 欧美bdsm另类| 中出人妻视频一区二区| 日日啪夜夜撸| 51国产日韩欧美| 在线天堂最新版资源| 嫩草影院新地址| 日日夜夜操网爽| 婷婷六月久久综合丁香| 特大巨黑吊av在线直播| 大型黄色视频在线免费观看| 在现免费观看毛片| av.在线天堂| 免费人成视频x8x8入口观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av日韩精品久久久久久密| 国产aⅴ精品一区二区三区波| 深夜a级毛片| 国产美女午夜福利| 成人av一区二区三区在线看| 男人和女人高潮做爰伦理| 久久久久久久久久成人| 十八禁网站免费在线| 欧洲精品卡2卡3卡4卡5卡区| 搡老岳熟女国产| 国产麻豆成人av免费视频| 麻豆成人av在线观看| 少妇裸体淫交视频免费看高清| 亚洲三级黄色毛片| 精品乱码久久久久久99久播| 国产在视频线在精品| 美女大奶头视频| 午夜激情欧美在线| 大型黄色视频在线免费观看| 能在线免费观看的黄片| 12—13女人毛片做爰片一| 日韩欧美精品免费久久| h日本视频在线播放| 我要搜黄色片| 国产伦在线观看视频一区| 有码 亚洲区| 久久国产精品人妻蜜桃| 色综合亚洲欧美另类图片| 久久草成人影院| 神马国产精品三级电影在线观看| 久99久视频精品免费| 国产熟女欧美一区二区| 国产成人影院久久av| 伦理电影大哥的女人| 超碰av人人做人人爽久久| 看片在线看免费视频| 97超视频在线观看视频| 国产又黄又爽又无遮挡在线| 久久久久久久久大av| 亚洲avbb在线观看| 国内揄拍国产精品人妻在线| 女的被弄到高潮叫床怎么办 | 国内久久婷婷六月综合欲色啪| 97人妻精品一区二区三区麻豆| bbb黄色大片| 精品一区二区三区视频在线观看免费| 18禁裸乳无遮挡免费网站照片| 看十八女毛片水多多多| 日韩av在线大香蕉| 亚洲最大成人中文| 女人十人毛片免费观看3o分钟| 国产成人一区二区在线| 少妇裸体淫交视频免费看高清| 人人妻人人看人人澡| 长腿黑丝高跟| 久久久精品大字幕| 长腿黑丝高跟| 成人欧美大片| 色综合站精品国产| 成人综合一区亚洲| 精品乱码久久久久久99久播| 日本一本二区三区精品| 日韩欧美三级三区| 久久久精品大字幕| 久久人人爽人人爽人人片va| 精品久久久噜噜| 国产精品免费一区二区三区在线| 国产久久久一区二区三区| 99国产极品粉嫩在线观看| 高清毛片免费观看视频网站| 国产白丝娇喘喷水9色精品| 免费大片18禁| 欧美人与善性xxx| 日韩强制内射视频| netflix在线观看网站| 色哟哟·www| 精品一区二区三区av网在线观看| 深夜精品福利| 国产成人a区在线观看| 国产精品免费一区二区三区在线| 男女那种视频在线观看| 欧美精品国产亚洲| 大又大粗又爽又黄少妇毛片口| 亚洲av免费在线观看| av黄色大香蕉| 我的女老师完整版在线观看| 精品一区二区三区av网在线观看| 少妇的逼好多水| 亚洲av成人精品一区久久| 国产精品国产三级国产av玫瑰| 欧美成人免费av一区二区三区| 国产真实伦视频高清在线观看 | 久久精品夜夜夜夜夜久久蜜豆| 人妻久久中文字幕网| 日韩欧美精品v在线| 国产69精品久久久久777片| 亚洲真实伦在线观看| avwww免费| 亚洲avbb在线观看| 免费av毛片视频| 直男gayav资源| 久久精品夜夜夜夜夜久久蜜豆| 欧美极品一区二区三区四区| 黄色女人牲交| 又黄又爽又刺激的免费视频.| 午夜a级毛片| 特大巨黑吊av在线直播| 婷婷六月久久综合丁香| 久久国产精品人妻蜜桃| 99在线人妻在线中文字幕| 日韩大尺度精品在线看网址| 色综合站精品国产| 亚洲无线观看免费| av视频在线观看入口| 黄色一级大片看看| 亚洲精品久久国产高清桃花| 国产一区二区三区视频了| 久久午夜亚洲精品久久| 中文字幕av成人在线电影| 亚洲成人久久爱视频| 97碰自拍视频| 在线看三级毛片| 2021天堂中文幕一二区在线观| 国产日本99.免费观看| 2021天堂中文幕一二区在线观| 国产视频内射| 精品人妻1区二区| 国产日本99.免费观看| 午夜影院日韩av| 欧美中文日本在线观看视频| 久久精品久久久久久噜噜老黄 | 国产精品三级大全| 亚洲国产日韩欧美精品在线观看| 日本爱情动作片www.在线观看 | 热99在线观看视频| 久久人妻av系列| 99热只有精品国产| 免费一级毛片在线播放高清视频| 黄色视频,在线免费观看| 国产伦一二天堂av在线观看| 两个人视频免费观看高清| 成人国产综合亚洲| 国模一区二区三区四区视频| 色哟哟哟哟哟哟| 国产淫片久久久久久久久| 一区二区三区高清视频在线| 亚洲欧美激情综合另类| 国产av在哪里看| 国产精品伦人一区二区| 国产蜜桃级精品一区二区三区| 国产 一区 欧美 日韩| 波多野结衣高清作品| 乱码一卡2卡4卡精品| 波多野结衣高清作品| 麻豆一二三区av精品| 综合色av麻豆| 亚洲,欧美,日韩| 国产精品三级大全| 精品久久久噜噜| 99在线人妻在线中文字幕| 大又大粗又爽又黄少妇毛片口| 亚洲美女黄片视频| 1024手机看黄色片| 男人舔奶头视频| 欧美中文日本在线观看视频| 国产伦精品一区二区三区视频9| 最近最新中文字幕大全电影3| 亚洲国产欧洲综合997久久,| netflix在线观看网站| 搡老妇女老女人老熟妇| 欧美日韩综合久久久久久 | 国产高清不卡午夜福利| 好男人在线观看高清免费视频| 大型黄色视频在线免费观看| 女生性感内裤真人,穿戴方法视频| 成人特级av手机在线观看| 国产综合懂色| 久久精品夜夜夜夜夜久久蜜豆| 99热这里只有是精品在线观看| 亚洲成人久久性| 亚洲欧美清纯卡通| 小说图片视频综合网站| 男女视频在线观看网站免费| 国语自产精品视频在线第100页| 亚洲国产精品成人综合色| 亚洲色图av天堂| 嫩草影院入口| 搞女人的毛片| 中文字幕熟女人妻在线| 联通29元200g的流量卡| 亚洲av.av天堂| 黄色丝袜av网址大全| 在线观看一区二区三区| 欧美极品一区二区三区四区| 午夜a级毛片| 色尼玛亚洲综合影院| 一进一出好大好爽视频| 国产综合懂色| 国产精品一区二区三区四区久久| 亚洲人成网站高清观看| 88av欧美| 啦啦啦韩国在线观看视频| 琪琪午夜伦伦电影理论片6080| 真人做人爱边吃奶动态| 亚洲精品国产成人久久av| 午夜日韩欧美国产| 一级毛片久久久久久久久女| 亚洲欧美日韩东京热| 成人午夜高清在线视频| 亚洲,欧美,日韩| 舔av片在线| 国产亚洲精品综合一区在线观看| 色哟哟·www| 免费一级毛片在线播放高清视频| 欧美一区二区国产精品久久精品| 看黄色毛片网站| 亚洲狠狠婷婷综合久久图片| 韩国av一区二区三区四区| 校园人妻丝袜中文字幕| 欧美不卡视频在线免费观看| 精品久久久久久久末码| av福利片在线观看| 男人的好看免费观看在线视频| 国产一级毛片七仙女欲春2| av女优亚洲男人天堂| 日韩欧美国产一区二区入口| 色视频www国产| 美女大奶头视频| 欧美高清性xxxxhd video| 久久久久久伊人网av| 久久久久久久亚洲中文字幕| 国产精品三级大全| 欧美日韩国产亚洲二区| 99国产精品一区二区蜜桃av| h日本视频在线播放| 亚洲精品久久国产高清桃花| 日本黄色视频三级网站网址| 在现免费观看毛片| 亚洲成人久久性| 国产中年淑女户外野战色| 日日啪夜夜撸| 久久久国产成人精品二区| av国产免费在线观看| 老司机午夜福利在线观看视频| 午夜福利欧美成人| 欧美另类亚洲清纯唯美| 美女大奶头视频| 亚洲av电影不卡..在线观看| 久久亚洲精品不卡| 日本在线视频免费播放| 久久久久久伊人网av| 99久久成人亚洲精品观看| 九九热线精品视视频播放| 国产乱人伦免费视频| 熟女人妻精品中文字幕| 久久精品国产自在天天线| 亚洲三级黄色毛片| 成人午夜高清在线视频| 91久久精品电影网| 蜜桃久久精品国产亚洲av| 国产精品三级大全| 成人国产综合亚洲| 日韩 亚洲 欧美在线| 久久精品夜夜夜夜夜久久蜜豆| 可以在线观看毛片的网站| 最近中文字幕高清免费大全6 | 国产黄a三级三级三级人| 少妇熟女aⅴ在线视频| 99久久精品国产国产毛片| 色在线成人网| 精品一区二区三区av网在线观看| 精品免费久久久久久久清纯| 欧美日韩黄片免| 干丝袜人妻中文字幕| 国产成年人精品一区二区| 一个人看视频在线观看www免费| 国内精品宾馆在线| av在线观看视频网站免费| 国产精品久久电影中文字幕| 成年女人永久免费观看视频| 免费av不卡在线播放| 欧美xxxx性猛交bbbb| 成人国产麻豆网| av专区在线播放| 91久久精品国产一区二区成人| 久久热精品热| 高清在线国产一区| 亚洲性夜色夜夜综合| 熟女电影av网| 国产精品伦人一区二区| 哪里可以看免费的av片| 国产乱人伦免费视频| 国产中年淑女户外野战色| 99热网站在线观看| 欧美+日韩+精品| 国产欧美日韩一区二区精品| 精华霜和精华液先用哪个| 国产av麻豆久久久久久久| 亚洲五月天丁香| av中文乱码字幕在线| 国产亚洲欧美98| 亚洲av美国av| 又黄又爽又免费观看的视频| 成人av在线播放网站| 午夜a级毛片| 欧美高清成人免费视频www| 嫁个100分男人电影在线观看| 精品人妻视频免费看| 联通29元200g的流量卡| 嫩草影院入口| 午夜影院日韩av| 啦啦啦观看免费观看视频高清| 免费人成视频x8x8入口观看| 欧美+亚洲+日韩+国产| 成人亚洲精品av一区二区| 九九热线精品视视频播放| 熟女电影av网| 欧美成人免费av一区二区三区| 琪琪午夜伦伦电影理论片6080| 在线a可以看的网站| 免费看av在线观看网站| 天堂√8在线中文| 国语自产精品视频在线第100页| 免费在线观看成人毛片| 丰满人妻一区二区三区视频av| 欧美激情国产日韩精品一区| 简卡轻食公司| 日日撸夜夜添| 亚洲最大成人中文| 在线免费观看不下载黄p国产 | 一本精品99久久精品77| 一级毛片久久久久久久久女| 国产高潮美女av| 精品久久久久久久久久久久久| 国产伦人伦偷精品视频| 在线观看av片永久免费下载| 国产综合懂色| 亚洲欧美日韩卡通动漫| 婷婷精品国产亚洲av在线| 亚洲在线观看片| 久久亚洲真实| 无遮挡黄片免费观看| 国产探花在线观看一区二区| 一区二区三区免费毛片| 亚洲美女黄片视频| 国产伦精品一区二区三区四那| 高清毛片免费观看视频网站| 国产一区二区在线观看日韩| 小蜜桃在线观看免费完整版高清| 国产亚洲av嫩草精品影院| 在线观看一区二区三区| av在线亚洲专区|