動(dòng)態(tài)脈搏信號(hào)的采集與處理

【作 者】張愛(ài)華，丑永新

1 蘭州理工大學(xué) 電氣工程與信息工程學(xué)院，甘肅，蘭州，730050

2 甘肅省工業(yè)過(guò)程先進(jìn)控制重點(diǎn)實(shí)驗(yàn)室，甘肅，蘭州，730050

近年來(lái)，隨著人們生活方式、飲食結(jié)構(gòu)的不斷改變，諸如高血壓、冠心病等疾病成為常見病，并且發(fā)病率不斷上升，特別是心血管疾病已成為威脅人類生命的主要疾病。為了避免和減輕這些疾病對(duì)人類健康的危害，如可穿戴式心電監(jiān)護(hù)儀，遠(yuǎn)程醫(yī)療監(jiān)護(hù)系統(tǒng)等應(yīng)運(yùn)而生。文獻(xiàn)[1]-[6]提到的脈搏信號(hào)檢測(cè)系統(tǒng)，將采集的人體生理信號(hào)作一定的處理，然后通過(guò)有線或無(wú)線通信的方式傳遞給PC機(jī)，由醫(yī)生診斷并給出相應(yīng)的治療措施。這樣給患者帶來(lái)許多方便，提高了醫(yī)生治病的靈活性和有效性。但由于心血管疾病發(fā)病的隨機(jī)性和快速性的特點(diǎn)，提高監(jiān)護(hù)系統(tǒng)的實(shí)時(shí)性就顯得十分重要。

脈搏信號(hào)蘊(yùn)含了豐富的人體生理病理信息，且用光電脈搏傳感器信號(hào)檢測(cè)方便[7]。本文以TMS320VC5509DSP為核心，研制包括脈搏信號(hào)采集電路在內(nèi)的動(dòng)態(tài)脈搏信號(hào)實(shí)時(shí)檢測(cè)與分析系統(tǒng)。針對(duì)動(dòng)態(tài)脈搏信號(hào)的特點(diǎn)，設(shè)計(jì)了整系數(shù)陷波和低通濾波器，提出了包絡(luò)濾波法。然后，通過(guò)優(yōu)化編程方法，在DSP處理器上實(shí)現(xiàn)動(dòng)態(tài)脈搏信號(hào)的實(shí)時(shí)濾波，對(duì)濾波后的信號(hào)進(jìn)行時(shí)域和頻域分析，計(jì)算脈搏信號(hào)的脈率、波形高度，并對(duì)脈搏信號(hào)進(jìn)行FFT變換，實(shí)現(xiàn)光電脈搏信號(hào)的實(shí)時(shí)檢測(cè)和簡(jiǎn)單處理。

1 整系數(shù)濾波器的設(shè)計(jì)

脈搏信號(hào)是隨著心臟的搏動(dòng)而產(chǎn)生的，一般比較微弱，特別是在不影響人體正?；顒?dòng)的條件下，采集的信號(hào)有許多來(lái)自外界的干擾。這些干擾主要有受呼吸和皮膚接觸滑動(dòng)影響產(chǎn)生的基線漂移，人體電位變化產(chǎn)生的肌電干擾，50 Hz交流電及其多次諧波引起的工頻干擾等[8]。這些干擾會(huì)使得脈搏采集系統(tǒng)的信噪比下降，甚至?xí)蜎](méi)有用的脈搏信號(hào)，僅靠信號(hào)采集系統(tǒng)硬件設(shè)計(jì)上的措施并不能完全消除這些干擾。因此，采取有效的手段消除或減少基線漂移、肌電干擾和工頻干擾，是識(shí)別和處理脈搏信號(hào)的前提，在脈搏檢測(cè)儀器中尤為重要。如文獻(xiàn)[9]-[13]，采用的IIR和FIR數(shù)字濾波等方法可以獲得相當(dāng)好的濾波效果，但這些濾波器系數(shù)一般為非整數(shù)，在實(shí)際應(yīng)用中，特別是實(shí)時(shí)信號(hào)處理場(chǎng)合，對(duì)信號(hào)的實(shí)時(shí)濾波速度有很大的影響。本文針對(duì)脈搏信號(hào)的特征，設(shè)計(jì)了整系數(shù)陷波、低通濾波器，并結(jié)合包絡(luò)法濾除動(dòng)態(tài)脈搏信號(hào)中的干擾信號(hào)。

1.1 陷波器的設(shè)計(jì)

為了去除脈搏信號(hào)中的工頻干擾，需要設(shè)計(jì)50 Hz的陷波器，同時(shí)為保證濾波前后信號(hào)不失真和便于DSP系統(tǒng)的實(shí)時(shí)濾波，該濾波器還有嚴(yán)格的線性相位且濾波計(jì)算量不宜太大的要求。因此，采用兩個(gè)具有相同傳輸延遲的線性相位濾波器相減得到陷波濾波，即從一個(gè)全通網(wǎng)絡(luò)中減去一個(gè)帶通網(wǎng)絡(luò)。帶通網(wǎng)絡(luò)可以用一個(gè)整系數(shù)的梳狀濾波器實(shí)現(xiàn)，全通網(wǎng)絡(luò)只具有延遲作用，起到抵消梳狀帶通網(wǎng)絡(luò)帶來(lái)的相位延遲作用[14]。系統(tǒng)的采樣頻率為400 Hz，則50陷波濾波器的傳遞函數(shù)和頻率響應(yīng)分別為：

式（1）中第一項(xiàng)為全通網(wǎng)絡(luò)，后一項(xiàng)為梳狀帶通濾波網(wǎng)絡(luò)；l、k為正整數(shù)，決定陷波器的阻帶寬度和通帶波紋。通過(guò)實(shí)驗(yàn)，選取l=12，k=2，則式（1）變?yōu)椋?/p>

顯然，式(3)前后兩項(xiàng)的純延時(shí)為96。由MATLAB仿真可得陷波濾波器的幅頻特性曲線如圖1(a)所示。由式(1)可求得該陷波濾波器的輸出方程為：

1.2 低通濾波器的設(shè)計(jì)

由于經(jīng)過(guò)陷波濾波器處理以后，脈搏信號(hào)中仍然存在較為嚴(yán)重肌電干擾，且其頻率范圍廣泛，不能從原始信號(hào)中根本的消除，只能通過(guò)低通濾波器減少其對(duì)脈搏信號(hào)的干擾。為此，采用零極點(diǎn)對(duì)消法設(shè)計(jì)簡(jiǎn)單的整系數(shù)低通濾波器。該濾波器的傳遞函數(shù)及頻率響應(yīng)如下：

由式(7)可以看出，此低通濾波器有M個(gè)零點(diǎn)和極點(diǎn)。其中，M個(gè)零點(diǎn)均勻的分布在單位圓上，在z=1處有一個(gè)極點(diǎn)，用來(lái)抵消第一個(gè)零點(diǎn)，其他的極點(diǎn)都集中在原點(diǎn)處，因而具有低通特性。整個(gè)單位圓圓周角取值范圍是0-2π，實(shí)際采樣頻率為fs，其對(duì)應(yīng)的圓周頻率為0～2 πfs。為了使得低通截止頻率相對(duì)較小，假設(shè)幅頻特性曲線第一次降為零的頻率為2πf(0)/fs，所有零點(diǎn)的個(gè)數(shù)為M=fs/f(0)。由于DSP系統(tǒng)采樣頻率為400 Hz，且必須保證M為整數(shù)，取f(0)=66.67 Hz，得到M=6，則傳遞函數(shù)為：

由MATLAB仿真可得低通濾波器的幅頻特性曲線如圖1(b)所示。其第一次降為零點(diǎn)時(shí)的截至頻率為66.7 Hz，在通帶范圍內(nèi)，濾波器具有線性相位，滿足濾波要求。由式(9)可求得該低通濾波器的輸出差分方程為：

圖1 濾波器的幅頻特性曲線Fig.1 Magnitude response of fi lters

1.3 包絡(luò)法去基線漂移

基線漂移屬超低頻信號(hào)，可以通過(guò)高通濾波器從原始的脈搏信號(hào)中消除。但在高頻濾波過(guò)程中，或多或少會(huì)使脈搏信號(hào)中的有效成份丟失。為此，提出一種新穎的濾波方法—包絡(luò)法，去除脈搏信號(hào)中的基線漂移。

如圖2所示，上、下兩包絡(luò)線是脈搏信號(hào)的峰值和谷值經(jīng)過(guò)多項(xiàng)式曲線擬合法擬合的曲線，中線為峰值和谷值包絡(luò)曲線的均值曲線，可以近似反映脈搏信號(hào)的基線漂移。用原始的脈搏信號(hào)減去均值曲線，即可去除脈搏信號(hào)中的基線漂移。

圖2 脈搏信號(hào)包絡(luò)示意圖Fig.2 Envelope diagram of pulse signal

圖2所示的脈搏包絡(luò)曲線是在PC機(jī)上仿真得到，在檢測(cè)出脈搏信號(hào)的峰值和谷值之后，需要用多項(xiàng)式曲線擬合法擬合出脈搏信號(hào)的峰值和谷值包絡(luò)。但由于DSP系統(tǒng)運(yùn)行速度的限制，很難用多項(xiàng)式曲線擬合法實(shí)時(shí)地?cái)M合出脈搏信號(hào)的包絡(luò)曲線。我們用相鄰兩峰值（谷值）之間的連線近似的代替擬合曲線，作為脈搏信號(hào)的包絡(luò)曲線，這樣極大地降低了算法的計(jì)算量，可以在DSP系統(tǒng)內(nèi)兩采樣點(diǎn)之間實(shí)時(shí)地實(shí)現(xiàn)包絡(luò)曲線的估算和基線漂移的去除。

2 脈搏信號(hào)的時(shí)域和頻域分析

2.1 脈搏信號(hào)的時(shí)域分析

時(shí)域分析主要是指在時(shí)域提取和分析脈搏特征，如幅值特征、時(shí)間特征等。本文主要分析脈搏信號(hào)的波高和脈率。

2.1.1 波高檢測(cè)

波高是指脈搏波波峰和相鄰波谷幅度的差值，要計(jì)算波高，首先要檢測(cè)脈搏波波峰和波谷。因?yàn)镈SP系統(tǒng)內(nèi)設(shè)置的數(shù)據(jù)緩沖區(qū)寬度為1024點(diǎn)，只能存儲(chǔ)1024點(diǎn)的脈搏數(shù)據(jù)，我們只需要在這1024點(diǎn)數(shù)據(jù)中求出脈搏波的波峰(波谷)即可。由于脈搏波波峰(波谷)為其鄰域內(nèi)最大值(最小值)，同時(shí)為了避免潮波和重搏波的影響，通過(guò)對(duì)大量脈搏波進(jìn)行反復(fù)分析，可以得到此鄰域的最佳寬度為100點(diǎn)。于是，在數(shù)據(jù)緩沖區(qū)內(nèi)，取100點(diǎn)寬度為一個(gè)固定窗，設(shè)為W[j]，然后將此窗固定在數(shù)據(jù)緩沖區(qū)內(nèi)，讓數(shù)據(jù)流過(guò)此窗，求出此窗內(nèi)數(shù)據(jù)的最大值(最小值)，再根據(jù)脈搏信號(hào)波峰(波谷)判決條件判斷此最大值(最小值)是否為波峰(波谷)。設(shè)脈搏信號(hào)數(shù)據(jù)序列為X[i]，其中i∈[0,1024]，由上面分析可得到在波峰鄰域內(nèi)波峰X[ip]的判決條件為：

式中max{W[j]}為固定窗內(nèi)數(shù)據(jù)最大值，ip為波峰位置。由于最大值鄰域?qū)挾鹊南拗?，ip的取值不能小于50，不能大于1153，由此決定j的取值上限和下限。j的取值決定著固定窗在數(shù)據(jù)緩沖區(qū)內(nèi)的位置，一旦 j值取定，就可以在數(shù)據(jù)緩沖區(qū)內(nèi)確定一個(gè)長(zhǎng)度為100點(diǎn)的固定窗。而k的取值就是為了限制j的取值，使其滿足要求，k可以為其定義域內(nèi)的任意整數(shù)。與波峰判決條件類似，可得波谷X[il]的判決條件為：

由上面方法求得的相鄰波峰與波谷幅度差便可獲得波形高度，即一個(gè)周期的波形高度是由此周期內(nèi)的波峰幅值減去其前面最鄰近的波谷幅值。

2.1.2 脈率計(jì)算

脈率是指每分鐘脈搏的次數(shù)，正常情況下與心率一致，健康成年人在安靜狀態(tài)下脈率為每分鐘60-100次。且脈率受諸如年齡、性別、體型等因素的影響，在一定的范圍內(nèi)波動(dòng)。而脈率的不同，其對(duì)應(yīng)的脈象也不同，反映著人體不同的生理和病理特征，故在臨床診斷中有必要對(duì)脈率進(jìn)行實(shí)時(shí)檢測(cè)。

檢測(cè)脈率的思路為：先檢測(cè)相鄰兩個(gè)脈搏波的波峰(波谷)值；然后求出兩個(gè)脈搏波峰(波谷)之間的數(shù)據(jù)點(diǎn)數(shù)(設(shè)為n)；由采樣率求出這n點(diǎn)所對(duì)應(yīng)的時(shí)間，即為脈搏信號(hào)的周期(設(shè)為T)；最后用60/T便可計(jì)算出實(shí)時(shí)脈率(設(shè)為N次/分)。具體檢測(cè)脈率的算法如下：

1) 脈搏信號(hào)峰值檢測(cè)并確定峰值位置。用式(11-1)(式(11-2))檢測(cè)出兩個(gè)連續(xù)周期波峰(波谷)及其位置，并求出波峰(波谷)之間的點(diǎn)數(shù)n。

2) 確定脈搏信號(hào)的周期T，計(jì)算脈率值。采樣頻率fs =400 Hz，則采樣周期為1/fs，故脈搏信號(hào)的周期為T=n/fs。根據(jù)定義，可求得脈率N=60/T=(60×fs)/n。

2.2 脈搏信號(hào)的頻域分析

雖然用時(shí)域法研究脈搏信息比較簡(jiǎn)單直觀，但用其在做脈象和脈勢(shì)等方面研究已經(jīng)不能完全滿足要求，有一定的局限性。脈搏信號(hào)有著自身獨(dú)特的頻域特征，這些頻域特征蘊(yùn)含著豐富的脈搏信息，求出脈搏信號(hào)的幅度譜十分必要。對(duì)信號(hào)進(jìn)行離散傅里葉變換（DFT），其基本變換公式為：

式（13）中，每項(xiàng)求和分別是原始序列的N/2個(gè)偶數(shù)點(diǎn)和N/2個(gè)奇數(shù)點(diǎn)的DFT。相比DFT，計(jì)算量極大的減少。可應(yīng)用式（13）對(duì)DSP系統(tǒng)數(shù)據(jù)緩沖區(qū)內(nèi)的128點(diǎn)脈搏數(shù)據(jù)進(jìn)行FFT變換，求出其頻譜圖。

3 基于DSP的動(dòng)態(tài)脈搏信號(hào)采集與處理

為了實(shí)時(shí)采集和處理動(dòng)態(tài)脈搏信號(hào)，以TI公司生產(chǎn)的TMS320VC5509 DSP芯片為核心，設(shè)計(jì)了動(dòng)態(tài)脈搏信號(hào)的實(shí)時(shí)檢測(cè)系統(tǒng)，如圖3所示。此系統(tǒng)主要由三個(gè)模塊組成：1）光電脈搏傳感器；2）脈搏信號(hào)調(diào)理電路；3）DSP實(shí)時(shí)檢測(cè)與處理系統(tǒng)。光電脈搏傳感器采集人體的生理脈搏信號(hào)，經(jīng)信號(hào)調(diào)理電路，被送入DSP系統(tǒng)進(jìn)一步做信號(hào)濾波、時(shí)域和頻域分析及結(jié)果顯示，實(shí)現(xiàn)脈搏信號(hào)的實(shí)時(shí)檢測(cè)和分析。

圖3 脈搏信號(hào)檢測(cè)系統(tǒng)結(jié)構(gòu)圖Fig.3 Architecture diagram of pulse signal detection system

3.1光電脈搏傳感器及脈搏信號(hào)的調(diào)理電路

光電脈搏傳感器主要完成生理脈搏信號(hào)向電信號(hào)的轉(zhuǎn)換。傳感器內(nèi)置一個(gè)近紅光（940 nm）發(fā)光二極管，用來(lái)提供光源；一個(gè)光敏二極管，用來(lái)感知光信號(hào)。由于動(dòng)脈血液對(duì)光的吸收量隨動(dòng)脈搏動(dòng)而變化，光敏二極管可將變化的光信號(hào)轉(zhuǎn)化為變化的電流信號(hào)，實(shí)現(xiàn)將指尖脈搏的波動(dòng)轉(zhuǎn)化為電流信號(hào)輸出。

脈搏信號(hào)的調(diào)理電路實(shí)現(xiàn)動(dòng)態(tài)脈搏信號(hào)的初步調(diào)理，先通過(guò)I/V轉(zhuǎn)換電路，將光電脈搏傳感器拾取的電流信號(hào)轉(zhuǎn)換為電壓信號(hào)；然后，通過(guò)一個(gè)截止頻率為0.16 Hz二階高通模擬濾波器，去除基線漂移；再由一個(gè)截止頻率為80 Hz的二階低通模擬濾波器濾波，降低肌電干擾。在濾波的同時(shí)，通過(guò)兩級(jí)放大器，將微弱的脈搏信號(hào)放大500倍，轉(zhuǎn)換為可供DSP內(nèi)ADC轉(zhuǎn)換器采樣的電壓信號(hào)。

3.2 DSP實(shí)時(shí)檢測(cè)與處理系統(tǒng)

脈搏動(dòng)態(tài)信號(hào)的檢測(cè)與處理系統(tǒng)其核心為TMS320VC5509DSP處理芯片，具有高處理速度(400MIPS)的定點(diǎn)內(nèi)核DSP芯片。自帶RAM(128 K×16 bit), 片上FLASH(4 M×16 bit)，2路10bit A/D轉(zhuǎn)換器，最大采樣率為12.5 KHz，采集電壓范圍-3V-+3V。該系統(tǒng)主要完成脈搏信號(hào)的AD轉(zhuǎn)換、整系數(shù)濾波、波高和脈率計(jì)算、頻譜計(jì)算，日期、時(shí)鐘和脈率顯示，并對(duì)脈率進(jìn)行監(jiān)控。當(dāng)脈率異常時(shí)，聲光報(bào)警。

圖4 系統(tǒng)程序流程圖Fig.4 Program fl ow chart

系統(tǒng)程序流程圖如圖4所示。設(shè)計(jì)思路為：采用DSP處理器自帶的定時(shí)器通過(guò)定時(shí)中斷的方式控制AD轉(zhuǎn)換器采集數(shù)據(jù)，每采集一個(gè)數(shù)據(jù)，就在定時(shí)中斷子程序內(nèi)完成一次數(shù)據(jù)處理。圖4(a)為程序主流程圖，主程序主要完成系統(tǒng)初始化、計(jì)算脈搏信號(hào)的頻譜、顯示數(shù)據(jù)處理結(jié)果。圖4(b)為AD轉(zhuǎn)換定時(shí)中斷子程序流程圖。本設(shè)計(jì)采取定時(shí)中斷方式控制AD轉(zhuǎn)換，采樣率為400 Hz。每次定時(shí)結(jié)束，定時(shí)器0便會(huì)向系統(tǒng)申請(qǐng)中斷，系統(tǒng)響應(yīng)中斷以后，開始運(yùn)行中斷子程序。中斷子程序主要完成工作為：首先啟動(dòng)AD轉(zhuǎn)換器開始轉(zhuǎn)換數(shù)據(jù)，讀取并存儲(chǔ)數(shù)據(jù)；然后對(duì)數(shù)據(jù)緩沖區(qū)的脈搏數(shù)據(jù)進(jìn)行陷波、低通和包絡(luò)濾波，檢測(cè)脈搏信號(hào)的峰值和谷值，計(jì)算波高和脈率，并判斷脈率是否異常，若脈率異常，發(fā)出聲光報(bào)警；最后返回主程序。

3.2 濾波算法在DSP上的編程簡(jiǎn)化

在實(shí)際應(yīng)用中，影響系統(tǒng)程序運(yùn)行速度除硬件設(shè)備的因素外，另一主要的因素就是算法內(nèi)乘法的運(yùn)算次數(shù)。雖然DSP處理器自帶硬件乘法器，但是采取有效的措施簡(jiǎn)化程序，減少運(yùn)算過(guò)程中乘法次數(shù)和循環(huán)次數(shù)，對(duì)于實(shí)現(xiàn)信號(hào)的實(shí)時(shí)處理和充分利用DSP內(nèi)部的硬件資源，有十分重要的意義。

如圖4(b)所示，在AD轉(zhuǎn)換的子流程圖中，要在兩個(gè)采樣點(diǎn)之間絕對(duì)實(shí)時(shí)地實(shí)現(xiàn)信號(hào)的陷波濾波、低通濾波、波高計(jì)算和脈率計(jì)算等，計(jì)算量非常大。為了解決這個(gè)問(wèn)題，我們對(duì)編程方法作如下改進(jìn)：

1）將乘法和除法運(yùn)算改為移位運(yùn)算。如式(6)所示，運(yùn)算中的乘法全為×2運(yùn)算，可優(yōu)化為左移運(yùn)算。

2）采用單點(diǎn)濾波。在DSP處理系統(tǒng)中，每采集一個(gè)數(shù)據(jù)點(diǎn)，就要對(duì)數(shù)據(jù)緩沖區(qū)內(nèi)的所有數(shù)據(jù)進(jìn)行移位，進(jìn)行數(shù)據(jù)的更新，然后再對(duì)數(shù)據(jù)進(jìn)行濾波。若采用全數(shù)據(jù)濾波，由式(3)和式(9)可知，每次濾波都需要對(duì)系統(tǒng)的群時(shí)延(206個(gè)數(shù)據(jù)點(diǎn))進(jìn)行初始化。為了解決這個(gè)問(wèn)題，每次在濾波結(jié)束時(shí)都要保留206個(gè)數(shù)據(jù)，用來(lái)進(jìn)行下次濾波的初始化，于是每次參與濾波的數(shù)據(jù)為1024+206=1230個(gè)。這不僅增加了乘法和循環(huán)的次數(shù)，而且還占用了多余的存儲(chǔ)空間。若采用單點(diǎn)濾波，即每采集一個(gè)數(shù)據(jù)點(diǎn)，只對(duì)此數(shù)據(jù)點(diǎn)進(jìn)行濾波運(yùn)算，如此所有的濾波過(guò)程只需要初始化群時(shí)延一次，同時(shí)降低了計(jì)算量和存儲(chǔ)空間的占有量。

4 結(jié)果

圖5(b)為經(jīng)過(guò)陷波和低通濾波以后的脈搏信號(hào)圖，對(duì)比圖5(a)原始脈搏信號(hào)可以看出，經(jīng)過(guò)陷波器和低通濾波器后，脈搏信號(hào)中的高頻干擾得到了有效的抑制。圖5(c)為經(jīng)包絡(luò)濾波后的脈搏信號(hào)，其幅值由400 mV～750 mV降為-158 mV～158 mV，基線漂移得到了有效地去除。 經(jīng)過(guò)濾波處理后的脈搏信號(hào)，基本上能滿足實(shí)時(shí)信號(hào)處理要求。圖6為將脈搏信號(hào)進(jìn)行FFT變換后的脈搏信號(hào)頻譜圖，由圖看出脈搏信號(hào)的低頻和高頻干擾得到了有效地抑制，方便脈搏信號(hào)的頻域特征分析。

圖5 脈搏信號(hào)預(yù)處理Fig.5 Pulse signal preprocessing

5 結(jié)論

系統(tǒng)運(yùn)行結(jié)果表明，所設(shè)計(jì)的整系數(shù)濾波和包絡(luò)濾波算法在保存原始有用信息的前提下，有效地抑制了動(dòng)態(tài)脈搏信號(hào)的噪聲干擾。在此基礎(chǔ)上，分別從時(shí)域和頻域?qū)λ杉拿}搏信號(hào)進(jìn)行了實(shí)時(shí)處理，計(jì)算出了脈搏信號(hào)的波形高度和脈率，以及FFT變換求得其頻譜圖。通過(guò)特征分析，實(shí)時(shí)顯示人體狀態(tài)和參數(shù)，并可實(shí)現(xiàn)異常狀態(tài)報(bào)警。該系統(tǒng)在健康監(jiān)護(hù)領(lǐng)域具有較強(qiáng)的實(shí)用性和應(yīng)用前景。

本課題組對(duì)脈搏信號(hào)中的高頻信號(hào)有了一定的研究。為了保存脈搏信號(hào)中的有用高頻信號(hào)， DSP處理系統(tǒng)采用400 Hz作為系統(tǒng)的采樣率，實(shí)時(shí)完成了脈搏信號(hào)的檢測(cè)和處理。本設(shè)計(jì)所提出的濾波算法和脈搏信號(hào)的處理算法，同樣可以在低采樣率的系統(tǒng)中實(shí)現(xiàn)動(dòng)態(tài)脈搏信號(hào)的實(shí)時(shí)檢測(cè)和處理。

圖6 脈搏信號(hào)頻譜圖Fig.6 Spectrum of pulse signal

[1] 于洋, 劉靜. 手機(jī)無(wú)線心電監(jiān)測(cè)技術(shù)系統(tǒng)的實(shí)現(xiàn)及性能測(cè)評(píng)[J].中國(guó)醫(yī)療器械雜志, 2010, 32(5): 391-395.

[2] 孟淑婷，王磊. 基于DSP的家庭健康監(jiān)護(hù)儀的設(shè)計(jì)[J]. 國(guó)際生物醫(yī)學(xué)工程雜志, 2010, 22(9): 280-282.

[3] Woosik Shin, Yong Dae Cha, Gilwon Yoon. ECG/PPG Integer signal processing for a ubiquitous health monitoring system[J].Journal of Medical Systems, 2009, 34: 891-898.

[4] Cheng Wen, Ming-Feng Yeh, Kuang-Chiung Chang, et al. Realtime ECG telemonitoring system design with mobile phone platform[J]. Measurement, 2007, 41(2008): 891-898.

[5] Sera fi m Tabakov, Ivo Iliev, Vessela Krasteva. Online digital fi lter and qrs detector applicable in low resource ECG Monitoring Systems[J]. Annals of Biomedical Engineering, 2008, 36(11):1805-1815.

[6] E Pinheiro, O Postolache. Non-Intrusive device for real-time circulatory system assessment with advanced signal processing capabilities [J]. Measurement Science Review, 2010, 10(5): 166-175.

[7] 張珣. 光電脈搏傳感器的設(shè)計(jì)與改進(jìn)[J]. 中國(guó)醫(yī)療器械雜志,2009, 33(5): 344-346.

[8] 楊穎飛. 強(qiáng)背景噪聲下的脈搏信號(hào)處理算法研究[D]. 西安：西安電子科技大學(xué), 2009.

[9] Jinseok Lee, Ki H Chon. An autoregressive model-based particle fi ltering algorithms for extraction of respiratory rates as high as 90 breaths per minute from pulse oximeter[J]. IEEE Transactions on Biomedical Engineering, 2010, 57(9): 2158-2167.

[10] Woosik Shin, Yong Dae Cha, Gilwon Yoon. ECG/PPG Integer Signal Processing for a Ubiquitous Health Monitoring System[J].Journal of Medical Systems, 2009, 34: 891-898.

[11] Kristjan Pilt, Kalju Meigas, Rain Ferenets, et al.Photoplethysmographic signal processing using adaptive sum comb filter for pulse delay measurement[J]. Estonian Journal of Engineering, 2010, 16(1): 78-94

[12] K Shafqat, D P Jones, R m langford, et al. Filtering techniques for the removal of ventilator artifact in oesophageal pulse oximetry[J].Med Bio Eng Comput, 2006, 44:729-737.

[13] S Poornachandra, N Kumaravel. A novel method for the elimination of power line frequency in ECG signal using hyper shrinkage function [J]. Digital Signal Processing, 2007, 18(2008): 116-126.

[14] 蘭瑞芬, 胡廣書. 高采樣率下簡(jiǎn)單整系數(shù)工頻陷波器的設(shè)計(jì)[J].航天醫(yī)學(xué)與醫(yī)學(xué)工程, 2008, 21(2): 152-156.

[15] 趙洪亮, 卜凡亮, 黃鶴松, 等. TMS320C5X應(yīng)用系統(tǒng)設(shè)計(jì)[M]. 北京: 北京航空航天大學(xué)出版社, 2008