關(guān)于脈搏波諧波產(chǎn)生機(jī)理、分布特性及相關(guān)心率檢測(cè)方法的研究

關(guān)鍵詞：脈搏波;諧波;頻譜;心率檢測(cè)

中圖分類號(hào)：TP301.6 文獻(xiàn)標(biāo)志碼：A

0 引言（Introduction）

脈診是中醫(yī)學(xué)一種獨(dú)特而歷史悠久的診斷方法，其臨床應(yīng)用在西漢時(shí)期就已經(jīng)發(fā)展成熟[1]。根據(jù)王唯工教授提出的氣血共振理論，脈搏波諧波分別與十二經(jīng)絡(luò)及其所屬臟腑相匹配，每一個(gè)諧波都有其相應(yīng)的中醫(yī)臨床診斷用途和意義，具體而言，前5個(gè)諧波所對(duì)應(yīng)的經(jīng)絡(luò)或臟腑依次為肝、腎、脾、肺、胃[2]。在一些動(dòng)物實(shí)驗(yàn)中發(fā)現(xiàn)，臟器功能的不同變化與特定的諧波有關(guān)[3]。在針灸治療時(shí)，徑向脈搏波的諧波分析可以先于一些血液動(dòng)力學(xué)指標(biāo)或參數(shù)出現(xiàn)相應(yīng)的變化[4]。研究脈搏波諧波產(chǎn)生機(jī)理及其基本分布特征具有重要的中西醫(yī)臨床應(yīng)用價(jià)值和意義。在實(shí)際應(yīng)用中，通常通過對(duì)脈搏波頻譜基頻的計(jì)算，實(shí)現(xiàn)心率參數(shù)的檢測(cè)，但這類方法的檢測(cè)效果不穩(wěn)定，因此本文通過實(shí)驗(yàn)提出了一種基于脈搏波諧波特性的更為穩(wěn)定的心率檢測(cè)方法，此方法也可用于類似的準(zhǔn)周期信號(hào)。

1 脈搏波諧波研究（Study of pulse wave harmonics）

根據(jù)傅里葉變換原理，包括脈搏波在內(nèi)的所有連續(xù)信號(hào)都可以分解為若干個(gè)不同頻率的正弦信號(hào)的線性組合，利用離散傅里葉變換對(duì)從連續(xù)信號(hào)中采樣處理得到的離散信號(hào)進(jìn)行傅里葉分析。假設(shè)某長(zhǎng)度為N 的離散脈搏波信號(hào)用序列x（n）表示，n=0，1，…，N -1，可以通過公式（1）離散傅里葉變換得到序列的頻譜序列：

由圖2（b）可以看出，復(fù)制并串接后的脈搏波信號(hào)頻譜在單周期頻譜相鄰的兩個(gè)頻率之間插入了零值，而偶數(shù)索引頻率處的幅值為圖1（b）所示的2倍。

為了進(jìn)一步研究脈搏波信號(hào)的變化規(guī)律，復(fù)制3個(gè)脈搏波的波形并將其串接在一起，得到理想脈搏波信號(hào)的原始波形及其頻譜如圖3所示。

由圖3（b）可以看出，復(fù)制3次脈搏波并串接所得信號(hào)的頻率分辨率變成原始譜頻率分辨率的3倍，而且頻譜相當(dāng)于在原始單周期頻譜的相鄰頻率之間插入了兩個(gè)零值。同時(shí)，幅值恰為原來(lái)的3倍，頻譜幅值與串接波形的個(gè)數(shù)成正比。以上現(xiàn)象給予我們啟發(fā)：當(dāng)復(fù)制更多數(shù)目的波形并進(jìn)行串接時(shí)，應(yīng)該存在與圖2（b）和圖3（b）相類似的頻譜分布規(guī)律。

實(shí)際上，與公式（5）的推導(dǎo)過程類似，將圖1（a）所示信號(hào)復(fù)制m 次再串接所得到的理想脈搏波信號(hào)，對(duì)其進(jìn)行離散傅里葉變換，可得到如公式（6）所示的推導(dǎo)過程。

通過以上實(shí)驗(yàn)和推導(dǎo)過程可以看出，在時(shí)域內(nèi)將波形復(fù)制m 次并進(jìn)行串接時(shí)，對(duì)于所得到的理想脈搏波序列而言，與所涉單一周期脈搏波頻譜相比，其頻譜的頻率分辨率有一定比例的提高，并且除各個(gè)幅度不為零的幅值及其相鄰頻率間隔皆變?yōu)閙 倍外，各個(gè)諧波分布皆保持相同的順序。也就是說(shuō)，新頻譜相當(dāng)于在原頻譜中每對(duì)相鄰頻率之間加了m-1個(gè)零，同時(shí)對(duì)應(yīng)諧波幅值也增加到m 倍。

1.2 真實(shí)脈搏波頻譜中諧波分布的特征

從以上實(shí)驗(yàn)和分析可以發(fā)現(xiàn)，通過波形復(fù)制和串接得到的理想脈搏波的頻譜呈現(xiàn)出規(guī)則的諧波分布規(guī)律，諧波的振幅和頻譜分辨率都與脈搏波的長(zhǎng)度有關(guān)。然而在真實(shí)情形下，脈搏波并不是規(guī)整地對(duì)單一波形進(jìn)行重復(fù)復(fù)制和串接，其最終得到的頻譜應(yīng)當(dāng)是許多個(gè)非常相似的串接脈搏波波形頻譜的疊加，并且相應(yīng)的諧波幅度所在頻率位置相對(duì)于彼此，一般都略有偏移。因此，隨著波形數(shù)目的增加和頻率分辨率的提高，在各個(gè)諧波附近彼此疊加進(jìn)而形成波群。雖然理想脈搏波信號(hào)的諧波幅值與所含波形的數(shù)目成正比，但是因?yàn)槊總€(gè)串接的脈搏波不再嚴(yán)格的相同，所以真實(shí)脈搏波頻譜中的諧波幅值并不會(huì)隨著波形數(shù)的增加而按比例地增加。為了考察真實(shí)脈搏波信號(hào)頻譜中的諧波分布規(guī)律，在原來(lái)截取的單個(gè)周期的脈搏波繼續(xù)沿著時(shí)間軸截取多個(gè)周期進(jìn)行實(shí)驗(yàn)，此時(shí)截取信號(hào)的第一個(gè)周期與本文“1.1”節(jié)中的單周期相同。其中，截取連續(xù)兩個(gè)周期的真實(shí)脈搏波信號(hào)及其頻譜如圖4所示。

對(duì)比圖4（b）和圖2（b）可以看出，真實(shí)脈搏波信號(hào)的頻譜分布特征與兩個(gè)周期的理想脈搏波的頻譜分布特征基本相似，即諧波的基頻分布的位置相對(duì)幅度的變化規(guī)律相一致。顯而易見，真實(shí)脈搏波的頻譜幅值在偶數(shù)索引頻率點(diǎn)處并不完全為零，這也表明脈搏波實(shí)際上不是一個(gè)理想的周期信號(hào)。

圖5和圖6分別顯示了繼續(xù)沿著時(shí)間軸截取連續(xù)3個(gè)周期和6個(gè)周期的真實(shí)脈搏波信號(hào)及其頻譜。

由圖5和圖6可以發(fā)現(xiàn)，除一些雜波的幅值接近于零外，真實(shí)脈搏波的頻譜幅值與理想脈搏波的頻譜幅值變化規(guī)律基本一致，但其相對(duì)于復(fù)制并串接后的理想脈搏波信號(hào)對(duì)應(yīng)的各諧波分量的幅度峰值略有減小，這是由于真實(shí)脈搏波具有明顯的準(zhǔn)周期特性。觀察圖6（b）發(fā)現(xiàn)，在0～1 Hz還存在較大幅值的分量，這是由于該脈搏波存在一定的基線漂移，觀察圖6（a）也可以明顯看出其信號(hào)呈現(xiàn)出先下降后上升的漂移趨勢(shì)。

真實(shí)脈搏波與通過復(fù)制相同數(shù)目的脈搏波波形并串接得到的理想脈搏波序列的頻譜中存在相似的諧波分量，這與前面的理論分析結(jié)果一致。考慮到脈搏可以看作由心臟的跳動(dòng)信號(hào)沿著時(shí)間線不斷綿延和準(zhǔn)周期重復(fù)，再經(jīng)過血管及器官的影響所形成的，這與上述實(shí)驗(yàn)中基于將信號(hào)進(jìn)行周期延拓后進(jìn)行傅里葉變換的諧波分析方法和過程完全相似。實(shí)驗(yàn)中采用傅里葉頻譜諧波分析方法所得結(jié)論在很大程度上適用于真實(shí)脈搏波信號(hào)的頻譜分析。

圖7對(duì)應(yīng)的提取到的脈搏波頻譜中諧波頻率的分布如表1所示。

表1中各個(gè)真實(shí)脈搏波信號(hào)的諧波數(shù)與諧波群中心頻率之間的線性關(guān)系如圖8所示。

從圖8中可以看出，脈搏波頻譜的諧波數(shù)與波群峰值對(duì)應(yīng)的頻率之間存在顯著的線性相關(guān)關(guān)系，說(shuō)明我們的推論是合理的。真實(shí)脈搏波的第1諧波是第一個(gè)波群的中心頻率，與理想脈搏波頻譜的第1諧波相對(duì)應(yīng)，而且第1諧波對(duì)應(yīng)了平均心率，同時(shí)數(shù)據(jù)擬合中線性關(guān)系的存在表明關(guān)于諧波次數(shù)所擬合直線的斜率也應(yīng)對(duì)應(yīng)平均心率。每個(gè)波群幅值峰值所在位置對(duì)應(yīng)于每個(gè)諧波頻率，而且相應(yīng)諧波分量的頻率是心跳頻率的整數(shù)倍。這也與王唯工提出的氣血共振理論的基本結(jié)論一致[5]。為了驗(yàn)證此結(jié)論的可靠性和適用性，本研究獲取858份采樣頻率為125 Hz的真實(shí)脈搏波信號(hào)，并求得每份數(shù)據(jù)的12個(gè)諧波群中心頻率與諧波數(shù)的皮爾森相關(guān)系數(shù)（PearsonCorrelation Coefficient）。皮爾森相關(guān)系數(shù)可以反映兩個(gè)變量的線性相關(guān)程度，其值為-1～1，其值越接近1，說(shuō)明兩個(gè)變量的線性相關(guān)程度越強(qiáng)。表2為858份脈搏波諧波群中心頻率與諧波數(shù)的皮爾森相關(guān)系數(shù)分布情況。從表2的數(shù)據(jù)可以看出，所有數(shù)據(jù)的皮爾森相關(guān)系數(shù)都大于0.990且計(jì)算出系數(shù)的平均值為0.998。該結(jié)果進(jìn)一步證明了上述結(jié)論具有普適性。

因此，對(duì)于真實(shí)的脈搏波，諧波的幅值會(huì)隨著信號(hào)長(zhǎng)度的增加（周期個(gè)數(shù)）而增大。但是，由于脈搏波具有準(zhǔn)周期性，每個(gè)周期對(duì)應(yīng)的諧波并非完全等同的，從而在頻域每個(gè)周期對(duì)應(yīng)的各個(gè)諧波數(shù)據(jù)也不完全相等，各周期對(duì)應(yīng)諧波疊加時(shí)彼此略有錯(cuò)位，進(jìn)而在信號(hào)頻譜中導(dǎo)致諧波群出現(xiàn)。通過對(duì)真實(shí)脈搏波的相對(duì)孤立波群幅度峰值進(jìn)行定位，證實(shí)了真實(shí)脈搏波頻譜中各個(gè)諧波所在位置體現(xiàn)為頻譜各個(gè)相對(duì)孤立波群幅度峰值所對(duì)應(yīng)的頻率，并且各個(gè)諧波幅度峰值對(duì)應(yīng)頻率是呈線性的。即第n 個(gè)波群幅度峰值對(duì)應(yīng)頻率約為第1波群幅度峰值對(duì)應(yīng)頻率（也就是基頻）的n 倍，這與理論推導(dǎo)得出的結(jié)論一致。

2 基于脈搏波諧波分布規(guī)律穩(wěn)定計(jì)算心率（Stableheart rate calculation based on the harmonicsdistribution of pulse）

靜息心率是冠心病的預(yù)前因素[6]，并且與中風(fēng)、猝死等非心血管等疾病高度相關(guān)[7]。本文針對(duì)脈搏波頻譜的諧波數(shù)與波群峰值位置對(duì)應(yīng)的頻率之間存在顯著的線性關(guān)系的研究，為相對(duì)穩(wěn)定檢測(cè)心率提供了新的可能。如前文所述，脈搏波信號(hào)頻譜第1諧波的頻率對(duì)應(yīng)平均心率，但由于脈搏波不是理想的周期信號(hào)，所以各心動(dòng)周期所對(duì)應(yīng)的時(shí)間間隔和相應(yīng)的脈搏波波形并不一致，并體現(xiàn)于其各個(gè)諧波的幅度和相位變化之中[8]，這使得基于第1諧波的心率檢測(cè)并不準(zhǔn)確。XU等[9]利用了心率信號(hào)的高次諧波估計(jì)心率和減少呼吸干擾。每個(gè)波群幅值的峰值所在位置對(duì)應(yīng)于每個(gè)諧波頻率，并且其所對(duì)應(yīng)諧波分量的頻率是心跳頻率的整數(shù)倍，那么可以利用每個(gè)波群幅值的峰值計(jì)算心率。基于本文“1.2”的研究成果，為了探究參與計(jì)算的脈搏波諧波數(shù)和心率的關(guān)系，本文針對(duì)5種基于諧波的心率計(jì)算方法進(jìn)行了實(shí)驗(yàn)比較，5種計(jì)算方法分別為第1諧波頻率、諧波線性擬合斜率、線性擬合斜率加上截距、諧波頻率加權(quán)平均值（Weighted Average of Harmonic Frequencies，WAoHF）和相鄰諧波頻率平均差（Mean Difference of AdjacentHarmonic Frequencies， MDoAHF）。其中，第1諧波頻率、多個(gè)諧波線性擬合斜率、線性擬合斜率加上截距所對(duì)應(yīng)的心率值分別由第1諧波頻率H、多個(gè)諧波擬合斜率a、線性擬合斜率a 加上截距b，即（a+b）乘以60得到，而諧波頻率加權(quán)平均HR 對(duì)應(yīng)心率數(shù)由公式（7）得到，其中H為第n 次諧波。

本研究截取采樣率為125 Hz、時(shí)長(zhǎng)為5 s，共計(jì)625個(gè)點(diǎn)的脈搏波進(jìn)行心率計(jì)算，對(duì)于每段信號(hào)計(jì)算不同諧波個(gè)數(shù)對(duì)應(yīng)的上述5種心率，總共計(jì)算960段信號(hào)，計(jì)算每種計(jì)算方法在不同諧波數(shù)參與下的心率，再以臨床測(cè)得的心跳數(shù)為標(biāo)準(zhǔn)值計(jì)算絕對(duì)平均誤差（Mean Absolute Error， MAE），計(jì)算結(jié)果如圖9所示。

通常，利用頻譜計(jì)算心率時(shí)，主要是通過傅里葉變換將信號(hào)分解為不同頻率的成分，從而找到主要的周期性成分作為心率，即第1諧波的頻率。但是，從圖9中可以看出，基于第1諧波頻率得到5 s平均心率誤差為4.526 46 BPM，而且在相較于第1諧波的計(jì)算方法，其他4種方法在3～6個(gè)諧波參與計(jì)算時(shí)的MAE明顯小于第1諧波主頻的結(jié)果，但當(dāng)諧波數(shù)目增加到6個(gè)之后，誤差會(huì)都大于第1諧波的心率計(jì)算結(jié)果。這可能是由諧波位置提取的準(zhǔn)確性降低導(dǎo)致的：對(duì)于時(shí)長(zhǎng)為5 s、采樣頻率為125 Hz的信號(hào)，其對(duì)應(yīng)的頻域分辨率為125/625=0.2 Hz，對(duì)應(yīng)的心率分辨率為12 BPM，分辨率不高的同時(shí)，高次諧波波群不明顯且易受噪聲影響，難以提取到準(zhǔn)確的諧波位置，最終導(dǎo)致心率計(jì)算的誤差增大。從圖9中還可以看出，采用4個(gè)或者5個(gè)諧波進(jìn)行計(jì)算得到的心率MAE是最小的，而前4個(gè)諧波頻率加權(quán)得到的心率MAE是最小的。除了擬合斜率方法，其他的算法都在第4次諧波之后開始出現(xiàn)MAE逐漸增大的情況。

將脈搏波信號(hào)時(shí)長(zhǎng)增加到10 s，實(shí)驗(yàn)結(jié)果如圖10所示。

10 s的信號(hào)頻率分辨率為0.1 Hz，對(duì)應(yīng)的心率分辨率為6 BPM。從圖10中可以看出，所有方法對(duì)應(yīng)的心率MAE都遠(yuǎn)小于5 s信號(hào)的心率MAE，說(shuō)明頻率分辨率對(duì)心率計(jì)算的準(zhǔn)確性有影響，理論上頻率分辨率越高，則計(jì)算越準(zhǔn)確。圖10中，第1諧波計(jì)算出的心率MAE為2.042 BPM，可以看出在利用前4 個(gè)諧波進(jìn)行心率的計(jì)算時(shí)，所獲得的平均絕對(duì)誤差依然最小，為1.011 BPM。并且，通過加權(quán)平均算法計(jì)算得到的心率MAE一直小于其他算法的心率MAE。由此可以得出，在利用頻譜計(jì)算短時(shí)脈搏波信號(hào)的心率時(shí)，利用前4個(gè)諧波進(jìn)行加權(quán)計(jì)算的心率是接近臨床真實(shí)心率的。此外，這種算法具有廣泛的應(yīng)用潛力，可以運(yùn)用于各種具有諧波性質(zhì)的信號(hào)中，同時(shí)它也為利用各類相關(guān)生理信號(hào)計(jì)算心率的研究提供了新的思路和方向。

3 結(jié)論（Conclusion）

對(duì)比理想脈搏波和真實(shí)脈搏波發(fā)現(xiàn)：隨著周期個(gè)數(shù)的增加，理想脈搏波的諧波幅值成比例地增加，而非諧波位置處則不斷增加零值。真實(shí)脈搏波各個(gè)周期相似但不完全相同，導(dǎo)致頻譜非諧波處插入相對(duì)應(yīng)周期個(gè)數(shù)的近零值，諧波處各周期對(duì)應(yīng)諧波彼此錯(cuò)位偏移形成諧波群。脈搏波頻譜中，各諧波所在位置體現(xiàn)為各波群幅度峰值所在位置，各個(gè)波群峰值對(duì)應(yīng)頻率的線性關(guān)系證明了該結(jié)論。進(jìn)一步從脈搏波頻譜出發(fā)對(duì)心率進(jìn)行計(jì)算，發(fā)現(xiàn)基于前4個(gè)諧波加權(quán)方法得到的心率更準(zhǔn)確和穩(wěn)定。

作者簡(jiǎn)介：

陳佳（1999-），男，碩士生。研究領(lǐng)域：醫(yī)學(xué)圖像與生物醫(yī)學(xué)信號(hào)處理。

王思研（1992-），女，碩士生。研究領(lǐng)域：醫(yī)學(xué)圖像與生物醫(yī)學(xué)信號(hào)處理。

顧金聯(lián)（2003-），男，本科生。研究領(lǐng)域：醫(yī)學(xué)圖像與生物醫(yī)學(xué)信號(hào)處理。

陳兆學(xué)（1975-），男，博士，副教授。研究領(lǐng)域：醫(yī)學(xué)圖像與生物醫(yī)學(xué)信號(hào)處理。本文通信作者。