崔艷斌，張 平

(1.長治醫(yī)學(xué)院生物醫(yī)學(xué)工程系，山西 長治 046000；2.第三軍醫(yī)大學(xué)大坪醫(yī)院野戰(zhàn)外科研究所，重慶 400042)

0 引言

心臟活動受到自主神經(jīng)系統(tǒng)的調(diào)控作用[1]，例如交感神經(jīng)的活性增強會釋放去甲腎上腺素，導(dǎo)致心率的增加；迷走神經(jīng)的活性增強會釋放乙酰膽堿，導(dǎo)致心率的降低，兩者相互依賴，共同維持人體生理環(huán)境的動態(tài)平衡。對于自主神經(jīng)系統(tǒng)功能的測定有很多方法，但是心率變異性作為唯一一個可以定量測定自主神經(jīng)系統(tǒng)活性的方法[2]，能夠測定自主神經(jīng)系統(tǒng)某一成分的變化情況，被廣泛用于醫(yī)療行業(yè)。心率變異性指的是連續(xù)心跳之間瞬時心率的微小變化，即相鄰RR間期的細微漲落。心率變異性反映了交感神經(jīng)、副交感神經(jīng)的張力以及兩者之間平衡的重要指標，與很多疾病尤其是心血管方面的疾病有著很大的關(guān)系，可以用來診斷甚至預(yù)防心血管疾病[3]，是預(yù)測心臟性猝死以及心律失常的一項極其重要的指標[4],對心率變異性的研究有著十分重要的臨床意義和科研意義。

心率變異性的分析方法有時域分析法和頻域分析法[5]。然而時域分析法丟掉了信號的時序信息,不夠全面，頻域分析法是研究史上的里程碑,它揭示了譜中高頻成份反映和呼吸相關(guān)的心臟迷走神經(jīng)的活動,低頻成份反映和血管舒縮系統(tǒng)的變化相關(guān)的心臟交感神經(jīng)活動或交感神經(jīng)和迷走神經(jīng)的共同活動[6]，故心率變異性的頻域指標能夠較好地反映自主神經(jīng)系統(tǒng)活性的變化情況。

在臨床上，5 min時長以內(nèi)的心率變異性分析稱為短時心率變異性。由于病人的不穩(wěn)定性無法采集長時間的數(shù)據(jù)進行分析，而隨著數(shù)據(jù)時長的改變，心率變異性頻域指標也會產(chǎn)生顯著差異[7]。對于心率變異性頻域指標的求取有Welch和Lomb-scargle兩種算法，本研究通過設(shè)置不同時長數(shù)據(jù)的對照組來分析心率變異性的頻域計算方法，得出可用于短時時長心率變異性頻域指標的求取方法。

1 方法

1.1 數(shù)據(jù)的獲取

共17名20～23歲的健康受試者參與了本研究的數(shù)據(jù)采集工作，實驗期間要求所有受試者以平躺姿勢正常呼吸10 min。同時采集受試者的心電信號，采樣率為200 Hz。

1.2 心電信號的處理

本研究中涉及的數(shù)據(jù)處理工作全部在Matlab 2015a平臺上完成。

心率變異性是基于RR間期進行分析的，故為提取正確的R波位置尤為重要。然而在心電信號的采集過程中存在噪聲及因為設(shè)備本身原因有基線漂移等問題存在，使得信號上下浮動和失真厲害，我們對其進行小波變換[8]來去除噪聲。由于引入的是隨機噪聲，小波變換去噪對于隨機信號效果很好，相比于傳統(tǒng)帶通濾波器，小波變換不僅可以濾除噪聲基線漂移，而且還能凸顯信號的特征部分，方便R波的提取，如圖1。然后利用閾值法對R波進行提取，如圖2。

圖1 心電信號的原始圖和小波變換去噪后的心電信號圖

圖2 R波提取圖，圖中小圓點為R波位置

我們對10 min時長的信號分別按照5 min、2 min、1 min的數(shù)據(jù)時長進行截取，然后分別求取兩種方法下的不同時長數(shù)據(jù)的心率變異性頻域指標nLF、nHF、LF/HF的值，并以10 min時長下心率變異性的頻域指標nLF、nHF、LF/HF的值作為標準進行對比，心率變異性頻域指標及其生理意義如表1。

表1 心率變異性頻域指標及其生理學(xué)意義

2 結(jié)果

利用Welch算法求取不同數(shù)據(jù)時長下的心率變異性頻域指標值如表2，利用Lomb-scargle算法求取不同數(shù)據(jù)時長下的心率變異性頻域指標值如表3。結(jié)果表明相比于10 min心率變異性頻域指標值，利用Welch算法求取的2 min時長下心率變異性指標值nLF和nHF均有顯著性差異(P<0.05)，在1 min時長下的心率變異性指標值nLF、nHF和LF/HF均有顯著差異(P<0.05)，而利用Lomb-scargle算法求取的心率變異性指標值不受數(shù)據(jù)時長的影響(P>0.05)。

表2 Welch算法求取心率變異性指標值(均值±標準差)

不同時長組的HRV指標分別與10min相比有統(tǒng)計學(xué)意義：*(P<0.05)表示兩者有顯著性差異。

表3 Lomb-Scargle算法求取心率變異性指標值

不同時長組的HRV指標分別與10 min相比有統(tǒng)計學(xué)意義：*(P<0.05)表示兩者有顯著性差異。

在本研究中我們使用配對T檢驗的方法對實驗組和對照組進行分析。所有的數(shù)據(jù)分析均使用SPSS(version22.0,Inc.Chicago,IL,USA)軟件進行計算，數(shù)據(jù)結(jié)果以均值±標準差的形式展現(xiàn)。同時定義P=0.05，當P<0.05表明兩組數(shù)據(jù)在統(tǒng)計學(xué)上有顯著差異，當P>0.05表明兩組數(shù)據(jù)在統(tǒng)計學(xué)上無顯著差異。

3 討論

本研究的結(jié)果表明，使用Welch算法求取的心率變異性頻域指標nLF、nHF、LF/HF的值會隨著數(shù)據(jù)長度的變化而改變，使用Lomb-scargle算法求取的心率變異性頻域指標nLF、nHF、LF/HF的值不會隨著數(shù)據(jù)長度的變化而改變。即使用Lomb-scargle算法求取短時心率變異性頻域指標是可行的，這樣在臨床上不必對病人進行長時間的監(jiān)測來計算心率變異性的頻域指標值，可以更加快速地發(fā)現(xiàn)發(fā)病原因。

由于Lomb-scargle算法是基于最小二乘法的原理，不需要對原始信號進行插值和重采樣，故可以避免頻譜的失真[9]。所以，當我們將其應(yīng)用于不同時長數(shù)據(jù)的分析時，結(jié)果能夠不受數(shù)據(jù)長度的影響，使得同一指標無統(tǒng)計學(xué)差異。

本研究還存在一些局限。第一，由于自身數(shù)據(jù)采集原因，沒有采集更長久的數(shù)據(jù)來進行對比。第二，實驗的受試者人數(shù)不是很多。在未來的研究中我們打算增加受試者人數(shù)以及監(jiān)測時間來使結(jié)果更有說服力。

綜上所述，利用Lomb-scargle算法可以用于求取短時心率變異性指標的頻域指標。當在臨床上應(yīng)用時，不需要對病人進行長時間的監(jiān)測來獲取數(shù)據(jù)，可以大大縮短對病人情況的監(jiān)測時間，同時結(jié)果準確可靠，具有一定的實際意義。