fast_apen1.m

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

load aaa.txt;%采样周期1/2000秒，1024点
s=aaa(1:1024);
t=0:1023;x=sin(2*390*pi*t/2000);
ss=s+x;
zzz=hmapdn(s);
zz=hmapdn(ss);
for i=1:10
  aaa=zzz(i,:,:); aa=zz(i,:,:); 
  iii(i)=fast_ApEn(aaa,0.25*std(aaa));
    ii(i)=fast_ApEn(aa,0.25*std(aa)); 
end
subplot(221)
plot(iii(1:10));ylabel('ApEn');xlabel('level (Group Ⅰ)');
subplot(222)
plot(ii(1:10));ylabel('ApEn');xlabel('level (Group Ⅱ)');