shuxuexingtai.m

数学形态学进行ECG 滤波处理算法,完全原代码。不是MATLAB 自带的函数。

        y(kl)=x(j+kl-1)-k(kl);
        if(y(kl)<min)
            min=y(kl);
        end
    end
    y4(j)=min;
    j=j+1;
end

% (x open s) Close s

x=y3;
j=1;
for i=5:n-0     
    max=-10000;
    for kl=1:5
        y(kl)=x(j+kl-1)+k(kl);
        if(y(kl)>max)
            max=y(kl);
        end
    end
    yy3(j)=max;
    j=j+1;
end

x=yy3;
j=1;
for i=5:n-0         
    min=10000;
    for kl=1:5
        y(kl)=x(j+kl-1)-k(kl);
        if(y(kl)<min)
            min=y(kl);
        end
    end
    y5(j)=min;
    j=j+1;
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% (x close s) open s
x=y4;
j=1;
for i=5:n-0     
    min=10000;
    for kl=1:5
        y(kl)=x(j+kl-1)-k(kl);
        if(y(kl)<min)
            min=y(kl);
        end
    end
    yy4(j)=min;
    j=j+1;
end

x=yy4;
j=1;
for i=5:n-0     
    max=-10000;
    for kl=1:5
        y(kl)=x(j+kl-1)+k(kl);
        if(y(kl)>max)
            max=y(kl);
        end
    end
    y6(j)=max;
    j=j+1;
end


y7=(y5+y6)/2;

% PVE(A) = A-(A o B) . B
y8=x0-y5;

subplot(2,1,1);
plot(x0);

subplot(2,1,2);
plot(y7);