HRV.rar_ECG_mother92p_neverezh_wavelet_小波 matlab

close all clear all tic % load ecg_1.mat load 69.mat % x0=block4(:,32); %x0=ecg2(1:14016); % x0=a1; point=length(x0); xx=-x0(1:point); figure(1) plot(xx); title('原始信号'); %intervel=260; for intervel=100:10:300 sp1=point/intervel; sp=ceil(sp1); s0=xx'; % % % N=9; % % % %wavelet denoise % % % % 用小波函数‘db5’对信号进行9层分解 % % % [c,l]=wavedec(data,N,'db8'); % % % c(1:l(1))=0; % % % c((sum(l(1:7))+1):length(c))=0; % % % s0=waverec(c,l,'db8'); % figure(2),plot(s0);%画出原始信号波形 lll=length(s0); % title('检出点');hold on; % % % Lo_D=[0.000,0.125,0.3750,0.3750,0.125,0.000]; % % % Hi_D=[-0.0061,-0.0869,-0.5798,0.5798,0.0869,0.0061]; % % % c=swt(s0,4,Lo_D,Hi_D); % % % d1=c(1,:); % % % d2=c(2,:); % % % d3=c(3,:); % % % d4=c(4,:); % % % %%d5(i)=c(5,i); % % % %%d6(i)=c(6,i); % % % figure(3); % % % plot(d4);%取尺度为2^4. % % % [RF,DF] = biorwavf('bior4.4'); %双正交样条小波滤波器，返回与指定双正交小波bior4.4相交联的两个尺度滤波器。rf为重构滤波器 df为分解滤波器 [Lo_D,Hi_D,Lo_R,Hi_R] = biorfilt(DF,RF);%双正交小波滤波器组 c=swt(s0,5,Lo_D,Hi_D);%平稳小波变换，lll191114猜测应该是滤波去噪 d1=c(1,:); d2=c(2,:); d3=c(3,:); d4=c(4,:); d5=c(5,:); %%d5(i)=c(5,i); %%d6(i)=c(6,i); % figure(3); % plot(d5);%250Hz取尺度为2^4.500Hz取尺度为2^5 a=0;b=0;m=zeros(sp,2);p=zeros(sp,2);q=1;x=zeros(1,sp);Q=zeros(1,sp);S=zeros(1,sp);P=zeros(1,sp);T=zeros(1,sp);%生成m*n全零矩阵 for n=1:sp max=0;min=0; for i=(intervel*(n-1)+3):(intervel*n+2) if i<=point if d3(i)>max max=d3(i);a=i;%a记录最大值标号 end if d3(i)<0 if abs(d3(i))>abs(min) min=d3(i);b=i; end end else break; end end %找出寻找区间的最大，最小值 if abs(b-a)<20 m(q,1)=a;m(q,2)=max; p(q,1)=b;p(q,2)=min; end %求极值阈值。 if q<3 s1=100;s2=-100; else a1=0;a2=0; for e=1:2 a1=a1+m(q-e,2);a2=a2+p(q-e,2); end %a1=a1/2;a2=a2/2; s1=1/5*(max+a1);s2=1/6*(min+a2);% end %求斜率阈值。 if q<4 if q<2 slope=0.05; else slope=0.0188; %3/8*0.05 end else slope=0; for f=(q-3):(q-1) xi=abs((m(f,2)-p(f,2))/(m(f,1)-p(f,1))); slope=slope+xi; end slope=1/8*slope/3;% end %if RR>1.5T if q>4 RR=0; for g=(q-3):(q-1) RR=RR+(x(1,g)-x(1,g-1)); end RR=RR/3; else RR=400;%%修改RR间期（因为采样率不一样） end if q>1&&x(1,q)>0 T=abs(x(1,q)-x(1,q-1)); if T>1.5*RR s1=.05*s1;s2=0.5*s2;slope=0.5*slope; end end if m(q,2)>s1&&p(q,2)<s2&&abs((m(q,2)-p(q,2))/(m(q,1)-p(q,1)))>slope %%%mm=0.2; if m(q,1)>p(q,1) for e=p(q,1):m(q,1) if d3(e)<=0&&d3(e+1)>=0 if abs(d3(e))<abs(d3(e+1)) x(1,q)=e; break; else x(1,q)=e+1; break; end end end else for e=m(q,1):p(q,1) if d3(e)>=0&&d3(e+1)<=0 if abs(d3(e))<abs(d3(e+1)) x(1,q)=e; break; else x(1,q)=e+1; break; end %%%if abs(d3(e))<mm %%% x(1,q)=e; %%% mm=abs(d3(e)); end end end if x(1,q)~=0 ma0=0;x2=x(1,q); for e=(x2-20):(x2+20) if e<=point if s0(e)>ma0 ma0=abs(s0(e)); x(1,q)=e; end else break; end end if q>1 %对有检重的进行修正。 if x(1,q)==x(1,q-1)||abs(m(q,1)-m(q-1,1))<100 %一分钟300次对应的最少间隔点数 q=q-1; end end figure(2); stem(x(1,q),s0(x(1,q)),'m.'); q=q+1; end % if x(1,q)~=0 end %满足阈值。 %-------------------------------------- end q=q-1; %for语句的end diffx=diff(x(1:q)); meanx=mean(diffx); stdx=std(diffx); figure(4) plot(diffx); end toc;