移动窗口中值滤波

% clc

clear

load('SignalMIT1.mat');

%设置窗函数

W=101;

for j=1:size(SignalMIT1)

ecg1=SignalMIT1(j,:);

%样本长度

L=length(ecg1);

%延拓ecg信号

for i=1:L+W

if i>=1 & i<=(W+1)/2

ecg2(i)=ecg1(1);

else if i>(W+1)/2 & i<=L+(W+1)/2

ecg2(i)=ecg1(i-(W+1)/2);

else if i>L+(W+1)/2 & i<=L+W-1

ecg2(i)=ecg1(L)

end

end

end

end

%窗口内ecg信号进行中值滤波

for i=1:L

BL(i)=median(ecg2(i:i+W-1));

end

%去除基线漂移

for i=1:L

ecg3(i)=ecg1(i)-BL(i);

end

SignalMIT2(j,:)=ecg3;

end

save('SignalMIT2.mat','SignalMIT2');