13QRS實(shí)驗(yàn)?zāi)縌RSQRSQRSQRSQRS實(shí)驗(yàn)原SS能計(jì)算心率并進(jìn)行心率變異分析，才能檢測(cè)ST段的參數(shù)和分析心電的其它細(xì)節(jié)信息。S波檢測(cè)包括QS少研究成果但是各種方法均有其不足之處因?yàn)樾碾娦盘?hào)波形的復(fù)雜性各種類(lèi)型噪聲()波的精確檢測(cè)有很大。QRS波檢測(cè)具有處QRS波的位置，很QRS波寬度和面積等參量。QRS波檢測(cè)可以靈活選擇調(diào)節(jié)各類(lèi)參數(shù)，對(duì)復(fù)雜情況進(jìn)行判濾波之后的信號(hào)經(jīng)過(guò)一些變換以提高QRS波的份量，進(jìn)而采用一系列閾值進(jìn)行P、T波的存在，若其濾波后超過(guò)其閾值便會(huì)產(chǎn)生假陽(yáng)性(FP，falsepositive)另外，當(dāng)心律失?；騋RS波幅度變小，閾值設(shè)置過(guò)高，會(huì)導(dǎo)致漏檢產(chǎn)生假幅度和寬度識(shí)別QRS復(fù)波如圖3-1所示他們實(shí)時(shí)算法包括以下環(huán)節(jié)：采樣的心電信號(hào)經(jīng)過(guò)數(shù)字帶通濾波器(5-15Hz)；微分后再平方，平方過(guò)程可突出微分后的頻響曲線(xiàn)的斜率，并有助于減少由于T波具有高頻分量而引起的假陽(yáng)性；再做移動(dòng)窗口(150ms)QRS波群的斜率，又包含QRS波群的寬度的信息，然后應(yīng)用兩種閾值進(jìn)行檢測(cè)，一是法還利重閾值法來(lái)尋找丟失的QRS波，以減少假，其方法是對(duì)上述每種閾值均設(shè)兩個(gè)閾值電平，其中一個(gè)是另一個(gè)的一半，當(dāng)發(fā)現(xiàn)可能丟失QRS波國(guó)MIT/BIH各種心律失常心電圖帶信號(hào)進(jìn)行實(shí)測(cè)檢驗(yàn)信號(hào)總計(jì)長(zhǎng)達(dá)24小時(shí)，48個(gè)，116173個(gè)心拍數(shù)，假陽(yáng)性為507，假為277，檢測(cè)的正確率為99.3%，其中波形形態(tài)變化較大者易導(dǎo)致假和假陽(yáng)性結(jié)果。

132132

實(shí)驗(yàn)內(nèi)QRS（一）QRSQRS(1z6H(z)(1z1y(n)2y(n1)y(n2)x(n)2x(n6)x(n11Hz5361z32Hlp(z)1

y(n)y(n1)x(n)x(n32155個(gè)采樣點(diǎn)

z16Hlp

y(n)x(n16)

1[y(n1)x(n)x(n15Hz161（提示：利用已知的傳遞函數(shù)實(shí)現(xiàn)濾波clear num=zeros(1,13);den=num(1)=1;num(7)=-2;num(13)=den(1)=1;den(2)=-2;den(3)=1;DATAafterLowfilter(num,den,OriginalData)/36;0 clearclearnum=zeros(1,33);den=num(1)=-1;num(17)=32;num(18)=-32;num(33)=1;den(1)=32;den(2)=-32;DATAafterHigh=filter(num,den,DATAafterLow);觀察一段采樣的ECG原始信號(hào)0 帶通濾波后信號(hào)0 QRS5~15HzP波，T50HZQRS（二）QRS為H(z)0.1(2z1z32z4num=zeros(1,5);den=zeros(1,5);num(1)=0.2;num(2)=0.1;num(4)=-0.1;num(5)=-DATAafterDif=filter(num,den,DATAafterHigh)*10;觀察采樣的ECG（三）信號(hào)逐點(diǎn)平方是一個(gè)非線(xiàn)性操作，方程為y(n)[x(n)]2DATAafterSqu=DATAafterDif.*DATAafterDif;fftDATAafterDif=fft(DATAafterDif);fftDATAafterSqu=fft(DATAafterSqu);n=(0:999)/999.*sfreq;title('平方后頻譜濾波后信號(hào)0 微分后信號(hào)0 （四）單靠RQRSQRS分抽取了R波斜率之外的特征。QRS復(fù)波寬度相近；QRSTime150;%nIntefloor(QRSTime/1000*sfreq);%forii=nInte:length(DATAafterSqu))=holdon;%將系統(tǒng)延時(shí)抵消0 實(shí)驗(yàn)結(jié)果分29采樣點(diǎn)系統(tǒng)延遲后，在上圖上顯示，可以看到積分QRSQRS復(fù)波稍寬。積分波形最高極值點(diǎn)即為R波波峰，下面的工作是要找到符合要求的積分曲線(xiàn)（五）QRST波、肌肉噪聲等的峰值。用QRS檢測(cè)算法的一組閾值是使信號(hào)峰值能被檢測(cè)到的那一組閾值。用于移動(dòng)SPKI＝0.125PEAKI+0.875SPKIPEAKI是信號(hào)峰值。NPKI＝0.125PEAKI+0.875NPKIPEAKI是噪聲峰值。THRESHOLDI2=0.5PEAKISPKI是信號(hào)峰值的實(shí)時(shí)估值；NPKI是噪聲峰值的實(shí)時(shí)估值；THRESHOLDI1是應(yīng)用的第一閾值；THRESHOLDI2值水平在信號(hào)首次分析期間就超過(guò)THRESHOLDI1，那么它是一個(gè)QRS峰。如QRS峰的信號(hào)峰值應(yīng)該超過(guò)THRESHOLDI2。QRS2SPKIPEAKImax(DATAmoveWinInte);初始化，使PEAKISPKI=PEAKI; PEAKL=0;NPKI=0; );%取大于最值0.05倍的曲線(xiàn)極大值QRSk=1;Rk=1;Nk=stepTime %step=floor(sfreq*stepTime/1000);%回檢步長(zhǎng)stepTimesTemp=1;STEPTIMES=5; forii=1:length(Pks)SPKI=0.125*PEAKI+0.875*SPKI;NPKI=0.125*PEAKL+0.875*NPKI; 測(cè)ifPks(ii)>THRESHOLDI1 =%%PksR(Rk)= locR(Rk)= PEAKI=Pks(ii);Rk=stepTimesTemp=QRSkQRSk+1;%SPKI=0.25*PEAKI+0.8*SPKI;%if %PEAKI=Pks(ii);QRSk=%%PksR(Rk)= locR(Rk)= PEAKI %確定為RRk=stepTimesTemp=1;

PEAKL %Nk=holdon;holdon;title('QRS起點(diǎn)標(biāo)示,紅色為R波波峰,未調(diào)整圖像');與原始的ECG綠色為QRS起點(diǎn)標(biāo)示,紅色為R波波峰,未調(diào)整圖像0 （六）次心跳的平均值。另一個(gè)均值RRAVERAGE2，是相鄰8次RR間隔均落入一定RRnRRAVERAGE10.125(RRn7RRRRnRRn是落入下述界限內(nèi)的RR間隔。RRLOWLIMIT92RRAVERAGE2RRHIGHLIMIT116RRAVERAGE2RR