劉冬冬，張玲，楊曉文，張博，武文芳

1.首都醫(yī)科大學(xué) 生物醫(yī)學(xué)工程學(xué)院，北京 100069；2.首都醫(yī)科大學(xué)附屬北京安定醫(yī)院 抑郁癥治療中心，北京100088；3.咸陽市中心醫(yī)院 睡眠疾患診療科，陜西 咸陽 712000；4.金紐合（北京）科技有限公司，北京 100011

基于經(jīng)驗(yàn)?zāi)Ｊ椒纸獾男姆务詈霞夹g(shù)在睡眠分析中的應(yīng)用

劉冬冬1，張玲2，楊曉文3，張博4，武文芳1

1.首都醫(yī)科大學(xué) 生物醫(yī)學(xué)工程學(xué)院，北京 100069；2.首都醫(yī)科大學(xué)附屬北京安定醫(yī)院 抑郁癥治療中心，北京100088；3.咸陽市中心醫(yī)院 睡眠疾患診療科，陜西 咸陽 712000；4.金紐合（北京）科技有限公司，北京 100011

目的探討基于經(jīng)驗(yàn)?zāi)Ｊ椒纸獾男姆务詈霞夹g(shù)在睡眠分析中的應(yīng)用。方法通過分析30例源自臨床多導(dǎo)睡眠監(jiān)測(cè)記錄的胸導(dǎo)心電信號(hào)，利用經(jīng)驗(yàn)?zāi)Ｊ椒纸夥椒ǐ@得瞬時(shí)頻率和瞬時(shí)相位信息，構(gòu)建心肺耦合圖譜。按照睡眠中的循環(huán)交替模式（CAP）分期方法，將睡眠過程劃分為CAP、non-CAP以及清醒/異相睡眠期，采取過零率（ZCR）衡量心肺耦合最大峰值的波動(dòng)程度，反映阻塞性睡眠呼吸暫停與低通氣綜合征（OSAHS）的嚴(yán)重程度。結(jié)果OSAHS患者譜圖分布頻帶集中于低頻區(qū)域，各時(shí)刻最大峰值波動(dòng)較小。通過對(duì)比人工和自動(dòng)分期結(jié)果，可以發(fā)現(xiàn)基于經(jīng)驗(yàn)?zāi)Ｊ椒纸獾男姆务詈霞夹g(shù)可以很精準(zhǔn)地區(qū)分睡眠過程的不同狀態(tài)。OSAHS患者和健康受試者在睡眠心肺耦合圖譜中的最大峰值波動(dòng)規(guī)律存在顯著差異，健康組和輕/中度組的ZCR值具有顯著性差異（P＜0.001）；輕/中度組與重度組的ZCR值存在顯著性差異（P＜0.001），因此耦合最大峰值和睡眠呼吸暫停低通氣指數(shù)均可以作為劃分OSAHS嚴(yán)重程度的指標(biāo)，且二者之間具有極強(qiáng)的負(fù)相關(guān)性（r=-0.77，P=5.8×10-18）。結(jié)論采用經(jīng)驗(yàn)?zāi)Ｊ椒纸饨Y(jié)合心肺耦合技術(shù)的方法，可提供可靠的睡眠微結(jié)構(gòu)以及睡眠呼吸障礙信息，其數(shù)據(jù)采集簡(jiǎn)單、容易執(zhí)行，在可穿戴健康管理以及臨床輔助診斷領(lǐng)域有巨大發(fā)展?jié)摿Α?/p>

睡眠呼吸暫停綜合征；心肺耦合；經(jīng)驗(yàn)?zāi)Ｊ椒纸猓凰叻制冢惠o助診斷

阻塞性睡眠呼吸暫停與低通氣綜合征（Obstructive Sleep Apnea-Hypopnea Syndrome，OSAHS）是一種最為常見的睡眠呼吸障礙，據(jù)國(guó)外調(diào)查數(shù)據(jù)顯示，成年男性發(fā)病率為4%，成年女性發(fā)病率為2%，我國(guó)的實(shí)際發(fā)病率高達(dá)7%～13%[1-3]。OSAHS與神經(jīng)系統(tǒng)疾病關(guān)系密切，是導(dǎo)致日間嗜睡、認(rèn)知功能障礙和某些神經(jīng)系統(tǒng)疾病的重要原因之一，是心腦血管病的重要危險(xiǎn)因素，同時(shí)可能與阿爾茨海默病發(fā)病密切相關(guān)[4-5]。及時(shí)診斷OSAHS并加以治療，可預(yù)防各種并發(fā)癥的發(fā)生，提高患者的生活質(zhì)量。

目前已有的睡眠呼吸暫停綜合征檢測(cè)算法都取得了較好的效果，但是還存在一些不足。使用多通道信號(hào)的算法，存在信號(hào)提取不方便、影響患者睡眠、檢測(cè)過程復(fù)雜等問題；使用單通道信號(hào)的算法存在準(zhǔn)確度不高、泛化能力弱、對(duì)信號(hào)質(zhì)量要求過高、需要人工矯正等問題。其中由于睡眠與自主神經(jīng)系統(tǒng)的調(diào)控有關(guān)，而心電信號(hào)蘊(yùn)含了與自主神經(jīng)系統(tǒng)活動(dòng)規(guī)律相關(guān)的信息，因此受到較多關(guān)注[6-7]，如Kawano等[8]運(yùn)用多元回歸分析法觀察睡眠和覺醒狀態(tài)的平均心率變化規(guī)律。da Silva等[9]通過計(jì)算心動(dòng)周期，心率變異率（Heart Rate Variability，HRV）高頻成分和射血前期來測(cè)量心率和自主神經(jīng)控制情況。Vigo等[10]采用小波變換、樣本嫡對(duì)HRV進(jìn)行分析。江朝暉等[11]運(yùn)用平均頻率、搏間自相關(guān)系數(shù)、LZ復(fù)雜度分析心率變異性。吳鋒等[12]利用譜分析和主成分分析方法挖掘心動(dòng)周期變異性的特征，通過Fisher分類準(zhǔn)則建立睡眠分期全自動(dòng)識(shí)別模型。佟光明等[13]用ROC方法確定HRV晝夜差異的敏感指標(biāo)及臨界值。

由于HRV受運(yùn)動(dòng)、血壓、生理節(jié)律等多種因素影響，得到的信號(hào)通常是不穩(wěn)定的、比較短而且包含大量噪音，使得傳統(tǒng)算法在生理研究和臨床應(yīng)用中有很大局限性。本研究采用心肺耦合（Cardiopulmonary Coupling，CPC）算法分析HRV和基于心電圖的呼吸提取信號(hào)兩者間的耦合關(guān)系，很大程度上避免了HRV不穩(wěn)定性的影響，利用經(jīng)驗(yàn)?zāi)Ｊ椒纸猓‥mpirical Mode Decomposition，EMD）將心電信號(hào)分解為只包含一種振動(dòng)模式，而沒有復(fù)雜疊加波的本征模態(tài)函數(shù)，計(jì)算出瞬時(shí)頻率和瞬時(shí)相位信息，提高CPC算法的時(shí)頻分辨率，通過與基于腦電的睡眠分期標(biāo)準(zhǔn)對(duì)比，驗(yàn)證了新方法在睡眠分期和OSAHS檢測(cè)中的有效性[14-15]。

1 材料與方法

1.1 一般資料

選取北京大學(xué)第一醫(yī)院呼吸內(nèi)科30例患者的多導(dǎo)睡眠監(jiān)測(cè)（Polysomnography，PSG）數(shù)據(jù)，其中包括采樣頻率為128 Hz的胸導(dǎo)聯(lián)心電記錄。數(shù)據(jù)采集使用的為澳大利亞Compumedics公司生產(chǎn)的Siesta 2監(jiān)測(cè)儀。入選的30例受試者均患有不同程度的呼吸暫停/低通氣綜合征，其中男性18例，女性12例，年齡19～75歲，平均年齡為46.5歲，標(biāo)準(zhǔn)差16.3歲。將受試者分為3組：健康組、輕/中度組和重度組，數(shù)據(jù)集中受試者的睡眠呼吸暫停-低通氣指數(shù)的范圍為1.5～88.3 h-1，均值為21.9 h-1，標(biāo)準(zhǔn)差為25.4 h-1，具體數(shù)值見表1。

表1 受試者的睡眠呼吸暫停低通氣指數(shù)值

數(shù)據(jù)排除標(biāo)準(zhǔn)如下：① 當(dāng)移除信號(hào)中的記錄誤差后，仍保持至少6 h的數(shù)據(jù)長(zhǎng)度；② 有80%以上可以用于執(zhí)行基于EMD的心肺耦合分析的信號(hào)數(shù)據(jù)。

1.2 方法

1.2.1 經(jīng)驗(yàn)?zāi)Ｊ椒纸?/p>

經(jīng)驗(yàn)?zāi)Ｊ椒纸獾倪^程大致可以描述為以下步驟：

（3）檢查這個(gè)新函數(shù)是否符合以下兩個(gè)條件：① 在整個(gè)信號(hào)中，過零點(diǎn)的個(gè)數(shù)和極值點(diǎn)的個(gè)數(shù)相等，或者數(shù)目相差不超過1個(gè)；② 在任意時(shí)間點(diǎn)，由上包絡(luò)線與下包絡(luò)線確定的平均包絡(luò)值為0。

圖1 固有模態(tài)函數(shù)提取過程

一般情況下，第一次分解得到的殘余量不會(huì)達(dá)到上述兩個(gè)標(biāo)準(zhǔn)，因此需要繼續(xù)令并重復(fù)上述(1)～(4)步驟，直到滿足設(shè)定的殘余量標(biāo)準(zhǔn)為止。

此時(shí)一個(gè)非線性、非平穩(wěn)的時(shí)間序列就被分解為一系列的IMF和一個(gè)殘余量的形式：

1.2.2 基于經(jīng)驗(yàn)?zāi)Ｊ椒纸獾男姆务詈霞夹g(shù)

（1）從心電信號(hào)中檢測(cè)出Q、R、S波發(fā)生時(shí)刻的幅值信息提取出呼吸相關(guān)信息（ECG-derived Respiratory Signal，EDR）[16]。

（2）得到R-R間期時(shí)間序列，通過窗寬為41個(gè)數(shù)據(jù)點(diǎn)的移動(dòng)平均濾波器去除由于R波檢測(cè)錯(cuò)誤而造成的異常點(diǎn)。由于心電信號(hào)的采樣頻率高達(dá)128 Hz，在實(shí)際的生理過程中人體的節(jié)律要遠(yuǎn)遠(yuǎn)低于這個(gè)值，因此為了加快運(yùn)算速度，文中將R-R間期信號(hào)和EDR信號(hào)進(jìn)行三次樣條插值亞采樣至2 Hz得到

（3）利用經(jīng)驗(yàn)?zāi)Ｊ椒纸鈱?duì)R-R間期信號(hào)和EDR信號(hào)進(jìn)行分解，分別將信號(hào)轉(zhuǎn)換為一系列的IMF和一個(gè)剩余量的形式：

其中P為柯西主值。于是可以分別構(gòu)建R-R序列和EDR序列的解析函數(shù)如下：

由于IMF所具備的特性使得希爾伯特變換后的瞬時(shí)頻率時(shí)時(shí)都有物理意義，因此得到N-N序列和EDR序列在希爾伯特變換后的瞬時(shí)頻率表示為：

（4）計(jì)算R-R序列與EDR序列間的相干及互功率譜。由于相干屬于統(tǒng)計(jì)特性，因此本文在計(jì)算過程中選擇16個(gè)數(shù)據(jù)點(diǎn)（8 s）作為一個(gè)窗進(jìn)行相干的統(tǒng)計(jì)學(xué)計(jì)算。在每個(gè)觀察窗中，首先計(jì)算窗內(nèi)瞬時(shí)頻率的平均值作為該窗的頻率。從R-R序列和EDR序列中找尋頻率匹配的兩個(gè)窗，例如，那么可以按照下面的公式首先計(jì)算互功率譜：

其中m和n表示在R-R序列中的第m個(gè)窗和EDR序列中的第n個(gè)窗。經(jīng)過計(jì)算，將時(shí)頻信息進(jìn)行匹配便可得到時(shí)間-頻率-互功率值的譜。

另一方面，相干是用于衡量?jī)蓚€(gè)信號(hào)間相位差一致性的參數(shù)，其計(jì)算公式如下：

利用互功率時(shí)頻譜和相干時(shí)頻譜所包含的信息，可以得到心肺耦合程度的定量結(jié)果，其計(jì)算公式如下：

1.3 統(tǒng)計(jì)學(xué)分析

采用SPSS軟件進(jìn)行統(tǒng)計(jì)學(xué)分析，利用方差分析（ANOVA）與基于Student-Newman-Kenls統(tǒng)計(jì)量的posthoc后驗(yàn)分析方法來統(tǒng)計(jì)不同分組的計(jì)量資料，以P＜0.001為差異有統(tǒng)計(jì)學(xué)意義。

2 結(jié)果

2.1 基于經(jīng)驗(yàn)?zāi)Ｊ椒纸獾男姆务詈蠄D譜

基于經(jīng)驗(yàn)?zāi)Ｊ椒纸獾男姆务詈蠄D譜按照縱軸方向可以劃分出3個(gè)頻率區(qū)間：高頻耦合區(qū)（0.15～0.4 Hz）、低頻耦合區(qū)（0.04～0.15 Hz）和超低頻耦合區(qū)（0.003～0.04 Hz）。OSAHS患者和正常受試者的心肺耦合圖譜對(duì)比見圖2，存在如下特征：① OSAHS患者譜圖分布頻帶集中于低頻區(qū)域，而正常受試者頻帶高低變化明顯；② OSAHS患者譜圖各時(shí)刻最大峰值波動(dòng)較小，而正常受試者各時(shí)刻最大峰值波動(dòng)明顯。

圖2 基于經(jīng)驗(yàn)?zāi)Ｊ椒纸獾男姆务詈蠄D譜

2.2 睡眠分期

作為PSG睡眠分期技術(shù)的延伸，睡眠中的循環(huán)交替模式（Cyclic Alternating Pattern，CAP）是一種評(píng)估睡眠質(zhì)量較好的補(bǔ)充方法。用CAP參數(shù)來衡量睡眠的不穩(wěn)定性，能提供PSG測(cè)量不能提供的微觀信息，有助于深入了解睡眠的微觀結(jié)構(gòu)。通常，CAP率的增加是睡眠質(zhì)量降低的一個(gè)顯著特征[17]。將心肺耦合圖譜中的高頻耦合、低頻耦合和極低頻耦合分別對(duì)應(yīng)于non-CAP、CAP和清醒/異相睡眠狀態(tài)。通過對(duì)比人工和自動(dòng)分期結(jié)果，可以發(fā)現(xiàn)基于經(jīng)驗(yàn)?zāi)Ｊ椒纸獾男姆务詈霞夹g(shù)可以很精準(zhǔn)地區(qū)分睡眠過程的不同狀態(tài)（圖3）。

圖3 人工與自動(dòng)睡眠分期結(jié)果對(duì)比

2.3 OSAHS嚴(yán)重程度評(píng)估

OSAHS患者和健康受試者在睡眠心肺耦合圖譜中的最大峰值波動(dòng)規(guī)律存在顯著差異，采取過零率（Zero Crossing Rate，ZCR）衡量心肺耦合最大峰值的波動(dòng)程度，進(jìn)而反映出OSAHS的嚴(yán)重程度。

健康組、輕/中度組和重度組的睡眠圖譜最大峰值的ZCR值見表2。健康組和輕/中度組的ZCR值具有顯著性差異（P＜0.001）；輕/中度組與重度組的ZCR值存在顯著性差異（P＜0.001），見圖4。因此耦合最大峰值和AHI均可以作為劃分OSAHS嚴(yán)重程度的指標(biāo)，且二者之間具有極強(qiáng)的負(fù)相關(guān)性（Spearman’s相關(guān)性分析）r=-0.77，P=5.8×10-18，見圖5。

表2 睡眠圖譜最大峰值的過零率

圖4 最大峰值過零率比較

圖5 最大峰值過零率與睡眠呼吸暫停-低通氣指數(shù)間的相關(guān)性分析

3 討論

目前臨床金標(biāo)準(zhǔn)所使用的PSG睡眠分期及睡眠疾病診斷方式，需要在患者身上接入多種生理信號(hào)的采集設(shè)備，不僅改變了受試者正常的睡眠環(huán)境，而且在診斷時(shí)存在個(gè)人差異，此外，基于形態(tài)學(xué)的R&K睡眠分期方式，目前也已有越來越多的研究者提出質(zhì)疑[18]?；诮?jīng)驗(yàn)?zāi)Ｊ椒纸獾男姆务詈霞夹g(shù)，信號(hào)數(shù)據(jù)獲取簡(jiǎn)便，在傳統(tǒng)HRV分析方法的基礎(chǔ)上，通過觀察耦合關(guān)系使計(jì)算結(jié)果更加穩(wěn)定，具有明確的生理意義；考慮到心電信號(hào)具有非線性特性，以經(jīng)驗(yàn)?zāi)Ｊ椒纸夥椒ㄗ鳛樾姆务詈铣潭扔?jì)算的前處理過程[19-20]；利用瞬時(shí)頻率和瞬時(shí)相位信息提高睡眠圖譜的時(shí)頻分辨率，對(duì)睡眠中微結(jié)構(gòu)的建立和保持睡眠的完整性具有至關(guān)重要的作用；對(duì)比于臨床CAP睡眠分期，結(jié)果有較高的一致性，作為PSG方法的一種延伸，可以更好地了解由疾病干擾引起睡眠變化的本質(zhì)及神經(jīng)生理機(jī)制。

另一方面，臨床中對(duì)于睡眠呼吸暫停綜合指標(biāo)AHI的獲取需要參考多種生理信號(hào)的變化，而從睡眠心肺耦合圖譜中提取得到的最大峰值過零率在評(píng)估OSAHS嚴(yán)重程度方面，與臨床常規(guī)指標(biāo)AHI存在很高的一致性，可以作為一種OSAHS診斷及治療效果評(píng)價(jià)的輔助指標(biāo)。

[1] Young T,Peppard PE,Gottlieb DJ．Epidemiology of obstructive sleep apnea:a population health perspective[J]．Am J Respir Crit Care Med,2002,165(9):1217-1239．

[2] Stradling JR,Davies RJ．Sleep 1:Obstructive sleep apnoea/hypopnoea syndrome:definitions, epidemiology,and natural history[J]．Tho rax, 2004,59(1):73-78．

[3] 趙忠新．臨床睡眠障礙學(xué)[M]．上海:第二軍醫(yī)大學(xué)出版社,2003, 16-17．

[4] 朱靜,張曉君．睡眠呼吸暫停綜合征認(rèn)知功能障礙研究新進(jìn)展[J]．國(guó)外醫(yī)學(xué)(耳鼻咽喉科學(xué)分冊(cè)),2004,28(5):273-275．

[5] Sateia MJ．Neuropsychological impairment and quality of life in obstructive sleep apnea[J]．Clin Chest Med,2003,24(2):249-259．

[6] Vanoli E,Adamson PB,Ba-Lin,et al．Heart rate variability during specific sleep stages．A comparison of healthy subjects with patients after myocardial infarction[J]．Circulation,1995,91(7):1918-1922．

[7] 江麗儀,吳效明．睡眠時(shí)相與心率變異性的關(guān)系研究[J]．生物醫(yī)學(xué)工程學(xué)雜志,2011,28(1):148-152．

[8] Kawano Y,Tamura A,Watanabe T,et al．Influence of the severity of obstructive sleep apnea on heart rate[J]．J Cardiol,2010,6(1):27-34．

[9] da Silva SP,Hulce VD,Backs R W．Effects of obstructive sleep apnea on autonomic cardiac control during sleep[J]．Sleep Breath,2009, 13(2):147-156．

[10] Vigo DE,Dominguez J,Guinjoan SM,et al．Nonlinear analysis of heart rate variability within independent frequency components during the sleep-wake cycle[J]．Auton Neurosci,2010,154(1-2):84-88．

[11] 江朝暉,李繼偉,馮煥清,等．R-R間期分析與睡眠分期[J]．生物醫(yī)學(xué)工程研究,2003,(3):17-20．

[12] 吳鋒,俞夢(mèng)孫,成奇明,等．基于譜分析方法提取心動(dòng)周期變異性中的睡眠結(jié)構(gòu)信息[J]．生物醫(yī)學(xué)工程學(xué)雜志,2004,21(2): 212-214．

[13] 佟光明,郭繼鴻,張兆國(guó),等．應(yīng)用心率變異性初篩睡眠呼吸暫停低通氣綜合征[J]．中國(guó)心臟起搏與心電生理雜志,2007, 21(1):44-46．

[14] Thomas RJ,Mietus JE,Peng CK,et al．An electrocardiogram-based technique to assess cardiopulmonary coupling during sleep[J]．Sleep, 28(9):1151-1161．

[15] Huang NE,Shen Z,Long SR,et al．The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]．P Roy Soc Lond B Bio,1998,454:903-995．

[16] Moody GB,Mark RG,Zoccola A,et al．Derivation of respiratory signals from multi-lead ECGs[J]．Computers in Cardiology, 1986,12:113-116．

[17] 彭小虎,王國(guó)鋒,劉軍,等．睡眠微觀結(jié)構(gòu)——CAP與睡眠質(zhì)量評(píng)估[J]．中國(guó)臨床心理學(xué)雜志,2013,21(6):920-923．

[18] Redline S,Kirchner HL,Quan SF,et al．The effects of age, sex, ethnicity, and sleep-disordered breathing on sleep architecture[J]．Arch Intern Med,2004,164(4):406-418．

[19] Liu D,Yang X,Wang G,et al．HHT based cardiopulmonary coupling analysis for sleep apnea detection[J]．Sleep Med,2012, 13(5):503-509．

[20] 祁富貴,李川濤,張華,等．一種非接觸呼吸暫停檢測(cè)技術(shù)的研究[J]．中國(guó)醫(yī)療設(shè)備,2014,29(12):24-27．

Application of the EMD-Based CPC Technique in Sleep Analysis

LIU Dong-dong1, ZHANG Ling2, YANG Xiao-wen3, ZHANG Bo4, WU Wen-fang1
1.School of Biomedical Engineering, Capital Medical University, Beijing 100069, China; 2.MDD Treatment Center, Beijing Anding Hospital, Capital Medical University, Beijing 100088, China; 3.Department of Sleep Respiration Diagnosis and Treatment, Xianyang Central Hospital, Xianyang Shanxi 712000, China; 4.Gold NOAH (Beijing) Tech Co., Ltd, Beijing 100011, China

ObjectiveTo explore application of the EMD (Empirical Mode Decomposition)-based CPC (Cardio-Pulmonary Coupling) technique in sleep analysis.MethodsThrough analysis of 30 cases of thoracic ECG signals recorded by PSG (Poly-Somno-Graphy), the instantaneous frequency and instantaneous phase were obtained with deployment of EMD so as to construct a CPC map. Then, CAP (Cyclic Alternating Pattern) was utilized to divide sleep into three stages: CAP Stage, Non-CAP Stage and Wake/REM (Rapid Eyes Movement) Stage. The waving degree of the maximum CPC peak was measured by ZCR (Zero Crossing Rate), which reflected the severity of OSAHS (Obstructive Sleep Apnea-Hypopnea Syndrome).ResultsThe frequency band of OSAHS patients’ map was distributed centralizedly in the low-frequency areas with small waving changes of the maximum peak at each time. Comparisons were made between manual staging and automatic staging, which revealed that EMD-based CPC could differentiate accurately between the different sleep statuses. Significant differences existed between the waving principles of the maximum peak in OSAHS Patient Group and Healthy Volunteer Group. ZCR values were significantly different between Slight/Middle OSAHS Patient Group and Healthy Volunteer Group (P＜0.001), and between Slight/Middle OSAHS Patient Group and Severe OSAHS Patient Group (P＜0.001). Therefore, the maximum coupling peak value and apneahypopnea Index could be used as indexes to identify the different severity of OSAHS patients. Moreover, strong negative correlation was seen between the two indexes (r=-0.77, P=5.8×10-18).ConclusionCombination of EMD and the CPC technique had proven its easy-to-operate features in data acquisition so as to provide reliable micro-structure and disorder information of sleep, which had huge development potentials in the fi elds of wearable health management and clinically-aided diagnosis.

obstructive sleep apnea-hypopnea syndrome; cardiopulmonary coupling; empirical mode decomposition; sleep staging; aided diagnosis

R54

A

10.3969/j.issn.1674-1633.2015.06.006

1674-1633(2015)06-0028-05

2015-04-13

首都醫(yī)科大學(xué)基礎(chǔ)-臨床科研合作基金（15JL15）；首都醫(yī)科大學(xué)青年教師科研啟動(dòng)基金（3500-114333050225）。

本文作者：劉冬冬，講師。

武文芳，副教授，專業(yè)方向：生理信號(hào)分析，睡眠醫(yī)學(xué)。

通訊作者郵箱：pqu2@163．com