基于非平穩(wěn)系統(tǒng)辨識(shí)的心音包絡(luò)自適應(yīng)分割

許春冬,周 靜,應(yīng)冬文,2,龍清華

(1.江西理工大學(xué) 信息工程學(xué)院,江西 贛州 341000;2.中國(guó)科學(xué)院聲學(xué)研究所 語(yǔ)言聲學(xué)與內(nèi)容理解重點(diǎn)實(shí)驗(yàn)室,北京 100190)

0 概述

目前,心血管疾病已成為危害人類(lèi)生命健康的嚴(yán)重疾病,心血管疾病死亡人數(shù)在全球疾病死亡人數(shù)中占據(jù)較大比例[1]。心音信號(hào)(Heart Sound Signal,HSS)分析是心血管疾病常用的輔助診斷方法之一,主要通過(guò)計(jì)算機(jī)輔助診斷分析系統(tǒng)的形式進(jìn)行應(yīng)用。心音信號(hào)的輔助診斷系統(tǒng)流程可劃分為分割與分類(lèi)兩個(gè)主要步驟[2-3],其中分割是心音信號(hào)分析與診斷的基礎(chǔ)與前提,分割結(jié)果對(duì)整個(gè)系統(tǒng)的分析結(jié)果具有重要影響。然而,心音信號(hào)相當(dāng)敏感,極易受到儀器噪聲、環(huán)境噪聲、肺音、呼吸道音、腸鳴音及檢測(cè)者自身活動(dòng)狀態(tài)等因素的影響。此外,各類(lèi)病理原因也會(huì)導(dǎo)致心音中出現(xiàn)心雜音[2-3]。這些干擾因素的存在使得分割任務(wù)變得十分困難,且實(shí)際中采集到的心音信號(hào)通常包含多種噪聲,導(dǎo)致其可分析性急劇下降,使得原本較弱的心音信號(hào)變成一個(gè)較復(fù)雜的聲學(xué)信號(hào)。因此,心音信號(hào)的有效分割逐漸成為目前研究的熱點(diǎn)之一。

隨著心音信號(hào)研究的深入,各種分割方法相繼被提出,主流的心音信號(hào)分割方法可分為[4-6]基于包絡(luò)的方法、基于特征提取的方法、基于隱馬爾科夫模型(Hidden Markov Model,HMM)的方法和基于機(jī)器學(xué)習(xí)的方法4類(lèi)。其中,基于包絡(luò)的方法因算法簡(jiǎn)單、實(shí)時(shí)性強(qiáng)等特點(diǎn)[6-7]得到廣泛研究與應(yīng)用,特別是在實(shí)時(shí)檢測(cè)需求下成為首選方法[8],而其他3類(lèi)分割方法的實(shí)時(shí)性相對(duì)較差[5,9-10]。本文借鑒包絡(luò)分割方法的思想,提出一種基于非平穩(wěn)性系統(tǒng)辨識(shí)的心音包絡(luò)提取方法,并將其應(yīng)用于心音閾值自適應(yīng)分割中。

1 基于包絡(luò)的分割方法

心音信號(hào)分割的首要任務(wù)是定位分割出基礎(chǔ)心音(Fundamental Heart Sound,FHS)信號(hào)s1(第一心音)和s2(第二心音),然后分割出收縮期(sys)和舒張期(dia),如圖1所示。基于包絡(luò)的分割方法主要包括預(yù)處理、包絡(luò)提取、閾值分割3個(gè)核心部分。

圖1 心音時(shí)域波形圖Fig.1 Oscillogram of heart sound time domain

預(yù)處理主要包括預(yù)加重、消除趨勢(shì)項(xiàng)和降噪,其中:預(yù)加重為高頻提升,通過(guò)一階FIR濾波器即可完成[2];對(duì)于趨勢(shì)項(xiàng),可采用最小二乘法擬合進(jìn)行消除[11];降噪中應(yīng)用較普遍的方法為基于小波變換的降噪方法,通過(guò)丟棄帶外小波系數(shù)、濾波帶內(nèi)小波系數(shù),實(shí)現(xiàn)噪聲抑制[12]。預(yù)處理后的心音信號(hào)特征更清晰,為包絡(luò)提取奠定了基礎(chǔ)。

心音信號(hào)包絡(luò)的提取方法主要包括香農(nóng)能量包絡(luò)法[3,5-6]、維奧拉積分包絡(luò)法[5]、希爾伯特變換包絡(luò)法[7,9]、樣條插值包絡(luò)法[6]以及經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition,EMD)包絡(luò)法[7,10]等。香農(nóng)能量包絡(luò)法在噪聲干擾較低的情況下能夠取得較好的效果,但若降噪后的噪聲干擾依舊較大,則提取的能量包絡(luò)特征會(huì)受到嚴(yán)重干擾。維奧拉積分包絡(luò)法能夠較好地提取出s1的特征,但對(duì)衰弱s2及心雜音的干擾信號(hào)包絡(luò)提取效果明顯下降。在希爾伯特變換法中,當(dāng)屬黃變換研究最深入,但該方法在通過(guò)極值點(diǎn)擬合包絡(luò)時(shí)通常存在一些固有缺陷,尤其是對(duì)于一些異常心音信號(hào)。樣條插值包絡(luò)法通過(guò)對(duì)極大值和極小值進(jìn)行插值獲得上包絡(luò)和下包絡(luò),可以直觀(guān)反映信號(hào)的幅值,但存在欠包絡(luò)的問(wèn)題,且實(shí)際中欠包絡(luò)問(wèn)題處理結(jié)果并不理想。EMD包絡(luò)法無(wú)需預(yù)先設(shè)定基函數(shù),自適應(yīng)能力較強(qiáng),但在處理非平穩(wěn)信號(hào)時(shí)存在模態(tài)混疊及信號(hào)失真等問(wèn)題。

分割過(guò)程主要包括閾值函數(shù)設(shè)計(jì)與分割成分判定兩部分。其中,閾值函數(shù)多數(shù)為雙閾值函數(shù)[7];任何一個(gè)心動(dòng)周期中均包含第一心音、收縮期、第二心音、舒張期4個(gè)部分,因此通過(guò)閾值分割判定基礎(chǔ)心音后,根據(jù)時(shí)域特征即可區(qū)分出收縮期和舒張期[6-7,9]。

2 非平穩(wěn)性系統(tǒng)辨識(shí)原理

系統(tǒng)辨識(shí)主要是根據(jù)輸入與輸出時(shí)間序列來(lái)描述系統(tǒng)行為的過(guò)程,用于估計(jì)與建立控制領(lǐng)域內(nèi)的系統(tǒng)模型[13],其采用的非平穩(wěn)性系統(tǒng)辨識(shí)模型衍生于波束形成等信號(hào)處理領(lǐng)域[13-14]。文獻(xiàn)[14]提出一種估計(jì)傳遞函數(shù)(Transfer Functions,TFs)比率的方法,用于衡量不同信道間的傳遞差異。首先,假設(shè)傳遞函數(shù)比所期望信號(hào)變化的更加緩慢,進(jìn)一步假設(shè)噪聲信號(hào)與期望信號(hào)的幅值比是緩慢變化,將時(shí)間軸劃分為一系列分析間隔,在每一個(gè)分析間隔期間內(nèi)假定TFs和噪聲信號(hào)均近似為平穩(wěn)。然后,將該分析時(shí)間間隔劃分為幀,使得期望信號(hào)在每幀內(nèi)看作不變。以第K幀的信號(hào)為例,得到:

(1)

(2)

(3)

(4)

由于Um(t,ejw)和Z1(t,ejw)一般是相關(guān)的,因此無(wú)法通過(guò)計(jì)算式(4)來(lái)直接獲得無(wú)偏估計(jì)。文獻(xiàn)[15]指出采用最小二乘法可求得超定方程組Hm(ejw)(m=2,3,…,M為一組獨(dú)立方程)的無(wú)偏估計(jì)。方程式(5)的解為方程式(6)。

(5)

(6)

其中:〈〉表示取均值運(yùn)算;Hm(ejw)為傳遞函數(shù)比率的估計(jì)結(jié)果,可用于衡量?jī)蓚€(gè)通道信號(hào)間的相關(guān)性,若兩個(gè)通道信號(hào)相關(guān)性強(qiáng),則Hm(ejw)值較大;若兩個(gè)通道信號(hào)相關(guān)性弱,則Hm(ejw)值較小。

3 基于非平穩(wěn)系統(tǒng)辨識(shí)的心音分割

3.1 基于系統(tǒng)辨識(shí)的包絡(luò)提取

雖然心音信號(hào)近似為準(zhǔn)周期性信號(hào),但其本質(zhì)上是一類(lèi)非平穩(wěn)信號(hào),具備混沌特性[16]。為較好地辨識(shí)出基礎(chǔ)心音信號(hào),需要提取更有效的特征包絡(luò),突出基礎(chǔ)心音特征,因此本文提出一種基于非平穩(wěn)系統(tǒng)辨識(shí)(Non-Stationary System Identification,NSSI)理論的包絡(luò)提取方法。

心音信號(hào)具有短時(shí)平穩(wěn)性,因此將時(shí)間段劃分得足夠小時(shí),相鄰兩個(gè)時(shí)間段內(nèi)的心音信號(hào)具有較大的相似性。而對(duì)噪聲和部分雜音而言,其非平穩(wěn)性特征更加明顯,在適當(dāng)長(zhǎng)度的時(shí)間段內(nèi),其相似性更低。根據(jù)非平穩(wěn)系統(tǒng)辨識(shí)原理,將原始心音段作為參考通道信號(hào),將后移一個(gè)時(shí)間段的信號(hào)作為另一個(gè)對(duì)比通道,通過(guò)非平穩(wěn)系統(tǒng)辨識(shí)來(lái)求解其相關(guān)性,從而提取得到心音信號(hào)的特征包絡(luò),具體步驟如下:

1)對(duì)采集的心音信號(hào)做分幀處理。在此過(guò)程中,幀長(zhǎng)wlen的選擇較為關(guān)鍵,其關(guān)系著心音與噪聲包絡(luò)特征的區(qū)分程度。文獻(xiàn)[17]指出,當(dāng)心音信號(hào)為10 ms～30 ms時(shí)具有短時(shí)平穩(wěn)性,因此選取的幀長(zhǎng)應(yīng)在此范圍內(nèi)。此外,幀移的選擇也需注意,幀移過(guò)小會(huì)加大計(jì)算量,幀移過(guò)大會(huì)使部分對(duì)應(yīng)幀間的相關(guān)度偏低,使得特征包絡(luò)不明顯,實(shí)際處理中通常取其為幀長(zhǎng)的一半[2,10]。本文按式(6)測(cè)試了多種類(lèi)型心音信號(hào)下的相關(guān)度,并給出部分結(jié)果,如圖2所示。由于噪聲及心音信號(hào)都具備非平穩(wěn)性,且噪聲具有隨機(jī)性,心音具有一定的隨機(jī)變異性,因此圖2中相關(guān)度值波動(dòng)較大屬于正常現(xiàn)象。在選擇幀長(zhǎng)時(shí),要選擇s1和s2相關(guān)度值差距較小且能明顯區(qū)分于噪聲的幀長(zhǎng)。圖2中用矩形框標(biāo)記了各情況下的可取幀長(zhǎng),為方便起見(jiàn),選取普適性更強(qiáng)的16 ms幀長(zhǎng)。

圖2 不同病況與幀長(zhǎng)下的相關(guān)度測(cè)試結(jié)果Fig.2 Test results of the correlation at different patient conditions and frame lengths

2)求取參考信號(hào)x1與對(duì)比信號(hào)x2。將分幀處理后的信號(hào)去掉最后一幀,得到參考信號(hào)x1;將分幀處理后的信號(hào)去掉第一幀,得到對(duì)比信號(hào)x2。

(7)

(8)

(9)

(10)

其中,xcorr表示求相關(guān),FFT表示快速傅里葉變換。

4)按照式(11)給出的非平穩(wěn)系統(tǒng)相關(guān)性辨識(shí)結(jié)果求得包絡(luò)特征Hx。

(11)

5)每幀的信號(hào)Hx值個(gè)數(shù)與幀長(zhǎng)相同,因此需要選擇一個(gè)代表性的值來(lái)表示該幀的包絡(luò)特征值,本文按照常規(guī)方法取升序排列50%的值,即取中值。由篩選的中值組成一組特征值構(gòu)成提取包絡(luò)。特征包絡(luò)Hx能夠較好地提取基礎(chǔ)心音信號(hào)特征,區(qū)分出基礎(chǔ)心音、噪聲和雜音。由于特征包絡(luò)的整體呈脈沖狀,若不做平滑處理,則會(huì)出現(xiàn)錯(cuò)誤分割現(xiàn)象,因此需要對(duì)包絡(luò)作進(jìn)一步處理。本文首先采用最值濾波來(lái)平滑特征包絡(luò)Hx,即取窗內(nèi)最大值為平滑值,平滑窗長(zhǎng)取為4(根據(jù)幀長(zhǎng)確定);然后采用鋸齒展寬處理,得到Hx_s特征包絡(luò)。

如圖3所示,本文將降噪后的心音信號(hào)圖與所提取的包絡(luò)圖進(jìn)行對(duì)比,并用矩形框分別標(biāo)記時(shí)域干擾音及其對(duì)應(yīng)的特征包絡(luò)。由圖3可以看出,提取的包絡(luò)能較好地體現(xiàn)基礎(chǔ)心音特征并湮沒(méi)噪聲,而香農(nóng)能量包絡(luò)法和希爾伯特變換包絡(luò)法卻無(wú)法湮沒(méi)難以抑制的噪聲,在多組心音信號(hào)測(cè)試中均存在此類(lèi)現(xiàn)象。由于良好的特征包絡(luò)將有助于分割出基礎(chǔ)心音信號(hào),因此包絡(luò)的有效性與普適性將通過(guò)分割結(jié)果進(jìn)一步驗(yàn)證。

圖3 特征包絡(luò)提取結(jié)果Fig.3 Extraction results of feature envelopes

3.2 自適應(yīng)閾值及分割點(diǎn)判定

在提取出包絡(luò)的基礎(chǔ)上進(jìn)行閾值選取分析。為避免人工設(shè)置閾值的不便,本文采用提取相關(guān)閾值參數(shù)來(lái)自適應(yīng)分配閾值。一般地,閾值是有用信號(hào)與噪聲的分界線(xiàn),因此所以首先選擇能夠反映信號(hào)中噪聲程度的參數(shù),其次閾值被用于檢測(cè)特征包絡(luò)是否滿(mǎn)足要求,所以閾值還應(yīng)與提取的特征包絡(luò)相關(guān)。此外,由于類(lèi)似基礎(chǔ)心音信號(hào)的雜音干擾或受試者心音異常,所以分割結(jié)果還需進(jìn)一步修正,以保證分割結(jié)果的正確性。整個(gè)過(guò)程具體如下:

1)選擇反映噪聲程度的參數(shù)——信噪比(Signal to Noise Ratio,SNR),但由于實(shí)際中的噪聲或者純凈信號(hào),而通常是不可預(yù)知的,因此本文采用降噪的方式來(lái)估計(jì)純凈信號(hào),獲取SNR。在包絡(luò)提取時(shí)已經(jīng)進(jìn)行小波降噪,但并不能將噪聲完全濾除,因此進(jìn)一步加大濾波參數(shù)為cora80%來(lái)抑制噪聲估計(jì)純凈心音,按式(12)計(jì)算SNR[12]。SNR估計(jì)值將用于提取閾值T1。

(12)

其中,s(n)為估計(jì)的純凈信號(hào),d(n)為噪聲,N為信號(hào)長(zhǎng)度。

(13)

其中,NH為包絡(luò)長(zhǎng)度。

3)按照選擇的參數(shù)建立閾值函數(shù)。在包絡(luò)分割中,常用的閾值函數(shù)形式為雙門(mén)限閾值,建立雙閾值函數(shù)如式(14)和式(15)所示:

(14)

(15)

若SNR為0,則:

(16)

其中,α和β為常數(shù),文獻(xiàn)[18]指出閾值系數(shù)可通過(guò)遺傳算法進(jìn)行全局優(yōu)化,實(shí)驗(yàn)得到α=1.4、β=1.6為最佳值。

4)根據(jù)閾值尋找分割點(diǎn)。分割點(diǎn)分為上分割點(diǎn)與下分割點(diǎn)[19],即沿時(shí)間軸取包絡(luò)與閾值T1的交點(diǎn)ui和包絡(luò)與閾值T2的交點(diǎn)rl的中點(diǎn)為分割點(diǎn),當(dāng)ui>rl時(shí)為下分割點(diǎn),當(dāng)ui

5)計(jì)算每組分割點(diǎn)間包絡(luò)所對(duì)應(yīng)的面積,對(duì)分割點(diǎn)進(jìn)行第一次篩檢[5,9]。由于所提取的基礎(chǔ)心音信號(hào)分割點(diǎn)間包絡(luò)所圍的面積要遠(yuǎn)大于噪聲信號(hào)分割點(diǎn)間的面積,因此可舍去所圍面積過(guò)小的分割點(diǎn),排除噪聲包絡(luò)的干擾。

6)通過(guò)時(shí)域檢查波峰對(duì)分割點(diǎn)進(jìn)行第二次篩檢。提取每組分割點(diǎn)間的最高波峰,沿時(shí)間軸可構(gòu)成一組波峰時(shí)間點(diǎn){pi|p1,)p2,…,pm},若波峰間的時(shí)間差滿(mǎn)足Δpi=pi+1-pi<0.1 s,則認(rèn)為pi+1或pi所對(duì)應(yīng)的分割點(diǎn)為錯(cuò)誤分割點(diǎn)。如果采集的受試者心音存在類(lèi)似基礎(chǔ)心音的心雜音、第三心音、第四心音,或者是奔馬律等異常類(lèi)型心音,又或者是首位存在不完整的基礎(chǔ)心音信號(hào),都極有可能導(dǎo)致分割點(diǎn)出錯(cuò)。而由于目前的包絡(luò)分割方法多數(shù)未考慮此類(lèi)情況,因此在實(shí)際應(yīng)用中容易加大分割誤差。本文根據(jù)心音信號(hào)的時(shí)域特征,提出相應(yīng)的誤分割點(diǎn)檢驗(yàn)方法,具體步驟如下:

(1)識(shí)別不完整的基礎(chǔ)心音,去除單個(gè)分割點(diǎn)。由于用于分割的心音信號(hào)開(kāi)始處和結(jié)束處可能存在不完整的基礎(chǔ)心音,會(huì)造成分割點(diǎn)判定結(jié)果存在單個(gè)的上分割點(diǎn)或下分割點(diǎn),單個(gè)分割點(diǎn)會(huì)影響υi和νl匹配,同時(shí)會(huì)影響起始基礎(chǔ)心音信號(hào)的識(shí)別,因此需要檢測(cè)出單個(gè)的υi和νl,并將其從檢測(cè)結(jié)果中去除。在不完整的基礎(chǔ)心音檢測(cè)中,將第一個(gè)上分割點(diǎn)與第一個(gè)下分割點(diǎn)進(jìn)行比較,若第一個(gè)上分割點(diǎn)位置在第一個(gè)下分割點(diǎn)位置的后面,則判定ν1為不完整的下分割點(diǎn),并利用類(lèi)似方法檢測(cè)末尾的分割點(diǎn),整個(gè)非完整基礎(chǔ)心音檢測(cè)及濾除過(guò)程可表示為:

(17)

其中,Qdel表示需要?jiǎng)h除的分割點(diǎn),υend表示最后一個(gè)上分割點(diǎn),νend表示最后一個(gè)下分割點(diǎn)。

(2)檢測(cè)心雜音或異常額外心音的分割點(diǎn)。實(shí)際中采集到的很多心雜音或異常額外心音與基礎(chǔ)心音非常類(lèi)似,這些非基礎(chǔ)心音成分會(huì)對(duì)分割點(diǎn)的判定造成干擾。文獻(xiàn)[20]指出收縮期時(shí)長(zhǎng)一般為160 ms～240 ms,在檢測(cè)誤檢分割點(diǎn)時(shí),取收縮期時(shí)長(zhǎng)最小值的一半(80 ms),由于基礎(chǔ)心音信號(hào)包絡(luò)值較非基礎(chǔ)心音信號(hào)包絡(luò)值更大,因此可采用如下判別方法:

(18)

當(dāng)誤檢點(diǎn)出現(xiàn)在舒張期的中間位置時(shí),上述檢測(cè)方法難以有效判定,因此需要進(jìn)一步檢測(cè)。文獻(xiàn)[20]指出心動(dòng)周期時(shí)長(zhǎng)一般為600 ms～1 000 ms,取心動(dòng)周期最短時(shí)長(zhǎng)作為又一判定參量,采用如下檢測(cè)方法:

(19)

7)提取時(shí)域特征,識(shí)別心音成分。舒張期持續(xù)時(shí)長(zhǎng)一般大于收縮期時(shí)長(zhǎng)[18-19],故根據(jù)時(shí)域時(shí)長(zhǎng)判定持續(xù)時(shí)間較長(zhǎng)的成分為舒張期,較短的成分為收縮期[20-21],按順序處于舒張期與收縮期之間的成分為s1,則剩余成分為s2。

3.3 分割流程及評(píng)價(jià)指標(biāo)

根據(jù)上文所述方法完成包絡(luò)提取及分割任務(wù),具體流程如圖4所示。

圖4 自適應(yīng)閾值分割流程Fig.4 Procedure of adaptive threshold segmentation

由于本文方法主要面向復(fù)雜心音信號(hào)實(shí)時(shí)分割處理,難以通過(guò)直觀(guān)的分割定位結(jié)果對(duì)分割算法進(jìn)行有效性檢驗(yàn),因此本文采用文獻(xiàn)[22]提出的4類(lèi)評(píng)價(jià)指標(biāo)作為檢驗(yàn)指標(biāo):

(20)

(21)

(22)

PS=run_time

(23)

其中,TPFHS為基礎(chǔ)心音信號(hào)檢出正確率,FPFHS為非基礎(chǔ)心音信號(hào)檢入錯(cuò)誤率,MPFHS為基礎(chǔ)心音未檢出率,TNFHS為正確檢入的基礎(chǔ)心音信號(hào)幀數(shù),FNFHS為錯(cuò)誤檢入的非基礎(chǔ)心音信號(hào)幀數(shù),MNFHS為未檢入的基礎(chǔ)心音信號(hào)幀數(shù),PS為算法運(yùn)行時(shí)間。TPFHS越高,FPFHS和TNFHS越低,代表分割精度越高;PS越小,代表實(shí)時(shí)性越好。

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

4.1 實(shí)驗(yàn)數(shù)據(jù)及其處理

實(shí)驗(yàn)數(shù)據(jù)由中國(guó)醫(yī)學(xué)科學(xué)院阜外醫(yī)院和南京郵電大學(xué)電子科學(xué)與工程學(xué)院提供。其中選用前面提供心音62條,后者提供心音43條,所有心音信號(hào)均降采樣至2 000 Hz。為模擬復(fù)雜環(huán)境,選用noiseX-92噪聲數(shù)據(jù)庫(kù)中的15類(lèi)不同領(lǐng)域的噪聲按5 dB、10 dB、15 dB這3種信噪比加入心音信號(hào)中構(gòu)成待處理復(fù)雜心音,噪聲疊加不超過(guò)兩種,選用的心音信號(hào)中還包括17類(lèi)異常心音信號(hào)。noiseX-92噪聲數(shù)據(jù)庫(kù)中的噪聲類(lèi)型具體為白噪聲、粉紅噪聲、高頻通道噪聲、飛機(jī)駕駛艙噪聲1、飛機(jī)駕駛艙噪聲2、廠(chǎng)房噪聲1、廠(chǎng)房噪聲2、類(lèi)語(yǔ)音噪聲、坦克噪聲、驅(qū)逐艦操作室噪音、驅(qū)逐艦機(jī)艙噪音、F16機(jī)艙噪聲、軍用車(chē)輛噪聲、機(jī)槍噪音和車(chē)內(nèi)噪聲。心音信號(hào)類(lèi)型及數(shù)量如表1所示。

表1 心音信號(hào)類(lèi)型及數(shù)量Table 1 Type and quantity of heart sound signal

105條心音信號(hào)經(jīng)處理后生成33 075條帶噪心音信號(hào),共計(jì)516 285個(gè)心動(dòng)周期,原始心音信號(hào)均在Cool Edit Pro下進(jìn)行手工標(biāo)注分割點(diǎn)。實(shí)驗(yàn)軟件為Matlab2017a,計(jì)算機(jī)CPU為Intel?CoreTMi7-8700k@3.70 GHz,RAM為16.0 GB,硬盤(pán)為1 TB。

4.2 結(jié)果分析

本文方法分割結(jié)果如圖5所示,直觀(guān)反映了其能夠較精準(zhǔn)地分割出基礎(chǔ)心音信號(hào),且不會(huì)將非基礎(chǔ)心音誤分割為基礎(chǔ)心音信號(hào)。本文方法與維奧拉積分包絡(luò)分割法[5]、改進(jìn)型希爾伯特-黃變換包絡(luò)雙閾值分割法[7]、基于心動(dòng)周期估計(jì)的心音包絡(luò)分割法[9]、基于周期與波峰的自適應(yīng)閾值包絡(luò)分割法[10]和基于個(gè)性化高斯混合建模的包絡(luò)分割法[4]進(jìn)行對(duì)比,其中基于個(gè)性化高斯混合建模的包絡(luò)分割法采用在線(xiàn)訓(xùn)練方式,無(wú)需通過(guò)不同的訓(xùn)練集進(jìn)行訓(xùn)練,能夠避免對(duì)比差異。

圖5 基于本文方法的心音信號(hào)分割結(jié)果Fig.5 Segmentation results of heart sound signal based on the proposed method

表2給出了6種分割方法的評(píng)價(jià)指標(biāo)結(jié)果,可以看出,基于個(gè)性化高斯混合建模的分割法檢出正確率TPFHS最高為91.49%,其次為本文方法的89.21%,其余方法的檢出準(zhǔn)確率相對(duì)較低,這是因?yàn)楸疚姆椒杀苊怆s音和誤檢干擾。此外,可以發(fā)現(xiàn)其他4種基于包絡(luò)特征的分割方法中基礎(chǔ)心音錯(cuò)誤檢入率FPFHS要高于未檢入率MPFHS,主要原因?yàn)樵诨A(chǔ)心音末尾處噪聲存在拖尾,容易使非基礎(chǔ)心音的噪聲部分被誤判為基礎(chǔ)心音幀,而本文方法不存在此類(lèi)問(wèn)題。

表2 6種方法的分割性能比較Table 2 Comparison of segmentation performance of six methods

另外,在算法耗時(shí)PS指標(biāo)評(píng)測(cè)中,采用的心音信號(hào)段為4.5 s,即每5個(gè)心動(dòng)周期為一段,獲取PS的平均值為最終結(jié)果。從表2可以看出,基于個(gè)性化高斯混合建模的包絡(luò)分割法的耗時(shí)為6.573 5 s,遠(yuǎn)高于其他分割方法,這是因?yàn)榛趥€(gè)性化高斯混合建模的包絡(luò)分割法采用在線(xiàn)訓(xùn)練方式,涉及到高斯混合建模、前向后向運(yùn)算、Viterbi解碼、誤差修正等計(jì)算[4],復(fù)雜度遠(yuǎn)高于其他基于包絡(luò)的雙閾值分割方法;基于改進(jìn)型希爾伯特-黃變換包絡(luò)的分割法先進(jìn)行經(jīng)驗(yàn)?zāi)B(tài)分解,然后做3次貝塞爾樣條插值運(yùn)算,同時(shí)對(duì)端點(diǎn)效應(yīng)進(jìn)行優(yōu)化[7],其運(yùn)算復(fù)雜度次于基于個(gè)性化高斯混合建模的包絡(luò)分割法;維奧拉積分包絡(luò)分割法在確定尺度與基礎(chǔ)心音的關(guān)系后直接提取包絡(luò)進(jìn)行分割[5],基于周期與波峰的自適應(yīng)閾值包絡(luò)分割法建立自相關(guān)估計(jì)周期與波峰的關(guān)系閾值后完成分割[10],本文方法在提取心音包絡(luò)后直接進(jìn)行分割與判別,復(fù)雜度均遠(yuǎn)低于前兩種方法。本文方法的耗時(shí)為0.105 2 s,能夠滿(mǎn)足實(shí)時(shí)分割要求,可為實(shí)時(shí)監(jiān)測(cè)的心音自動(dòng)分析技術(shù)提供參考。

5 結(jié)束語(yǔ)

基于包絡(luò)的心音信號(hào)分割方法因算法簡(jiǎn)單、時(shí)效性強(qiáng)等特點(diǎn)而得到廣泛應(yīng)用,其中的特征包絡(luò)提取與閾值處理算法設(shè)計(jì)為關(guān)鍵步驟。本文提出一種基于非平穩(wěn)系統(tǒng)辨識(shí)的心音特征包絡(luò)提取方法,采用遺傳算法設(shè)計(jì)自適應(yīng)雙閾值函數(shù)實(shí)現(xiàn)閾值分割。在近似基礎(chǔ)心音的心雜音處理上,根據(jù)時(shí)域特征進(jìn)行分割點(diǎn)篩查,識(shí)別并去除誤檢分割點(diǎn)和單獨(dú)分割點(diǎn),從而得到正確的基礎(chǔ)心音分割點(diǎn)。實(shí)驗(yàn)結(jié)果表明,本文方法相比維奧拉積分包絡(luò)分割法、改進(jìn)型希爾伯特-黃變換包絡(luò)雙閾值分割法、基于心動(dòng)周期估計(jì)的心音包絡(luò)分割法和基于周期與波峰的自適應(yīng)閾值分割法分割精度更高,相比基于個(gè)性化高斯混合建模的包絡(luò)分割法實(shí)時(shí)性更強(qiáng),提取的包絡(luò)能夠更好地反映基礎(chǔ)心音信號(hào)特征。但本文方法對(duì)于人體運(yùn)動(dòng)狀態(tài)下的心音信號(hào)分割精度不理想,后續(xù)將對(duì)此做進(jìn)一步研究。