qrs波检测程序.docxVIP

clear all;close all;points=1000; level=4; sr=360; %读入ECG信号%load ecgdata.mat;load ecgdata.txt;%ecgdata=mydata;plot(ecgdata(1:points));grid on;axis tight;axis([0,points,-2,5]);title(ECG信号);swa=zeros(4,points);swd=zeros(4,points);signal=ecgdata(0*1000+1:1*1000);%算小波系数和尺度系数for i=1:points-3 swa(1,i+3)=1/4*signal(i+3-2^0*0)+3/4*signal(i+3-2^0*1)+3/4*signal(i+3-2^0*2)+1/4*signal(i+3-2^0*3); swd(1,i+3)=-1/4*signal(i+3-2^0*0)-3/4*signal(i+3-2^0*1)+3/4*signal(i+3-2^0*2)+1/4*signal(i+3-2^0*3);endj=2;while j=level for i=1:points-24 swa(j,i+24)=1/4*swa(j-1,i+24-2^(j-1)*0)+3/4*swa(j-1,i+24-2^(j-1)*1)+3/4*swa(j-1,i+24-2^(j-1)*2)+1/4*swa(j-1,i+24-2^(j-1)*3); swd(j,i+24)=-1/4*swa(j-1,i+24-2^(j-1)*0)-3/4*swa(j-1,i+24-2^(j-1)*1)+3/4*swa(j-1,i+24-2^(j-1)*2)+1/4*swa(j-1,i+24-2^(j-1)*3); end j=j+1;end%画出原信号和尺度系数，小波系数figure;subplot(level,1,1); plot(ecgdata(1:points)); grid on;axis tight;title(ECG信号及其在j=1,2,3,4尺度下的尺度系数及小波系数);for i=1:level subplot(level+1,2,2*(i)+1); plot(swa(i,:)); axis tight;grid on;xlabel(time); ylabel(strcat(a ,num2str(i))); subplot(level+1,2,2*(i)+2); plot(swd(i,:)); axis tight;grid on; ylabel(strcat(d ,num2str(i)));end%画出原图及小波系数figure;subplot(level,1,1); plot(real(ecgdata(1:points)),b); grid on;axis tight;title(ECG信号及其在j=1,2,3,4尺度下的小波系数);for i=1:level subplot(level+1,1,i+1); plot(swd(i,:),b); axis tight;grid on; ylabel(strcat(d ,num2str(i)));end%**************************************求正负极大值对*****************************************%ddw=zeros(size(swd));pddw=ddw;nddw=ddw;%小波系数的大于0的点posw=swd.*(swd0);%斜率大于0pdw=((posw(:,1:points-1)-posw(:,2:points))0);%正极大值点pddw(:,2:points-1)=((pdw(:,1:points-2)-pdw(:,2:points-1))0);%小波系数小于0的点negw=swd.*(swd0);ndw=((negw(:,1:points-1)-negw(:,2:points))0);%负极大值点nddw(:,2:points-1)=((ndw(:,1:points-2)-ndw(:,2:points-1))0);%或运算ddw=pddw|nddw;ddw(:,1)=1;ddw(:,points)=1;%求出极值点的值,其他点置0wpeak=ddw.*swd;wpeak(:,1)=wpeak(:,1)+1e-10;wpeak(:,points)=wpeak(:,points)+1e-10;%画出各尺度下极值点figure;for i=1:level subplot(level,1,i); plot(wpeak(i,:)); a