2392004年7月中國(guó)醫(yī)學(xué)物理學(xué)雜志J.,2004第21卷第4期Chinese Journal of Medical Physicsvol.21. No.4心音的分析何飛燕(東莞理工學(xué)校，廣東東莞523006)摘要:自上世紀(jì)70年代以來(lái),由于計(jì)算機(jī)和現(xiàn)代數(shù)字信號(hào)處理技術(shù)的發(fā)展,生物醫(yī)學(xué)工作者對(duì)心音進(jìn)行了大量的分析研究，相繼運(yùn)用了譜分析、時(shí)頻分析、小波變換和最優(yōu)匹配法等方法。譜分析不適合對(duì)非平穩(wěn)隨機(jī)心音信號(hào)的分析。國(guó)內(nèi)外工作者通過(guò)加核函數(shù)和采用模糊函數(shù)等方法來(lái)減少時(shí)頻分布中的交叉項(xiàng)。具有自適應(yīng)性時(shí)一頻窗的小波變換可以得到更能反映心音病例特征的信息。最優(yōu)匹配法是一種沒有交又項(xiàng)的時(shí)頻分析方法。關(guān)鍵詞:心音;時(shí)須分析;小波變換; MP中圖分類號(hào): R540.4文獻(xiàn)識(shí)別碼: A文章編號(hào): 1005-202(2004)04 -0239-02The Analysis on Heart SoundsHE Fei-yen(Dongguan Science and Technical School of Guang Dong Prorince, Donggruan Guangdong 523006, China)Abstract: In this paper, the siuation and the destination of the heart sound study were introduced, at the same time, some methodswhich were used to research heart sound by these biomedical researchers since the 70's on the last century were introduced, theyincluded the spectral analysis, the distribution of time frequency. the wavelet translation, the malching pursuit, and so on. The spectralanalysis is not suit for analyzing the typical non- stationary stochastic signal. These researchers on the world have reduced the Cro8s-terms on the distribution of time- frequency by nuclear function or fuzzy function. The wavelet translation with the adapive time-frequency window could gain the characters that more rllet the symploms of heart sound. The matching pursuit is the analysis oftime-frequeney without the cross-terms.Key words: the heart Bound; the distribution of time - frequency; the wavelet trunslation; the matching pursuit感區(qū)以內(nèi)的信息外,還記錄了大量10 Hz以下的信息叫。但是，前言單純心音圖是無(wú)法顯示心臟雜音的強(qiáng)度分級(jí)、雜音的性質(zhì)、頻近年來(lái),隨著物質(zhì)生活水平的提高,心血管疾病在我國(guó)星率和許多柔軟的高音調(diào)雜音。上升趨勢(shì)。有資料表明，單純冠心病發(fā)病率在我國(guó)男性增加70年代后，與其它檢查手段相比，心音圖的發(fā)展相對(duì)緩26%,女性增加19% ,并且對(duì)國(guó)民健康構(gòu)成了極大的威脅,其死慢,主要原因有:(1)心音圖是- -種典型的非平穩(wěn)隨機(jī)信號(hào),所亡率已上升為各類疾病之首。在發(fā)達(dá)國(guó)家中心血管疾病也是以不能用確定的數(shù)學(xué)函數(shù)來(lái)描述;(2)對(duì)心音的產(chǎn)生機(jī)理仍存在各種死亡原因之最。有爭(zhēng)議的說(shuō)法,這使得心音信號(hào)的特征與生理病理的結(jié)合不緊心音信號(hào)是人體最重要的生理信號(hào)之一,它含有關(guān)于心臟密,往往會(huì)有一病多癥,一癥多病的情況;(3)缺乏對(duì)心音信號(hào)各個(gè)部分(如心房、心室、大血管、心血管及各個(gè)瓣膜)功能狀態(tài)系統(tǒng)的、定量的分析;(4)沒有標(biāo)準(zhǔn)化的心音采集儀器。的大量的生理病理信息。盡管心臟聽診可以診斷許多心臟疾90年代以來(lái)，由于計(jì)算機(jī)技術(shù)和現(xiàn)代數(shù)字信號(hào)處理技術(shù)病。但是,人耳聽力敏感局限(102Hz以上)和聽診者的臨床經(jīng)驗(yàn)的發(fā)展,大大地提高了對(duì)傳統(tǒng)心音的分析能力。因此,國(guó)內(nèi)外生在很大程度上影響了聽診的準(zhǔn)確度。計(jì)算機(jī)技術(shù)和現(xiàn)代數(shù)字信物醫(yī)學(xué)界的人士對(duì)心音信號(hào)進(jìn)行了較全面的、系統(tǒng)的研究。心號(hào)處理技術(shù)的發(fā)展,為研究心音并最終更準(zhǔn)確地診斷心血管疾音信號(hào)的分析與研究主要在以下幾個(gè)方面:(1)對(duì)S,和S.的生病提供了必要的手段。.理病理研究;(2)對(duì)人工心臟瓣膜的無(wú)創(chuàng)傷檢測(cè);(3)對(duì)心音微弱成分(第三心音和第四心音)的分析研究;(4)分析心臟雜音1心音基本發(fā)展及分析方法概況,的頻率變化規(guī)律;(5)從一個(gè)心動(dòng)周期中定位提取心音成分;(6)1.1研究心音的發(fā)展過(guò)程對(duì)心音傳導(dǎo)機(jī)制建模。自上世紀(jì)60年代以來(lái),心音圖(phonocardiogram ,PCG)的出現(xiàn)在分析所用的PCG信號(hào)，是通過(guò)心音傳感器把人體心現(xiàn)曾為心血管疾病的診斷提供過(guò)幫助,心音圖除了記錄聽力敏音轉(zhuǎn)化為電信號(hào)，再通過(guò)A/D采集卡量化后記錄于計(jì)算機(jī)中的數(shù)字信號(hào),因此,心音圖有了數(shù)字心音圖之稱。運(yùn)用現(xiàn)代數(shù)字收稿日期:2003-11-28作者簡(jiǎn)介:何飛燕(1969-),女,講師,工學(xué)碩士。主要從事電子通訊和行深,言號(hào)劉|中國(guó)煤化工均生理及病理關(guān)系進(jìn)二的基礎(chǔ)研究和臨床診信息系統(tǒng)專業(yè)的教學(xué)與研究。電話: (762000-288;斷上都YHCNMHG疾病中,用現(xiàn)代數(shù)字E- -mail: lghfy@21cn.com。-240-中國(guó)醫(yī)學(xué)物理學(xué)雜志第21卷第4期_ 2004年7月信號(hào)處理技術(shù)分析心音具有很大的潛力。(1)1.2心音的分析研究方法πOTV品 n!(n+1)!自上世紀(jì)60年代以來(lái)，生物醫(yī)學(xué)界對(duì)心音進(jìn)行了大量的J是第一類一階函數(shù),而∞0為尺度因子。研究,按其所運(yùn)用的方法和研究的目的不同,主要有下面四種。(2)采用模糊函數(shù):模糊函數(shù)是雷達(dá)、通信領(lǐng)域中設(shè)計(jì)分析1.2.1譜分析信號(hào)時(shí)常用的函數(shù),它的定義尚未統(tǒng)一。鄭玲等"應(yīng)用式(2)定上世紀(jì)90年代以前，心音分析工作者主要采用FFT、AR義的模糊函數(shù)來(lái)減少心音時(shí)頻分布的交叉項(xiàng)。和ARMA等譜分析,而且其研究對(duì)象又局限在狗、鴿和豬等動(dòng)A5(6,)站L_s(+2)》"(一號(hào))“*山(2)物的心音。采用這些方法對(duì)心音從心腔到胸表的傳導(dǎo)特性進(jìn)行建模和心音的共振頻率。研究表明:心胸聲傳系統(tǒng)還是一個(gè)時(shí)時(shí)頻分析能夠反映心音的動(dòng)態(tài)變化特性,但是其時(shí)頻窗是變模型。Durand 等人田,基于信號(hào)與系統(tǒng)的TFR確立了如圖1固定的,如果其窗口長(zhǎng)度取的過(guò)長(zhǎng),盡管可提高頻率分辨率,但所示的建立了參量模型，該模型能夠反映時(shí)變的心胸聲傳系時(shí)間分辨率變得較低;反之,時(shí)間分辨率提高,頻率分辨率變充。低。圖1中x(t)表示心內(nèi)心音,n(1)表示輸出背景噪音,0()是1.2.3小波分布心胸聲傳系統(tǒng)的輸出y(l)胸表記錄的心音,h(t,r)是心胸聲傳由于小波變換中的時(shí)一頻窗具有自適應(yīng)性,所以可以得到系統(tǒng)的時(shí)變頻率響應(yīng)。更能反映心音病例特征的信息。心音信號(hào)的小波變換譜在不同尺度上的分布相應(yīng)于心臟不同部分的振動(dòng)信息。尺度較小的部x(1)h0,n]- 0()分相應(yīng)于信號(hào)的高頻成分，可以表示心臟瓣膜的振動(dòng)信息、血流本身及其對(duì)其它部分相互作用的信息。而尺度較大的部分相應(yīng)于信號(hào)的低頻部分,這是心肌與大血管等心臟的質(zhì)量較大的部分在瓣膜血流和肌肉的作用下產(chǎn)生的振動(dòng)。團(tuán)1心胸聲傳系統(tǒng)因此它在時(shí)域和頻域同時(shí)具有良好的局部化性質(zhì),從而可Fig,1 The modeling of heart-thorax acoustic system以把分析的重點(diǎn)聚焦到任意的細(xì)節(jié),被譽(yù)為數(shù)學(xué)“顯微鏡"。近無(wú)論是采用FFT,還是采用參數(shù)模型(LPC、AR和ARMA)年來(lái),國(guó)內(nèi)生物醫(yī)學(xué)工作者7應(yīng)用連續(xù)小波變換對(duì)心音的頻率的諧分析,都是基于平穩(wěn)信號(hào)的分析方法,這樣難以反映出心范圍進(jìn)行估計(jì)和對(duì)心音做三維分析。音信號(hào)的動(dòng)態(tài)變化過(guò)程。如果采用參數(shù)模型,例如LPC模型，但是在分析信號(hào)時(shí),在有些情況下,選用不同的小波基函則應(yīng)該正確選擇LPC的階次,數(shù)據(jù)長(zhǎng)度等。而且,在對(duì)心音信數(shù),會(huì)得出截然不同的結(jié)果。選取合適的小波基,會(huì)使計(jì)算出的號(hào)進(jìn)行譜分析時(shí)應(yīng)該針對(duì)不同的分析目的,采用相應(yīng)的參數(shù)選小波系 數(shù)能夠最好地反映信號(hào)的特征。取性能最好的算法。功率譜分析的另一個(gè)缺點(diǎn)是:功率譜只能12.4 MP(the matching pursuit)反應(yīng)幅頻特性,而忽視了往往蘊(yùn)含著大量信息的相頻特性。實(shí)現(xiàn)MP方法的具體過(guò)程是，首先針對(duì)與心音信號(hào)特點(diǎn)建但這并不意昧譜分析就不再用在心音的研究中,九十年代立一個(gè)有確定頻率和已知衰減系數(shù)的信號(hào)冗余庫(kù),再將實(shí)際心后,SavaHP等間就采用政進(jìn)的前后向預(yù)算的Prory方法對(duì)用機(jī)音信號(hào)與該信息庫(kù)中的每個(gè)基信號(hào)進(jìn)行匹配。首先匹配系數(shù)最械瓣膜修復(fù)后的心臟心音的譜分布進(jìn)行了研究。大的實(shí)際心音信號(hào)的第一個(gè)分量，然后將實(shí)際信號(hào)減去該分1.2.2時(shí)頻分布(TFRs)量,用余信號(hào)在冗余庫(kù)中重新匹配以獲取第二個(gè)分量,依此類由于心臟搏動(dòng)過(guò)程是一個(gè)動(dòng)態(tài)過(guò)程,因此,要獲得心臟的推,直到余量少到- -定指標(biāo)為止。依此方法得到的分量可完全血液動(dòng)力學(xué)特征,就需要研究心音信號(hào)頻率成分隨時(shí)間的變化消除交叉項(xiàng)。規(guī)律。1998年,Xuan Zhang等人采用MP(matching pursuil)對(duì)心時(shí)一頻分析是研究非平穩(wěn)隨機(jī)信號(hào)的有效工具。一般來(lái)音進(jìn)行分析和重構(gòu)。設(shè)法將聽力敏感區(qū)以外的心音轉(zhuǎn)到敏感區(qū)說(shuō),時(shí)一頻分析方法具有很強(qiáng)的能量聚集作用,不需要知道信內(nèi),以提高醫(yī)生診斷準(zhǔn)確度。王威廉等人運(yùn)用該方法對(duì)正常的號(hào)頻率隨時(shí)間的確定關(guān)系，只要信噪比足夠高,通過(guò)時(shí)一頻分第一心音和不正常的第一心音作了系統(tǒng)的分析及識(shí)別分類。此析方法就可以在時(shí)間-頻率平面上得到信號(hào)的時(shí)間頻率關(guān)系。外,還有其他分析方法來(lái)研究心音。2001年,inpping Xu等人閉但是,時(shí)一頻分布有一個(gè)嚴(yán)重的不足即存在交叉項(xiàng)。交叉應(yīng)用nonlinear tansient chip signal model 把第二心音的兩個(gè)項(xiàng)是二次型或雙線性的時(shí)一頻分布的固有結(jié)果,它們來(lái)自多分主要成分分離出來(lái)(A2和P2)。量信號(hào)中不同信號(hào)分量之間的交叉作用。交叉項(xiàng)對(duì)心音信號(hào)特綜上所述,心音研究的方法基本上是隨現(xiàn)代數(shù)字信號(hào)處理征頻率的提取將產(chǎn)生影響。技術(shù)的發(fā)展而發(fā)展。為減少交叉項(xiàng)，國(guó)內(nèi)外的生物醫(yī)學(xué)工作者作了大量工作，2總結(jié)取得了一定的研究成果。歸結(jié)起來(lái)有以下幾種:(1)加核Guo 等人叫應(yīng)用Bessel核的TFRs對(duì)第一心音和中國(guó)煤化工信號(hào)間的交叉項(xiàng)的影心臟雜音分析時(shí)可得到比較好的時(shí)-頻分布。他在時(shí)頻計(jì)算式響,以:MHC N MH G征參量。隨著數(shù)字信中加了如式(1)的函數(shù):號(hào)處理技術(shù)的提高,通過(guò)心音分析工作者的努力,心音的研究(下轉(zhuǎn)第243頁(yè))中國(guó)醫(yī)學(xué)物理學(xué)雜志第21卷第4期2004年7月2.3工作原理染色標(biāo)本數(shù)較多,可達(dá)10片左右;④操作步驟簡(jiǎn)單,染色質(zhì)量首先,根據(jù)標(biāo)本染色的需要,在儀器的染色程序編程器上清晰;⑤獨(dú)特的廢液回收系統(tǒng)具有綠色環(huán)保功能符合國(guó)際潮設(shè)定其染色自動(dòng)程序(其中包括染液種類的選定、沖洗時(shí)間的流;設(shè)定、環(huán)境溫度的選擇以及染色步驟的設(shè)計(jì)等)。另外,儀器提(3)玻片標(biāo)本自動(dòng)染色機(jī)是一-種適用于生物學(xué)、微生物學(xué)、供了兩種標(biāo)準(zhǔn)程序:標(biāo)準(zhǔn)程序1和標(biāo)準(zhǔn)程序2。免疫學(xué)病理學(xué)以及臨床病理輔助診斷等學(xué)科的標(biāo)本切片染色具體的工作方式為:在設(shè)定的染色自動(dòng)程序的控制下,步的全新自動(dòng)化儀器設(shè)備。特別是采用單片機(jī)編程技術(shù)來(lái)實(shí)現(xiàn)進(jìn)電動(dòng)機(jī)首先驅(qū)動(dòng)放置在托架導(dǎo)軌上的玻片標(biāo)本托架作水平玻片標(biāo)本的自動(dòng)染色、沖洗功能、自動(dòng)控時(shí)和控溫等功能,是現(xiàn)運(yùn)動(dòng),使其托架.上的玻片標(biāo)本正對(duì)著所需染液的滴管;然后,在代醫(yī)學(xué)實(shí)驗(yàn)設(shè)備發(fā)展的必然趨勢(shì);繼電器J-Js的控制下完成其托架上所有玻片標(biāo)本的第一種染(4)該產(chǎn)品可用于醫(yī)學(xué)生物標(biāo)本的革蘭氏染色法、HE染色液的滴染操作過(guò)程;經(jīng)過(guò)一定的染色時(shí)間后 ,啟動(dòng)繼電器J。進(jìn)法以及相關(guān)的染色法,而且該自動(dòng)染色機(jī)已經(jīng)獲得了國(guó)家知識(shí)行玻片標(biāo)本的清水沖洗，從而完成標(biāo)本的第-種染液的染色;產(chǎn) 權(quán)局頒發(fā)的專利證書(專利號(hào):ZL 02 2400.7)。再經(jīng)過(guò)一定的時(shí)間間隔,儀器重復(fù)上述步驟,從而完成標(biāo)本的第二、第三.--、第n種染液的染色,直至玻片標(biāo)本的全部染色。參考文獻(xiàn):同時(shí)，顯示器將用數(shù)字和代碼的形式實(shí)時(shí)顯示染色的全部過(guò)[]傅揚(yáng)烈.單片機(jī)原理與應(yīng)用教程[M]北京:電子工業(yè)出版社，2002,圣。[2]姜寶港.智能家用電器原理與維修[M].北京:機(jī)械工業(yè)出版社，3結(jié)論2002,5.[3]甘平,等用電荷耦合提高pH計(jì)測(cè)最準(zhǔn)確度[]傳感器技術(shù), 2001,(1)據(jù)調(diào)查目前我國(guó)大部分實(shí)驗(yàn)室、醫(yī)院以及醫(yī)學(xué)科研所20(9):49-51.對(duì)玻片標(biāo)本的染色都主要是采用手工操作,普遍存在著染色質(zhì)[4] Berelowits J T, Poluta M A J. Design and implementation of PC-based system for acquisition of a neonatal vital signs database[].量不高、工序繁瑣、耗費(fèi)時(shí)間和染液等缺點(diǎn)。所以,國(guó)內(nèi)急需適Autonedica, 194,161):97-107.合中國(guó)實(shí)際國(guó)情的醫(yī)用玻片標(biāo)本自動(dòng)染色機(jī)。據(jù)《中國(guó)醫(yī)療器[s]甘平,等心電模擬信號(hào)發(fā)生器的設(shè)計(jì)與制作[J],重慶大學(xué)學(xué)報(bào)(自械產(chǎn)業(yè)年鑒>提供的資料顯示,國(guó)外醫(yī)用染色機(jī)的售價(jià)普遍都然科學(xué) 版), 2001,24:129-130.在25萬(wàn)元人民幣左右，國(guó)內(nèi)廠家生產(chǎn)的組織染色機(jī)的售價(jià)也[6]田源,等.誘發(fā)電位儀電極接觸電阻檢刮系統(tǒng)的研制方法[辦中國(guó)醫(yī)學(xué)物理學(xué)雜志, 203.202)115116.6在4萬(wàn)元左右。而本自動(dòng)染色機(jī)的價(jià)格大約只有國(guó)內(nèi)組織染色[7]黃力字,等.基于單片機(jī)的血壓監(jiān)護(hù)儀的研制[J].中國(guó)醫(yī)療器被雜機(jī)價(jià)格的1/20,而且國(guó)內(nèi)目前還未發(fā)現(xiàn)有此同類產(chǎn)品。因此，志，20264):253-255.在市場(chǎng)上的競(jìng)爭(zhēng)力是顯易而見的:[8]杜曉蘭,等.呼吸、體溫兩參數(shù)檢測(cè)電路的設(shè)計(jì)I生物醫(yī)學(xué)工程雜(2)該全自動(dòng)染色機(jī)具有以下特點(diǎn):①全自動(dòng)程序化操作: .志, 2001,18(4):538~540.②體積小、重量輕、易于搬動(dòng)和產(chǎn)品運(yùn)行環(huán)境容易滿足;③-次[9]湯黎明,等.精密神經(jīng)閾值刺激儀的研制小醫(yī)療衛(wèi)生裝備，2001,22(6):12~14.(上接第240頁(yè))將取得更多的成果。.[5] Wang WL, Guo Z, Yang J, et al. analysis of the first heart soundusing the matching pursuit method[J]. Med. Bio. Eng, Comput,2001,396):644- 648.[6)] Sava HP， Edward IT, MoDonell. Spectral compsition of heart[1]鄭玲,宋長(zhǎng)宏高小格,等采用模棚函數(shù)消除維格納分布的交叉項(xiàng)及sounds before and after mechanical heart valve implantation using其在心音信號(hào)分析中的應(yīng)用[I.航天醫(yī)學(xué)與醫(yī)學(xué)工程,1995,8(1):a roodified forward -backward prony's meathod [J]. IEE Trans.15-22.Bioned. Eng, 196.437):734.[2] Durand LG, Pibard P. Digial Sigal Proessing of the Ponocar[7]李橋，趙玲,那慶余等應(yīng)用小波變換進(jìn)行心音三維時(shí)顎分析的研究diogum:Review of the Most Recent Advancements[], Critical Re[JI.中國(guó)醫(yī)學(xué)物理學(xué)雜志0110-11221views in Biomedical Engneg199.23/4):163-219.[8]胡曉,王威廉.基于小波變換第一心音頻率檢測(cè)川生物醫(yī)學(xué)工程學(xué)[3] Zhang X Durand LG. Time -frequency scaling tnsformation of雜志2003.(3).the phonoeardiogram based of the matching pursuit method [J].9] 胡曉,王威廉.基于小波變換心音定位用生物醫(yī)學(xué)工程與臨床, .IEEE Trans. Biomed. Eng1998.458>92-979.2002./4)[4] Xu JP, Durand LG, Pibern P. Extnction of the Aorie and Pul-中國(guó)煤化工ime-frequey dsibuiosmonary Components of the seond Heart Sound Using a Nonlinersel Kermel[] IEEE Trans.Transient Chirp Signal Model[J]. IEEE Trans Biomed. Eng, 2001,YHCNMHG48(3)277-283.