雷 鳴 郭興明* 肖守中，2

1(重慶大學(xué)生物工程學(xué)院，重慶 400044)

2(重慶博精醫(yī)學(xué)信息研究所 400044)

引言

近年來(lái)，心血管疾病的發(fā)病率和死亡率越來(lái)越高，因此心臟系統(tǒng)疾病的防治和診斷成為醫(yī)學(xué)界面臨的重要問(wèn)題。2005年全世界心臟學(xué)界在關(guān)注女性心臟健康這個(gè)主題上達(dá)成共識(shí)：女性的第一主要死因是心臟病[1]，但其原因和機(jī)理尚待研究和查明。心音信號(hào)是人體重要的生理信號(hào)之一[2]，是心臟及心血管系統(tǒng)機(jī)械運(yùn)動(dòng)狀況的反映，其中包含著心臟各部分的生理信息和病理信息。因此，心音信號(hào)的分析和處理在心臟病的臨床診斷中具有重要意義。

孕婦作為女性的特殊群體，其身體狀況更應(yīng)受到人們的關(guān)注。前期基于心力-心音關(guān)系的研究表明[3-4]，孕婦在妊娠期間動(dòng)用了心臟儲(chǔ)備，心臟負(fù)擔(dān)增大，心臟疾病的出現(xiàn)概率比正常人更大。臨床試驗(yàn)表明[5]，心音是無(wú)創(chuàng)檢測(cè)心臟病的有效手段，當(dāng)心血管疾病尚未發(fā)展到足以產(chǎn)生某種病變以前，心音中就會(huì)出現(xiàn)關(guān)于心臟各個(gè)部分(如心房、心室等)功能狀態(tài)的有助于診斷的重要病理信息。所以，可以利用心音信號(hào)特征來(lái)輔助診斷孕婦的心臟疾病。

最近幾十年興起的分形理論是非線性科學(xué)的一門(mén)分支，是描述自然界和非線性系統(tǒng)中不光滑和不規(guī)則的幾何形體的有效工具。國(guó)外相關(guān)學(xué)者用非線性分形理論中的Higuchi維數(shù)作為腦電數(shù)據(jù)的特征值，分析腦部受傷前后的腦皮層電信號(hào)，并得出分形維數(shù)值能夠較好地反映腦部警覺(jué)性變化的結(jié)論[6]。國(guó)內(nèi)相關(guān)學(xué)者通過(guò)非線性分形方法對(duì)肌電信號(hào)進(jìn)行分析和分類(lèi)處理，證明了分形維數(shù)值可以反映肌肉疲勞狀態(tài)的特征變化[7]。相關(guān)研究也指出了分形理論可以應(yīng)用到心電信號(hào)QRS波群的識(shí)別，為臨床診斷心率不齊提供了有價(jià)值的特征信息[8]。本研究將分形理論引入到心音信號(hào)，不僅為心音信號(hào)的分析與識(shí)別提供了新的手段和方法，而且也為分形理論的應(yīng)用進(jìn)行了有益的拓展。

1 心音信號(hào)的分形特征

心音信號(hào)是一種典型的非平穩(wěn)信號(hào)，具有比較復(fù)雜的時(shí)變性特征，其研究方法分為線性分析和非線性分析。線性分析是傳統(tǒng)的分析方法，其基本思想是：盡管心音信號(hào)是一種時(shí)變、非平穩(wěn)的過(guò)程，但在一小段時(shí)間里，可以將其近似地看成平穩(wěn)過(guò)程。線性分析方法包括時(shí)域分析、頻域分析和時(shí)頻分析。然而，心音信號(hào)的非平穩(wěn)特性決定了采用非線性方法更合適。目前，心音信號(hào)的研究主要集中在進(jìn)行定性的分析，缺少系統(tǒng)的定量分析，影響了心音的臨床應(yīng)用。

分形理論是非線性學(xué)科的一個(gè)分支，在信號(hào)處理方面有著廣泛的應(yīng)用。分形具有兩個(gè)典型的性質(zhì)：一是自相似性，是指局部(部分)與整體或另一局部在形態(tài)、功能、信息等方面具有某種相似性，適當(dāng)?shù)胤糯蠡蚩s小分形對(duì)象的幾何尺寸，整個(gè)結(jié)構(gòu)并不改變;二是無(wú)標(biāo)度性，是指無(wú)論測(cè)量的尺度如何改變，所測(cè)量對(duì)象本身的特性均不發(fā)生變化[9-10]。研究分形理論的關(guān)鍵就是分形維數(shù)值的提取，心音信號(hào)的分形維數(shù)值不是整數(shù)，而是分?jǐn)?shù)。這是因?yàn)樾囊粜盘?hào)是一種非平穩(wěn)、復(fù)雜的非線性生理信號(hào)，有無(wú)窮的細(xì)節(jié)，所以其分形維數(shù)值在一維直線和二維平面之間，又由于信號(hào)沒(méi)有填滿一個(gè)有限的平面，所以其維數(shù)不能等于2，因此它的測(cè)量維數(shù)必定在1～2之間。計(jì)算分形維數(shù)的方法有多種，如霍斯多夫維數(shù)、信息維數(shù)、關(guān)聯(lián)維數(shù)、相似維數(shù)、計(jì)盒維數(shù)等。

本研究采用計(jì)盒維數(shù)分析健康孕婦和妊娠期疾病孕婦的心音信號(hào)，并將其與健康未孕婦女的心音信號(hào)進(jìn)行比較，以證明計(jì)盒維數(shù)值對(duì)于健康孕婦、妊娠期疾病孕婦和健康未孕婦女3組樣本的心音信號(hào)有著顯著的區(qū)分能力，為初步探索其在診斷孕婦心音信號(hào)是否異常方面提供有價(jià)值的信息。

2 計(jì)盒維數(shù)算法

計(jì)盒維數(shù)又稱盒維數(shù)(box dimension，DimB)，由于其數(shù)學(xué)計(jì)算及其經(jīng)驗(yàn)估計(jì)相對(duì)容易，所以是應(yīng)用較廣泛的維數(shù)之一[11-13]。

考慮Rn中的 δ坐標(biāo)網(wǎng)立方體，即[m1δ，(m1+1)δ]× … × [mnδ，(mn+1)δ]形式的立方體，其中m1，…，mn都是整數(shù)。顯而易見(jiàn)，“立方體”在 R1中表示為區(qū)間，在R2中表示為正方形。設(shè)N'δ(F)是δ網(wǎng)立方體與 F相交的個(gè)數(shù)，顯然這是 N'δ(F)個(gè)直徑為 δ的覆蓋 F的集類(lèi)，因此有 Nδ■n(F)≤N＇δ(F)，如果＜1，則

計(jì)盒維數(shù)這個(gè)形式的定義在實(shí)際中有廣泛的應(yīng)用，為計(jì)算一個(gè)平面集F的計(jì)盒維數(shù)，可以構(gòu)造一些邊長(zhǎng)為δ的正方形(或稱為盒子)，然后計(jì)算不同δ值的“盒子”和 F相交的個(gè)數(shù) Nδ(F)。這個(gè)維數(shù)即是當(dāng)δ→0時(shí)Nδ(F)增加的對(duì)數(shù)數(shù)率，或者可以由函數(shù)lg Nδ(F)相對(duì)于lgδ圖的斜率值來(lái)估計(jì)。

3 實(shí)驗(yàn)對(duì)象和方法

3.1 實(shí)驗(yàn)對(duì)象

本研究選擇了378例自愿受試者：妊娠期疾病孕婦64例(研究組)，年齡在 25～37歲，(26.78±9.35)歲，其中患有重度妊娠子癇24例，患有妊娠期高血壓20例，患有妊娠期肝內(nèi)膽汁淤積癥20例;健康孕婦153例(對(duì)照組)，年齡在 25～37歲，(25.49±8.27)歲;在校教師 161例(健康未孕婦女對(duì)照組)，年齡在 25 ～37 歲，(28.53 ±9.63)歲。各組年齡差異均無(wú)統(tǒng)計(jì)學(xué)意義，且所有健康孕婦和未孕婦女受試者無(wú)心血管病史。測(cè)試于2008年7月分別在重慶醫(yī)科大學(xué)附屬一院、鄭州大學(xué)第三附屬醫(yī)院及重慶大學(xué)心臟儲(chǔ)備研究實(shí)驗(yàn)室進(jìn)行，當(dāng)時(shí)受試者均得到知情權(quán)，并簽署了知情同意書(shū)。

3.2 方法

3.2.1 儀器

本研究應(yīng)用博精醫(yī)學(xué)信息研究所研制的“運(yùn)動(dòng)心力監(jiān)測(cè)儀”(ECCM，專(zhuān)利號(hào)01256971.2，第一代產(chǎn)品注冊(cè)證號(hào)：渝藥管械(試)字99第220007)進(jìn)行現(xiàn)場(chǎng)測(cè)試。該系統(tǒng)的硬件組成包括心音脈搏傳感器、CCM信號(hào)采集器、計(jì)算機(jī);心音信號(hào)處理軟件運(yùn)用Visual Basic語(yǔ)言編寫(xiě)，主要由心音信號(hào)顯示模塊、指標(biāo)提取模塊、數(shù)據(jù)處理模塊等組成，運(yùn)行環(huán)境為Windows 98/2000/XP操作系統(tǒng)。

3.2.2 測(cè)試方法

設(shè)備在正常心音的整個(gè)頻率范圍內(nèi)(30～800 Hz)有平坦的頻率響應(yīng)。心音信號(hào)通過(guò)聲卡記錄，采樣頻率為11025 Hz。受試者在安靜狀態(tài)下平靜呼吸，心音脈搏傳感器置于心前區(qū)心尖搏動(dòng)最明顯處，記錄心音圖。操作上的要點(diǎn)是受試者的體位、傳感器放置部位、檢查者對(duì)傳感器的施壓情況應(yīng)該盡量一致，從而盡可能減少人為誤差。

3.2.3 統(tǒng)計(jì)學(xué)方法

組間均值比較采用 t檢驗(yàn)，P＜0.05為差異有統(tǒng)計(jì)學(xué)意義。

4 實(shí)驗(yàn)結(jié)果

根據(jù)非線性分形理論計(jì)盒維數(shù)算法，利用Matlab軟件編程，計(jì)算3組受試者的分形維數(shù)值。分形維數(shù)值的分布：健康未孕婦女組如圖1所示，健康孕婦組如圖2所示，妊娠期疾病孕婦組如圖3所示。分布趨勢(shì)線亦如圖中所示。

由圖1可以看出，該組心音信號(hào)的分形維數(shù)值的均值在1.31左右，所有分形維數(shù)值在該值附近上下波動(dòng)，其中最高值在 1.42左右，最低值在 1.24左右。

由圖2可以看出，該組心音信號(hào)的分形維數(shù)值的均值在1.27左右，所有分形維數(shù)值在該值附近上下波動(dòng)，其中最高值在 1.32左右，最低值在 1.20左右。

由圖3可以看出，該組心音信號(hào)的分形維數(shù)值的均值接近1.39，所有分形維數(shù)值在該值附近上下波動(dòng)，其中最高值在 1.49左右，最低值在 1.29左右。

圖1 健康未孕婦女組分形維數(shù)值Fig.1 FD value of healthy non-pregnant women

圖2 健康孕婦組分形維數(shù)值Fig.2 FD value of healthy pregnant women

圖3 妊娠期疾病孕婦組分形維數(shù)值Fig.3 FD value of pathological pregnant women

根據(jù)所得測(cè)試結(jié)果，利用SPSS軟件計(jì)算各組的均值和方差，并采用 t檢驗(yàn)對(duì)組間均值進(jìn)行比較。統(tǒng)計(jì)結(jié)果表明健康孕婦的分形維數(shù)值(1.27±0.29)與對(duì)照組健康未孕婦女(1.31 ±0.41)相比降低，顯著性差異具有統(tǒng)計(jì)學(xué)意義(P＜0.05)。妊娠期疾病孕婦(1.38±0.05)與對(duì)照組健康未孕婦女(1.31±0.41)相比升高，顯著性差異具有統(tǒng)計(jì)學(xué)意義(P ＜0.05)，見(jiàn)表 1。

表1 3組受試者分形維數(shù)值的均值數(shù)據(jù)比較Tab.1 Comparison of FD average value

5 討論和結(jié)論

本研究選取3組樣本來(lái)分析各組分形維數(shù)值，數(shù)據(jù)結(jié)果顯示健康孕婦最低，妊娠期疾病孕婦最高。一般來(lái)說(shuō)，正常的心音信號(hào)包括第一心音和第二心音，無(wú)其他外加音和心雜音，而病理心音是由于心臟或者血管病變時(shí)，心音的強(qiáng)度和性質(zhì)發(fā)生改變而產(chǎn)生，從而導(dǎo)致了心音信號(hào)中出現(xiàn)了額外的心音成分，這就使得病理心音信號(hào)的復(fù)雜程度高于正常心音。

由表1可知，健康孕婦心音信號(hào)的分形維數(shù)值較健康未孕婦女降低，分形維數(shù)可以定量地描述心音信號(hào)的不規(guī)則程度和復(fù)雜程度，因此與健康未孕婦女相比，健康孕婦心音信號(hào)的復(fù)雜程度更低。這可能是由于孕婦妊娠期間心臟負(fù)擔(dān)進(jìn)行性增加，發(fā)生代償性生理變化，導(dǎo)致心臟應(yīng)激狀態(tài)適應(yīng)性降低。患有重度妊娠子癇等的妊娠期疾病孕婦，其分形維數(shù)值最高，這并不表明妊娠期疾病孕婦的心臟應(yīng)激狀態(tài)適應(yīng)性增強(qiáng)，而主要是由于心音信號(hào)中增加了病理成分，復(fù)雜程度增大。

對(duì)于健康孕婦和妊娠期疾病孕婦兩組樣本結(jié)果，其分形維數(shù)值具有統(tǒng)計(jì)學(xué)差異，妊娠期疾病孕婦的分形維數(shù)值較健康孕婦顯著升高，因此采取分形理論來(lái)分析孕婦的心音信號(hào)是可行的。將分形維數(shù)作為心音信號(hào)的特征值，其值大小可以明顯地區(qū)分健康孕婦和妊娠期疾病孕婦，所以該方法能為臨床診斷孕婦心音信號(hào)是否異常、保障孕期婦女健康提供重要依據(jù)。

孕婦的心臟安全性是產(chǎn)科醫(yī)生和婦幼保健工作的一個(gè)重要問(wèn)題。孕婦是一個(gè)特殊的群體，大多數(shù)孕婦的懷孕和分娩對(duì)心臟的影響都處在生理范圍內(nèi)，由于各種因素的影響，少數(shù)孕婦在長(zhǎng)達(dá)9個(gè)月的懷孕和分娩過(guò)程中在心臟負(fù)擔(dān)進(jìn)行性增加的作用下，會(huì)出現(xiàn)亞健康狀態(tài)，或者導(dǎo)致病理性變化。妊娠期疾病的診斷更要求具有無(wú)創(chuàng)性和安全性，本研究初步將分形理論應(yīng)用于孕婦的心音信號(hào)分析方面，保證了診斷方法的無(wú)創(chuàng)性和安全性。在后續(xù)的研究中，還需要進(jìn)一步做多中心、大樣本的調(diào)研，探索分形理論在孕婦心音信號(hào)中的應(yīng)用。

[1]He J，Gu D，Wu X，et al.Major causes of death among men and women in China[J].New England Journal of Medicine，2005，353(11)：1124-1134.

[2]陳劍，郭興明，肖守中.心音信號(hào)識(shí)別的意義及其方法的研究[J].國(guó)外醫(yī)學(xué)生物醫(yī)學(xué)工程分冊(cè)，2004，27(2)：87-89.

[3]Xiao Shouzhong，Bu Bin，Ge Ping，et al.An evaluation the merit of the method for simultaneous measuring of cardiac contractility and heart rate for athletes[J].Journal of Science and Medicine in Sport，2004;7(4)：39.

[4]王琳，郭興明，成敏，等.孕期和產(chǎn)程心臟儲(chǔ)備無(wú)創(chuàng)檢測(cè)和評(píng)估方法的初步研究[J].第三軍醫(yī)大學(xué)學(xué)報(bào)，2006，28(4)：360-362.

[5]Gabriel KM.Heart disease diagnosis and therapy：a practical approach[M].London：Williams and Wilkins，1996.

[6]Sladana S.Surrogate data test for nonlinearity of the rat cerebellar electrocorticogram in the model of brain injury [J].Signal Processing，2010，(90)：3015-3025.

[7]任焱，楊基海，尹少華，等.局部肌疲勞過(guò)程與表面肌電信號(hào)分維值變化關(guān)系的研究[J].生物醫(yī)學(xué)工程研究，2005，23(4)：215-217.

[8]Lin Chiahung，Du Yichun. Fractal QRS-complexes pattern recognition for imperative cardiac arrhythmias[J].Digital Signal Processing，2010，(10)：1274-1285.

[9]許可，孟濬.基于分形理論的非線性信號(hào)噪聲污染程度檢測(cè)[J].機(jī)電工程，2008，25(9)：4-7.

[10]Ning Xinbao， Bian Chunhua， Wang Jun， et al. Research progress in nonlinear analysis of heart electric activities[J].Chinese Science Bulletin，2006，51(4)：385-393.

[11]孫博文，張本祥.上證指數(shù)的分形結(jié)構(gòu)及盒維數(shù)測(cè)量的研究[J].哈爾濱理工大學(xué)學(xué)報(bào)，2001，6(6)：32-33.

[12]趙健，雷蕾，蒲小勤.分形理論及其在信號(hào)處理中的應(yīng)用[M].北京：清華大學(xué)出版社，2008.

[13]王堅(jiān)，田社平，張亞利，等.腸道壓力數(shù)據(jù)的盒維數(shù)分形分析[J].中國(guó)生物醫(yī)學(xué)工程學(xué)報(bào).2007，26(3)：389-393.