matlab求振动信号信息熵？

题主给出的用matlab求振动信号信息熵的代码有下列几个问题：

1、shannon=-sum(x.*log2(x))，该代码不能写在for循环中，因为sum是求和命令，其要求x为数值向量

2、循环语句中 x=b(i)/sum(b) ，该代码最后的得到结果是一个x数值，并没有得到一消兆组数值向量，所以该循环应这样来写

for i=1:3072

x(ⅰ)=b(i)/sum(b)

end

3、为了帮你理解，你可运行下列代码，就可以发现，茄庆他们的区别颤桥握。

A=1:10

for i=1:10,x=A(i)/sum(A)end,x

for i=1:10,x(i)=A(i)/sum(A)end,x

内容：排列熵是一种基于复杂性量度的非线指森迟性动力学参唯李数,能够快速、有春梁效地反映系统的特征.目前在心电信号检测研究中还少有此类方法的应用,通过对算法进行研究,并应用于心室病症的检测中,利用MIT-BIH数据库,对心室病症心...

如何计算一个信号的多尺度熵 Multiscale Entropy

1、多尺度熵，是从样本熵演变过来的

2、样本熵是基于近似熵(ApEn)的一种用于度量时间序列复杂性的改进方法；

SampEn has two advantages over ApEn: data length independence and a relatively trouble-free implementation.

3、首先理解一下什么是近似熵：

（1）我有个疑问，时间t怎么定义啊；明明k的个片段都用上了，还怎么区分不同的时刻啊。

（2）会不会有一种可能，整个时间序列可以得到一个近似熵的值。也就是一个值。

（3）给一个时间序列后，如果这个时间序列是一维的，你只能得到一个值。

（4）如果这个时间序列是多维的，你可以把某个时刻，这些不同维度的值，组织在一起，计算得到一个熵。然后你就能得到熵随时间的演化了；

---