第04课 matlab心电信号R波检测 - 图文(3)

%定位R波 rvalue=[]; for i=1:thrlen if sigmax(i,1)>thr rvalue=[rvalue;sigmax(i,2)]; end; end; rvalue_1=rvalue; %排除误检，如果相邻两个极大值间距小于0.4，则去掉幅度较小的一个 lenvalue=length(rvalue); i=2; while i<=lenvalue if (rvalue(i)-rvalue(i-1))*rate<0.4 if yabs(rvalue(i))>yabs(rvalue(i-1)) rvalue(i-1)=[]; else rvalue(i)=[]; end; lenvalue=length(rvalue); i=i-1; end; i=i+1; end; lenvalue=length(rvalue); %在原信号上精确校准 for i=1:lenvalue if (wtsig1(rvalue(i))>0) k=(rvalue(i)-5):(rvalue(i)+5); [a,b]=max(sig(k)); rvalue(i)=rvalue(i)-6+b; else k=(rvalue(i)-5):(rvalue(i)+5); [a,b]=min(sig(k)); rvalue(i)=rvalue(i)-6+b; end; end; %打印纠正及校准前后的R波信号 figure(2); subplot(2,1,1),plot(1:lensig,wtsig1,rvalue_1,wtsig1(rvalue_1),'r.'); title('NOISYECG纠正及校准前的R波信号'); subplot(2,1,2),plot(1:lensig,sig,rvalue,sig(rvalue),'r.'); title('NOISYECG纠正及校准后的R波信号'); 结论（画出要求的图形）: ECG-R波检测所获得的结果如下图3-1与图4-1所示： 图3-1 图4-1 而NOISYECG-R波检测所获得的结果如下图5-1与图6-1所示： 图5-1 图6-1 总结： 实习报告分数： 指导教师：