sampenc_data.m

样本熵 结合样本熵和谐波小波分解

close all;
clear all;
clc
load('original signals.mat');
for q=1:27
y_wear=wear(q,:);
y_sharp=sharp(q,:);
zz=hmapdn(y_sharp);
zzz=hmapdn(y_wear);
m=2;
r=0.1;
for i=1:13
aa=zz(i,:,:);
aa=aa(end-1023:end);
sd = std(aa);
aa = (aa-mean(aa))/sd;
aaa=zzz(i,:,:);
aaa=aaa(end-1023:end);
sd = std(aaa);
aaa = (aa-mean(aaa))/sd;
[e1,A,B]=sampenc(aa,m,r);
[e2,A,B]=sampenc(aaa,m,r);
sampenc_sharp(q,i)=e1(end);
sampenc_wear(q,i)=e2(end);
end
end
a=1:13;
for q=1:27
figure
plot(a,sampenc_sharp(q,:),a,sampenc_wear(q,:),'r');
ylabel('样本熵','fontsize',18);
xlabel('层数(Blue_sharp&Rad_wear)','fontsize',18);
set(gca,'FontSize',13, 'FontName','Symbol');
end
save('Sampenc_data.mat','sampenc_sharp','sampenc_wear');