matlab - Matlab 中向量的谱熵和谱能量

Question

我将使用谱能量和谱熵作为基于窗口的时间序列数据的特征。但是，我对在线使用的公式有点困惑，尤其是关于特殊熵。

我使用了来自 Matlab 的熵，但这不适用于时间序列数据。它只是让我对一切都为零。 http://www.mathworks.nl/help/images/ref/entropy.html此版本的测试窗口熵结果 = 0

然后我用的就是这个版本。 http://www.mathworks.com/matlabcentral/fileexchange/28692-entropy此版本的测试窗口熵结果 = 4.3219

在将 imhist 应用于数据窗口后，我还尝试了 -sum(p.*log2(p)) (p = imhist(aw1(:));)。从联机帮助中获得此信息。
此版本的测试窗口熵结果 = 0.0369

他们都报告了不同的值(value)。

对于频谱能量，我使用的是 fft 系数的平方和。总和(绝对值(fft(数据窗口))。^2)

任何人都可以给我任何建议，哪个是正确的版本？

Answer 1

对于谱熵，步骤是:

• 计算功率谱密度(PSD)
• 规范化 PSD

• 计算熵 −∑(P)log2(P)，其中 P = PSD

  P=sum(abs(fft(data-window)).^2)
  %Normalization
  d=P(:);
  d=d/sum(d+ 1e-12);
  
  %Entropy Calculation
  logd = log2(d + 1e-12);
  Entropy(inc) = -sum(d.*logd)/log2(length(d));