介紹了一種基于虛擬儀器的電腦胎兒心電圖系統。該系統由兩部分組成：以PIC單片機為核心硬件的數據采集裝置和PC平臺的處理系統。 
　　硬件電路功耗低，直接采用端口供電。軟件平臺采用可視化語言LabView6.1編程實現，由于使用其大量成熟軟件模塊，大大提高了開發效率。介紹了LabView實現的匹配濾波器等算法，并改進了匹配濾波的模板更新。臨床試用取得了滿意效果。

　　1 前言
　　虛擬儀器是機技術與儀器相結合產生的一種新型儀器模塊，它通常是由計算機、模塊化的功能硬件與帶有數據處理、過程控制的可視化軟件構成。本系統使用NI公司出品的LabView6.1圖形語言編程實現。LabView采用圖標代替傳統的文字代碼實現編程，圖標代表的程序模塊都是NI公司開發的成熟程序模塊，因此，采用LabView編程可以大大提高系統開發效率。
　　胎兒監護及圍產期監護是產科中一項重要的監護措施。常用的胎兒監護方法有：超聲多普勒監護、胎心音監護和胎心電監護。其中，胎兒心電監護因其從母體腹部提取胎兒心電信號，對母體和胎兒*，特別是能同時觀察胎兒心電圖，其優點明顯，是產科常用的診斷方法。

　　雖然腹部心電信號拾取方便，但是，母體腹部信號復雜，信噪比低，特別是一般母體心電信號比胎兒心電信號大10-20倍，而且胎兒心電約有10％-30％與母體心電重合，形成復雜的QRS波群，并且胎兒心電的頻譜與母體心電的頻譜重疊。這些因素給胎兒心電的檢測帶來了比較大的困難。學者們探索了一系列的信號處理方法解決這個，包括：匹配濾波法[1]、自適應濾波法[2]等。本系統采用比較常用的匹配濾波法。
　　2  系統的實現
　　2.1系統概要
　　本系統由硬件數據采集器和普通PC機兩部分組成，通過RS-232口連接。PC機上的應用程序采用圖形化編程語言LabView6.1實現。硬件數據采集器采樣獲得的數據以串行通信的方式發送給PC機，PC機上的軟件通過RS-232接口模塊讀入波形數據，并送預處理模塊，工頻陷波和抑制基線漂移后，進行匹配濾波去除母親心電的干擾，zui終在顯示器上顯示胎兒心電信號的波形。另外，PC機也可以串行通信方式控制硬件數據采集器的工作狀態，如調整增益、時間常數等。
　　2.2  硬件信號采集系統
　　系統硬件由前置放大、工頻陷波、主放大器、單片機系統等幾部分組成。前置放大器將導聯傳人的信號適當放大后，送工頻陷波器陷波，以防止可能出現的過強工頻干擾阻塞系統放大通道。隨后，信號經過主放大器和末級放大器幾萬倍放大之后，單片機控制的A／D對其進行400×12 bits／s采樣，數據經RS-232接口送PC機處理。硬件部分是PIC單片機控制的高精度心電信號采集系統。該系統內部采用12位高精度串行輸出A／D轉換器AD7895。采樣得到的數據通過單片機處理后將按照下面的幀格式通過串口發給PC機。系統采樣率為400Hz，數據傳輸速率為19.2kbps。該系統以低功耗為主要設計思想，其工作電流小于3mA，可以采用串口供電[3]，無須外接電源。
　　其中，12個X表示12位A／D采樣得到的數據；00、11為標志位，00代表高6位，11代表低6位。另外，硬件部分也接收PC機的控制字。
　　2.3軟件部分
　　軟件部分分為RS-232接口模塊、信號預處理、匹配濾波算法3部分，其中匹配濾波又包括相關檢測和波形對消2個部分。
　　2.3.1  RS-232接口模塊
　　（1）串口供電。由于硬件數據采集器工作電流比較小，本系統使用RS-232中空閑的握手信號線制成電源，使其工作因此，在主程序打開串口讀人數據前，必須對硬件數據采集器加電，使其處于工作狀態。在LabView中采用Serial Line Ctrl.Vi模塊，將RTS和DTR置相應電子；（2）串口讀寫。由于硬件數據采集器在整個工作期間，以19.2kbps的速率連續向外發送數據，為了數據處理方便，采用了Serial Read With Time-out FT.vi模塊，以保證每次從串口緩存中讀出等量偶數個數據。用Serial Port Write.vi模塊發送控制字給硬件模塊。控制字的發送采用了LabView6.1中新增的事件驅動[4]功能。
　　2.3.2信號預處理
　　主要功能是濾除工頻干擾和抑制基線漂移。
　　由于硬件系統采用串口供電，采集的信號中常混入50Hz的工頻干擾，必須加以抑制。抑制工頻干擾的可以從硬件和軟件2個方面來實現。在硬件方面，工頻干擾是一種共模信號，可以設計高質量的前置放大器和良好的屏蔽以及的電源模塊減弱工頻干擾對系統的。在軟件方面，使用50Hz的數字陷波的方法。本系統的數字陷波器采用Lab-View6.1 Signal Processing Toolbox中的Equi-Ripple Band-Stop PtByPt.vi控件實現。
　　基線漂移的原因有多種，其中由于運動造成體表電極與皮膚之間細微的滑動和硬件元件熱效應造成參數變化是主要原因。本系統采用非線性高通濾波的方法抑制基線漂移。具體方法是將腹部電信號通過中值濾波得到基線信號，并與原腹部信號相減。本系統采用的是Median Filter PtByPt.vi控件。顯示了實際處理效果，A為原始腹部信號；B為經過中值濾波后的基線信號;C為抑制基線漂移后的腹部信號。
　　2.4  匹配濾波算法
　　2.4.1  數據隊列的生成
　　隊列是一種常用的數據結構，本系統使用這種數據結構進行匹配濾波。但是LabView6.1已有的隊列控件queue，并不能滿足使用要求，必須重新創建一個。這里采用與C語言中類似的作法，創建一個N個元素的數組變量main Buffer，將數據不斷從尾部移動到首部。是用LabView編程實現的長對為1000的雙精度浮點數隊列，與C語言實現不同之處在于隊列中數據移動的次數是1000，而不是999，這是因為LabView初始化與 
C語言不同所致。在LabView中，如果數組某個單元沒有被附值，那么，該單元不被初始化。
　　因此，如果N=999，則第998個單元將被初始化，第999個單元因為是讀，而沒有被初始化，而致使后面的新數據輸入操作無效。
　　2.4.2  相關檢測算法
　　信號與檢測模板的互相關運算,其中：S表不經過抑制基線操移和去除工頻干擾的腹部電信號；T表示系統內部的檢測模板。經過互相關運算后就可以進行R峰檢測了。互相關運算采用Signal Processing Toolbox中的ColssCorrelation.vi控件。
　　R峰檢測實際上是對互相關運算的結果進行峰值檢測，以求得母親心電R峰的位置。用LabView系統自帶的峰值檢測控件Threshold Peak Deicer PtByPt.vi進行檢測。
　　2.4.3  波形對消算法
　　波形對消算法的基本思想是把當前母親心電QRS波群與先前求得的QRS波群模板進行加權平均，出新的QRS波群模板，并與當前QRS波群相減，濾除母親心電。臨床實驗獨立于主程序之外的定時中斷服務子程序完成顯示功能，程序設計思路簡單清晰、可讀性強。系統所用的24C01C操作簡單，3個管腳外，其余管腳接地即可。
　　3  系統軟件
　　系統軟件由主程序、顯示數字中斷服務子程序、功率控制中斷服務子程序、計時中斷服務子程序以及故障處理子程序組成。人眼視覺暫留時間為20ms，每位顯示時間間隔必須小于20ms。顯示中斷每次中斷間隔9ms，輪流切換顯示3位LED，每位數字的顯示時間間隔為18ms（<20ms），可以達到穩定的視覺效果。計時中斷服務子程序包括正計時和，其中又分為以min為單位計時和以s為單位計時。故障處理子程序完成的功能是關閉系統中斷、關閉加熱輸出、根據采樣溫度判斷并顯示故障的類型（短路、斷路以及溫度超標）、發出報警以及等待用戶確認故障信息后程序自動復位等功能。為了編程方便，本系統具體編制了按鍵檢測確認子程序、從D轉換子程序、電壓＼溫度轉換子程序、讀EEPROM以及寫EEPROM等子程序。
　　4  結束語
　　自動煎藥機以其突出的*性受到了各大與制藥廠家的歡迎，當然它也存在一些需要解決的，即：如何解決先煎、后下問題，增加濃縮功能等。這些問題有待于在進一步的實踐中不斷完善。過程中發現傳統的模板更新[1]對母親心電QRS波群變化的跟蹤能力不強，某些情況下不能*濾除母親心電QRS波群，本系統改進了該方法。系統初始化時建立一個容量為N的模板庫，將檢測到的母親心電QRS波群存入該庫中，替換zui初入庫的母親心電QRS波群，即構建一個母親心電QRS波群隊列，隨后按照新高舊低的優先級原則進行加權平均。各母親心電QRS波群的權系數。
　　用信號流圖的形式表示該算法。展示了系統對腹部電信號進行匹配濾波處理后的效果，A為從腹部提取的電信號；B為改進模板更新方法后的胎兒心電信號；C為采用傳統模板更新方法的胎兒心電信號，*表示母親心電的殘余。
　　5  結論    
　　本系統硬件無須另配電源，系統結構更簡單。軟件采用LabView6.1進行圖形化編程，直接使用NI公司提供的多種成熟的軟件模塊，因而大大節省系統開發時間和成本。本系統具有信號基線校正功能，改進了傳統匹配濾波模板更新的方法，經過近20例臨床測試，均取得良好效果。本系統在胎兒心電波群的識別上仍有改進的空間。