引言
散度熵（Diversity Entropy，翻译可能不准确哈，请谅解）于2021年发表在中科院1区期刊 IEEE Transactions on Industrial Informatics（影响因子突破12）。散度熵首次提出并运用于机械设备的状态识别和故障诊断，从利用熵值表征不同的状态。而实际上，振动信号的本质是时间序列，那么散度熵可用于其他时间序列中，如：电能质量数据、风速、功率、声音、温度、交通、水流等等，您能想到的时间序列皆有可能。目前，基于散度熵的研究还不多。
基于熵值理论的分析方法可以不经过数据的分解或变换，直接度量数据的复杂度，可完成不同时间序列的分类。状态识别和故障诊断本质是一个分类问题，解决这类问题主要有两步：1.利用熵的方法从信号中提取信息，即信号的熵值为故障特征；2.基于提取的熵值，利用机器学习方法作为分类器进行识别。
本期推出散度熵及其5种多尺度熵结合多种机器学习方法进行状态识别和故障诊断。
散度熵（Diversity Entropy）多尺度散度熵（Multiscale Diversity Entropy）复合多尺度散度熵（Composite Multiscale Diversity Entropy）精细复合多尺度散度熵（Refined Composite Multiscale Diversity Entropy）时移多尺度散度熵（Time-shift Multiscale Diversity Entropy）层次散度熵（Hierarchical Multiscale Diversity Entropy）
散度熵的理论请参考原始论文：
X. Wang, S. Si and Y. Li, "Multiscale Diversity Entropy: A Novel Dynamical Measure for Fault Diagnosis of Rotating Machinery," in IEEE Transactions on Industrial Informatics, vol. 17, no. 8, pp. 5419-5429, Aug. 2021, doi: 10.1109/TII.2020.3022369.
5种多尺度的理论知识在中、英文期刊都有很多，可以参考其他论文，本期不在赘述。
02. 实操
数据集
本文使用故障识别和诊断任务中经典的凯斯西储大学CWRU数据集。随机选择了5类状态（当然可以加入更多的状态），选用DE端的振动数据，暂且命名为正常，故障1，故障2，故障3，故障4。设置每个样本长度为1024，每类状态共100个样本，总样本大小为：500*1024。即样本数*样本长度。
样本分割设置了重叠个数，可以随意选择重叠个数进行样本划分。非重叠即为重叠个数为0的情况，具体的原理如下：
非重叠：
重叠：
特征提取
计算每个样本的熵值（以上6种熵可随意切换），形成500*n的特征空间，将特征输入到后续的分类器中。对于散度熵而言，n=1，对于多尺度熵而言，n=你设置的时间尺度。
分类器
将上述特征空间以一定的比例，划分了训练集、测试集，可非常方便地输入到任何分类器中。
暂时集成了3种常用的分类器（可随意切换，也可自行添加其他分类器）：
-
极限学习机（Extreme Learning Machine，ELM）
-
支持向量机（Support Vector Machine, SVM）
-
决策树（Decision Tree，DT）
为了消除随机性带来的影响，使结果更具说服力，模型运行了多次，每次运行都随机选择新的训练集和测试集，以平均结果作为最后的结果。
诊断流程
故障识别、诊断流程图大致如下：
TSMDivEn-ELM结果展示：
不同状态对比
混淆矩阵及分类结果图
特征值 和 多次运行结果保存在mat中。可以方便使用matlab，python调用和查看这里考虑快速验证代码和算法，仅运行了2次，可根据实际需要修改。可以看到各状态的特征类间分离明显，类内高度聚合。
友情提示
1.理论上，该代码可适用于所有时间序列分类数据集，仅需将数据格式转换成：样本数N x 特征个数M。因此，可以浅显的理解：多分类问题实质是一个数据格式转换过程；
2.理论上，该代码可替换的部分只有：数据加载和参数设置。好处是：代码修改量极少，复用率高；
3.如有matlab基础语法和代码bug问题，请多看资料，百度或者GPT解决；
4.如有其他专业疑问，请多阅读相关文献。
Matlab代码下载
微信搜索并关注-优化算法侠，或扫描下方二维码关注，以算法名字搜索历史文章即可下载。
点击链接跳转：
cec2022测试函使用教程及matlab代码免费下载
绘制cec2017/018/2019/2020/2021/2022函数的三维图像教程，SO EASY！
175种群智能优化算法python库
求解cec测试函数-matlab
解决12工程设计优化问题-matlab
求解11种cec测试函数-python
解决12种工程设计优化问题-python
用于改进所有优化算法：21种混沌映射方法-混沌初始化（附matlab代码）
沙场大点兵：24种信号分解方法（附matlab代码）
沙场大点兵：27种一维数据转换成二维图像的方法-matlab代码沙场大点兵：27种一维数据转换成二维图像的方法-matlab代码
https://mp.weixin.qq.com/s?__biz=MzkxMDQ5MDk4Ng==&mid=2247486260&idx=1&sn=81b1970cb89364c0289ccdfb403e5388&chksm=c12be731f65c6e273a85456326b503b7f35d9f035405050932ff1926e0b1bfa8076b1bc2d1f2&token=25423484&lang=zh_CN#rd
3759
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
