侯麗敏，施曉宇，童 超，劉煥成

(上海大學(xué)通信與信息工程學(xué)院，上海 200444)

0 引 言

睡眠呼吸暫停低通氣綜合征(Sleep Apnea Hypopnea Syndrome,SAHS)是一種患病率高，具有潛在危險(xiǎn)的疾病[1]。在我國(guó)估計(jì)患病率達(dá)3.5%，約有五千萬(wàn)的SAHS患者[2]。隨著醫(yī)療輔助診斷和居家醫(yī)療的發(fā)展，便捷的鼾癥檢測(cè)儀成為國(guó)內(nèi)外的研究熱點(diǎn)[3]。與醫(yī)院臨床診斷使用的多導(dǎo)睡眠監(jiān)測(cè)儀(Polysomnography,PSG)不同的是，便攜式鼾癥自動(dòng)診斷系統(tǒng)僅通過(guò)采集含有豐富病理信息的鼾聲信號(hào)來(lái)分析患者的鼾聲聲學(xué)特性，從而實(shí)現(xiàn)對(duì)SAHS的輔助診斷[4-5]。

鼾聲是SAHS最為明顯的癥狀，易于采集[6]?；l是鼾聲聲學(xué)特性研究中基本且重要的特征之一，已被大量應(yīng)用在鼾聲的檢測(cè)和輔助診斷中[7]。文獻(xiàn)[8]用基頻以及基頻的相關(guān)特征對(duì)整夜鼾聲的錄音自動(dòng)檢測(cè)鼾聲和非鼾聲。文獻(xiàn)[9]用鼾聲響度、基頻和頻譜等多個(gè)特征，基于貝葉斯分類(lèi)器構(gòu)建的核密度模型區(qū)分SAHS患者和單純打鼾者。文獻(xiàn)[10]中研究者計(jì)算鼾聲的基頻密度等特征參數(shù)，用多元回歸法判斷SAHS病癥的嚴(yán)重程度。

大部分研究中用的是歐美人鼾聲的錄音，作為不同的人種，中國(guó)漢族人與西方人的身體外形存在較大差異，本文針對(duì)漢族人的鼾聲錄音做了詳細(xì)的分析研究。用改進(jìn)的線性預(yù)測(cè)殘差倒譜法提取鼾聲的基頻，對(duì)93名受試者的鼾聲做了分析。實(shí)驗(yàn)結(jié)果表明，鼾聲基頻分布特性與SAHS的嚴(yán)重程度存在著密切的關(guān)聯(lián)。SAHS嚴(yán)重程度不同，鼾聲的分布不同，基頻軌跡的抖動(dòng)大小也不同。本研究對(duì)互聯(lián)網(wǎng)時(shí)代的居家醫(yī)療和篩查有廣泛的應(yīng)用價(jià)值。

1 鼾聲產(chǎn)生機(jī)理和基頻計(jì)算

1.1 鼾聲的發(fā)聲原理

打鼾者的上氣道中某個(gè)部位或多個(gè)部位比正常人狹窄。在清醒狀態(tài)下上氣道肌肉收縮使氣道通暢。但是在夜間睡眠狀態(tài)時(shí)，神經(jīng)興奮性下降，肌肉松弛，上氣道坍塌，阻礙氣流的通過(guò)。當(dāng)氣流累積到一定程度就會(huì)沖開(kāi)阻塞部位，引起局部振動(dòng)，從而導(dǎo)致鼾聲的產(chǎn)生。圖1為鼾聲產(chǎn)生機(jī)理示意圖。鼾聲來(lái)源于阻塞位置振動(dòng)的準(zhǔn)周期信號(hào)和呼吸氣流引起的噪聲信號(hào)疊加，共鳴系統(tǒng)指上氣道或部分上氣道，因?yàn)樽枞课豢赡軐⑸蠚獾婪指舫蓛刹糠帧?/p>

圖1 鼾聲的產(chǎn)生機(jī)理圖Fig.1 Mechanism diagram of snore production

每個(gè)受試者整夜打鼾中既有正常鼾聲又有發(fā)生睡眠呼吸暫停之后引發(fā)的SAHS鼾聲，這里簡(jiǎn)稱為簡(jiǎn)單鼾聲(simple snores,SIMP)和SAHS鼾聲，波形如圖2所示。圖2(a)中鼾聲節(jié)奏均勻，每個(gè)完整鼾聲片段幅值相當(dāng)，持續(xù)時(shí)間約3 s，這樣的鼾聲盡管聲量較大，但沒(méi)有出現(xiàn)呼吸暫停，也不會(huì)伴有血氧下降，是非SAHS的，圖2(b)是圖2(a)中紅色虛線框出的單個(gè)簡(jiǎn)單鼾聲片段，稱為簡(jiǎn)單鼾聲；而圖2(c)的鼾聲分布大部分無(wú)規(guī)則，幅值有大有小，持續(xù)時(shí)間有長(zhǎng)有短。其中有一段持續(xù)時(shí)間大于10 s的呼吸暫停，用天藍(lán)色虛線標(biāo)注開(kāi)始時(shí)間和結(jié)束時(shí)間，圖2(c)中紅色虛線代表SAHS鼾聲，圖2(d)是圖2(c)的局部時(shí)域波形圖，其中綠色線標(biāo)注出吸氣段，黑色線標(biāo)注出呼氣段，兩者之間是過(guò)渡段。

無(wú)論SIMP鼾聲還是SAHS鼾聲，一個(gè)完整的打鼾事件包含三個(gè)部分：吸氣段、呼氣段和他們之間的過(guò)渡段。在吸氣段，氣流努力沖破阻塞部分，造成阻塞部位的振動(dòng)，振動(dòng)頻率則為鼾聲的基頻。這段時(shí)間的鼾聲信號(hào)大部分具有準(zhǔn)周期性，與濁音類(lèi)似，幅度高。在過(guò)渡段，進(jìn)行體內(nèi)的氣流交換， 此時(shí)的聲音信號(hào)微弱，類(lèi)似于靜音段。在呼氣段，氣道一般相對(duì)通暢，此時(shí)氣流噪聲起主要作用，聲音波形的幅度小，只有少數(shù)具有準(zhǔn)周期性。因此，本文研究分析以吸氣段為主的鼾聲的聲學(xué)特性。

圖2 鼾聲時(shí)域波形圖Fig.2 The time domain waveforms of snoring sounds

1.2 實(shí)驗(yàn)數(shù)據(jù)

錄音環(huán)境為上海市第六人民醫(yī)院耳鼻喉科睡眠實(shí)驗(yàn)室。錄音設(shè)備：Dell Inspiration 570，聲卡型號(hào)：Creative Audigy 4 Value，麥克風(fēng)：Sony ECMC10，采樣頻率：8 kHz，采樣精度：16 bit。PSG型號(hào)為偉康A(chǔ)lice4。鼾聲錄音與PSG診斷同步進(jìn)行[5]。

研究對(duì)象：由上海市第六人民醫(yī)院耳鼻喉科提供的接收PSG診斷的受試者。每人的PSG報(bào)告給出了呼吸暫停低通氣指數(shù)(Apnea Hypopnea Index,AHI)等相關(guān)診斷結(jié)果。AHI的含義為平均每小時(shí)呼吸暫停低通氣的次數(shù)，單位：次/h。AHI值對(duì)應(yīng)著SAHS嚴(yán)重程度，共分為四種不同的級(jí)別：AHI＜5為單純打鼾者(Non-SAHS，N)，5≤AHI≤15為輕度SAHS患者(Mild-SAHS，L)，15＜AHI≤30為中度SAHS患者(Moderate-SAHS,M)，AHI＞30為重度SAHS患者(Severe-SAHS,S)。

本文訓(xùn)練數(shù)據(jù)集包含單純打鼾N型10人，輕度L型23人，中度M型24人，重度S型36人，共計(jì)93人，受試者的信息如表1所示。對(duì)照著PSG人工切出簡(jiǎn)單鼾聲和SAHS鼾聲，共計(jì)3 995個(gè)片段，如表1中最后兩列。其中，簡(jiǎn)單鼾聲共2 168個(gè)片段，SAHS鼾聲共1 847個(gè)片段。

表1 訓(xùn)練數(shù)據(jù)信息Table 1 General information of training data

1.3 基頻計(jì)算

為了消除上氣道共鳴系統(tǒng)的影響，更好地提取聲源信號(hào)，采用兩次線性預(yù)測(cè)法，用逆濾波得到信號(hào)的殘差，即聲源信號(hào)[11]。然后對(duì)預(yù)測(cè)誤差信號(hào)進(jìn)行倒譜計(jì)算和峰值檢測(cè)，從倒譜域中尋找最大值，最大值對(duì)應(yīng)的位置即為基音周期。具體的計(jì)算基頻過(guò)程如圖3所示。鼾聲信號(hào)s(n)經(jīng)預(yù)處理(分幀、端點(diǎn)檢測(cè))、線性預(yù)測(cè)編碼(Linear Prediction Coding,LPC)逆濾波后得到初始?xì)埐钚盘?hào)e1(n)，經(jīng)過(guò)截止頻率為1 kHz的低通濾波器后再經(jīng)LPC逆濾波后得到聲源信號(hào)e2(n)，進(jìn)行快速傅里葉變換(Fast Fourier Transform,FFT)，得到頻譜信號(hào)Y(ω)，對(duì)其幅值取對(duì)數(shù)后再進(jìn)行傅里葉逆變換得到倒譜信號(hào)c(n)，最后進(jìn)行峰值檢測(cè)即可得到基頻值F0。這種兩次線性預(yù)測(cè)倒譜法，簡(jiǎn)稱DLPCE(Double Linear Prediction And Cepstrum,DLPCE)法。

圖3 DLPCE計(jì)算基頻流程圖Fig.3 Flowchart of pitch calculation by DLPCE

為了說(shuō)明DLPCE法對(duì)鼾聲基頻計(jì)算的準(zhǔn)確性以及確定合適的幀長(zhǎng)，選取20個(gè)準(zhǔn)周期性明顯的鼾聲片段，分別用傳統(tǒng)的倒譜法(Cepstrim Method,CEP)、線性預(yù)測(cè)倒譜法(Linear Prediction and Cepstrum,LPCE)以及DLPCE法計(jì)算鼾聲的基頻。

計(jì)算鼾聲基頻時(shí)，有兩種錯(cuò)誤類(lèi)型：一是受聲道響應(yīng)的影響，峰值檢測(cè)出的最大值在倒譜原點(diǎn)附近，導(dǎo)致得出的基頻值大于實(shí)際值，此類(lèi)錯(cuò)誤稱為甲類(lèi)錯(cuò)誤。第二種是受倒譜二次諧波的影響，估算結(jié)果為實(shí)際基音周期的二倍，導(dǎo)致基頻值是實(shí)際基頻值的一半，即半倍頻。此類(lèi)錯(cuò)誤稱為乙類(lèi)錯(cuò)誤。

鼾聲的變化速率不及語(yǔ)音快，短時(shí)分析幀長(zhǎng)需重新設(shè)置。設(shè)置幀長(zhǎng)分別為40、60、80、100 ms和120 ms，幀移設(shè)為幀長(zhǎng)的30%，分別統(tǒng)計(jì)三種算法不同幀長(zhǎng)計(jì)算出基頻軌跡中錯(cuò)誤的幀數(shù)占所有鼾聲片段幀數(shù)的比例，統(tǒng)計(jì)結(jié)果如圖4所示。

圖4 估計(jì)基頻的錯(cuò)誤比例Fig.4 Error proportion of pitch estimation

由圖4可見(jiàn)：DLPCE法的甲類(lèi)錯(cuò)誤小于CEP法和LPCE法。這是由于使用DLPCE法計(jì)算基頻時(shí)，兩次LPC處理明顯消弱了上氣道的影響；三種算法的乙類(lèi)錯(cuò)誤相當(dāng)。因此，利用DLPCE法計(jì)算鼾聲的基頻要優(yōu)于傳統(tǒng)的LPCE法和CEP法。

幀長(zhǎng)的選取對(duì)錯(cuò)誤率的影響很大，如圖4所示。當(dāng)幀長(zhǎng)選取過(guò)短時(shí)，會(huì)由于幀內(nèi)周期信號(hào)太少導(dǎo)致基頻判斷不準(zhǔn)確，當(dāng)幀長(zhǎng)選取過(guò)長(zhǎng)時(shí)，又會(huì)由于幀內(nèi)周期信號(hào)過(guò)多而導(dǎo)致倒譜二次諧波影響加重，使乙類(lèi)錯(cuò)誤大幅增加。選取幀長(zhǎng)為80 ms更適合于鼾聲基頻的分析。以下的分析幀長(zhǎng)選取為80 ms。

2 鼾聲基頻統(tǒng)計(jì)特性

2.1 鼾聲分類(lèi)

要研究鼾聲的基頻，首先要對(duì)鼾聲的時(shí)域波形有清楚的認(rèn)識(shí)。從鼾聲的時(shí)域波形是否有準(zhǔn)周期的角度，本文將其分為三類(lèi)。時(shí)域波形有一個(gè)或一個(gè)以上較明顯的準(zhǔn)周期為Ⅰ型，時(shí)域波形中一部分有準(zhǔn)周期、一部分無(wú)準(zhǔn)周期為Ⅱ型，時(shí)域波形無(wú)明顯準(zhǔn)周期為Ⅲ型。下面分別針對(duì)這三類(lèi)鼾聲進(jìn)行分析。Ⅰ型鼾聲為時(shí)域波形有一個(gè)或多個(gè)較明顯的準(zhǔn)周期的鼾聲，如圖5(a)所示，其計(jì)算所得的基頻軌跡如圖5(b)所示。Ⅱ型鼾聲為時(shí)域波形一部分有準(zhǔn)周期、一部分無(wú)準(zhǔn)周期的鼾聲，如圖5(c)所示，基頻軌跡如圖5(d)所示。后部分的鼾聲的時(shí)域波形雜亂，沒(méi)有明顯的重復(fù)性。Ⅲ型鼾聲為時(shí)域波形無(wú)明顯的準(zhǔn)周期的鼾聲，如圖5(e)所示，基頻軌跡如圖5(f)所示。該類(lèi)型鼾聲的時(shí)域波形有較多的快變信息，時(shí)域波形沒(méi)有明顯的重復(fù)性，只有隨機(jī)性較強(qiáng)的紊亂波形。

圖5 三種類(lèi)型的鼾聲的波形和基頻軌跡Fig.5 Waveforms and pitch contours of three types of snores

2.2 三類(lèi)鼾聲的分布

用1.3中的兩次LPCE法估算三類(lèi)鼾聲片段的基頻和基頻軌跡。Ⅰ型鼾聲計(jì)算出的基頻軌跡較為準(zhǔn)確，甲類(lèi)和乙類(lèi)錯(cuò)誤甚少，如圖5(b)所示。而對(duì)于Ⅱ型和Ⅲ型鼾聲，計(jì)算其無(wú)準(zhǔn)周期的部分時(shí)，由于信號(hào)本身無(wú)周期，導(dǎo)致計(jì)算出的基頻軌跡抖動(dòng)幅度大且頻繁，如圖5(d)和圖5(f)所示。

根據(jù)估計(jì)的鼾聲基頻軌跡的特點(diǎn)，將基頻軌跡中連續(xù)4幀基頻值相差不超過(guò)10 Hz的情況判為存在準(zhǔn)周期。若存在準(zhǔn)周期的幀數(shù)占鼾聲片段總幀數(shù)的占比超過(guò)70%，則判為Ⅰ型鼾聲，若存在準(zhǔn)周期且占比未超過(guò)70%，則判為Ⅱ型鼾聲，若不存在準(zhǔn)周期，則判為Ⅲ型鼾聲。統(tǒng)計(jì)出不同嚴(yán)重程度患者的SIMP和SAHS鼾聲中三種類(lèi)型鼾聲的分布情況，如圖6所示。圖6(a)～6(d)表示不同嚴(yán)重程度鼾癥中SIMP鼾聲包含三種類(lèi)型鼾聲的比例。圖6(e)～6(h)，表示不同嚴(yán)重程度鼾癥中SAHS鼾聲包含三種類(lèi)型鼾聲的比例。

統(tǒng)計(jì)結(jié)果表明：?jiǎn)渭兇蝼洼p度SAHS患者的鼾聲以I型為主，圖6(a)、6(b)中I型占80%左右，圖6(e)、6(f)中占70%左右，Ⅱ型、Ⅲ型鼾聲占比很少。而中度和重度SAHS患者的鼾聲中Ⅱ型、Ⅲ型鼾聲占比明顯增多，M-SIMP和S-SIMP鼾聲中占比增大到30%左右，M-SAHS和S-SAHS鼾聲中占比增大到45%左右，如圖6(c)、6(d)和圖6(g)、6(h)所示。Ⅲ型鼾聲在SAHS鼾聲中的占比比在SIMP鼾聲中的占比也增加明顯，如圖6(e)、6(f)和圖6(g)、6(h)所示，在M-SAHS和S-SAHS鼾聲中占比增大到25%左右。

圖6 三類(lèi)鼾聲在不同SAHS嚴(yán)重程度患者中的占比Fig.6 The proportions of the three types of snores in patients with different severities of SAHS

這些分布表明：SIMP鼾聲的準(zhǔn)周期性相對(duì)于SAHS鼾聲更明顯；隨著SAHS病癥的加重，無(wú)論SIMP鼾聲還是SAHS鼾聲，更多鼾聲的準(zhǔn)周期性消失，而噪聲信號(hào)增多。

2.3 利用基頻特征診斷SAHS病癥

記基頻差分值dPi為鼾聲的基頻軌跡中，前一幀基頻值Pi減后一幀基頻值Pi+1的絕對(duì)值，即：

由于SAHS病理特性較為復(fù)雜，僅靠鼾聲的基頻特征無(wú)法對(duì)病癥的嚴(yán)重程度有效區(qū)分，因此本文在提取基頻的差分值dP作為特征的同時(shí)，提取Mel頻率倒譜系數(shù)(Mel-Frequency Cepstral Coefficients,MFCC)特征作為實(shí)驗(yàn)的特征，記為dPMFCC。

本文選取與表1不重疊的另外120位受試者的整夜錄音作為測(cè)試數(shù)據(jù)。測(cè)試組受試者的信息見(jiàn)表2。

表2 測(cè)試數(shù)據(jù)信息及測(cè)試結(jié)果Table 2 General information of test data and the test results

實(shí)驗(yàn)采用高斯混合模型(Gaussian Mixture Model,GMM)對(duì)表1中8類(lèi)不同嚴(yán)重程度的SIMP和SAHS鼾聲的特征分別建模。對(duì)測(cè)試者的整晚錄音自動(dòng)檢測(cè)出鼾聲片段，用鼾聲特有的節(jié)律做進(jìn)一步約束，得到候選的呼吸事件，再通過(guò)高斯模型分析判斷是否為呼吸事件。按照式(2)估算出每個(gè)人的AHIdPMFCC值(次/h)[5]。

本文方法獲得的AHIdPMFCC值與PSG 診斷結(jié)果AHIPSG值對(duì)比效果如圖7所示，圖中中心綠色實(shí)線表示同一性，綠色虛線表示95% 的置信區(qū)間。根據(jù)診斷一致性定義：若AHIPSG與AHIdPMFCC兩者均大于40，則為一致；若AHIPSG＜40，兩者的差即AHIPSG-AHIdPMFCC＜10，為一致；若兩者的差＞10，為欠估；若兩者的差＜-10，為過(guò)估[9]。本文方法與黃金標(biāo)準(zhǔn)PSG的一致性為90.8%(109/120)，欠估率為4.2%，過(guò)估率為5%。說(shuō)明本文方法與臨床標(biāo)準(zhǔn)有較好的一致性。

為說(shuō)明基頻特征對(duì)提升SAHS病癥診斷嚴(yán)重程度的有效性，本文利用dPMFCC與單純使用MFCC特征進(jìn)行對(duì)比實(shí)驗(yàn)，實(shí)驗(yàn)結(jié)果見(jiàn)表2中的MFCC正確率與dPMFCC正確率兩欄，正確率的得出是由估計(jì)正確的人數(shù)與該嚴(yán)重程度人數(shù)相除取得的。可以發(fā)現(xiàn)，基頻特征的引入有效提升了輕度患者的診斷效果，正確率從60%上升為76%。結(jié)果表明，利用dPMFCC作為特征判斷鼾癥嚴(yán)重程度的準(zhǔn)確率為85.8%(103/120)，優(yōu)于將MFCC作為特征的80%(96/120)的準(zhǔn)確率。

3 分析及討論

圖7 AHIPSG與AHIdPMFCC對(duì)比圖Fig.7 Comparative chart of AHIPSG and AHIdPMFCC

本文按照時(shí)域波形的準(zhǔn)周期性把鼾聲劃分為三類(lèi)，與文獻(xiàn)[12]中結(jié)果是一致的。文獻(xiàn)發(fā)現(xiàn)軟腭產(chǎn)生的鼾聲是由一些重復(fù)出現(xiàn)且有相似波形的脈沖信號(hào)組成的，而舌根產(chǎn)生的鼾聲比較密集、雜亂，沒(méi)有明顯的重復(fù)性。因此，上氣道阻塞的部位不同，引起的鼾聲信號(hào)時(shí)域波形也不同。文獻(xiàn)[13]發(fā)現(xiàn)單純打鼾者鼾聲頻域曲線表現(xiàn)出明顯的基頻-諧波結(jié)構(gòu)，而SAHS患者鼾聲的時(shí)域曲線表現(xiàn)為多個(gè)振幅、間隔不規(guī)則的復(fù)合波，頻域曲線中無(wú)明顯的基頻-諧波結(jié)構(gòu)。文獻(xiàn)[14]用Praat語(yǔ)音分析軟件對(duì)鼾聲的基頻值做了估計(jì)，單純打鼾者有明顯的基頻，而阻塞性睡眠呼吸暫停低通氣綜合征(Obstructive Sleep Apnea Hypopnea Syndrome,OSAHS)組的鼾聲的基頻分布離散，部分OSAHS鼾聲基頻不能測(cè)出。本文觀察的鼾聲時(shí)域波形和基頻與文獻(xiàn)[13-14]是一致的。

文獻(xiàn)[15]采用中心削波自相關(guān)法計(jì)算鼾聲基頻，計(jì)算出鼾聲的基音周期的范圍為8～40 ms(25～125 Hz)。文獻(xiàn)[16]采用倒譜法計(jì)算得到的基頻范圍為10～35 ms(28～100 Hz)，由于鼾聲的時(shí)變節(jié)奏沒(méi)有語(yǔ)音快，對(duì)鼾聲短時(shí)分析幀長(zhǎng)都大于語(yǔ)音信號(hào)的幀長(zhǎng)，本文也如此。本文的基頻值分布范圍為27～286 Hz，大于這兩篇文獻(xiàn)的結(jié)果。這可能與人種有關(guān)，亞洲人體型小于歐洲人，上氣道結(jié)構(gòu)也會(huì)正比于體型?；l值與體積或重量成反比。

來(lái)自日本的研究文獻(xiàn)[17]對(duì)74名患者進(jìn)行測(cè)試，運(yùn)用上呼吸道測(cè)壓的方法測(cè)得鼾聲產(chǎn)生的振動(dòng)部位有：軟腭振動(dòng)、扁桃體/舌根振動(dòng)、軟腭和扁桃體/舌根混合振動(dòng)、咽喉振動(dòng)，用FFT分析儀計(jì)算出軟腭類(lèi)型鼾聲的基頻為(102.8±34.9)Hz，扁桃體/舌根類(lèi)型鼾聲的基頻為(331.7±144.8)Hz，軟腭和扁桃體/舌根混合類(lèi)型鼾聲的基頻為(115.7±58.9)Hz，咽喉類(lèi)型鼾聲的基頻為250 Hz左右。本文對(duì)107名受試者的五千多個(gè)鼾聲片段做了基頻計(jì)算，基頻分布的范圍27～286 Hz，與文獻(xiàn)[18]的結(jié)果是一致的。

文獻(xiàn)[18]用睡眠鼻內(nèi)窺鏡檢測(cè)法對(duì)54名成年打鼾者的鼾聲發(fā)聲部位進(jìn)行監(jiān)測(cè)，觀察到打鼾時(shí)咽部振動(dòng)的主要形式有：①軟腭或懸雍垂撲動(dòng)，②扁桃體振動(dòng)，③會(huì)厭軟骨振動(dòng)，④舌根振動(dòng)。其中單純舌根振動(dòng)占8%，單純會(huì)厭振動(dòng)占2%，扁桃體、會(huì)厭、舌根共同振動(dòng)產(chǎn)生的鼾聲占20%，而軟腭或軟腭與其他部位共同振動(dòng)占70%。圖6中鼾聲的基頻分布與文獻(xiàn)[18]的結(jié)果一致，說(shuō)明中度和重度SAHS患者的鼾聲基頻包含的軟腭或軟腭與其它部位共同發(fā)生阻塞的占70%左右(Ⅰ、Ⅱ型)，舌根及扁桃體、會(huì)厭、舌根共同發(fā)生阻塞占30%左右(Ⅲ型)。

本文使用的dPMFCC特征聯(lián)合聲源特征中的基頻特征和上氣道特征中的MFCC特征對(duì)鼾癥嚴(yán)重程度實(shí)現(xiàn)有效分類(lèi)。對(duì)單純打鼾型和重度鼾癥判斷的準(zhǔn)確率為100%，而對(duì)輕度和中度鼾癥判斷的正確率不夠理想，可能是這兩類(lèi)患者的鼾癥的阻塞部位和上氣道狀態(tài)不夠穩(wěn)定，導(dǎo)致SIMP和SAHS鼾聲的基頻和MFCC特征的差異性不夠凸顯，從而引起較多的誤判。

4 結(jié) 論

本文對(duì)鼾聲的基頻分布做了較為詳細(xì)的統(tǒng)計(jì)分析。鼾聲基頻值主要取決于上氣道發(fā)生阻塞的部位和程度，鼾聲基頻值的分布為27～286 Hz。鼾聲基頻抖動(dòng)及相關(guān)特性可以一定程度上反映出SAHS病癥嚴(yán)重程度。本文將基頻特征與MFCC特征結(jié)合對(duì)打鼾者的SAHS嚴(yán)重程度進(jìn)行估計(jì)，正確率為85.5%。鼾聲基頻抖動(dòng)及相關(guān)特性可作為鑒別SAHS病癥嚴(yán)重程度的參數(shù)之一。本文的研究成果推動(dòng)了利用鼾聲進(jìn)行SAHS診斷的居家健康醫(yī)療的應(yīng)用。

致謝感謝上海交通大學(xué)附屬上海第六人民醫(yī)院耳鼻喉科的支持。