hmapdn.m

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

function A=hmapdn(f)
N=length(f);
n=round(log(N)/log(2));
a=hwtdn(f);
for j=1:2^(n-2)
  A(1,j)=abs(a(1))^2;
  A(n+1,j)=abs(a(2^(n-1)+1))^2;
end
for j=2:n
  for k=1:2^(j-2)
     for m=1:2^(n-j)
       A(j,(k-1)*2^(n-j)+m)=abs(a(2^(j-2)+k))^2+abs(a(N+2-2^(j-2)-k))^2;
     end
end
end