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

    一種基于血氧飽和度與鼾聲檢測的呼吸暫停監(jiān)測系統(tǒng)設(shè)計(jì)

    2021-03-17 09:27:24謝運(yùn)泉鄔楊波謝建軍
    關(guān)鍵詞:鼾聲脈搏血氧

    謝運(yùn)泉, 李 宏, 鄔楊波, 謝建軍

    一種基于血氧飽和度與鼾聲檢測的呼吸暫停監(jiān)測系統(tǒng)設(shè)計(jì)

    謝運(yùn)泉, 李 宏, 鄔楊波, 謝建軍*

    (寧波大學(xué) 信息科學(xué)與工程學(xué)院, 浙江 寧波 315211)

    長期打鼾嚴(yán)重的人往往都患有呼吸暫停綜合征, 該癥狀對人的身體健康有著極大的危害. 針對用鼾聲實(shí)現(xiàn)呼吸暫停檢測技術(shù)準(zhǔn)確度不夠的問題, 設(shè)計(jì)了一種基于血氧飽和度與鼾聲檢測相結(jié)合的呼吸暫停監(jiān)測系統(tǒng), 利用微控制器結(jié)合外部的血氧檢測電路以及鼾聲檢測電路實(shí)現(xiàn)血氧飽和度的實(shí)時(shí)監(jiān)測與鼾聲片段的識別, 當(dāng)監(jiān)測到用戶睡眠時(shí)正在打鼾同時(shí)血氧飽和度在一定程度上下降, 系統(tǒng)便可以識別出該打鼾者的呼吸暫停癥狀. 而對于患有呼吸暫停癥狀的打鼾者, 后續(xù)可以對其進(jìn)行止鼾, 從而一定程度上可以緩解呼吸暫停癥狀. 整個(gè)系統(tǒng)在準(zhǔn)確實(shí)時(shí)的血氧飽和度監(jiān)測以及鼾聲識別基礎(chǔ)上, 可以較為方便準(zhǔn)確地識別出打鼾者的呼吸暫停癥狀, 為呼吸暫停綜合征的實(shí)時(shí)監(jiān)測提供了一種新的方法.

    呼吸暫停綜合征; 血氧飽和度; 鼾聲檢測; 止鼾

    睡眠呼吸暫停綜合征(Sleep Apnea Syndrome, SAS)是一種常見的睡眠障礙性疾病, 嚴(yán)重影響人們的生活質(zhì)量以及身體健康. SAS最常見的臨床表現(xiàn)便是打鼾, 相比一般的打鼾, SAS的特有癥狀是呼吸暫停和低通氣事件頻繁發(fā)生. 所謂呼吸暫停是指夜間睡眠時(shí)呼吸氣流完全停止10s以上, 而低通氣是指睡眠期間的呼吸氣流相對一般基礎(chǔ)水平減少50%以上, 并伴有血氧飽和度減少4%以上[1]. 血氧飽和度下降會產(chǎn)生更加嚴(yán)重的危害, 如高血壓、低氧血癥、冠心病等, 部分心血管疾病患者甚至?xí)霈F(xiàn)夜間猝死[2-3]. 根據(jù)相關(guān)調(diào)查結(jié)果, 全球每天有3000多人發(fā)生與SAS相關(guān)的夜間死亡[4], 且在中老年肥胖患者中更易出現(xiàn)[5]. 因此, 對于打鼾者所患的SAS有必要進(jìn)行監(jiān)測和治療.

    多導(dǎo)睡眠儀(Polysomnography, PSG)是呼吸暫停癥狀的主要檢測設(shè)備, 但其操作復(fù)雜, 價(jià)格昂貴, 給患者帶來的負(fù)擔(dān)較重, 不易在廣大人群中推廣[6]. 通過動態(tài)腦電監(jiān)測來進(jìn)行SAS檢測的方法應(yīng)運(yùn)而生, 但該方法在檢測實(shí)施過程中需要專業(yè)醫(yī)護(hù)人員進(jìn)行實(shí)時(shí)操作, 不易在家庭中推廣[7]. 通過血氧飽和度來獨(dú)立判斷SAS的發(fā)生可以作為一種方便有效的SAS初篩手段, 但單一的血氧飽和度監(jiān)測會低估患病的嚴(yán)重程度[8]. 目前, 基于動態(tài)心電圖機(jī)監(jiān)測的心率變異性(Heart Rate Variability, HRV)檢測法是替代PSG檢測的主要方法之一, 但該檢測方法需要病人長期佩戴心電電極, 容易造成皮膚過敏[6]. 通過直接分析用戶鼾聲信號來進(jìn)行SAS檢測與篩選是一種簡單且常見的方法[9], 但該方法不易準(zhǔn)確區(qū)分單純打鼾者和患有呼吸暫停癥狀的打鼾者, 存在一定的誤差, 準(zhǔn)確度有待提高.

    針對已有SAS檢測方法準(zhǔn)確度低的問題, 結(jié)合血氧飽和度與呼吸暫停癥狀的密切關(guān)系[10], 本文提出了一種通過鼾聲檢測區(qū)分鼾聲與非鼾聲, 再通過監(jiān)測其血氧飽和度的變化來檢測呼吸暫停癥狀的方法, 該方法可以有效地識別出打鼾者的呼吸暫停癥狀, 明顯提高檢測SAS的準(zhǔn)確度.

    1 血氧與鼾聲檢測原理

    1.1 血氧檢測原理

    血氧飽和度(SpO2)用于表示血液中的血氧含量, 定義為氧合血紅蛋白(HbO2)容量占全部血紅蛋白(HbO2+Hb)容量的百分比:

    式中:為血氧飽和度,為血紅蛋白容量.

    氧合血紅蛋白與還原血紅蛋白對波長為940和660nm的入射光吸收系數(shù)差異最大[11-13], 這兩種光可以作為血氧飽和度測量的光源. 本文采用透射式方法, 如圖1所示, 發(fā)射端是波長1=660 nm的紅光(Red)和波長2=940nm的紅外光(IRed)發(fā)光二極管, 接收端是光電接收管, 光透過手指被接收管接收轉(zhuǎn)換為電信號輸出. 根據(jù)朗伯-比爾定律推導(dǎo)出血氧飽和度的標(biāo)定公式:

    式中:、為常數(shù);AC、DC分別為脈搏信號的交流分量、直流分量. 為了提高測量精度, 對式(2)取泰勒級數(shù)展開的前兩項(xiàng), 得出實(shí)際設(shè)計(jì)應(yīng)用中的血氧飽和度標(biāo)定公式:

    式中,為常數(shù). 由于動脈血液會隨著心臟的搏動而產(chǎn)生節(jié)律性變化, 使得光透過手指末端時(shí)光強(qiáng)改變, 而其他組織相對靜止, 最后只需要獲取光經(jīng)過人體組織后光強(qiáng)的變化量(AC)與最大值(DC)的比值, 便可以計(jì)算出血氧飽和度.

    圖1 透射式血氧檢測原理

    1.2 鼾聲檢測原理

    鼾聲的發(fā)聲機(jī)制與語音信號類似, 都是在聲道中形成的, 不過鼾聲是由咽部振動而不是聲帶振動產(chǎn)生, 且產(chǎn)生在吸氣階段. 可以使用語音信號的常用特征來作為研究鼾聲信號的依據(jù), 如梅爾頻率倒譜系數(shù)(Mel Frequency Cepstral Coefficents, MFCC)、共振峰等.

    目前, 對于語音信號的研究已較為成熟, 但對于鼾聲的識別研究相對較少. 對于鼾聲的檢測與識別, 主要分為監(jiān)督與非監(jiān)督兩種方法, 兩者的區(qū)別在于是否需要訓(xùn)練數(shù)據(jù). 本文采用監(jiān)督學(xué)習(xí)的方法, 利用訓(xùn)練數(shù)據(jù)來識別出鼾聲. 首先通過麥克風(fēng)檢測聲音信號并采集到微控制器中, 然后進(jìn)行信號處理, 主要包括預(yù)處理、端點(diǎn)檢測、特征提取以及鼾聲識別. 預(yù)處理過程主要是對采集的語音數(shù)據(jù)進(jìn)行預(yù)加重、分幀、加窗, 是為了消除人類發(fā)聲機(jī)制等因素給語音信號分析帶來的影響. 端點(diǎn)檢測主要是檢測出有聲段, 去除睡眠時(shí)的寂靜段. 特征提取過程主要是提取與人耳感知特性相關(guān)的MFCC. 最后采用機(jī)器學(xué)習(xí)算法K近鄰(K-Nearest Neighbor, KNN)來識別出鼾聲.

    2 基于血氧飽和度與鼾聲檢測的呼吸暫停監(jiān)測與識別

    打鼾嚴(yán)重的人通常都患有睡眠呼吸暫停綜合征. 在基于聲學(xué)的分析中, 鼾聲對區(qū)分單純打鼾以及呼吸暫停癥狀具有很大的潛力. 分析整晚的鼾聲信號特征, 可以初步判定是否患有呼吸暫停癥狀. 但研究發(fā)現(xiàn), 呼吸暫停患者的正常打鼾時(shí)間間隔不定且振幅不穩(wěn)定[1], 通過分析鼾聲信號特征并不能準(zhǔn)確地檢測出呼吸暫停癥狀, 但患有呼吸暫停癥狀的打鼾者睡眠時(shí)血氧飽和度會一定程度上降低, 可以對睡眠時(shí)的打鼾者監(jiān)測其血氧飽和度的下降程度來判斷是否患有呼吸暫停癥狀. 然而, 對于長期患有心血管疾病的患者, 在睡眠時(shí)也會因?yàn)楹粑豁槙硨?dǎo)致血氧飽和度下降, 單純監(jiān)測血氧飽和度并不能準(zhǔn)確地檢測出呼吸暫停癥狀. 本文的呼吸暫停監(jiān)測系統(tǒng)主要是監(jiān)測睡眠時(shí)打鼾者血氧飽和度的下降程度, 識別出打鼾者的呼吸暫停癥狀, 無需從鼾聲信號特征上進(jìn)行分析. 因?yàn)槌霈F(xiàn)呼吸暫停的最基礎(chǔ)癥狀便是打鼾, 系統(tǒng)首先需要檢測出用戶睡眠時(shí)是否在打鼾. 對于打鼾者, 監(jiān)測其血氧飽和度的變化, 如果血氧飽和度明顯下降, 則識別出用戶是一名呼吸暫停癥狀患者, 后續(xù)可以對用戶進(jìn)行止鼾, 一定程度上緩解用戶的呼吸暫停癥狀.

    3 系統(tǒng)硬件設(shè)計(jì)

    整個(gè)系統(tǒng)主要由血氧檢測電路、鼾聲采集電路及處理器外圍相關(guān)電路3部分組成. 3個(gè)部分協(xié)同工作, 當(dāng)系統(tǒng)檢測到睡眠用戶打鼾的同時(shí)血氧飽和度在一定范圍內(nèi)下降, 便識別該打鼾者是呼吸暫停癥狀患者. 系統(tǒng)硬件框圖如圖2所示.

    圖2 系統(tǒng)硬件框圖

    3.1 血氧檢測電路設(shè)計(jì)

    血氧檢測電路以STM32F407單片機(jī)為核心, 由指夾式脈搏血氧探頭、光源驅(qū)動電路以及信號調(diào)理電路構(gòu)成. 血氧探頭采集人體手指末端的光電容積脈搏波(Photoplethysmographic, PPG), 由于探頭輸出的原始脈搏波信號微弱, 且干擾較多, 對其進(jìn)行前置放大、信號分離以及濾波等信號調(diào)理電路的處理, 再經(jīng)過模數(shù)轉(zhuǎn)換, 由單片機(jī)采集原始的脈搏波信號, 處理器內(nèi)部對PPG信號處理, 顯示實(shí)時(shí)的血氧飽和度值.

    3.1.1光源驅(qū)動電路設(shè)計(jì)

    光源驅(qū)動電路主要是驅(qū)動血氧探頭正常工作, 如圖3所示, 為了使血氧探頭內(nèi)部的紅光與紅外光發(fā)光二極管交替點(diǎn)亮[14], 利用H橋電路, 通過四路PWM波對信號進(jìn)行分時(shí)驅(qū)動.

    圖3 光源驅(qū)動電路

    工作時(shí)紅光與紅外光交替點(diǎn)亮. Q1、Q4導(dǎo)通, Q2、Q3截止, 紅光發(fā)光二極管點(diǎn)亮, 流過的電流為15mA, 滿足驅(qū)動紅光二極管的要求; Q2、Q3導(dǎo)通, Q1、Q4截止, 紅外發(fā)光二極管點(diǎn)亮, 流過的電流為17mA, 滿足紅外發(fā)光二極管工作的要求. 電路通過分時(shí)交替發(fā)光來驅(qū)動血氧探頭正常工作.

    3.1.2前置放大及信號分離電路設(shè)計(jì)

    由血氧探頭產(chǎn)生的原始脈搏波信號微弱, 且極易受到干擾, 為了提升電路的抗干擾能力, 需要加入前置放大電路, 而前置放大電路的性能會對電路抗干擾能力產(chǎn)生重大影響, 因此在選擇前置放大器時(shí)要求運(yùn)放具有高共模抑制比、低噪聲電壓、高信噪比等特點(diǎn)[15]. 經(jīng)對比, 本文采用INA128儀表運(yùn)算放大器. P1和P2是血氧探頭正常工作時(shí)的輸出端, 與前置放大電路相連, 經(jīng)前置放大電路后輸出的是紅光與紅外光混合信號, 需要將兩路信號進(jìn)行分離以便后續(xù)處理. 本文采用模擬電子開關(guān)CD4066芯片來實(shí)現(xiàn), 其內(nèi)部由4個(gè)雙向電子開關(guān)組成. 需要注意的是, 為了將兩路波長信號成功分離, 該芯片的分離控制信號必須與光源驅(qū)動電路的導(dǎo)通時(shí)序保持一致. 具體電路如圖4所示.

    圖4 前置放大及信號分離電路

    在INA128的1腳與8腳之間外接一個(gè)電阻7, 可以對前置放大電路的增益進(jìn)行設(shè)置, 公式如下,

    3.1.3濾波電路設(shè)計(jì)

    濾波電路主要是為了得到計(jì)算血氧飽和度的交流分量以及直流分量, 兩路信號各自處理. 以紅外光信號為例, 將信號通過一個(gè)截止頻率為31.2 Hz的低通濾波器, 濾去信號中的雜波, 此時(shí)的信號可以看成是原信號的直流分量. 而交流分量需要通過高通濾波器來獲得, 脈搏波頻率很低, 將高通濾波器截止頻率設(shè)為0.159Hz, 信號經(jīng)高通濾波器后直流分量被濾除, 得到的便是交流分量, 信號所呈現(xiàn)的形態(tài)便是生理學(xué)上的脈搏波信號. 具體電路如圖5所示.

    圖5 濾波電路

    濾波器的截止頻率為

    當(dāng)8=9=,1=2=時(shí), 式(5)可化為

    設(shè)置8、9為51kΩ,10、11為1MΩ,1、2為0.1μF,3、4為1μF, 可得低通濾波器、高通濾波器的截止頻率分別為31.2、0.159Hz.

    通過高通濾波器后的信號即為原始脈搏波信號, 但是信號十分微弱, 需進(jìn)一步放大才能進(jìn)行后續(xù)處理. 同時(shí), 需要一個(gè)加法器將信號抬升到0V以上, 以便于控制器的采集, 保證信號采集完整.

    3.2 鼾聲采集電路設(shè)計(jì)

    鼾聲檢測的電路主要由鼾聲采集電路以及主控芯片STM32F407組成. 采用MAX9814作為鼾聲采集電路, 如圖6所示, MAX9814是一款低成本高性能麥克風(fēng)放大器, 具有自動增益控制(AGC), 能夠設(shè)置40、50、60dB 3種增益, 電源電壓范圍2.7~5.5V, 內(nèi)置2V低噪聲麥克風(fēng)偏置. 麥克風(fēng)用于采集聲音, 輸出(OUTPUT)接入控制器.

    圖6 鼾聲采集電路

    通過鼾聲采集電路采集睡眠時(shí)的聲音, 通過微控制器對采集聲音進(jìn)行相關(guān)算法處理, 從而檢測出鼾聲.

    4 血氧及鼾聲信號處理算法設(shè)計(jì)

    4.1 PPG信號去噪及峰值點(diǎn)提取

    監(jiān)測血氧飽和度最重要的是對脈搏波進(jìn)行消噪以及峰值點(diǎn)提取, 兩者的結(jié)果會影響血氧飽和度監(jiān)測的準(zhǔn)確度. 對脈搏波消噪算法的研究主要針對的是基線漂移等低頻噪聲. 環(huán)境干擾等高頻噪聲可以通過低通濾波器濾除, 而基線漂移會改變信號的形態(tài)和幅值, 影響血氧飽和度計(jì)算的準(zhǔn)確度, 需要對其濾除. 本文利用小波變換方法予以濾除, 小波變換在處理非平穩(wěn)低頻信號時(shí)優(yōu)勢較大, 很適合文中脈搏波基線漂移的濾除. 基于小波變換的脈搏波基線漂移濾除步驟如下:

    (1)小波基選取. 根據(jù)各個(gè)小波基的形態(tài), 通過實(shí)驗(yàn)比較, 選取db4作為基線校準(zhǔn)的小波基.

    (2)分解層數(shù)選取. 根據(jù)采樣率、采樣點(diǎn)數(shù)、脈搏波低頻特性, 通過實(shí)驗(yàn)選取5層分解最為合適.

    (3)噪聲估計(jì). 利用5層分解得到第5層上的近似分量來估計(jì)PPG信號中的低頻噪聲.

    (4)信號重構(gòu). 從PPG信號中直接去除第5層的近似分量進(jìn)行重構(gòu), 得出去除基線漂移的信號.

    通過以上步驟, 基本消除了原始脈搏波的基線漂移, 如圖7所示.

    圖7 原始波形及去除基線漂移后波形

    對于脈搏波峰值點(diǎn)的提取, 主要是提取脈搏波的波峰以及波谷, 根據(jù)波峰與波谷的差值得出計(jì)算血氧飽和度的交流分量. 本文利用微分閾值法來提取脈搏波的波峰波谷, 步驟如下:

    (1)獲取一段原始脈搏波去噪后的數(shù)據(jù), 對其進(jìn)行一階差分.

    (2)找出差分結(jié)果中的最大值, 根據(jù)設(shè)定一個(gè)閾值, 閾值系數(shù)為, 即閾值為, 通過實(shí)驗(yàn)測得取0.9最為合適.

    (3)在差分結(jié)果中尋找大于閾值的點(diǎn), 可以得出此點(diǎn)左邊第一個(gè)過零點(diǎn)以及右邊第一個(gè)過零點(diǎn), 對應(yīng)于原始數(shù)據(jù)中便是需要提取的波谷與波峰點(diǎn).

    通過以上步驟, 可以準(zhǔn)確地識別出脈搏波的波峰波谷點(diǎn), 如圖8所示.

    圖8 脈搏波波峰波谷識別

    差分結(jié)果中每個(gè)周期內(nèi)的最大值左右兩邊的第一個(gè)過零點(diǎn)對應(yīng)于原始脈搏波的波谷與波峰.利用微分閾值法, 提取脈搏波的波峰與波谷, 波峰與波谷的均值為計(jì)算血氧飽和度的直流分量, 差值為交流分量. 得到PPG信號的直流分量以及交流分量, 經(jīng)定標(biāo)后便可以計(jì)算出血氧飽和度.

    通過小波變換消除原始脈搏波的基線漂移, 微分閾值法提取脈搏波主波的波峰波谷, 經(jīng)數(shù)據(jù)處理后, 便可以實(shí)時(shí)地監(jiān)測血氧飽和度.

    4.2 鼾聲檢測算法設(shè)計(jì)

    目前語音信號處理技術(shù)較為成熟, 但對于鼾聲檢測與識別卻很少, 鼾聲信號的檢測是指對傳感器采集的語音數(shù)據(jù)進(jìn)行分析, 檢測和識別出其中的鼾聲. 根據(jù)語音信號的處理過程, 對鼾聲信號進(jìn)行處理主要分為4個(gè)步驟: 預(yù)處理、端點(diǎn)檢測、特征提取以及鼾聲識別.

    4.2.1預(yù)處理

    對語音信號進(jìn)行預(yù)處理, 以消除因人類發(fā)聲機(jī)制給語音分析帶來的影響, 預(yù)處理過程主要分為預(yù)加重和分幀加窗, 預(yù)加重就是對語音信號的高頻部分進(jìn)行增強(qiáng)補(bǔ)償, 實(shí)現(xiàn)方式如下:

    語音信號在長時(shí)間上是一個(gè)非平穩(wěn)過程, 但是在短時(shí)間內(nèi)(10~30ms)近似為一個(gè)穩(wěn)定過程, 即短時(shí)相對穩(wěn)定, 可以將語音信號分為一小段一小段來處理, 即進(jìn)行分幀操作. 而分幀操作通常是通過窗函數(shù)來截取, 本文采用漢明窗進(jìn)行截取, 漢明窗具有更好的低通特性, 能更好地保留鼾聲信號.

    4.2.2端點(diǎn)檢測

    從傳感器采集到的用戶睡眠聲音, 除鼾聲以外, 還可能存在呼吸聲、關(guān)門聲、空調(diào)聲等噪聲, 這些都屬于有聲段, 而環(huán)境寂靜聲屬于無聲段, 需要很好地區(qū)分出有聲段與無聲段, 檢測出有聲段便是端點(diǎn)檢測的目的. 本文采用較為簡單的雙門限法進(jìn)行檢測. 基于短時(shí)能量與短時(shí)過零率的雙門限法如圖9所示. 通過設(shè)置短時(shí)能量的門限th, 初步判斷出有聲段范圍在1左邊以及2右邊, 在1左邊和2右邊搜索. 設(shè)置短時(shí)過零率的門限th, 即可判斷出有聲段的范圍在1和2之間. 對每幀數(shù)據(jù)的短時(shí)能量和短時(shí)過零率判斷, 一旦滿足其中一個(gè)條件就被視為有聲段. 需要注意的是, 文中需要以環(huán)境聲的短時(shí)能量和短時(shí)過零率作為參考來設(shè)定門限閾值.

    圖9 雙門限法示意

    4.2.3特征提取

    端點(diǎn)檢測之后, 獲取語音信號的有聲段, 有聲段中包含鼾聲以及非鼾聲. 從有聲段中找出鼾聲片段, 首先需要做的便是特征提取, 選擇一個(gè)特征來區(qū)分這些聲音. 本文提取的是MFCC, MFCC可以更好地反映鼾聲與非鼾聲片段的差異[16], 通過將語音信號的線性頻率轉(zhuǎn)為人耳感知的梅爾頻率, 就可以更加貼合人耳對語音信號的聽覺感知機(jī)制. 特征提取步驟如圖10所示.

    圖10 提取MFCC基本步驟

    4.2.4鼾聲識別

    提取了特征參數(shù), 就可以對有聲段語音數(shù)據(jù)進(jìn)行分類, 也就有了將鼾聲信號識別出的手段. 鼾聲識別方法主要有隱馬爾科夫模型、K近鄰算法以及K均值算法等[17]. 對比3種方法, 并考慮到本系統(tǒng)需要實(shí)時(shí)監(jiān)測呼吸暫停癥狀, 采用K近鄰(KNN)較為簡單, 其他兩種方法不能做到實(shí)時(shí)監(jiān)測. K近鄰算法的實(shí)現(xiàn)如圖11所示.

    圖11 K近鄰算法

    K近鄰算法應(yīng)用于鼾聲識別首先需要確定訓(xùn)練集樣本, 且每個(gè)樣本的特征參數(shù)以及所屬類別需要給出, 為了準(zhǔn)確識別, 訓(xùn)練集樣本需要多樣性, 除了鼾聲, 還需要電視聲音、呼吸聲、關(guān)門聲等, 樣本種類越多識別的準(zhǔn)確率越高. 確定完訓(xùn)練集樣本, 便可以對新樣本特征參數(shù)與訓(xùn)練集樣本特征參數(shù)進(jìn)行分類識別, 具體步驟是首先測試新樣本與訓(xùn)練集樣本特征參數(shù)的歐氏距離, 然后按距離排序, 選取與新樣本特征參數(shù)距離最近的個(gè)訓(xùn)練集樣本, 確定最近的個(gè)訓(xùn)練集樣本中類別出現(xiàn)最多的樣本, 并將其作為新樣本的預(yù)測類別, 這樣便可以將每一幀的語音數(shù)據(jù)類別識別出來, 從而識別出鼾聲片段.

    5 系統(tǒng)軟件設(shè)計(jì)

    系統(tǒng)軟件設(shè)計(jì)主要包括血氧飽和度檢測軟件設(shè)計(jì)以及鼾聲檢測軟件設(shè)計(jì), 整個(gè)系統(tǒng)軟件設(shè)計(jì)流程如圖12所示.

    圖12 系統(tǒng)軟件流程

    系統(tǒng)實(shí)現(xiàn)的功能是先檢測出鼾聲, 再對打鼾者進(jìn)行血氧飽和度的監(jiān)測, 血氧飽和度低于臨界值, 系統(tǒng)識別出該打鼾者患有呼吸暫停癥狀. 各個(gè)部分協(xié)同工作, 對血氧飽和度正常的單純打鼾者, 系統(tǒng)便不會識別出呼吸暫停癥狀.

    6 系統(tǒng)測試結(jié)果

    6.1 血氧飽和度測試結(jié)果

    10名受試者保持靜坐姿勢10s, 佩戴血氧探頭進(jìn)行5次測量, 得出平均值. 為了模擬睡眠時(shí)缺氧, 采用屏息的方法使血氧飽和度下降, 屏息20s, 佩戴血氧探頭進(jìn)行5次測量, 得出平均值. 結(jié)果對比見表1. 受實(shí)驗(yàn)條件限制, 采用華為手環(huán)FDS-B19作為標(biāo)定儀器, 對本文設(shè)計(jì)的系統(tǒng)和手環(huán)測量的血氧飽和度進(jìn)行對比, 兩者具有很好的一致性. 由于個(gè)體的差異, 誤差會有一定范圍的波動, 但總體而言, 本文系統(tǒng)的測量誤差在正常范圍內(nèi), 測量較為精確.

    表1 血氧飽和度對比測試實(shí)驗(yàn)結(jié)果 %

    6.2 鼾聲識別測試結(jié)果

    在實(shí)驗(yàn)室環(huán)境中, 利用真實(shí)睡眠環(huán)境下的含有1min鼾聲片段的錄音進(jìn)行鼾聲識別測試. 而對于鼾聲的標(biāo)定, 采用人工對錄音進(jìn)行鼾聲數(shù)量統(tǒng)計(jì), 可以得出最終的鼾聲識別準(zhǔn)確率. 具體測試結(jié)果見表2.

    表2 鼾聲識別測試結(jié)果

    從采集的5份錄音樣本的識別率來看, 鼾聲識別達(dá)到了較好的水平. 但本文的重點(diǎn)并不是識別鼾聲數(shù), 而是能準(zhǔn)確識別出用戶打鼾. 由于運(yùn)行在微控制器上, 不能采用較為復(fù)雜的算法, 而且處理器處理速度有限, 訓(xùn)練集樣本也不能太多, 這樣會導(dǎo)致最終的識別率略低于一些復(fù)雜算法的識別結(jié)果. 總體而言, 本文系統(tǒng)已經(jīng)取得了較好的識別效果, 滿足通常的需求.

    6.3 伴有呼吸暫停癥狀的鼾聲識別結(jié)果

    為了模擬實(shí)際睡眠環(huán)境, 受試者佩戴血氧探頭進(jìn)行實(shí)時(shí)的血氧飽和度監(jiān)測, 同時(shí)利用含有鼾聲片段的錄音以及非鼾聲片段錄音進(jìn)行鼾聲識別, 測試結(jié)果見表3.

    表3 伴有呼吸暫停的鼾聲患者測試結(jié)果

    根據(jù)鼾聲檢測結(jié)果以及血氧飽和度監(jiān)測結(jié)果, 當(dāng)系統(tǒng)檢測到用戶睡眠時(shí)在打鼾, 同時(shí)血氧飽和度下降, 系統(tǒng)便可以自動識別出該打鼾者患有呼吸暫停癥狀. 而對于未打鼾的用戶, 當(dāng)監(jiān)測到其血氧飽和度下降時(shí), 可能是患有心血管或者高血壓等疾病, 這些由其他疾病引起的血氧飽和度下降并不是本文研究的重點(diǎn), 系統(tǒng)不作識別.

    7 結(jié)論

    本文提出了一種基于血氧飽和度及鼾聲檢測的呼吸暫停監(jiān)測系統(tǒng), 將血氧飽和度與睡眠時(shí)的鼾聲相結(jié)合, 對傳統(tǒng)的用鼾聲來檢測呼吸暫停癥狀方法進(jìn)行了改進(jìn), 增加了血氧飽和度的監(jiān)測, 通過血氧飽和度來檢測SAS的準(zhǔn)確率相比于單純鼾聲檢測提高了2%左右. 采集用戶睡眠時(shí)的語音數(shù)據(jù), 識別出其中的鼾聲部分, 便可以檢測出用戶的打鼾狀態(tài). 對正在打鼾的用戶, 同時(shí)監(jiān)測其血氧飽和度下降的程度來識別呼吸暫停癥狀. 識別出用戶患有呼吸暫停癥狀后可以對其進(jìn)行一定的刺激, 從而在一定程度上緩解呼吸暫停癥狀. 該系統(tǒng)可以實(shí)時(shí)地監(jiān)測血氧飽和度以及鼾聲, 工作穩(wěn)定, 檢測方便, 無創(chuàng)傷, 具有一定的應(yīng)用前景.

    [1] 楊賀. 基于AHI指數(shù)的OSAHS篩查系統(tǒng)的研究與實(shí)現(xiàn)[D]. 大連: 大連理工大學(xué), 2015.

    [2] van Maanen J P, de Vries N. Long-term effectiveness and compliance of positional therapy with the sleep position trainer in the treatment of positional obstructive sleep apnea syndrome[J]. Sleep, 2014, 37(7):1209-1215.

    [3] 王輝, 施戀, 王玉梅. 51例睡眠打鼾患者應(yīng)用多導(dǎo)睡眠監(jiān)測系統(tǒng)監(jiān)測的臨床分析[J]. 中國醫(yī)學(xué)工程, 2019, 22(12):55.

    [4] Farishta F, Hadi M S, Srinivas V. Obstructive sleep apnea hypopnea syndrome (OSAHS): A case review[J]. International Journal of Pharmaceutics, 2018, 3(3):15-16.

    [6] Tuncer A T, Bilgen M. A viable snore detection system: Hardware and software implementations[J]. Turkish Journal of Electrical Engineering & computer sciences, 2020, 28(2):1059-1069.

    [7] Hassan A R, Haque M A. An expert system for automated identification of obstructive sleep apnea from single-lead ECG using random under sampling boosting[J]. Neuro- computing, 2017, 235:122-130.

    [8] Pathinarupothi R K, Dhara P J, Rangan E S, et al. Single sensor techniques for sleep apnea diagnosis using deep learning[C]. 2017 IEEE International Conference on Healthcare Informatics, 2017:524-529.

    [9] Pevernagie D, Aarts R M, de Meyer M. The acoustics of snoring[J]. Sleep Medicine Reviews, 2010, 14(2):131-144.

    [10] Liu Y, Huang W, Zou J, et al. Sea level nocturnal minimal oxygen saturation can accurately detect the presence of obstructive sleep apnea in a population with high pretest probability[J]. Sleep and Breathing, 2020, 10(10):1007- 1010.

    [11] 黎圣峰, 龐宇, 高小鵬. 便捷式血氧信號檢測裝置設(shè)計(jì)[J]. 傳感器與微系統(tǒng), 2017, 36(3):111-113.

    [12] 王倩, 丁虹霞, 張?jiān)? EMD在消除光電容積脈搏波中運(yùn)動偽跡的應(yīng)用[J]. 先進(jìn)技術(shù)研究通報(bào), 2010, 4(6): 2-6.

    [13] 陳帝良, 方震, 趙湛. 集成多生理參數(shù)監(jiān)測的終端設(shè)計(jì)[J]. 傳感器與微系統(tǒng), 2016, 35(1):95-97.

    [14] 盧蒙蒙. 實(shí)現(xiàn)與血氧光譜分離的脈搏色素譜測量方法研究[D]. 長春: 吉林大學(xué), 2018.

    [15] 康華光, 陳大欽, 張林. 電子技術(shù)基礎(chǔ)模擬部分[M]. 5版. 北京: 高等教育出版社, 2005:228-266.

    [16] Jiang Y, Peng J, Zhang X. Automatic snoring sounds detection from sleep sounds based on deep learning[J]. Physical and Engineering Sciences in Medicine, 2020, 43(2):679-689.

    [17] Vesperini F, Romeo L, Principi E. Convolutional re- current neural networks and acoustic data augmentation for snore detection[M]//Esposito A, Faundez-Zanuy M, Morabito F C, et al. Neural Approaches to Dynamics of Signal Exchanges. Singapore: Springer Singapore, 2019: 35-46.

    (責(zé)任編輯 韓 超)

    Design of an apnea monitoring system by combining blood oxygen saturation and snoring detection

    XIE Yunquan, LI Hong, WU Yangbo, XIE Jianjun*

    ( Faculty of Electrical Engineering and Computer Science, Ningbo University, Ningbo 315211, China )

    People with chronic snoring often suffer from obstructive sleep apnea syndrome (SAS), which is extremely harmful to their health. To tackle the issue of inaccuracy with the conventional technique of using snoring for apnea detection, an apnea monitoring system based on the combination of blood oxygen saturation and snoring detection is proposed. A microcontroller is used in combination with external blood oxygen detection circuit and snoring detection circuit to realize real-time monitoring of blood oxygen saturation and identification of snoring spells. When the system finds that the user is snoring while the blood oxygen saturation decreases to a certain level, it can recognize the apnea symptoms of the snorer. For snorers suffering from apnea symptoms, it can help them stop further snoring, so as to mitinate the apnea symptoms to a certain extent. On the basis of accurate and real-time monitoring of blood oxygen saturation and snoring recognition, the whole system can detect the apnea symptoms of snorers more conveniently and more accurately, which provides a new method for real-time monitoring of apnea syndrome.

    SAS;blood oxygen saturation;snoring detection;stop snoring

    TP311

    A

    1001-5132(2021)02-0009-08

    2020?10?13.

    寧波大學(xué)學(xué)報(bào)(理工版)網(wǎng)址: http://journallg.nbu.edu.cn/

    寧波市自然科學(xué)基金(2018A610401); 浙江省基礎(chǔ)公益研究計(jì)劃(LGF19B050003).

    謝運(yùn)泉(1994-), 男, 安徽壽縣人, 在讀碩士研究生, 主要研究方向: 嵌入式系統(tǒng)設(shè)計(jì). E-mail: 83672311@qq.com

    謝建軍(1974-), 男, 浙江余姚人, 工程師, 主要研究方向: 電路與系統(tǒng). E-mail: xiejianjun@nbu.edu.cn

    猜你喜歡
    鼾聲脈搏血氧
    用心感受狗狗的脈搏
    金橋(2018年2期)2018-12-06 09:30:14
    智能血氧飽和度監(jiān)測系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)
    《今日評說》評今日浙江 說時(shí)代脈搏
    基于51系列單片機(jī)的穿戴式脈搏儀的設(shè)計(jì)
    電子制作(2017年14期)2017-12-18 07:07:50
    基于STM32血氧心率檢測儀的研制
    近紅外可穿戴設(shè)備中脈搏波的呼吸率檢測
    基于血氧模擬前端AFE4490的無創(chuàng)血氧測量模塊設(shè)計(jì)
    如雷鼾聲驚醒“名校夢”,15歲少年不堪重負(fù)臥軌自殺
    爸爸的鼾聲
    乳腺血氧與多普勒聯(lián)合檢查在乳腺癌早期篩查中的應(yīng)用
    国产爱豆传媒在线观看| 久久久久久久久久黄片| 精品一区二区三区人妻视频| 一区二区三区高清视频在线| 日本成人三级电影网站| av在线蜜桃| 成人毛片a级毛片在线播放| 99久久久亚洲精品蜜臀av| 91久久精品电影网| 人妻制服诱惑在线中文字幕| 啦啦啦啦在线视频资源| 国产精品一二三区在线看| 麻豆av噜噜一区二区三区| 亚洲精品456在线播放app| 一本精品99久久精品77| 精品国内亚洲2022精品成人| 精品久久久久久久人妻蜜臀av| 国产精品一区二区三区四区免费观看| 熟妇人妻久久中文字幕3abv| 国产精品1区2区在线观看.| 亚洲一区二区三区色噜噜| 中文字幕av成人在线电影| 亚洲自拍偷在线| 久久精品国产自在天天线| 91麻豆精品激情在线观看国产| 亚洲欧洲日产国产| 国产高清有码在线观看视频| 一边摸一边抽搐一进一小说| 国产精品一区www在线观看| 简卡轻食公司| 干丝袜人妻中文字幕| 国产 一区 欧美 日韩| 小蜜桃在线观看免费完整版高清| 日韩国内少妇激情av| 亚洲中文字幕一区二区三区有码在线看| 91久久精品电影网| 日韩中字成人| 日韩制服骚丝袜av| 国产国拍精品亚洲av在线观看| 国产精品国产三级国产av玫瑰| 一级二级三级毛片免费看| 成人特级av手机在线观看| 69人妻影院| 亚洲成人精品中文字幕电影| 国产免费一级a男人的天堂| 少妇人妻一区二区三区视频| 中文精品一卡2卡3卡4更新| 最近视频中文字幕2019在线8| 午夜亚洲福利在线播放| 一本一本综合久久| 极品教师在线视频| 国产色婷婷99| 国内精品一区二区在线观看| 成人特级黄色片久久久久久久| 国产精品久久久久久亚洲av鲁大| 日韩中字成人| 国产国拍精品亚洲av在线观看| 日本免费a在线| 麻豆国产av国片精品| 婷婷色av中文字幕| 国产乱人偷精品视频| 国产伦在线观看视频一区| 婷婷精品国产亚洲av| ponron亚洲| 国产精品av视频在线免费观看| 淫秽高清视频在线观看| 毛片一级片免费看久久久久| 亚洲欧美中文字幕日韩二区| 精品99又大又爽又粗少妇毛片| 可以在线观看毛片的网站| 一级黄色大片毛片| 99久久久亚洲精品蜜臀av| 日本免费一区二区三区高清不卡| 12—13女人毛片做爰片一| 欧美一区二区精品小视频在线| 91av网一区二区| 国产精品,欧美在线| 国产亚洲91精品色在线| 欧美性猛交╳xxx乱大交人| 我要搜黄色片| 国产黄片视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品久久国产高清桃花| 一本久久中文字幕| 国产中年淑女户外野战色| 国产美女午夜福利| 欧美另类亚洲清纯唯美| 成人美女网站在线观看视频| 亚洲精华国产精华液的使用体验 | 亚洲一区二区三区色噜噜| 神马国产精品三级电影在线观看| 亚洲精品成人久久久久久| 色吧在线观看| a级毛片a级免费在线| 日日干狠狠操夜夜爽| 99久久精品国产国产毛片| 免费在线观看成人毛片| 日本撒尿小便嘘嘘汇集6| 国产av一区在线观看免费| 97超碰精品成人国产| 午夜福利在线观看免费完整高清在 | av天堂中文字幕网| ponron亚洲| 久久精品国产亚洲av香蕉五月| 综合色av麻豆| 在现免费观看毛片| 欧美高清性xxxxhd video| 久久久久久久久中文| 国产白丝娇喘喷水9色精品| 22中文网久久字幕| 最新中文字幕久久久久| 国产视频首页在线观看| 免费搜索国产男女视频| 在线免费十八禁| 成人二区视频| 亚洲人成网站高清观看| 成人永久免费在线观看视频| 亚洲第一区二区三区不卡| 黄片无遮挡物在线观看| 97超视频在线观看视频| 免费大片18禁| 成人性生交大片免费视频hd| 午夜久久久久精精品| 久久久久久久久中文| 高清午夜精品一区二区三区 | 少妇人妻精品综合一区二区 | 国产高清不卡午夜福利| 午夜福利在线在线| 99热只有精品国产| 我要搜黄色片| 亚洲欧美日韩高清专用| 国内久久婷婷六月综合欲色啪| 欧美变态另类bdsm刘玥| 九草在线视频观看| 午夜福利在线在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产高清视频在线观看网站| 好男人视频免费观看在线| 亚洲国产欧美人成| 身体一侧抽搐| 少妇被粗大猛烈的视频| 国产精品日韩av在线免费观看| 欧美在线一区亚洲| 蜜桃亚洲精品一区二区三区| 黄色欧美视频在线观看| 99热6这里只有精品| 日韩成人av中文字幕在线观看| 国产成人a∨麻豆精品| 深夜a级毛片| 国产男人的电影天堂91| 91在线精品国自产拍蜜月| 99在线视频只有这里精品首页| 亚洲久久久久久中文字幕| 国产激情偷乱视频一区二区| 插逼视频在线观看| 亚洲自拍偷在线| 亚洲av成人av| 男女视频在线观看网站免费| 在线免费观看不下载黄p国产| 十八禁国产超污无遮挡网站| 内射极品少妇av片p| 国产高潮美女av| 久久久久久久久大av| 女人十人毛片免费观看3o分钟| 精品无人区乱码1区二区| 亚洲av成人精品一区久久| 网址你懂的国产日韩在线| 亚洲av中文av极速乱| 美女黄网站色视频| 男女视频在线观看网站免费| 欧美3d第一页| kizo精华| 男的添女的下面高潮视频| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品日韩av片在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久九九热精品免费| 国产一级毛片在线| 国产激情偷乱视频一区二区| 一进一出抽搐gif免费好疼| av免费在线看不卡| 欧美激情在线99| 国产69精品久久久久777片| 人人妻人人看人人澡| 国内精品美女久久久久久| 97超视频在线观看视频| 夜夜爽天天搞| 成人二区视频| 91精品一卡2卡3卡4卡| 日韩人妻高清精品专区| 天堂√8在线中文| 久久久久网色| 国产精品蜜桃在线观看 | 婷婷色综合大香蕉| 男女边吃奶边做爰视频| 男人和女人高潮做爰伦理| 国产乱人偷精品视频| 少妇高潮的动态图| 哪个播放器可以免费观看大片| 欧美一区二区亚洲| 成人av在线播放网站| 99精品在免费线老司机午夜| 插逼视频在线观看| 国产亚洲91精品色在线| 亚洲欧美日韩高清专用| 久久久a久久爽久久v久久| 精品不卡国产一区二区三区| www.色视频.com| 最近中文字幕高清免费大全6| 欧美丝袜亚洲另类| 欧美潮喷喷水| 成人毛片a级毛片在线播放| 麻豆国产97在线/欧美| 给我免费播放毛片高清在线观看| 老师上课跳d突然被开到最大视频| 麻豆久久精品国产亚洲av| 色综合亚洲欧美另类图片| 91在线精品国自产拍蜜月| 女的被弄到高潮叫床怎么办| 国内精品一区二区在线观看| 26uuu在线亚洲综合色| 青春草亚洲视频在线观看| 亚洲欧美中文字幕日韩二区| 亚洲三级黄色毛片| 亚洲中文字幕一区二区三区有码在线看| 菩萨蛮人人尽说江南好唐韦庄 | 久久精品综合一区二区三区| 美女 人体艺术 gogo| 欧美变态另类bdsm刘玥| 日产精品乱码卡一卡2卡三| 最近视频中文字幕2019在线8| 美女高潮的动态| 欧美日韩国产亚洲二区| eeuss影院久久| 岛国在线免费视频观看| 九九爱精品视频在线观看| 人体艺术视频欧美日本| 精品午夜福利在线看| 在线免费十八禁| 色5月婷婷丁香| 看非洲黑人一级黄片| 好男人视频免费观看在线| 国产 一区 欧美 日韩| 麻豆久久精品国产亚洲av| 麻豆国产97在线/欧美| 秋霞在线观看毛片| 97超碰精品成人国产| 中文欧美无线码| 亚洲欧美日韩卡通动漫| 麻豆av噜噜一区二区三区| 特级一级黄色大片| 欧美+亚洲+日韩+国产| 少妇裸体淫交视频免费看高清| 午夜激情福利司机影院| 午夜视频国产福利| 欧美又色又爽又黄视频| 听说在线观看完整版免费高清| 欧美成人免费av一区二区三区| 人人妻人人澡人人爽人人夜夜 | 久久人人爽人人片av| 18+在线观看网站| 国产三级在线视频| 我要搜黄色片| 老司机福利观看| 在线天堂最新版资源| 床上黄色一级片| 欧美丝袜亚洲另类| 婷婷色综合大香蕉| 深爱激情五月婷婷| 免费在线观看成人毛片| 欧美人与善性xxx| 成人国产麻豆网| 99久久精品国产国产毛片| 午夜精品一区二区三区免费看| ponron亚洲| 亚洲av一区综合| 一卡2卡三卡四卡精品乱码亚洲| 亚洲色图av天堂| 日日撸夜夜添| 好男人在线观看高清免费视频| 日本一本二区三区精品| 一进一出抽搐gif免费好疼| 亚洲无线在线观看| 久久精品久久久久久久性| 久久精品国产亚洲av涩爱 | 国产单亲对白刺激| 色综合色国产| www.色视频.com| 亚洲欧美日韩高清专用| 免费av观看视频| 久久久午夜欧美精品| 一夜夜www| 精品人妻一区二区三区麻豆| 国产亚洲精品久久久久久毛片| 波多野结衣高清作品| 非洲黑人性xxxx精品又粗又长| 尾随美女入室| 久久精品久久久久久噜噜老黄 | 免费观看的影片在线观看| 99精品在免费线老司机午夜| 亚洲精品乱码久久久久久按摩| 中文字幕熟女人妻在线| 国产精品久久电影中文字幕| av在线亚洲专区| 国产乱人视频| 欧美性感艳星| 亚洲图色成人| 寂寞人妻少妇视频99o| 亚洲av一区综合| 老师上课跳d突然被开到最大视频| 不卡一级毛片| 中国国产av一级| 欧美高清性xxxxhd video| 国产精品久久久久久久电影| a级毛片a级免费在线| 精品久久久久久成人av| 久久精品91蜜桃| 波多野结衣高清作品| 日韩一区二区三区影片| 亚洲经典国产精华液单| 哪里可以看免费的av片| 国产精品伦人一区二区| 国产亚洲精品久久久久久毛片| 国产精品日韩av在线免费观看| 男人狂女人下面高潮的视频| 精品不卡国产一区二区三区| 亚洲欧洲日产国产| 人妻久久中文字幕网| 亚洲人成网站在线播放欧美日韩| 亚洲国产欧美人成| 久久精品久久久久久久性| 白带黄色成豆腐渣| 少妇的逼好多水| 又爽又黄无遮挡网站| 欧美三级亚洲精品| a级毛片免费高清观看在线播放| 欧美成人精品欧美一级黄| 夜夜看夜夜爽夜夜摸| 国产精品精品国产色婷婷| 国产乱人偷精品视频| 日韩一本色道免费dvd| 夜夜爽天天搞| 可以在线观看毛片的网站| 国产一区二区在线观看日韩| 最近手机中文字幕大全| 18+在线观看网站| 少妇人妻一区二区三区视频| 久久久久九九精品影院| kizo精华| 亚洲欧美成人综合另类久久久 | 99在线视频只有这里精品首页| 精品少妇黑人巨大在线播放 | 哪里可以看免费的av片| 性欧美人与动物交配| 欧美日韩综合久久久久久| 美女xxoo啪啪120秒动态图| 亚洲国产精品久久男人天堂| 91午夜精品亚洲一区二区三区| 夜夜爽天天搞| 欧美高清成人免费视频www| 精品一区二区免费观看| 99热这里只有是精品在线观看| av福利片在线观看| 欧美性感艳星| 日本一本二区三区精品| 亚洲精品日韩av片在线观看| 国产淫片久久久久久久久| 波野结衣二区三区在线| 国产一区二区激情短视频| 美女高潮的动态| 成人特级黄色片久久久久久久| 欧美一区二区亚洲| 午夜久久久久精精品| 91av网一区二区| 国产麻豆成人av免费视频| 国产精品一区二区三区四区免费观看| 深夜精品福利| 一级av片app| 久久热精品热| 午夜福利在线观看吧| 日韩国内少妇激情av| 成年女人永久免费观看视频| 欧美成人免费av一区二区三区| 日韩欧美 国产精品| 亚洲高清免费不卡视频| 国产伦理片在线播放av一区 | 51国产日韩欧美| 亚洲成人av在线免费| 久久精品综合一区二区三区| 精品人妻一区二区三区麻豆| 久久99蜜桃精品久久| 好男人在线观看高清免费视频| 久久鲁丝午夜福利片| 久久久精品94久久精品| 日韩欧美国产在线观看| 欧美又色又爽又黄视频| 18禁裸乳无遮挡免费网站照片| 国产午夜精品论理片| 亚洲国产欧美在线一区| 亚洲一区高清亚洲精品| 成人特级黄色片久久久久久久| 最近中文字幕高清免费大全6| 欧美最黄视频在线播放免费| 村上凉子中文字幕在线| 97热精品久久久久久| 天堂av国产一区二区熟女人妻| 久久久久久伊人网av| 如何舔出高潮| 久久欧美精品欧美久久欧美| 日本一二三区视频观看| 两性午夜刺激爽爽歪歪视频在线观看| 色综合站精品国产| 级片在线观看| 日韩中字成人| 亚洲电影在线观看av| 日日干狠狠操夜夜爽| 精品久久久久久久人妻蜜臀av| 国产精品一及| 又爽又黄无遮挡网站| 22中文网久久字幕| 午夜免费激情av| 久久鲁丝午夜福利片| 国产精品三级大全| 老师上课跳d突然被开到最大视频| 此物有八面人人有两片| 99在线视频只有这里精品首页| 黄色日韩在线| 午夜亚洲福利在线播放| 国产激情偷乱视频一区二区| 99久久中文字幕三级久久日本| 老司机福利观看| 亚洲精品日韩在线中文字幕 | 三级男女做爰猛烈吃奶摸视频| 国产精品女同一区二区软件| 亚洲在线观看片| 热99在线观看视频| 婷婷色av中文字幕| 久久久午夜欧美精品| 在线免费十八禁| 在线天堂最新版资源| 亚洲av第一区精品v没综合| 日韩在线高清观看一区二区三区| 久久热精品热| 精品久久久久久久末码| 中国美白少妇内射xxxbb| 天堂√8在线中文| 直男gayav资源| 亚洲成人中文字幕在线播放| 99热6这里只有精品| 1000部很黄的大片| 亚洲欧美清纯卡通| 国产精品1区2区在线观看.| 亚洲av一区综合| 在线播放无遮挡| 狂野欧美激情性xxxx在线观看| 男女那种视频在线观看| 精品一区二区免费观看| 精品99又大又爽又粗少妇毛片| www日本黄色视频网| 久久99热6这里只有精品| 午夜福利高清视频| 国产成人影院久久av| 国产午夜精品一二区理论片| 99在线人妻在线中文字幕| 成人美女网站在线观看视频| 亚洲av.av天堂| 高清毛片免费看| 男女下面进入的视频免费午夜| 成人鲁丝片一二三区免费| 亚洲av中文av极速乱| 亚洲成人久久爱视频| 12—13女人毛片做爰片一| 亚洲最大成人手机在线| 日本黄色视频三级网站网址| 一本久久中文字幕| 一个人看的www免费观看视频| 女人十人毛片免费观看3o分钟| 久久精品国产亚洲av涩爱 | 亚洲在线自拍视频| 国产精品国产三级国产av玫瑰| 丝袜喷水一区| 国产精品免费一区二区三区在线| 亚洲第一电影网av| 久久久久久伊人网av| 国产高清三级在线| 97超视频在线观看视频| 观看免费一级毛片| 久久中文看片网| 国产精品蜜桃在线观看 | 一个人观看的视频www高清免费观看| 欧美日韩国产亚洲二区| 美女xxoo啪啪120秒动态图| 亚洲国产日韩欧美精品在线观看| 国国产精品蜜臀av免费| 日韩视频在线欧美| 97超视频在线观看视频| 国产av一区在线观看免费| 久久99热这里只有精品18| 日本在线视频免费播放| 精品一区二区免费观看| 国产在视频线在精品| 搡女人真爽免费视频火全软件| 晚上一个人看的免费电影| 亚洲精品成人久久久久久| 校园春色视频在线观看| 色噜噜av男人的天堂激情| 亚洲av第一区精品v没综合| 不卡视频在线观看欧美| 一个人免费在线观看电影| 伊人久久精品亚洲午夜| а√天堂www在线а√下载| 熟妇人妻久久中文字幕3abv| 一个人免费在线观看电影| 免费看光身美女| 综合色丁香网| 欧美日本视频| 听说在线观看完整版免费高清| 亚洲av中文字字幕乱码综合| 久久久久久国产a免费观看| 国产伦精品一区二区三区四那| 色哟哟·www| 日本在线视频免费播放| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲欧美日韩高清专用| 国产精品综合久久久久久久免费| 精品一区二区三区视频在线| 我的老师免费观看完整版| 日韩高清综合在线| 国产精品一区二区在线观看99 | 网址你懂的国产日韩在线| 久久6这里有精品| 成人亚洲欧美一区二区av| 最近手机中文字幕大全| av.在线天堂| 99在线视频只有这里精品首页| 午夜福利高清视频| 韩国av在线不卡| 亚洲精品日韩av片在线观看| 九九爱精品视频在线观看| 99久久中文字幕三级久久日本| 国产综合懂色| 国产成年人精品一区二区| 成人欧美大片| 噜噜噜噜噜久久久久久91| 精品少妇黑人巨大在线播放 | 色播亚洲综合网| 国产黄片美女视频| 一边摸一边抽搐一进一小说| 插逼视频在线观看| 99在线视频只有这里精品首页| 久久热精品热| 少妇熟女欧美另类| 亚洲自拍偷在线| 日本在线视频免费播放| 久久久久久久久久黄片| 啦啦啦韩国在线观看视频| 亚洲最大成人av| 欧美激情久久久久久爽电影| 夜夜爽天天搞| 国产熟女欧美一区二区| av在线播放精品| 国产熟女欧美一区二区| 男女视频在线观看网站免费| 看免费成人av毛片| 国模一区二区三区四区视频| 欧美三级亚洲精品| 给我免费播放毛片高清在线观看| 日韩欧美 国产精品| 少妇熟女欧美另类| 欧美日韩精品成人综合77777| 国产精品久久电影中文字幕| 亚洲乱码一区二区免费版| 校园春色视频在线观看| 久久婷婷人人爽人人干人人爱| 99久国产av精品| 菩萨蛮人人尽说江南好唐韦庄 | 国产成人aa在线观看| 日本与韩国留学比较| 欧美另类亚洲清纯唯美| 亚洲综合色惰| 成人一区二区视频在线观看| 成人漫画全彩无遮挡| 国产毛片a区久久久久| 午夜老司机福利剧场| 永久网站在线| a级毛色黄片| 久久精品久久久久久噜噜老黄 | 国产色婷婷99| АⅤ资源中文在线天堂| 国产精品人妻久久久久久| or卡值多少钱| 特级一级黄色大片| 午夜福利在线观看吧| 国内精品美女久久久久久| 秋霞在线观看毛片| 亚洲欧美日韩无卡精品| 色尼玛亚洲综合影院| 欧美精品国产亚洲| 尾随美女入室| 99久久无色码亚洲精品果冻| 人人妻人人澡人人爽人人夜夜 | 国产精品乱码一区二三区的特点| 精品久久久久久久末码| 国产女主播在线喷水免费视频网站 | 美女高潮的动态| 99久久精品热视频| 国产av麻豆久久久久久久| 日本免费a在线| 插逼视频在线观看| 三级毛片av免费| 少妇人妻一区二区三区视频| 韩国av在线不卡| 在线观看66精品国产| 国产精品嫩草影院av在线观看| 日本免费一区二区三区高清不卡|