刁愛芹

基于Labview的心電信號關(guān)聯(lián)維數(shù)的快速實現(xiàn)*

刁愛芹①

目的：研究基于Labview的心電信號關(guān)聯(lián)維數(shù)的快速實現(xiàn)方法。方法：以PCI-6023E板卡作為硬件，以圖形化編程語言LabVIEW為開發(fā)軟件，實現(xiàn)心電信號的采集、顯示和處理，并通過軟件編程實現(xiàn)心電信號關(guān)聯(lián)維數(shù)的快速計算。結(jié)果：對冠心病患者和健康人的心電信號關(guān)聯(lián)維數(shù)進行分析計算，可以快速計算出關(guān)聯(lián)維數(shù)，冠心病患者的關(guān)聯(lián)維數(shù)要明顯低于正常人。結(jié)論：基于Labview的心電信號關(guān)聯(lián)維數(shù)系統(tǒng)，實現(xiàn)心電信號的關(guān)聯(lián)維數(shù)的快速計算，解決從心電信號中提取心率變異(HRV)信息較難以及提取到有用信息不足的問題，為臨床診斷和治療提供重要的參考信息。

LabVIEW；心電信號；關(guān)聯(lián)維數(shù)；心率變異

心臟是人體重要的器官之一，通過心肌細胞的收縮按照一定的節(jié)奏向人體全身供血，受神經(jīng)、血壓、激素等多方面因素的制約，因此是一個復(fù)雜的非線性生物系統(tǒng)，而現(xiàn)代醫(yī)學(xué)對于心電信號的分析仍以頻譜分析和小波分析為代表的線性方法為主。

20世紀80年代以來，對心臟電活動非線性現(xiàn)象的研究有了長足的進步，但多數(shù)研究集中在心率變異信號(heart rate variability，HRV)方面[1]。由于對HRV信號進行非線性研究，需要采集足夠長(數(shù)小時)的信號，并且要保證研究的系統(tǒng)處于長時間的平穩(wěn)狀態(tài)，實踐中往往難以實現(xiàn)。因此，近年來對心電圖(electrocardiogram，ECG)信號的非線性研究[2-4]得到了重視，并在臨床上證實了許多重要的臨床綜合癥，如心臟猝死、心臟衰竭、急性心肌梗塞等，由于自主神經(jīng)功能的改變而擾動心率控制，使心電圖信號的關(guān)聯(lián)維數(shù)發(fā)生變化。

分形和混沌理論是當今非線性科學(xué)一個重要的方面，特別適合研究各種復(fù)雜非線性系統(tǒng)，其復(fù)雜程度可以用分數(shù)維來描述。關(guān)聯(lián)維作為分數(shù)維的一種，能較準確的反映系統(tǒng)的動態(tài)結(jié)構(gòu)，但是分形和混沌理論需要復(fù)雜的數(shù)學(xué)計算，心電信號的關(guān)聯(lián)維數(shù)對諸多心臟病臨床綜合癥具有重要的臨床參考價值[5-6]。傳統(tǒng)的ECG機無法對采集到的ECG信號完成關(guān)聯(lián)維的計算，因此，本研究以LabVIEW為開發(fā)軟件[7]、以NI公司的PCI-6023E板卡作為硬件，實現(xiàn)心電信號的采集、顯示和處理。LabVIEW軟件結(jié)合了圖形化編程方式的高效和靈活，其提供的大量圖形化的算法函數(shù)為關(guān)聯(lián)維數(shù)的快速計算提供了方便。

1 硬件設(shè)計

采用按照臨床實踐確立的12導(dǎo)聯(lián)電極，安放在人體的相應(yīng)部位來采集心電信號，ECG機作為信號調(diào)理設(shè)備，可以對電極采集到的信號進行放大和濾波，消除信號中的工頻干擾和肌電干擾，并將心電信號轉(zhuǎn)換成標準電信號[8-9]。通過ECG機上的CRO接口和DAQ采集板卡連接，將ECG機輸出的標準電信號轉(zhuǎn)換成數(shù)字信號輸入計算機[10]。系統(tǒng)采用PCI-6023E板卡進行數(shù)據(jù)采集，該板卡可以實現(xiàn)16單通道輸入和8對差分輸入，實現(xiàn)多個信號的同時采集，為了減小誤差，提高信號精度，系統(tǒng)采用差分輸入方式[11-12]。計算機通過信號分析軟件將DAQ采集板卡采集到的信息進行處理，并形象地顯示出來，其硬件組成如圖1所示[13-14]。

圖1 系統(tǒng)硬件組成框圖

圖2 系統(tǒng)軟件流程圖

2 軟件設(shè)計

采用NI公司的LabVIEW軟件作為開發(fā)平臺，此軟件采用圖形化的編程方式，可較簡便直觀的進行軟件開發(fā)，同時由于軟件中包含了豐富的功能函數(shù)庫，極大縮短系統(tǒng)的開發(fā)時間，并降低軟件的開發(fā)難度。系統(tǒng)軟件流程如圖2所示。

2.1 信號采集

數(shù)據(jù)采集模塊程序的作用是對ECG機調(diào)理后的信號進行連續(xù)采集，程序中AI Config VI函數(shù)用來配置采樣時的各個參數(shù)，如通道配置、緩沖區(qū)配置等[15-17]。AI Start VI函數(shù)用來指定程序中用于開始采樣的參數(shù)，如掃描次數(shù)配置(-1表示采集一個緩沖區(qū)的數(shù)據(jù))、采樣率配置等。AI Read VI設(shè)定從采樣緩沖區(qū)讀取數(shù)據(jù)的參數(shù)，如每次讀取的數(shù)據(jù)量設(shè)定，scan backlog用來顯示緩沖區(qū)中剩余的數(shù)據(jù)量。采集完成后調(diào)用AI Clear VI函數(shù)將task ID清除，從而釋放采集占用的空間(如圖3所示)。

圖3 數(shù)據(jù)采集程序圖

2.2 小波降噪

數(shù)據(jù)采集的過程中會受到各種噪聲信號的干擾，如呼吸波(低頻)、市電干擾(中頻)及高頻噪聲[18-19]。這些噪聲會嚴重的影響非線性參數(shù)的計算，因此在關(guān)聯(lián)維數(shù)的計算之前必須能同時消除這些干擾信號。傳統(tǒng)的濾波方法可以有效的過濾噪聲，但是也會衰減實際信號造成測量的誤差，小波包濾波技術(shù)可以很好的解決上述問題[20-22]。

MATLAB軟件附帶的小波工具箱中含有多種常用的小波和小波包，并且可以方便靈活的實現(xiàn)小波分析。因此，本研究通過在LabVIEW中調(diào)用MATLAB Script節(jié)點來實現(xiàn)小波降噪。LabVIEW與MATLAB Script節(jié)點間參數(shù)的傳遞可以通過添加輸入、輸出接口完成，這樣參數(shù)就可以通過LabVIEW與MATLAB接口傳遞給MATLAB應(yīng)用程序，完成信號的分析和處理，最后再回傳給LabVIEW。本系統(tǒng)前端12導(dǎo)聯(lián)放大器的頻響范圍為0.5～1000 Hz，選用小波基6.8，軟閾值去噪，設(shè)計濾波的頻響為2～250 Hz。對比發(fā)現(xiàn)降噪后的心電信號在保證幅值和形狀基本不變的基礎(chǔ)上，各種干擾噪聲信號得到了極大的削弱。降噪前后的信號如圖4所示。

圖4 小波降噪前后的結(jié)果對比圖

2.3 關(guān)聯(lián)維數(shù)計算

2.3.1 相空間的重構(gòu)

相空間的重構(gòu)一般采用延時法，設(shè)系統(tǒng)嵌入的維數(shù)為k，延時時間為τ，則狀態(tài)空間向量即為公式1：

xi=(xi，xi+τ，…，xi+(k-1)τ) i=1，2，…，N N=n-kt+τ(1)

Xi對應(yīng)k維相空間的一個點，N個點構(gòu)成相空間中一條動力學(xué)軌道。將心電信號第一次經(jīng)過零值時所對應(yīng)的t選為延時時間。相空間重構(gòu)的程序如圖5所示，其中內(nèi)層For循環(huán)構(gòu)造相空間中的一個點；外層For循環(huán)結(jié)構(gòu)構(gòu)造完整的相空間；循環(huán)體外面的程序用來計算N和延時時間τ。

圖5 相空間的重構(gòu)程序圖

表1 冠心病患者和正常人關(guān)聯(lián)維數(shù)比較

2.3.2 關(guān)聯(lián)維數(shù)的計算

計算關(guān)聯(lián)維數(shù)的步驟如下。

(1)首先給定一個距離數(shù)值ε，計算符號函數(shù)θ(α)，即為公式2：

式中θ為重構(gòu)的相空間中距離小于ε的點對數(shù)。

(2)關(guān)聯(lián)積分的計算為公式3：

式中C(ε)表示距離小于ε的點對在整個點對中所占有的比例，當ε值在一定范圍內(nèi)時，應(yīng)滿足公式4：

(3)計算關(guān)聯(lián)維數(shù)，對公式4的兩邊分別取對數(shù)得公式5：

式中d即為關(guān)聯(lián)維數(shù)。

通常關(guān)聯(lián)維數(shù)的大小是通過求取標度曲線lnε～lnC(ε)的直線部分的斜率來確定。通過LabVIEW的“線性擬合”函數(shù)來擬合直線，然后通過讀取擬合直線的斜率來確定關(guān)聯(lián)維數(shù)，其具體程序如圖6所示。

圖6 關(guān)聯(lián)維數(shù)計算程序圖

3 關(guān)聯(lián)維數(shù)分析

系統(tǒng)對冠心病患者和正常人的12導(dǎo)聯(lián)ECG信號進行采集、處理和計算，得到的關(guān)聯(lián)維數(shù)見表1。

比較發(fā)現(xiàn)，不論正常人還是冠心病患者，從不同導(dǎo)聯(lián)信號分析得到的關(guān)聯(lián)維數(shù)各不相同，但從同一導(dǎo)聯(lián)信號得到的關(guān)聯(lián)維數(shù)來看，冠心病患者的關(guān)聯(lián)維數(shù)要明顯低于正常人，顯示了在心肌缺血情況下，心臟的復(fù)雜程度有所降低。因此，通過本系統(tǒng)對ECG信號計算得到的關(guān)聯(lián)維數(shù)可以為心臟病患者提供有價值的臨床診斷信息。

4 結(jié)論

本系統(tǒng)在傳統(tǒng)ECG機的基礎(chǔ)上，以LabVIEW為軟件開發(fā)平臺，以DAQ采集板卡為硬件支持，實現(xiàn)了人體心電信號的采集、處理以及關(guān)聯(lián)維數(shù)的計算。實驗證明，系統(tǒng)不僅可以快速計算出心電信號的關(guān)聯(lián)維數(shù)，而且計算得到的數(shù)據(jù)真實可靠，解決了之前從心電信號中提取HRV信息較難，以及提取到有用信息不足的問題，為臨床診斷和治療提供了重要的參考信息。

[1]吳旭輝，卞春華，范虎偉，等.HRV信號關(guān)聯(lián)維計算中參數(shù)選取的研究[J].北京生物醫(yī)學(xué)工程，2009，28(4)：358-362.

[2]白建軍，山方，郭盟.心電模型中的混沌特性研究[J].電腦知識與技術(shù)，2012，8(7)：1164-1167.

[3]寧新寶，卞春華，王俊，等.心臟電活動過程的非線性分析[J].科學(xué)通報，2006，51(7)：764-771.

[4]張濤，章偉睿，楊永青，等.正常心電信號及室性心動過速和室顫快速識別方法的初步研究[J].中國醫(yī)學(xué)裝備，2017，14(7)：37-40.

[5]溫曉通，孟麗艷，朱勁松，等.一種非線性時間序列的關(guān)聯(lián)維快速算法[J].北京師范大學(xué)學(xué)報(自然科學(xué)版)，2005，15(4)：358-361.

[6]謝川.基于混沌關(guān)聯(lián)維特征提取的大數(shù)據(jù)聚類算法[J].計算機科學(xué)，2016，43(6)：229-232.

[7]李江全.LabVIEW虛擬儀器從入門到測控應(yīng)用130例[M].北京：電子工業(yè)出版社，2013：302-308.

[8]于杰，李川勇，賈林壯.基于LabVIEW的心電信號采集系統(tǒng)的設(shè)計[J].生物醫(yī)學(xué)工程與臨床，2001，5(3)131-133.

[9]朱芳，李淮江，姜恩華，等.心電信號預(yù)處理數(shù)字濾波器的設(shè)計與分析[J].淮北師范大學(xué)學(xué)報：自然科學(xué)版，2017，38(1)：36-40.

[10]行鴻彥，黃敏松.基于LabVIEW虛擬儀器的心電信號采集系統(tǒng)的設(shè)計[J].電子測量與儀器學(xué)報，2008，22(5)：107-111.

[11]何伶俐.基于LabVIEW的虛擬心電監(jiān)護儀系統(tǒng)設(shè)計[J].醫(yī)療衛(wèi)生裝備，2017，38(2)：24-27.

[12]王宗剛，潘崢嶸，祝燎.基于LabVIEW的心電信號采集系統(tǒng)的設(shè)計與實現(xiàn)[J].工業(yè)儀表與自動化裝置，2011(3)：40-46.

[13]鞏萍，夏文艷，倪紅艷.基于虛擬儀器技術(shù)的心電信號采集與分析系統(tǒng)[J].醫(yī)療衛(wèi)生裝備，2013，34(9)：17-19.

[14]曼世超，盧會國，盧勇.基于LabVIEW的心電信號檢測平臺的設(shè)計[J].成都信息工程學(xué)院學(xué)報，2016，31(3)：255-260.

[14]吳旭輝，卞春華，范虎偉，等.HRV信號關(guān)聯(lián)維計算中參數(shù)選取的研究[J].北京生物醫(yī)學(xué)工程，2009，28(4)：358-362.

[15]黃進文.基于LabVIEW的ECG的混沌特性研究[J].保山學(xué)院學(xué)報，2013，32(2)：16-18.

[16]胡廣平，周華強.基于小波變換的心電采集分析系統(tǒng)設(shè)計與實現(xiàn)[J].計算機測量與控制，2014，22(3)：854-856.

[17]李濤，劉延武.基于labview的心電信號采集和分析[J].醫(yī)療設(shè)備信息，2007，22(6)：15-17.

[18]王芳，季忠，彭承琳.基于雙樹復(fù)小波變換的心電信號去噪研究[J].儀器儀表學(xué)報，2013，34(5)：1161-1166.

[19]王吉鳴，呂穎瑩，董晗，等.一種濾除高采樣心電工頻干擾的改進算法[J].中國生物醫(yī)學(xué)工程學(xué)報，2016，35(6)：744-748.

[20]韋高.基于小波變換的心電信號去噪法[J].中國醫(yī)學(xué)裝備，2010，7(8)：9-11.

[21]鄭敏敏，高小榕，謝海鶴.心電信號小波去噪的改進算法研究[J].中國生物醫(yī)學(xué)工程學(xué)報，2017，36(1)：114-118.

[22]劉書俊，李生林，將明，等.一種基于加權(quán)平均的改進型小波閾值降噪算法[J].化工自動化及儀表，2017(3)：239-242.

Fast-implementation of correlation dimension of ECG signal base on LabVIEW

/DIAO Ai-qin

Objective:To research a fast-implementation method of correlation dimension of ECG signal base on LabVIEW.Methods:In this paper, acquisition ,display and processing of ECG signal were implemented through applied PCI-6023E as hardware and applied graphical programming language LabVIEW as development software. And then through software programming to realize fast calculation for correlation dimension of ECG signal.Results:This method was used to analyze and calculate the correlation dimension of ECG between patients with coronary heart disease and healthy people. And the experiment indicated that it can fast and reliably calculate the correlation dimension, and the correlation dimension of patients with coronary heart disease was obviously lower than that of healthy people.Conclusion:The correlation dimension system of ECG signal based on LabVIEW can realize the fast calculation of correlation dimension of ECG signal, and solve the two problems that heart rate variability (HRV) information are difficultly extracted from ECG signal and the extracted useful information are inadequate. Therefore, it provides important reference information for clinical diagnosis and treatment.

LabVIEW; ECG signal; Correlation dimension; Heart rate variability

School of Medical Technology, Taizhou Polytechnical College, Taizhou 225300,China.

刁愛芹,女,(1986- ),碩士，講師。泰州職業(yè)技術(shù)學(xué)院醫(yī)學(xué)技術(shù)學(xué)院，研究方向：心血管信號轉(zhuǎn)導(dǎo)。

1672-8270(2017)12-0023-04

R197.39

A

10.3969/J.ISSN.1672-8270.2017.12.007

泰州市科技社會發(fā)展項目(TS201523)“雌激素通過抑制Integrin?1/FAK/p38 MAPK信號通路減輕機械牽張誘導(dǎo)的心肌細胞肥大的研究”；泰州職業(yè)技術(shù)學(xué)院科研重點項目(TZXY-15-8)“大豆異黃酮下調(diào)TLR4/MyD88介導(dǎo)的信號通路保護心肌肥厚”

①泰州職業(yè)技術(shù)學(xué)院醫(yī)學(xué)技術(shù)學(xué)院 江蘇 泰州 225300

//China Medical Equipment,2017,14(12):23-26.

2017-06-22