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

    基于DCG心電信號的R波檢測算法

    2014-06-05 09:50:14張英濤黃劍華李明達
    關(guān)鍵詞:波峰心電電信號

    張英濤,黃劍華,李明達,宋 濤

    (1. 哈爾濱工業(yè)大學(xué)計算機科學(xué)與技術(shù)學(xué)院,哈爾濱 150001;2. 哈爾濱醫(yī)科大學(xué)附屬第一醫(yī)院心內(nèi)科,哈爾濱 150001)

    基于DCG心電信號的R波檢測算法

    張英濤1,黃劍華1,李明達1,宋 濤2

    (1. 哈爾濱工業(yè)大學(xué)計算機科學(xué)與技術(shù)學(xué)院,哈爾濱 150001;2. 哈爾濱醫(yī)科大學(xué)附屬第一醫(yī)院心內(nèi)科,哈爾濱 150001)

    提出了一種基于動態(tài)心電圖(DCG)心電信號的R波檢測的新算法.該算法以平均雙向斜率和相對高度為主要特征,能夠快捷準確地檢測R波.通過對MIT-BIH Long-Term ECG數(shù)據(jù)庫以及哈爾濱醫(yī)科大學(xué)醫(yī)院提供的Holter記錄進行R波檢測,驗證了該算法的可行性.針對檢測結(jié)果分別與最大值雙重搜索技術(shù)和差分運算方法(DOM)的檢測結(jié)果進行了比較,發(fā)現(xiàn)所提算法的識別正確率(98.3%)高于以上兩種算法(95.2%和90.7%).

    動態(tài)心電圖心電信號;平均雙向斜率;相對高度;美國麻省理工心律失常數(shù)據(jù)庫

    R波的準確定位對于心臟性猝死的預(yù)警診斷有著積極的臨床意義[1].動態(tài)心電圖(dynamic electrocardiography,DCG)作為一種診斷心律失常的輔助工具,以其無創(chuàng)性而被廣泛應(yīng)用于臨床[2].心率變異性分析[3]、心率減速力分析[4]、竇性心率震蕩分析[5]以及T波電交替分析[6]是幾種基于DCG的有效預(yù)警心律失常的方法,而這幾種分析方法都依賴于準確的R波檢測.

    非實時的R波檢測方法包括小波分析法和神經(jīng)網(wǎng)絡(luò)等方法,雖然檢測效果較好,但運算量大、原理復(fù)雜,不易重現(xiàn);而實時的R波檢測方法,其檢測效率高、速度快、原理簡單、易于重現(xiàn),已經(jīng)成為近幾年的研究熱點.

    文獻[7-8]提出了兩種實時的R波檢測算法——最大值雙重搜索技術(shù)和差分運算方法,兩種算法的共同點是利用心電信號中R波斜率大的特點,以相鄰點之間的斜率為主要特征對R波進行檢測.上述算法對噪聲干擾比較敏感,在檢測R波之前需要對原始信號進行基線漂移校正,對于采樣頻率較低的DCG心電信號的處理效果不夠理想.

    筆者提出了一種基于DCG心電信號的R波檢測新算法,利用了R波波峰斜率大、幅度高的特點,提取了平均雙向斜率和相對高度兩項特征,對MIT-BIH Long-Term ECG數(shù)據(jù)庫中的全部7個樣本實現(xiàn)了R波波峰點坐標的定位.通過與另外兩種實時的R波檢測算法的對比,驗證了本文所提算法的可行性.

    1 R波檢測原理

    一個正常心電圖的每個心動周期由P波、Q波、R波、S波、T波以及U波組成(如圖1所示).其中特征最為顯著的是R波,它充分反映了心室收縮時心臟的電行為[9],因此解決心動周期提取問題的關(guān)鍵就是R波的準確定位.

    觀察發(fā)現(xiàn),QRS波群中的QR波和RS波斜率高、幅度大.利用這兩個特點,對心電信號這兩方面的特征進行提取可以突顯R波波峰的位置,從而達到準確定位R波的目的.

    圖1 正常心電圖的一個心動周期Fig.1 A cardiac cycle of normal DCG

    1.1 斜率步長的確定

    在對心電信號的處理中,定義斜率步長為用于求取斜率的兩采樣點間的采樣點間隔.

    設(shè)()N i為心電圖中采樣點i的縱坐標(i點對應(yīng)的電壓值),f為心電信號的采樣頻率.采樣點i和采樣點j之間的斜率為

    由于采樣頻率f是一個定值,因此相鄰采樣點間的時間間隔相等,于是式(1)可以簡化為

    在式(2)中,采樣點i和j的斜率步長就是ij-.

    因為DCG心電信號的采樣頻率較小,如果斜率步長設(shè)置得不夠合理,容易誤將偽R波(如圖2(a)所示)或畸形T波(如圖2(b)所示)誤識為R波.為了防止這兩種情況的發(fā)生,需要對斜率步長進行合理設(shè)定.

    本文所采用的MIT-BIH Long-Term ECG數(shù)據(jù)庫所有心電信號的采樣頻率均為128,Hz,文獻[10]表明,一般QRS波的時間是0.06,s左右,也就是說QR斜率步長與RS斜率步長的和大約在7個采樣點左右(0.06÷(1÷128)=7.68),為突顯QR斜率和RS斜率高的同時保持R波的對稱性,QR坡的斜率步長和RS坡的斜率步長統(tǒng)一取3.

    圖2 MIT-BIH Long-Term ECG數(shù)據(jù)庫中容易誤識的波形Fig.2Waveforms easy to be misdetected in the MIT-BIH Long-Term ECG database

    1.2 平均雙向斜率

    為突顯出R波兩側(cè)斜率高的特性,本文定義了前向斜率、后向斜率以及平均雙向斜率3個概念.

    設(shè)t為斜率步長.對采樣點i而言,點i與點it-之間的斜率1()k i為點i的前向斜率,即

    而點i與點it+之間的斜率的相反數(shù)2()k i為點i的后向斜率,即

    點i的前向斜率1()k i和后向斜率2()k i的平均值()m i為點i的平均雙向斜率,即

    將式(3)和式(4)代入式(5),得

    由于對每個點i而言t都相等,因此在計算平均雙向斜率時,式(6)可簡化為

    如圖3所示,在平均雙向斜率序列中,R波波峰的位置更加突出,識別起來更加容易.

    圖3 MIT-BIH Long-Term ECG數(shù)據(jù)庫中14046號心電記錄Fig.3 The 14046th waveforms in the MIT-BIH Long-Term ECG database

    1.3 第1個R波波峰的確定

    定位第1個R波波峰是定位全部R波前需要解決的首要問題.只有確定了第1個R波波峰的位置才能獲得R波形態(tài)特征的閾值,進而迭代地把R波檢測進行下去.

    第1個R波波峰的定位依賴于初始平均雙向斜率閾值的確定:首先在信號開頭的至少含有2個R波的信號段中求取最大平均雙向斜率x;然后選擇一個合適的閾值r(如0.6),以xr為平均雙向斜率的閾值重新搜索這一段信號,可以定位第1個最大平均雙向斜率的位置loc.在該點附近(oc5l-~oc5l+)搜索到的最大峰值top(如式(8)所示)就是第1個R波波峰的位置.

    1.4 相對高度

    為突顯R波幅度特性的同時又能克服某些具有基線漂移(如圖4(a)和4(b)所示)的信號在R波檢測中的干擾,本文定義前向相對高度、后向相對高度以及相對高度的概念.

    對采樣點i而言,點i與點it-之間的縱坐標之差1()h i為點i的前向相對高度,即

    點i與點it+之間的縱坐標之差2()h i為點i的后向相對高度,即

    點i的前向相對高度1()h i與后向相對高度2()h i的平均值()H i為點i的相對高度,即

    將式(9)和式(10)代入式(11),得

    圖4 MIT-BIH Long-Term ECG數(shù)據(jù)庫中14046號和15814號心電記錄Fig.4The 14046th and 15814th waveforms in the MITBIH Long-Term ECG database

    如圖4(b)和4(c)所示,相對高度序列不僅突出了R波波峰的位置,而且有效解決了基線漂移的干擾.

    2 R波檢測流程

    本文算法的流程大致可分為3部分:雙向斜率和相對高度的計算;第1個R波波峰的定位;迭代地依次定位R波.圖5給出了該算法的算法流程,圖6給出了基于平均雙向斜率和相對高度的R波檢測算法以及相關(guān)過程的偽代碼.

    圖5 R波檢測程序流程Fig.5 Flow chart of R-wave detection

    2.1 掃描定位

    以當(dāng)前已經(jīng)檢測出的R波波峰a的平均雙向斜率()M a和相對高度()H a為標準,選擇合適的閾值d、rd和ru,在一定的范圍(mind~maxd)內(nèi)搜索,滿足下列條件的采樣點i即為下一個R波波峰.

    (1) 控制平均雙向斜率的條件為

    (2) 控制相對高度的條件為

    (3) 判斷是否是波峰的條件為

    經(jīng)實驗整理,可以參考的閾值為d=0.5,dr=0.5,ur=3,dmin=40,dmax=640.

    2.2 偽R波排除與R波波峰疑似點定位

    尋找下一個R波波峰疑似點是在R波定位中需要解決的另一問題.當(dāng)定位的“R波”是由噪聲等干擾產(chǎn)生的偽R波(如圖7所示)或在搜索范圍內(nèi)未能定位出R波時,就需要識別并排除偽R波,繼而定位后續(xù)信號中的R波波峰疑似點,以便于再次迭代搜索時進行篩選和判斷.

    圖6 R波檢測算法Fig.6 Algorithm of R-wave detection

    本文的算法分為兩種情況進行處理.

    (1) 搜索時出現(xiàn)異常波峰(如圖7所示).

    若H( i)>urH( a),需立即排除該點,并將i+1~i+dmax范圍內(nèi)的平均雙向斜率最大值點s1(如式(16)所示)作為R波波峰疑似點,并以該點為起點開始下一輪迭代.

    (2) 在搜索范圍內(nèi)未能搜索到符合閾值的波峰.

    在a+1~a+dmax范圍內(nèi)未能搜索到滿足式(16)的s1作為R波波峰疑似點,應(yīng)將a+dmax+1~a+2dmax范圍內(nèi)的平均雙向斜率最大值點s2(如式(17)所示)作為R波波峰疑似點,并以該點為起點開始下一輪迭代.

    圖7 含有異常波峰的心電記錄Fig.7 Waveforms with abnormal peak

    在每輪迭代之前需要對該疑似點進行式(13)、式(14)和式(15)3項條件的審核以避免誤識的情況發(fā)生,若未通過審核,則判定該疑似點為偽R波波峰并再次搜索,直到搜索到通過審核的疑似點為止.

    3 實驗結(jié)果分析

    針對MIT-BIH Long-Term ECG數(shù)據(jù)庫中的全部7個心電信號記錄——10446、14134、14149、14157、14172、14184和15814,以每個心電信號第1導(dǎo)的前50個R波波峰坐標的實際位置為標準,運用本文算法以及文獻[7]和文獻[8]算法進行測試,檢測結(jié)果如表1所示.

    正確率AR的計算式為

    式中:TP為心拍數(shù);FP為漏檢數(shù);NP為誤檢數(shù).

    實驗中所采用的MIT-BIH Long-Term ECG數(shù)據(jù)庫的全部7個樣本涵蓋了DCG心電信號中可能出現(xiàn)的各種波形.

    15814號樣本存在基線漂移和畸形T波(如圖2和圖4所示);14134和14172號樣本存在一定的噪聲干擾(如圖8和圖9所示);14157號樣本存在倒置的T波(如圖10所示).

    圖8 MIT-BIH Long-Term ECG數(shù)據(jù)庫中14134號心電記錄Fig.8The 14134th waveforms in MIT-BIH Long-Term ECG database

    圖9 MIT-BIH Long-Term ECG數(shù)據(jù)庫中14172號心電記錄Fig.9The 14172nd waveforms in MIT-BIH Long-Term ECG database

    表1 對MIT-BIH Long-Term ECG數(shù)據(jù)的R波檢測結(jié)果Tab.1 Comparison of R-wave detection in MIT-BIH Long-Term ECG database

    圖10 MIT-BIH Long-Term ECG數(shù)據(jù)庫中14157號心電記錄Fig.10 The 14157th waveforms of the MIT-BIH Long-Term ECG database

    在這些波形的影響下,文獻[7]的誤檢個數(shù)較多,而文獻[8]的漏檢和誤檢個數(shù)都較多(如圖11所示,其中的星狀點為各算法檢測出的R波波峰點,圖11(b)和11(c)中的黑點為漏檢的波峰位置).這主要是由于兩種算法在求取斜率時只考慮了相鄰采樣點間的單向斜率,而本文算法求取的是平均雙向斜率,加上相對高度的協(xié)調(diào)控制,能夠排除各種波形的干擾,達到準確定位的目的.

    在檢測精度方面,將本文算法檢測出的R波位置與心電信號中R波的實際位置相對照,其誤差不大于1個采樣點(其中約99.7%能準確定位),具有相當(dāng)高的定位精度.

    為進一步驗證本文算法在多組心電信號數(shù)據(jù)中檢測的穩(wěn)定性,針對哈爾濱醫(yī)科大學(xué)第一附屬醫(yī)院心內(nèi)科提供的24例心肌梗死后患者的Holter記錄(采樣頻率為125,Hz),同樣以每個心電信號第1導(dǎo)的前50個R波波峰坐標的實際位置為標準,運用本文算法以及文獻[7]和文獻[8]算法進行測試,檢測結(jié)果如表2所示.

    檢測結(jié)果表明在多組Holter數(shù)據(jù)的測試中,本文算法依然可以保持很高的準確率,而另兩種方法則略顯遜色.這說明本算法既具有準確性高的特點,同時具有很強的穩(wěn)定性.

    表2 對實際臨床采集Holter數(shù)據(jù)的R波檢測結(jié)果Tab.2 Comparison of R-wave detection of Holter data of clinical collection

    圖11 檢測結(jié)果對比Fig.11 Comparison of detection results

    4 結(jié) 論

    (1) 利用平均雙向斜率和相對高度雙重特征進行定位的策略,使得該算法對信號噪聲的敏感性較低,可以自動排除因信號噪聲產(chǎn)生的偽R波波峰.因此在檢測R波之前不用處理基線漂移和消噪,較好地保持了原信號波形的細節(jié),提高了心電信號處理的效率.

    (2) 具有很好的穩(wěn)定性,能夠準確識別R波波峰的位置,克服了文獻[7]和文獻[8]中難以避免的漏識和誤識的情況.

    (3) 算法易于重現(xiàn),為心率變異性分析、心率減速力分析提供了準確的數(shù)據(jù)特征.

    [1] 鐘杭美. 動態(tài)心電圖預(yù)警心臟性猝死[J]. 臨床心電學(xué)雜志,2007,16(5):352-355.

    Zhong Hangmei. Dynamic electrocardiogram warning sudden cardiac death[J]. Journal of Clinical Electrocardiology,2007,16(5):352-355(in Chinese).

    [2] 栗 浩,錢春麗,丁彥利. 動態(tài)心電圖臨床應(yīng)用的現(xiàn)狀與進展[J]. 中原醫(yī)刊,2008,35(6):86-87.

    Li Hao,Qian Chunli,Ding Yanli. Status and progress of dynamic electrocardiogram clinical application[J]. Central Plains Medical Journal,2008,35(6):86-87(in Chinese).

    [3] Lombardi F. Chaos theory,heart rate variability,and arrhythmic mortality[J]. Circulation,2000,101:8-10.

    [4] Bauer Axel,Kantelhardt Jan W,Barthel Petra,et al. Deceleration capacity of heart rate as a predictor of mortality after myocardial infarction:Cohort study[J]. Lancet,2006,367(9523):1674-1681.

    [5] Schmidt Georg,Malik Marek,Barthel Petra,et al. Heart-rate turbulence after ventricular premature beats as a predictor of mortality after acute myocardial infarction[J]. Lancet,1999,353(9162):1390-1396.

    [6] Verrier Richard L,Klingenheben Thomas,Malik Marek,et al. Microvolt T-wave alternans:Physiological basis,methods of measurement,and clinical utility—Consensus guideline by international society for Holter and noninvasive electrocardiology[J]. Journal of the American College of Cardiology,2011,58(13):1309-1324.

    [7] 張開滋,劉海祥,吳 杰,等. 臨床心電信息學(xué)[M].長沙:湖南科學(xué)技術(shù)出版社,2002.

    Zhang Kaizi,Liu Haixiang,Wu Jie,et al. Clinical Electrocardiographic Informatics[M]. Changsha:Hunan Science and Technology Press,2002(in Chinese).

    [8] Yeh Y,Wang W. QRS complexes detection for ECG signal:The difference operation method[J]. Comput Meth Prog Bio,2008,91(3):245-254.

    [9] 楊 鵬,田愛英,郭 欣. 基于差分斜率法檢測QRS波算法[J]. 南京理工大學(xué)學(xué)報,2009,33(增):128-132.

    Yang Peng,Tian Aiying,Guo Xin. Detection of QRS wave based on difference -slope method[J]. Journal of Nanjing University of Science and Technology,2009,33(Suppl):128-132(in Chinese).

    [10] Saurabh Pal,Madhuchhanda Mitra. Detection of ECG characteristic points using multiresolution wavelet analysis based selective coefficient method[J]. Measurement,2010,43(2):255-261.

    (責(zé)任編輯:金順愛)

    Novel R-Wave Detection Algorithm of DCG Signal

    Zhang Yingtao1,Huang Jianhua1,Li Mingda1,Song Tao2
    (1. School of Computer Science and Engineering,Harbin Institute of Technology,Harbin 150001,China;2. Department of Cardiology,F(xiàn)irst Affiliated Hospital of Harbin Medical University,Harbin 150001,China)

    A novel R-wave detection algorithm of DCG signal is proposed. The algorithm takes the average two-way slope and relative height as its main characteristics. It can detect the R-wave fast and accurately. The feasibility of the algorithm is verified by MIT-BIH Long-Term ECG database and the Holter records from FAHHMU. Experimental results indicate that this algorithm has a much more precise detection rate(98.3%)than the maximum double-searching technology and the difference operation method(DOM)(95.2% and 90.7%).

    DCG signal;average two-way slope;relative height;MIT-BIH Long-Term ECG database

    TP391.4

    A

    0493-2137(2014)01-0074-07

    10.11784/tdxbz201206014

    2012-06-09;

    2012-11-06.

    國家自然科學(xué)基金資助項目(61100097).

    張英濤(1975— ),女,博士,副教授.

    張英濤,yingtao@hit.edu.cn.

    猜你喜歡
    波峰心電電信號
    基于聯(lián)合聚類分析的單通道腹部心電信號的胎心率提取
    心電向量圖診斷高血壓病左心室異常的臨床應(yīng)用
    作用于直立堤墻與樁柱的波峰高度分析計算
    基于非接觸式電極的心電監(jiān)測系統(tǒng)
    電子制作(2019年19期)2019-11-23 08:41:40
    穿戴式心電:發(fā)展歷程、核心技術(shù)與未來挑戰(zhàn)
    基于Code Composer Studio3.3完成對心電信號的去噪
    科技傳播(2019年24期)2019-06-15 09:29:28
    更正啟事
    基于隨機森林的航天器電信號多分類識別方法
    兒童標準12導(dǎo)聯(lián)T波峰末間期的分析
    Dynamic Loads and Wake Prediction for Large Wind Turbines Based on Free Wake Method
    国产成人欧美| 久久精品成人免费网站| 首页视频小说图片口味搜索| xxxhd国产人妻xxx| 建设人人有责人人尽责人人享有的| 天天躁狠狠躁夜夜躁狠狠躁| 国产有黄有色有爽视频| 亚洲免费av在线视频| 免费一级毛片在线播放高清视频 | 捣出白浆h1v1| 免费人妻精品一区二区三区视频| 国产aⅴ精品一区二区三区波| 亚洲专区字幕在线| 午夜免费成人在线视频| videosex国产| 国产欧美日韩综合在线一区二区| 大片免费播放器 马上看| 国产无遮挡羞羞视频在线观看| 人妻 亚洲 视频| 免费女性裸体啪啪无遮挡网站| 国产一区二区在线观看av| av国产精品久久久久影院| 欧美日韩黄片免| 久久久欧美国产精品| 欧美精品一区二区大全| 亚洲国产中文字幕在线视频| 欧美日韩av久久| 美女午夜性视频免费| 香蕉丝袜av| av片东京热男人的天堂| 欧美国产精品va在线观看不卡| 国产男靠女视频免费网站| 午夜福利免费观看在线| 在线av久久热| www日本在线高清视频| 免费黄频网站在线观看国产| 99精国产麻豆久久婷婷| 欧美黄色淫秽网站| 一本综合久久免费| 巨乳人妻的诱惑在线观看| 久久人妻福利社区极品人妻图片| 一本色道久久久久久精品综合| 欧美日韩国产mv在线观看视频| 韩国精品一区二区三区| 亚洲性夜色夜夜综合| 一区二区av电影网| 久久精品国产综合久久久| 丝瓜视频免费看黄片| 国产精品九九99| 久久香蕉激情| 午夜福利影视在线免费观看| 亚洲美女黄片视频| 国产福利在线免费观看视频| 美女扒开内裤让男人捅视频| 欧美精品亚洲一区二区| 18禁国产床啪视频网站| 99国产综合亚洲精品| 99香蕉大伊视频| 日日摸夜夜添夜夜添小说| 国产欧美日韩一区二区三区在线| 久久久水蜜桃国产精品网| 成人黄色视频免费在线看| 黄色a级毛片大全视频| 国产无遮挡羞羞视频在线观看| 91精品三级在线观看| 如日韩欧美国产精品一区二区三区| 国产日韩一区二区三区精品不卡| 亚洲精品中文字幕在线视频| 成人亚洲精品一区在线观看| 欧美一级毛片孕妇| 欧美成人午夜精品| 19禁男女啪啪无遮挡网站| 精品国产乱码久久久久久小说| 久久青草综合色| 狠狠精品人妻久久久久久综合| 国产在线免费精品| 热99re8久久精品国产| 美女福利国产在线| 丁香六月天网| 一级,二级,三级黄色视频| av网站在线播放免费| 欧美激情久久久久久爽电影 | 色视频在线一区二区三区| 999久久久精品免费观看国产| 欧美黑人欧美精品刺激| 性少妇av在线| 黄色视频在线播放观看不卡| 亚洲熟妇熟女久久| 国产97色在线日韩免费| 一级a爱视频在线免费观看| 他把我摸到了高潮在线观看 | 99国产极品粉嫩在线观看| 日韩视频一区二区在线观看| 成年女人毛片免费观看观看9 | 美女视频免费永久观看网站| 俄罗斯特黄特色一大片| 午夜久久久在线观看| 国产亚洲欧美在线一区二区| 亚洲第一欧美日韩一区二区三区 | 俄罗斯特黄特色一大片| 丰满人妻熟妇乱又伦精品不卡| 在线观看66精品国产| 欧美激情 高清一区二区三区| 欧美精品一区二区免费开放| 亚洲国产欧美网| 免费黄频网站在线观看国产| 一区二区三区精品91| 高潮久久久久久久久久久不卡| 老汉色∧v一级毛片| 亚洲精品粉嫩美女一区| 精品欧美一区二区三区在线| 天堂俺去俺来也www色官网| 捣出白浆h1v1| 夜夜夜夜夜久久久久| 免费日韩欧美在线观看| 久久精品国产亚洲av香蕉五月 | 久久热在线av| 亚洲伊人色综图| 中文字幕另类日韩欧美亚洲嫩草| 成人影院久久| 亚洲精品在线观看二区| 交换朋友夫妻互换小说| 日本av手机在线免费观看| 好男人电影高清在线观看| 91av网站免费观看| 亚洲情色 制服丝袜| 每晚都被弄得嗷嗷叫到高潮| 性色av乱码一区二区三区2| 亚洲精品国产一区二区精华液| 精品少妇一区二区三区视频日本电影| 五月开心婷婷网| e午夜精品久久久久久久| 欧美激情久久久久久爽电影 | 久久精品国产综合久久久| 巨乳人妻的诱惑在线观看| 免费看十八禁软件| 王馨瑶露胸无遮挡在线观看| 精品一区二区三卡| 久久性视频一级片| 高清欧美精品videossex| 久久精品aⅴ一区二区三区四区| 十八禁人妻一区二区| 啦啦啦中文免费视频观看日本| 国产精品影院久久| av网站免费在线观看视频| 久久热在线av| 日韩欧美国产一区二区入口| 1024视频免费在线观看| 看免费av毛片| 免费少妇av软件| 又大又爽又粗| 咕卡用的链子| 国产精品国产av在线观看| 99精品久久久久人妻精品| 国产在视频线精品| 宅男免费午夜| 成人国产av品久久久| 狠狠狠狠99中文字幕| 正在播放国产对白刺激| 极品少妇高潮喷水抽搐| 日韩人妻精品一区2区三区| 怎么达到女性高潮| 成年动漫av网址| 欧美成狂野欧美在线观看| 99九九在线精品视频| h视频一区二区三区| 如日韩欧美国产精品一区二区三区| 亚洲av成人一区二区三| a级片在线免费高清观看视频| 日本av免费视频播放| 王馨瑶露胸无遮挡在线观看| 久久av网站| 久久这里只有精品19| 国产精品国产高清国产av | 欧美人与性动交α欧美精品济南到| 欧美乱码精品一区二区三区| 一二三四社区在线视频社区8| 菩萨蛮人人尽说江南好唐韦庄| 女性被躁到高潮视频| 亚洲国产看品久久| 国产精品av久久久久免费| 极品人妻少妇av视频| 中文字幕精品免费在线观看视频| 欧美老熟妇乱子伦牲交| 少妇的丰满在线观看| 日韩一区二区三区影片| 日韩欧美免费精品| 国产精品电影一区二区三区 | 99国产精品一区二区蜜桃av | 国产欧美日韩一区二区三区在线| 免费看十八禁软件| 99国产精品一区二区三区| 国产男女内射视频| 欧美+亚洲+日韩+国产| 人成视频在线观看免费观看| 国产成人精品久久二区二区免费| 亚洲男人天堂网一区| av片东京热男人的天堂| 一区二区三区激情视频| 大型av网站在线播放| 国产99久久九九免费精品| 亚洲国产毛片av蜜桃av| 18禁黄网站禁片午夜丰满| 亚洲午夜精品一区,二区,三区| 欧美黑人欧美精品刺激| √禁漫天堂资源中文www| 午夜视频精品福利| 久久精品成人免费网站| 免费日韩欧美在线观看| 亚洲人成77777在线视频| 国产精品 欧美亚洲| 成人国产一区最新在线观看| 久久久久久久久久久久大奶| 成人影院久久| 精品一区二区三区视频在线观看免费 | 国精品久久久久久国模美| 亚洲欧美日韩另类电影网站| 国产成人欧美| 亚洲精品国产精品久久久不卡| 久久午夜亚洲精品久久| 黑人猛操日本美女一级片| 国产精品一区二区精品视频观看| 一个人免费看片子| 久久国产精品大桥未久av| 日本一区二区免费在线视频| 美女福利国产在线| 12—13女人毛片做爰片一| 国产不卡av网站在线观看| 中国美女看黄片| 国产精品一区二区在线不卡| 国产成+人综合+亚洲专区| 日韩制服丝袜自拍偷拍| 亚洲精品粉嫩美女一区| 一区二区av电影网| 中文字幕精品免费在线观看视频| 精品亚洲成a人片在线观看| 久久精品成人免费网站| a级毛片黄视频| 视频区欧美日本亚洲| 国产片内射在线| 极品少妇高潮喷水抽搐| 一进一出抽搐动态| 亚洲精品国产区一区二| 视频在线观看一区二区三区| 国产精品国产高清国产av | 18禁黄网站禁片午夜丰满| 一二三四在线观看免费中文在| 国产av精品麻豆| 精品第一国产精品| 久久毛片免费看一区二区三区| av片东京热男人的天堂| 久久人人97超碰香蕉20202| 亚洲av国产av综合av卡| 美国免费a级毛片| 亚洲成av片中文字幕在线观看| 精品欧美一区二区三区在线| 嫩草影视91久久| 桃红色精品国产亚洲av| 757午夜福利合集在线观看| 成人免费观看视频高清| 免费在线观看日本一区| 宅男免费午夜| 在线播放国产精品三级| 久久狼人影院| 亚洲成a人片在线一区二区| 婷婷成人精品国产| 91国产中文字幕| 一本综合久久免费| 天天添夜夜摸| 少妇猛男粗大的猛烈进出视频| 亚洲精品久久成人aⅴ小说| 操美女的视频在线观看| 日韩制服丝袜自拍偷拍| 在线观看免费午夜福利视频| 999精品在线视频| 亚洲国产欧美一区二区综合| 欧美性长视频在线观看| 俄罗斯特黄特色一大片| 亚洲午夜理论影院| videosex国产| 真人做人爱边吃奶动态| 国产在线免费精品| 窝窝影院91人妻| 亚洲第一青青草原| 制服人妻中文乱码| 精品一品国产午夜福利视频| 亚洲综合色网址| 亚洲少妇的诱惑av| 色婷婷av一区二区三区视频| 国产在线一区二区三区精| 青草久久国产| 午夜两性在线视频| a在线观看视频网站| 美女国产高潮福利片在线看| 性色av乱码一区二区三区2| 精品少妇内射三级| 视频区欧美日本亚洲| 午夜免费鲁丝| 亚洲精品美女久久久久99蜜臀| 成人三级做爰电影| 国产不卡一卡二| 日韩大片免费观看网站| 精品亚洲成a人片在线观看| 精品国产乱码久久久久久小说| 日本黄色视频三级网站网址 | 亚洲avbb在线观看| 国产又爽黄色视频| 成人三级做爰电影| av在线播放免费不卡| 日本vs欧美在线观看视频| 国产欧美日韩一区二区三区在线| 欧美日韩福利视频一区二区| 悠悠久久av| 大香蕉久久成人网| 成人特级黄色片久久久久久久 | 美国免费a级毛片| 亚洲av日韩精品久久久久久密| 久久99热这里只频精品6学生| 国产精品国产av在线观看| 美女扒开内裤让男人捅视频| 成年动漫av网址| 一本色道久久久久久精品综合| 精品亚洲乱码少妇综合久久| 操美女的视频在线观看| 精品亚洲成国产av| 女人精品久久久久毛片| 99国产极品粉嫩在线观看| 国产欧美日韩综合在线一区二区| 天天躁夜夜躁狠狠躁躁| 热99re8久久精品国产| 久久久久久人人人人人| 一进一出抽搐动态| 久久精品国产a三级三级三级| 69精品国产乱码久久久| 999久久久国产精品视频| 国产欧美日韩一区二区精品| 在线 av 中文字幕| 黄色视频不卡| 国产一区二区 视频在线| 丝袜喷水一区| 两性夫妻黄色片| 日日摸夜夜添夜夜添小说| 精品一品国产午夜福利视频| 无人区码免费观看不卡 | 99久久精品国产亚洲精品| 自线自在国产av| 欧美亚洲 丝袜 人妻 在线| 久热爱精品视频在线9| 久久久久久久久久久久大奶| 亚洲专区国产一区二区| 丝袜人妻中文字幕| 热re99久久精品国产66热6| 久久久久国内视频| 亚洲精品在线美女| 美女主播在线视频| 999久久久国产精品视频| 亚洲人成电影免费在线| 成人手机av| 97人妻天天添夜夜摸| 国产福利在线免费观看视频| 午夜激情av网站| 久久久久久久精品吃奶| 亚洲专区国产一区二区| 亚洲五月婷婷丁香| 每晚都被弄得嗷嗷叫到高潮| 久久久国产成人免费| 欧美 日韩 精品 国产| 男人操女人黄网站| 亚洲专区中文字幕在线| 极品人妻少妇av视频| 99国产极品粉嫩在线观看| 极品人妻少妇av视频| 亚洲av成人不卡在线观看播放网| 亚洲免费av在线视频| 在线观看人妻少妇| 老司机午夜福利在线观看视频 | 亚洲伊人久久精品综合| 新久久久久国产一级毛片| 国产区一区二久久| 亚洲全国av大片| 一边摸一边抽搐一进一小说 | 欧美黑人欧美精品刺激| 国产精品电影一区二区三区 | 亚洲国产av新网站| 少妇粗大呻吟视频| 99久久国产精品久久久| 999久久久国产精品视频| 久久精品亚洲av国产电影网| 久久精品亚洲精品国产色婷小说| 另类亚洲欧美激情| 一本色道久久久久久精品综合| 女人高潮潮喷娇喘18禁视频| 久久精品国产99精品国产亚洲性色 | 国产一区二区三区在线臀色熟女 | 啦啦啦在线免费观看视频4| 岛国毛片在线播放| netflix在线观看网站| 久久中文字幕人妻熟女| 一边摸一边抽搐一进一出视频| 夜夜爽天天搞| 黑人巨大精品欧美一区二区mp4| 超色免费av| 热re99久久精品国产66热6| 麻豆国产av国片精品| 纯流量卡能插随身wifi吗| 99热国产这里只有精品6| 丝瓜视频免费看黄片| av片东京热男人的天堂| 国产不卡一卡二| 色尼玛亚洲综合影院| 亚洲成人国产一区在线观看| 亚洲色图综合在线观看| 国产免费现黄频在线看| 久久精品国产综合久久久| 免费女性裸体啪啪无遮挡网站| www日本在线高清视频| 电影成人av| 又紧又爽又黄一区二区| 欧美成狂野欧美在线观看| 中文字幕制服av| 久久久久国内视频| 我要看黄色一级片免费的| 女人久久www免费人成看片| 欧美大码av| 午夜福利欧美成人| 国产极品粉嫩免费观看在线| 又紧又爽又黄一区二区| 欧美激情高清一区二区三区| 精品少妇久久久久久888优播| 日韩精品免费视频一区二区三区| 精品一区二区三区四区五区乱码| 妹子高潮喷水视频| a级毛片在线看网站| 欧美 亚洲 国产 日韩一| 五月开心婷婷网| 国产在线精品亚洲第一网站| 日韩精品免费视频一区二区三区| 亚洲情色 制服丝袜| 久久国产精品男人的天堂亚洲| 菩萨蛮人人尽说江南好唐韦庄| 天天躁狠狠躁夜夜躁狠狠躁| av一本久久久久| 亚洲精品美女久久久久99蜜臀| 国产成人免费无遮挡视频| 老司机影院毛片| 久久精品亚洲熟妇少妇任你| 亚洲精品自拍成人| 亚洲国产成人一精品久久久| 夜夜夜夜夜久久久久| 国产精品香港三级国产av潘金莲| 久久热在线av| 亚洲人成77777在线视频| 亚洲va日本ⅴa欧美va伊人久久| 日韩人妻精品一区2区三区| 热re99久久精品国产66热6| 国产国语露脸激情在线看| 色在线成人网| 国产精品 欧美亚洲| 中文字幕人妻丝袜一区二区| 天天躁日日躁夜夜躁夜夜| 大型av网站在线播放| 一边摸一边做爽爽视频免费| 国产国语露脸激情在线看| 国产一区有黄有色的免费视频| 一级毛片电影观看| 国产精品九九99| 亚洲自偷自拍图片 自拍| 男女边摸边吃奶| 悠悠久久av| 久久久久久亚洲精品国产蜜桃av| 久久国产精品男人的天堂亚洲| 搡老乐熟女国产| av有码第一页| 一区二区三区乱码不卡18| 欧美在线一区亚洲| 美女视频免费永久观看网站| 一本—道久久a久久精品蜜桃钙片| 一区二区日韩欧美中文字幕| 黄色怎么调成土黄色| 五月开心婷婷网| 热99国产精品久久久久久7| 国产成人影院久久av| 国产欧美日韩一区二区三区在线| 欧美一级毛片孕妇| 男女床上黄色一级片免费看| 午夜福利影视在线免费观看| 一边摸一边抽搐一进一小说 | 久久久久精品人妻al黑| 精品久久久精品久久久| 中亚洲国语对白在线视频| 国产成人欧美| 亚洲第一欧美日韩一区二区三区 | 亚洲国产中文字幕在线视频| 99在线人妻在线中文字幕 | 波多野结衣av一区二区av| 国产高清videossex| 亚洲精品国产精品久久久不卡| 99国产精品一区二区三区| 嫁个100分男人电影在线观看| 日韩制服丝袜自拍偷拍| 夫妻午夜视频| 亚洲av日韩精品久久久久久密| 久久久久网色| 别揉我奶头~嗯~啊~动态视频| 在线观看人妻少妇| 两个人看的免费小视频| 欧美日韩亚洲综合一区二区三区_| 97人妻天天添夜夜摸| 美女主播在线视频| av在线播放免费不卡| 久久99热这里只频精品6学生| 狠狠婷婷综合久久久久久88av| 国产亚洲av高清不卡| 国产麻豆69| 午夜福利一区二区在线看| 精品国内亚洲2022精品成人 | 免费人妻精品一区二区三区视频| 极品少妇高潮喷水抽搐| 91大片在线观看| 成在线人永久免费视频| 国产精品1区2区在线观看. | 成人三级做爰电影| 老汉色av国产亚洲站长工具| 五月天丁香电影| 免费观看a级毛片全部| 搡老岳熟女国产| 黄色怎么调成土黄色| 18禁美女被吸乳视频| 国产精品久久久久久精品古装| 首页视频小说图片口味搜索| 美女高潮喷水抽搐中文字幕| 超碰成人久久| 69精品国产乱码久久久| avwww免费| 丁香六月欧美| a级毛片在线看网站| 久久人人97超碰香蕉20202| 国产亚洲午夜精品一区二区久久| 国产精品九九99| 欧美日韩精品网址| 久久久精品国产亚洲av高清涩受| 久久久欧美国产精品| 欧美日韩国产mv在线观看视频| 十八禁人妻一区二区| videos熟女内射| 黑人欧美特级aaaaaa片| 香蕉丝袜av| 欧美日韩黄片免| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲熟妇熟女久久| 水蜜桃什么品种好| 岛国在线观看网站| 69精品国产乱码久久久| 成年动漫av网址| 成人精品一区二区免费| 久久精品成人免费网站| 亚洲成人免费av在线播放| 精品第一国产精品| 亚洲av第一区精品v没综合| 悠悠久久av| 国产成人啪精品午夜网站| 91成年电影在线观看| 麻豆成人av在线观看| 国产97色在线日韩免费| www.999成人在线观看| 免费日韩欧美在线观看| 日本一区二区免费在线视频| 少妇被粗大的猛进出69影院| 搡老熟女国产l中国老女人| 欧美精品人与动牲交sv欧美| 亚洲av成人不卡在线观看播放网| 成人18禁高潮啪啪吃奶动态图| 男女之事视频高清在线观看| 老熟妇仑乱视频hdxx| www.精华液| 成在线人永久免费视频| 超碰97精品在线观看| 亚洲国产欧美在线一区| 日韩熟女老妇一区二区性免费视频| 亚洲欧美精品综合一区二区三区| 日本vs欧美在线观看视频| 欧美黑人欧美精品刺激| 色老头精品视频在线观看| 午夜老司机福利片| 黄片小视频在线播放| 757午夜福利合集在线观看| 国产国语露脸激情在线看| 国产一卡二卡三卡精品| 国产精品久久久久成人av| 欧美日韩福利视频一区二区| 久久久欧美国产精品| 国产精品久久久久成人av| 在线观看免费午夜福利视频| 亚洲视频免费观看视频| 超碰成人久久| 无人区码免费观看不卡 | 啦啦啦视频在线资源免费观看| 久久精品91无色码中文字幕| 波多野结衣av一区二区av| 啦啦啦在线免费观看视频4| 1024视频免费在线观看| 欧美精品高潮呻吟av久久| 国产1区2区3区精品| 高清毛片免费观看视频网站 | 欧美日韩精品网址| 母亲3免费完整高清在线观看| 九色亚洲精品在线播放| 欧美日韩一级在线毛片| 日本欧美视频一区| 91大片在线观看| 国产有黄有色有爽视频| 纯流量卡能插随身wifi吗| 成人18禁高潮啪啪吃奶动态图| 亚洲性夜色夜夜综合| 欧美 日韩 精品 国产|