sampencandapen.m

来自「样本熵 结合样本熵和谐波小波分解」· M 代码 · 共 33 行

close all;
clear all;
clc;
m=2;
i=0;
a=1:0.003:25;
rr=rand(1,8001);
for r=1/100:1/100:100/100
    i=i+1;
    y=0.9;
    mix2=1+(1-y)*(sin(2*pi*a/12))+y*rr;
    mix2=mix2(end-499:end);
    sd = std(mix2);
    mix2 = (mix2-mean(mix2))/sd;
    [e2,A,B]=sampenc23(mix2,m,r);
    sampenc_mix2(i)=(e2(end));
end
for r=1/100:1/100:100/100
    i=i+1;
    y=0.9;
    mix2=1+(1-y)*(sin(2*pi*a/12))+y*rr;
    mix2=mix2(end-499:end);
    [e1,AA,BB]=fast_ApEn(mix2,r);
    fast_ApEn_mix2(i)=(e1(end));
end
b=1:100;
    plot(b,fast_ApEn_mix2(101:200),b,sampenc_mix2,'-r');
    set(gca,'Xscale','log');
xlabel('logistic parameter','fontsize',20)
ylabel('sample entropy','fontsize',20)
set(gca,'XTickLabel',{'0.01','0.1','1'});
set(gca,'FontName','Symbol');
set(gca,'FontSize',18);