第04课 matlab心电信号R波检测

《生物医学信号处理》实习报告学生姓名：学号：实验室名称：到的MIT-BIH数据）波段及特征点所代表的生理学意义将在下一节中进行较为详细的"缩所算法首先对心电信号作连续小波变换并对信号按照心动周期进行具体算法实现步骤如下：面上以灰度图的形式表示。中部。极值线。EMD分解:分解得到的第一、二、三本征模函数分量直接相加重构得到一个小波基的选取算法递归计算出信号的小波变换。分小波检测算法需通过检测小波模极大值对的过零点来定位信号falsepositive)结果；另外，当心律失常或QRS波幅度变小，阈值设置falsenegative)结果。由于固定阈编写的源程序：具体程序如下：clearall;clc;z=textread('ECG.txt');input=ECG(1:256);rate=ECG(100);sig=input;lensig=length(sig);wtsig1=cwt(sig,6,'mexh');lenwtsig1=length(wtsig1);wtsig1(lenwtsig1-20:lenwtsig1)=0;y=wtsig1;yabs=abs(y);%?sigtemp=y;siglen=length(y);sigmax=[];fori=1:siglen-2if))sigmax=[sigmax;abs(sigtemp(i+1)),i+1];end;end;figure(1);subplot(2,1,1),plot(sig);title('ECG原信号');subplot(2,1,2),plot(wtsig1);title('ECG变换信号');thrtemp=sort(sigmax);thrlen=length(sigmax);thr=0;fori=(thrlen-7):thrlenthr=thr+thrtemp(i);end;值点的平均值zerotemp=sort(y);zerovalue=0;zerovalue=zerovalue+zerotemp(i);end;zerovalue=zerovalue/100;%最小幅度平均值，对消幅度，100thr=(thrmax-zerovalue)*0.3;%最大、最小幅度的差值的30%为判%定位R波rvalue=[];fori=1:thrlenifsigmax(i,1)>thrrvalue=[rvalue;sigmax(i,2)];end;end;rvalue_1=rvalue;lenvalue=length(rvalue);i=2;whilei<=lenvalueif(rvalue(i)-rvalue(i-1))*rate<0.4ifyabs(rvalue(i))>yabs(rvalue(i-1))elseend;lenvalue=length(rvalue);i=i-1;end;i=i+1;end;lenvalue=length(rvalue);fori=1:lenvalueif(wtsig1(rvalue(i))>0)k=(rvalue(i)-5):(rvalue(i)+5);[a,b]=max(sig(k));rvalue(i)=rvalue(i)-6+b;elsek=(rvalue(i)-5):(rvalue(i)+5);[a,b]=min(sig(k));rvalue(i)=rvalue(i)-6+b;end;end;figure(2);subplot(2,1,1),plot(1:lensig,wtsig1,rvalue_1,wtsig1(rtitle('ECG纠正及校准前的R波信号');subplot(2,1,2),plot(1:lensig,sig,rvalue,sig(rvalue),'title('ECG纠正及校准后的R波信号');NOISYECG-R波检测：clearall;clc;z=textread('NOISYECG.txt');input=ECG(1:256);rate=ECG(100);sig=input;lensig=length(sig);wtsig1=cwt(sig,6,'mexh');lenwtsig1=length(wtsig1);wtsig1(lenwtsig1-20:lenwtsig1)=0;y=wtsig1;yabs=abs(y);%?sigtemp=y;siglen=length(y);sigmax=[];fori=1:siglen-2if))sigmax=[sigmax;abs(sigtemp(i+1)),i+1];end;end;figure(1);subplot(2,1,1),plot(sig);title('NOISYECG原信号');subplot(2,1,2),plot(wtsig1);title('NOISYECG变换信号');thrtemp=sort(sigmax);thrlen=length(sigmax);thr=0;fori=(thrlen-7):thrlenthr=thr+thrtemp(i);end;值点的平均值zerotemp=sort(y);zerovalue=0;zerovalue=zerovalue+zerotemp(i);end;zerovalue=zerovalue/100;%最小幅度平均值，对消幅度，100thr=(thrmax-zerovalue)*0.3;%最大、最小幅度的差值的30%为判%定位R波rvalue=[];fori=1:thrlenifsigmax(i,1)>thrrvalue=[rvalue;sigmax(i,2)];end;end;rvalue_1=rvalue;lenvalue=length(rvalue);i=2;whilei<=lenvalueif(rvalue(i)-rvalue(i-1))*rate<0.4ifyabs(rvalue(i))>yabs(rvalue(i-1))elseend;lenvalue=length(rvalue);i=i-1;end;i=i+1;end;lenvalue=length(rvalue);fori=1:lenvalueif(wtsig1(rvalue(i))>0)k=(rvalue(i)-5):(rvalue(i)+5);[a,b]=max(sig(k));rvalue(i)=rvalue(i)-6+b;elsek=(rvalue(i)-5):(rvalue(i)+5);[a,b]=min(sig(k));rvalue(i)=rvalue(i)-6+b;end;end;figure(2);subplot(2,1,