[发明专利]一种测量随机数据非指数性的熵方法

说明书

技术领域

本发明属于随机数据挖掘领域，具体涉及一种测量随机数据非指数性的熵方法。

背景技术

目前，非指数性随机数据已广泛存在于水文、网络、生化、金融和土建等复杂系统。这类随机数据的最大特点是其概率分布或尾部分布不服从传统的指数分布或指数衰减。而随机数据的非指数统计性质很大程度上决定着整个系统的稳定性。因此，测量随机数据的非指数性是十分必要的。

现有测量随机数据非指数性的方法中，经验累积分布的尾部估计法是最常用的方法。该方法在双对数坐标下，描绘经验累积分布的右尾。如其右尾近似为一条直线，则被测量的随机数据具有非指数性。该方法简单且直观，但仅是一种定性的方法，而且样本的大小直接影响其精度。熵是一种测量随机系统不确定性的有力工具。目前，香农熵已应用于风险管理、肿瘤诊断和信息通信等领域。但是，香农熵依赖随机数据直方图分组区间的大小或需要计算复杂的积分，计算成本大且精度较低。

在国内外，已有多项专利技术涉及熵方法，如专利CN103886186A“一种确定钢管混凝土承载力设计误差分布的熵方法”，通过结合稳定分布的累计分布函数计算累计熵，选择误差分布的概型；CN102622649A“一种基于信息熵的改进进化多目标优化方法”，结合高斯函数，利用传统的香农熵进行目标优化；US20140189792“Method and system for electronic content storage and retrieval using Galois Fields and information entropy on cloud computing networks”，同样采用传统的香农熵用于云网络计算。

因此，需要一种新的测量随机数据非指数性的方法以解决上述问题。

发明内容

本发明的目的是针对现有技术中测量随机数据非指数性的缺陷，提供一种测量随机数据非指数性的熵方法。

为实现上述发明目的，本发明测量随机数据非指数性的熵方法可采用如下技术方案：

一种测量随机数据非指数性的熵方法，包括如下步骤：

1)、获取分析对象的随机数据；

2)、利用对数矩估计法计算步骤1)得到的随机数据对应的Mittag-Leffler分布的参数，Mittag-Leffler分布包括以下两个参数：

α^=π3(σ2-π2/12)]]>

其中，是Mittag-Leffler分布的稳定指数，是Mittag-Leffler分布的尺度参数，μ为随机数据X的对数变换log(X)的均值，σ为log(X)的标准差，为欧拉常数；

3)、计算随机数据的实际熵；

4)、利用步骤2)得到的随机数据对应的Mittag-Leffler分布计算随机数据的预测熵；

5)、结合步骤3)的实际熵和步骤4)的预测熵，计算相对熵指标，其中，相对熵指标R为预测熵与实际熵比的绝对值；

6)、当相对熵指标R小于等于0.01时，则表明随机数据具有非指数性。

更进一步的，步骤3)中利用分数阶矩法和/或对数矩法计算随机数据的实际熵。更进一步的，分数阶矩法通过下式计算随机数据的实际熵：