陳浩 聶永琦 嵇秀娟

摘 要：為研究冠心病與脈象圖的相關(guān)性，分別提取冠心病患者與健康者的脈象信號，但人體脈象信號是一種信噪比較低的非平穩(wěn)隨機(jī)信號，因此采用改進(jìn)小波變換，利用Mallatb對脈搏波信號去噪，采用自定的極值法進(jìn)行冠心病患者與健康者的脈搏提取，并根據(jù)特征值的正態(tài)性檢驗和方差齊性檢驗結(jié)果，分別采用非參數(shù)檢驗，參數(shù)檢驗驗證兩組的差異性。

關(guān)鍵詞：脈搏信號 改進(jìn)小波變換去噪 特征值提取 非參數(shù)檢驗

引 言

脈診起源于中國，是中醫(yī)診斷最主要的輔助方法，但是受到醫(yī)生的一些主觀判斷和自身知識和經(jīng)驗的影響，診斷有可能出現(xiàn)一些誤差，因此就需要與科學(xué)技術(shù)結(jié)合起來使之成為一種簡潔、準(zhǔn)確、高效的科學(xué)手段。

冠心病作為全球死亡率最高的疾病之一，已經(jīng)引起國內(nèi)外的廣泛關(guān)注。冠心病在中醫(yī)上又為胸痹，而其中的一些病理特征又可以通過中醫(yī)的“切”感知到脈象而獲得，與脈診有著一定的聯(lián)系。中醫(yī)作為中華文化的瑰寶，也擁有許多脈象診斷的方法。脈象，即指脈搏的快慢、強(qiáng)弱、深淺的情況。脈象是人體機(jī)能狀態(tài)的客觀反應(yīng)之一，而它的脈搏波也能反映出很多重要的生意病理信息。脈診現(xiàn)在也隨著現(xiàn)代科技的發(fā)展不斷走向客觀化，與計算機(jī)技術(shù)建立了聯(lián)系。

但是由于脈象的生理變化是由人體機(jī)能狀態(tài)的變化引起的，不僅受心臟的調(diào)劑，還受到各動脈和人體各種活動的影響，特征的提取就變的很困難。而脈搏的形態(tài)等特征是反映人體各項生理活動最好的方法，現(xiàn)在脈搏特征提取也受到廣泛關(guān)注，但各種分析方法還不成熟，這次我們采用自定的極值法來提取脈象信號的特征。

1 脈搏信號的采集

本次脈搏信號采集系統(tǒng)為自主研發(fā)的CVFD-I型心血管系統(tǒng)狀態(tài)監(jiān)測儀;采集對象是30名中年健康者，30名老年患者和30名冠心病患者，采集過程在安靜的專用試驗室中完成，每位受試者在試驗前一天無飲酒、熬夜等情況，且在試驗開始前仰臥靜息10分鐘。采集過程中受試者采取仰臥姿勢，保持靜息狀態(tài)，平穩(wěn)呼吸，不能說話或移動身體，避免人為因素引起的信號改變。信號的采樣頻率為1000Hz，采集時間為5-10分鐘。

2 脈搏特征分析與去噪

2.1 脈搏特征分析

對于時域來說，正常脈搏波形主要分為主波，重搏前波，重博波，主波與重博前波之間是潮波前谷，重博波與重博前波之間是降中峽：

其中A點為主波波峰，B點為降中峽，C點為重搏波，但A與B之間的重博前波與潮波前谷幾乎消失，這是因為正常脈搏發(fā)生變異，其因素有很多，比如氣候，年齡，體質(zhì)等，這些都會引起脈象的變化，不同的人群在不同的環(huán)境下都會出現(xiàn)特殊的脈象，并且在采集過程中，所有脈象都有這種現(xiàn)象，因此，這不屬于病脈，而是類似于緩脈與滑脈的正常脈搏波形[1]。

h1：主波波峰相對值，反映動脈內(nèi)壓力與容積的最大值。主動脈壓力的最高點，此處是主波，是波形圖的基線至主波峰頂?shù)囊粭l上升曲線，構(gòu)成主波的上升支，反映心室快速射血，動脈壓迅速上升，管壁突然擴(kuò)張，其上升速度主要與心輸出量、心室射血速度、動脈阻力和管壁彈性有關(guān)，可用上升支斜率來表示。冠心病人的主動脈壓力小，造成缺血，相對于正常人來說h1較小。

h2：降中峽相對值，其受動脈血管外周阻力的大小和主動脈瓣關(guān)閉的功能的影響。當(dāng)血管外周阻力降低時，h4的高度也會隨之降低，反之增加。

h3：表示重搏波的幅值。重博波是主動脈M超聲曲線圖，受主動脈壁彈性影響，反應(yīng)主動脈硬化后，舒縮能力下降，即稱彈性降低。冠心病人的主動脈粥樣硬化，動脈壁的彈性降低，收縮能力下降。

h4：重搏波與降中峽波幅值的差值。它反映主動脈瓣的功能狀況、血管彈性和血流流動狀態(tài)。因反流的血液的沖擊，主動脈瓣突然關(guān)閉，反流的血液撞擊在驟然關(guān)閉的主動脈瓣上而被彈回，使主動脈壓再次稍有上升，動脈敢比亦隨之稍有擴(kuò)張。冠心病又稱冠狀動脈性心臟病，它是以導(dǎo)致缺氧的冠狀動脈病，使冠狀動脈粥樣硬化，血流速度慢，血管彈性小，h4較小

h4/h1：重搏波和降中峽幅度差值與主幅度波的比值。與大動脈的順應(yīng)性、左心室的射血功能有關(guān)，也受外周阻力的影響。當(dāng)外周阻力降低時，高度降低，所以增大，h5/h1隨之增大。當(dāng)遇到冠心病人時，冠心病人圖像的不明顯，所以冠心病人的h5/h1趨近于0，而健康人的數(shù)據(jù)明顯。

w/t：主波上1/3高度處的寬度值與其脈搏周期的比值，相當(dāng)于動脈內(nèi)持續(xù)高壓時間在脈動周期中所占比例。主要反映動脈管壁彈性與外周阻力大小。另設(shè)脈搏信號為P（t），對于K值[2]，為平均動脈壓，K值代表脈搏波形面積的變化程度，當(dāng)人體的脈搏波形受到血管外周阻力變化和血管壁硬化程度影響時，K值可較為直觀的反映出其中存在的差異，大小也會在0.350-0.50之間相應(yīng)變化。

2.2 脈搏信號的去噪

人體脈博信號是一種50HZ以下的低頻信號，容易受到機(jī)械、人體等方面因素的影響，從而產(chǎn)生基線漂移，工頻干擾等誤差，因此，對脈搏信號的去噪十分重要。

而傳統(tǒng)的軟、硬閾值法存在缺陷，因而參考文獻(xiàn)[3]改進(jìn)的小波系數(shù)：

式中j為小波分解尺度;為正實數(shù)，且>1。對于標(biāo)準(zhǔn)差為的高斯白噪聲，它所對應(yīng)的小波系數(shù)在每個尺度上是均勻的，并且隨著尺度的增長，的幅值有所減小，考慮到閾值與分解尺度的關(guān)系，取閾值。

進(jìn)行小波變換去噪后，用matlab中的detrend函數(shù)去除基線漂移，并結(jié)合‘constant選項修正零漂移，再將信號在通入通帶截止頻率20阻帶截止頻率30，邊帶區(qū)衰減DB數(shù)設(shè)置為0.1截止區(qū)衰減DB數(shù)設(shè)置為15自定義的低通濾波函數(shù)lowp平滑曲線。

從上圖中可以看出，原脈搏信號中含有較多的高頻噪聲干擾，出現(xiàn)了基線漂移，對比去噪后的脈搏信號如圖1.2.2所示，基線漂移與噪聲明顯去除。另外為方便數(shù)據(jù)處理，將脈搏的始點預(yù)設(shè)為坐標(biāo)軸原點。

3 脈搏信號特征提取

健康者與冠心病患者的h1，h4，w具有較大差異，這是因為冠心病人的主動脈壓力小，造成缺血，相對于正常人來說h1較小。冠心病是導(dǎo)致心肌缺氧，缺血的冠狀動脈病，使冠狀動脈粥樣硬化，血流速度慢，血管彈性小，因此h4，w/t具有明顯差異。所以，本文提取h1，h4，h4/h1，w/t這四個值，并采用北京工業(yè)大學(xué)羅志昌教授提出的脈搏波形特征量K值，建立特征指標(biāo)。

在時域特征值提取過程中，采用自定的極值點算法便可得出h1，h2，h3，h4，w等值，再經(jīng)過部分運(yùn)算處理后便可得到h4/h1，w/t，k值。

3 特征值的統(tǒng)計檢驗

對上述提取的特征值分別進(jìn)行正態(tài)性檢驗，對服從正太分布的h1w/tK進(jìn)行方差齊性檢驗再對h1，w/t，k進(jìn)行t'檢驗，主要結(jié)果：方差齊性檢驗結(jié)果，若P>0.05，則選擇方差齊的t檢驗結(jié)果;否則取方差不齊的t。

本例結(jié)果：除h1顯著性雙側(cè)P值0.693大于0.05外，w/t，K值均小于0.05，分別為0.02，0.00，因此健康與冠心病兩組的特征值h1無統(tǒng)計學(xué)意義，w/t與K值差別有統(tǒng)計學(xué)意義，可認(rèn)為兩組的w/t與K值不同。

對不服從正太分布的特征值h2 h4 h4/h1進(jìn)行方差齊性檢驗，顯著性均小于0.05，因此進(jìn)行H檢驗結(jié)果如下：

由結(jié)果可知兩組的漸進(jìn)顯著性h2>0.05，h3<0.05，h4/h1<0.05，因此可認(rèn)為h4，h4/h1兩者的總體分布不同，具有差異性。

4.討論

脈搏是人體血液循環(huán)系統(tǒng)的產(chǎn)物，其特征與人體的也肌收縮的力量、外周的阻力以及動脈壁的彈性等因素有直接的關(guān)系，也是反映人體氣血病癥的重要生理參數(shù)。所以，對于脈搏波特征的研究和分析使得我們更加清楚自己的身體健康情況。

本文通過提取冠心病患者與健康者脈搏特征值，對建立冠心病與健康的區(qū)分指標(biāo)有著重要意義，本文通過提取特征值，驗證了兩組K值，w/t，h4，h4/h1具有差異性可作為區(qū)分指標(biāo)。

參考文獻(xiàn)

[1] 焦琪玉.脈像特征信號特征值提取與分類識別[J].長春理工大學(xué)學(xué)報.2015.

[2] 羅志昌，張松，楊文鳴.脈搏波波形特征信息的研究[J].北京工業(yè)大學(xué)學(xué)報，1996，22（1）：71-79.

[3] 洪文學(xué)，陳媛媛，張璋. 基于小波變換的脈搏信號去噪[J].電子測量技術(shù)2009，32（9）3.