專利名稱：自組織脈搏傳感器中的濾波融合方法
技術(shù)領(lǐng)域：
本發(fā)明涉及信號(hào)處理技術(shù)領(lǐng)域，尤其涉及一種脈搏波智能檢測(cè)裝置中的濾波融合方法。
背景技術(shù)：
現(xiàn)有的脈搏波檢測(cè)裝置在測(cè)量病人脈搏時(shí)要求被測(cè)者處于平靜狀態(tài)，得出脈搏波 用于醫(yī)學(xué)上的輔助診斷。但是，在動(dòng)態(tài)測(cè)量脈搏波的情況下，現(xiàn)有的脈搏波檢測(cè)裝置往往精 度不能滿足要求，且醫(yī)學(xué)參考價(jià)值不大。在需要對(duì)被檢者進(jìn)行全天候監(jiān)測(cè)的情況下，靜態(tài)的 脈搏波檢測(cè)裝置不能取得很好的效果。已有的多路脈搏檢測(cè)裝置，在脈搏波后續(xù)的處理過程中沒有按照中醫(yī)中的切脈理 論去考慮合適的融合處理方法，比如寸關(guān)尺三部脈波形各異，還有不同的運(yùn)動(dòng)負(fù)荷對(duì)脈搏 造成的干擾如何消除。導(dǎo)致測(cè)量出來的脈搏波缺乏魯棒性。

發(fā)明內(nèi)容
本發(fā)明的目的在于針對(duì)現(xiàn)有技術(shù)的不足，提供一種應(yīng)用于自組織脈搏傳感器中的 脈搏波的濾波融合方法。一種自組織脈搏傳感器中的濾波融合方法，包括以下步驟(1)用小波分析的方法進(jìn)行數(shù)字濾波；(2)經(jīng)過濾噪處理之后的9路脈搏信號(hào)的融合處理。進(jìn)一步地，所述步驟(1)具體如下(A)選擇小波基選擇sym8作為小波基波；(B)選擇尺度在尺度9上進(jìn)行小波分解；(C)閾值去噪采用軟、硬閾值折衷的算法進(jìn)行閾值去噪；(D)重構(gòu)脈搏波。所述步驟(2)具體如下(a)因子分析可將原始9個(gè)脈搏通道數(shù)據(jù)序列劃分為3個(gè)脈搏波數(shù)據(jù)組，每個(gè)數(shù) 據(jù)組由在某個(gè)因子上的載荷較大，而在其它2個(gè)因子上的載荷均較小的幾個(gè)時(shí)間序列所組 成。采用如下簡(jiǎn)便的方法進(jìn)行空間降維在由3個(gè)因子決定的3個(gè)變量組中，分別取且僅取 一個(gè)典型變量，組成降維之后的變量集合；(b)選取主因子脈搏波通道選擇第1因子和第2因子上載荷最大的脈搏波通道 為兩路主因子脈搏波通道；(c)輸出融合脈搏波對(duì)兩路脈搏波信號(hào)做加權(quán)融合，最終輸出理想脈搏波信號(hào)。本發(fā)明的有益效果是本發(fā)明在不限制被測(cè)者正?；顒?dòng)的情況下，對(duì)其脈搏波進(jìn) 行24小時(shí)實(shí)時(shí)測(cè)量，并輸出接近真實(shí)的脈搏波形曲線，用于后續(xù)分析。本發(fā)明采用自組織 傳感器網(wǎng)絡(luò)技術(shù)，以點(diǎn)陣形式把脈搏波傳感器排列在手腕內(nèi)側(cè)。在佩戴很短的時(shí)間內(nèi)完成 傳感器網(wǎng)絡(luò)的自適應(yīng)調(diào)整，輸出連續(xù)真實(shí)的脈搏波形曲線。


圖1為本發(fā)明的硬件結(jié)構(gòu)圖；圖2為本發(fā)明的脈搏波分析處理流程圖；圖3為本發(fā)明的原始脈搏波經(jīng)過小波分解之后的各頻段疊加圖；圖4為本發(fā)明的在MATLAB環(huán)境中采用SWT ( 一維平穩(wěn)小波降噪)模擬小波降噪處 理的過程圖；圖5為本發(fā)明的脈搏波變量在3個(gè)因子下的聚類圖；圖6為本發(fā)明的原始脈搏波變量在3個(gè)因子上的載荷圖；圖7為本發(fā)明的A類自相關(guān)度較大的脈搏波；圖8為本發(fā)明的B類自相關(guān)度較大的脈搏波；圖9為本發(fā)明的最終處理結(jié)束后輸出脈搏波形圖。
具體實(shí)施例方式本發(fā)明自組織脈搏傳感器中的濾波融合方法通過在手腕內(nèi)側(cè)橈動(dòng)脈處放置多個(gè) 壓電式傳感器點(diǎn)陣得到多路脈搏波曲線，其中每一路單點(diǎn)傳感器可單獨(dú)測(cè)得該點(diǎn)處脈搏 波。多路脈搏波曲線在信息融合中心進(jìn)行數(shù)字濾波和因子分析。首先，采用小波處理方法進(jìn) 行數(shù)字濾波，能有效地消除了脈搏波中混入的基線漂移，肌電干擾，工頻干擾等噪聲信號(hào)。 然后，通過因子分析對(duì)濾波后的多路脈搏信號(hào)做降維處理，去除其中嚴(yán)重失真的信號(hào)通道， 同時(shí)優(yōu)化選擇兩路脈搏波信號(hào)，再計(jì)算其連接權(quán)重系數(shù)，用于抵消不同運(yùn)動(dòng)負(fù)荷下的脈搏 波畸變，最終融合為一路脈搏波輸出。采用本發(fā)明中的方法可以實(shí)時(shí)調(diào)整傳感器點(diǎn)陣的輸 出權(quán)重，輸出一路最接近真實(shí)的脈搏信號(hào)。再通過無線通信模塊與后續(xù)處理模塊接合。實(shí) 現(xiàn)實(shí)時(shí)的脈搏信號(hào)監(jiān)測(cè)分析。具體來說，本發(fā)明采用手表鏈?zhǔn)酵髱?，在手腕?nèi)側(cè)橈動(dòng)脈處放置3行3列壓電式傳 感器點(diǎn)陣。在靠近脈搏待測(cè)區(qū)域下用可充氣氣墊或彈力柱裝置使傳感器充分接觸橈動(dòng)脈處 的脈搏。每一路單點(diǎn)傳感器可單獨(dú)測(cè)得該點(diǎn)處脈搏波。9路脈搏波曲線在信息融合中心進(jìn) 行融合判斷，實(shí)時(shí)調(diào)整傳感器點(diǎn)陣的輸出權(quán)重。輸出一路最接近真實(shí)的脈搏信號(hào)，通過無線 通信模塊與后續(xù)處理模塊接合。本發(fā)明的檢測(cè)方法采用如下技術(shù)方案如圖1所示，脈搏傳感器網(wǎng)絡(luò)通過濾波電路對(duì)采集到的9路脈搏波信號(hào)進(jìn)行預(yù)處 理，預(yù)處理之后的信號(hào)經(jīng)過A/D變換送入信息融合中心(DSP模塊)，在這里對(duì)9路脈搏信號(hào) 進(jìn)行數(shù)字濾波和融合處理。第一步用小波分析的方法進(jìn)行數(shù)字濾波。在數(shù)字濾波過程中，本發(fā)明中采用小波分解和小波包分解的降噪方法，來克服脈 搏波中混入的基線漂移，肌電干擾，工頻干擾等噪聲信號(hào)。1、小波濾噪的步驟1. 1選擇小波基同一信號(hào)，選取不同的小波函數(shù)進(jìn)行處理，將得到不同的效果，所以小波函數(shù)的選 取顯得尤為重要。為此，我們選取小波從以下四個(gè)方面著手1、支集長度，表征了當(dāng)時(shí)間或頻率趨向于無窮大時(shí)，尺度函數(shù)與小波函數(shù)從一個(gè)有限值趨向于零的速度。2、對(duì)稱性，在信 號(hào)處理中對(duì)避免相移有非常重要的作用。3、消失矩階數(shù)在數(shù)據(jù)壓縮中有非常重要的作用。 4、正則性，對(duì)信號(hào)重構(gòu)及獲得較好的平滑效果十分有用。對(duì)脈搏信號(hào)濾波來說，選擇支集長度較短的小波可提高處理的實(shí)時(shí)性；選取對(duì)稱 性的小波可滿足相移為基本線性，使脈搏信號(hào)不失真；選取正則性的小波可使重構(gòu)以后的 信號(hào)比較平滑。滿足以上三個(gè)條件的小波有Daubechies小波，Symmetry小波和Coiflet小 波。選取db3、sym8、coif4小波進(jìn)行反復(fù)仿真，結(jié)果表明“sym8”小波濾波最適合脈搏波的 濾波處理，因此，選擇sym8作為小波基波。1. 2選擇尺度如圖3，對(duì)原始脈搏波信號(hào)在尺度9上的分解波形可以看出，信號(hào)中的直流分量及 趨勢(shì)項(xiàng)明顯地顯現(xiàn)在較大的尺度上，尺度8，尺度9。由于漂移信號(hào)主要為超低頻信號(hào)分量， 考慮到可以在小波重構(gòu)的過程中用來消除基線漂移，故以sym8為小波基對(duì)脈搏波在尺度9 上進(jìn)行小波分解。便于后面重構(gòu)的時(shí)候消除大尺度上的基線漂移。1. 3閾值去噪如圖4，經(jīng)過基線校正后的脈搏波信號(hào)進(jìn)入一維平穩(wěn)小波降噪處理階段。本發(fā)明中 選擇sym8基本小波，尺度為8。然后進(jìn)行小波分解。肌電干擾噪聲主要分布在低尺度上。 采用軟閾值法對(duì)低尺度(level 1 3)上的小波系數(shù)進(jìn)行大幅度的衰減進(jìn)而去除它與脈搏 波信號(hào)的頻帶相重疊的部分。即可達(dá)到消除高頻噪聲的目的。為消除肌電干擾噪聲的影響， 對(duì)4 8層尺度上的系數(shù)，采用軟、硬閾值折衷的算法進(jìn)行閾值去噪。保證重構(gòu)信號(hào)與真實(shí) 信號(hào)的最小方差。1.4重構(gòu)脈搏波去掉噪聲后的細(xì)節(jié)信號(hào)和平滑信號(hào)可以用來重構(gòu)信號(hào)，為了消除基線偏移，只要 在小波變換重構(gòu)的過程中，將該尺度下的分量置零，就可以得到去除了直流及緩變趨勢(shì)分 量的合成信號(hào)。因此，將含噪脈搏信號(hào)利用sym8小波進(jìn)行9層小波分解，得到各層小波系數(shù)。去 除N大于等于9的小波分量。在本研究中脈搏信號(hào)的采樣頻率為200Hz，sym8小波分解在尺度9下逼近信號(hào)的 頻率和功率極低，因此，原始脈搏波信號(hào)的低頻信號(hào)的主要成分在經(jīng)過基線矯正后不受影 響。在matlab小波工具箱中按照上述步驟演示濾波過程(如圖4)(1)啟動(dòng)SWT并裝載信號(hào)采用Matlab(The MathWorks，Inc.)作為數(shù)據(jù)處理軟件。版本為R2009a。在 Wavelet Toolbox Main Menu窗口選擇SWT De-noising 1_D按鈕，裝載要分析原始脈搏 波數(shù)據(jù)文件。(2)對(duì)信號(hào)進(jìn)行一維平穩(wěn)小波變換選擇基本小波為sym8，尺度為8。經(jīng)過小波分解，如圖4左邊一列圖形。顯示分解后的低頻和高頻信號(hào)。(3)利用平穩(wěn)小波變換來消噪選擇固定軟閾值，并且按照1.3中的方法調(diào)整各層次的閾值。點(diǎn)擊降噪命令。結(jié)果如(圖4)右下列波形所示。最終輸出波形見(圖4)右上波形。第二步經(jīng)過濾噪處理之后的9路脈搏信號(hào)的融合處理。通過因子分析(Factor Analysis)對(duì)濾噪后的數(shù)據(jù)進(jìn)行空間降維，目的是消除多通道數(shù)據(jù)的空間線性相關(guān)性和剔除無用信息的干擾。最終確定最優(yōu)的傳感器網(wǎng)絡(luò)模型。2. 1因子分析設(shè)定因子數(shù)目為3，采用promax算子進(jìn)行因子旋轉(zhuǎn)。表1為原始9個(gè)變量在3個(gè) 因子上的載荷(Loading)。本發(fā)明中對(duì)脈搏波信號(hào)進(jìn)行因子分析，可將原始9個(gè)脈搏通道數(shù)據(jù)序列劃分為3 個(gè)脈搏波數(shù)據(jù)組，每個(gè)數(shù)據(jù)組由在某個(gè)因子上的載荷較大，而在其它2個(gè)因子上的載荷均 較小的幾個(gè)時(shí)間序列所組成。采用如下簡(jiǎn)便的方法進(jìn)行空間降維在由3個(gè)因子決定的3 個(gè)變量組中，分別取且僅取一個(gè)典型變量，組成降維之后的變量集合。其中，典型變量的數(shù) 學(xué)意義是在所有原始變量中，該變量在特定因子上具有最高的載荷(Loading)。三個(gè)主因 子的物理意義是第1因子上承載的信息最接近真實(shí)脈搏波；第2因子上承載的信息最接 近運(yùn)動(dòng)干擾脈搏信號(hào)；第3因子上承載的信息最接近測(cè)量噪聲。在SPSS中進(jìn)行因子分析，步驟如下(1)啟動(dòng)SPSS并裝載信號(hào)SPSS軟件版本為SPSS16. 0。導(dǎo)入9路脈搏波的時(shí)間序列數(shù)據(jù)。(2)選擇因子分析工具開始進(jìn)行因子分析選擇顯示因子聚類圖；設(shè)定因子數(shù)目為3 ；選擇提取因子方法為PCA ；旋轉(zhuǎn)方法為 Promax。開始分析。(3)輸出分析結(jié)果因子分析結(jié)束后導(dǎo)出分析結(jié)果，見圖5，圖6。2. 2選取主因子脈搏波通道脈搏波信號(hào)因子分析之后產(chǎn)生3列主因子脈搏波信號(hào)。把第1因子和第2因子上 載荷(Loading)最大的脈搏波通道，標(biāo)記為A路和B路脈搏波信號(hào)。去除第3因子上對(duì)應(yīng) 的脈搏波采集通道，對(duì)應(yīng)的是被嚴(yán)重失真的脈搏波通道。其中，A路信號(hào)最能反映真實(shí)脈搏 信息；B路信號(hào)為參考脈搏信號(hào)，來自無脈搏處的皮膚波動(dòng)。用于對(duì)不同運(yùn)動(dòng)情況下對(duì)A路 脈搏信號(hào)去除擾動(dòng)。A路脈搏波信號(hào)與B路脈搏波信號(hào)以合適的權(quán)重融合輸出最終的脈搏 波信號(hào)。同時(shí)，A路脈搏波信號(hào)和B路脈搏波信號(hào)獲得傳感器陣列中的最大鏈接權(quán)重。2. 3輸出融合脈搏波如圖5所示。在3個(gè)因子上9路變量被分成3類。通過各個(gè)變量在因子上的載荷 得分(圖6)把脈搏波分成3類。其中脈搏波變量V1，V5，V6，V7自相關(guān)度較大，聚成一類， 計(jì)算其分形維數(shù)，屬于標(biāo)準(zhǔn)脈搏波的分形維數(shù)區(qū)域內(nèi)，故定義為A類脈搏波(如圖7)。去除 嚴(yán)重失真的一路脈搏波。剩下的脈搏波變量￥2，￥4，￥8，￥3(見圖8)屬于B類脈搏波。分別選擇A類和B類脈搏波中因子載荷(Loading)最大的脈搏波通道，標(biāo)記為A 路和B路脈搏波信號(hào)。對(duì)兩路脈搏波信號(hào)做加權(quán)融合。B路的權(quán)重Ic1在-0. 3到-0. 5之間 選擇，A路的權(quán)重k2在1到0. 8之間選擇。輸出脈搏波序列F為A = Ii1Xa^k2Xbi其中，fi為t = i時(shí)刻的脈搏波值，a,和bi分別為A路脈搏波和B路脈搏波t =i時(shí)刻的脈搏波值。最終輸出理想脈搏波信號(hào)(圖9)。當(dāng)兩路輸出脈搏波通道被選定之后，相應(yīng)的對(duì)脈搏波傳感器網(wǎng)絡(luò) SN(sensornetwork)的輸出權(quán)重進(jìn)行調(diào)整。對(duì)應(yīng)A路和B路脈搏波信號(hào)的傳感器測(cè)量點(diǎn)得 到權(quán)重值為1。其它測(cè)量點(diǎn)權(quán)重置零，即只接收來自主因子通道的脈搏波信號(hào)。以減小后續(xù)的運(yùn)算量。當(dāng)脈搏傳感器的佩戴位置改變時(shí)，傳感器網(wǎng)絡(luò)可以根據(jù)上述方法自適應(yīng)的調(diào)整 連接權(quán)重，迅速找到最合適的脈搏采集節(jié)點(diǎn)并輸出結(jié)果。
權(quán)利要求
一種自組織脈搏傳感器中的濾波融合方法，其特征在于，包括以下步驟(1)用小波分析的方法進(jìn)行數(shù)字濾波。(2)經(jīng)過濾噪處理之后的9路脈搏信號(hào)的融合處理。
2.根據(jù)權(quán)利要求1所述自組織脈搏傳感器中的濾波融合方法，其特征在于，所述步驟(1)具體如下(A)選擇小波基選擇sym8作為小波基波。(B)選擇尺度在尺度9上進(jìn)行小波分解。(C)閾值去噪采用軟、硬閾值折衷的算法進(jìn)行閾值去噪。(D)重構(gòu)脈搏波。
3.根據(jù)權(quán)利要求1所述自組織脈搏傳感器中的濾波融合方法，其特征在于，所述步驟(2)具體如下(a)因子分析可將原始9個(gè)脈搏通道數(shù)據(jù)序列劃分為3個(gè)脈搏波數(shù)據(jù)組，每個(gè)數(shù)據(jù)組 由在某個(gè)因子上的載荷較大，而在其它2個(gè)因子上的載荷均較小的幾個(gè)時(shí)間序列所組成。 采用如下簡(jiǎn)便的方法進(jìn)行空間降維在由3個(gè)因子決定的3個(gè)變量組中，分別取且僅取一個(gè) 典型變量，組成降維之后的變量集合。(b)選取主因子脈搏波通道選擇第1因子和第2因子上載荷最大的脈搏波通道為兩 路主因子脈搏波通道。(c)輸出融合脈搏波對(duì)兩路脈搏波信號(hào)做加權(quán)融合，最終輸出理想脈搏波信號(hào)。
全文摘要
本發(fā)明公開了一種應(yīng)用于自組織脈搏傳感器中的濾波融合方法。此方法通過在手腕內(nèi)側(cè)橈動(dòng)脈處放置多個(gè)壓電式傳感器點(diǎn)陣(每個(gè)單點(diǎn)傳感器可單獨(dú)測(cè)得脈搏波)獲得多路脈搏波曲線。曲線在信息融合中心進(jìn)行數(shù)字濾波和因子分析首先，數(shù)字濾波采用小波處理方法，能有效消除脈搏波中混入的如基線漂移、肌電干擾、工頻干擾等噪聲信號(hào)；然后，因子分析對(duì)濾波后的多路脈搏信號(hào)做降維處理，去除其中嚴(yán)重失真的信號(hào)通道，同時(shí)優(yōu)化選擇兩路脈搏波信號(hào)，最終融合為一路脈搏波輸出。采用本發(fā)明的方法可以實(shí)時(shí)調(diào)整傳感器點(diǎn)陣的輸出權(quán)重，從而得到最接近真實(shí)脈搏信號(hào)的一路脈搏信號(hào)，然后通過無線通信模塊與后續(xù)處理模塊接合，最終實(shí)現(xiàn)實(shí)時(shí)的脈搏信號(hào)監(jiān)測(cè)與分析。
文檔編號(hào)A61B5/02GK101803911SQ201010139528
公開日2010年8月18日 申請(qǐng)日期2010年4月2日 優(yōu)先權(quán)日2010年4月2日
發(fā)明者孟濬, 王磊, 陳嘯, 黃小靜, 黃德樣 申請(qǐng)人:浙江大學(xué)