[转]关联维数(G-P算法)
关联维数不仅是我们在相空间重构过程中求解嵌入维数的一种重要方法，而且在机械故障诊断中发挥了重要的作用。关联维数计算的方法目前最主要的就是GP算法了，虽然目前也提出来一些其他的方法，但大都是基于GP算法的改进，譬如在GP算法前加入迭代奇异值降噪，减少噪声的影响，用小波变换方法计算关联维数，逐步迭代方法计算关联维数等等。我这里还是以GP算法为主要的介绍目标。分为几个内容：理论基础、GP算法求解步骤、程序中的参数理解和设置。
1 GP算法理论基础
设点X1,X2,...,XN为相空间内吸引子上的点，用Br(Xi)表示以参考点Xi为中心、半径为r的球形盒子，盒子的形状对维数的计算不会产生影响，盒子Br(Xi)的概率测度见图中的公式（1）
由此可得关联维数的定义式（2）。
由式（5）可知，当r->0时，关联积分C(r)和r之间存在标度关系，因此做出LnC(r)对Lnr的变化图，曲线斜率k即为关联维数D2。
2 GP算法的计算步骤
（1）利用时间序列X1,X2,...,XN，先给定一个较小的嵌入维数m0，重构相空间，得到新的序列{Yi}
（2）计算关联积分C(r)
（3）对于r的某个取值范围，吸引子的维数d与累积分布函数C(r)应满足对数线性关系，即d(m)=LnC(r)/Lnr，从而可用最小二乘拟合得到对应于m0的关联维数估计d(m0)
（4）增加嵌入维数m0，重新计算步骤（2）和（3），知道相应的维数估计值d(m)不再随着m的增加而在一定误差范围内不变为止。
————以上参考自吕金虎《混沌时间序列及其应用》
3 参数选取
一般计算过程中，可选定时间延迟、嵌入维数，做出LnC(r)~Lnr的曲线，选取线性度较好的一段作线性拟合，计算出拟合直线的斜率值，即得到关联维数。
（1）时间延迟——可根据自相关方法》椒扑愕玫
（2）嵌入维数，可先选定一个较小区域的m，看计算结果而定，如果曲线收敛较好，则进行拟合求得关联维数，如果不行，则继续定义m的区域
（3）r的取值，这个取值应限定在重构相空间中各点间最小距离与最大距离之间，即dmindmax，这点在GP程序中已有体现，那么程序中输入的ss是什么呢？这个ss其实就是r的序列的长度，一般来说，r的长度不能少于20。
r的取值有好几种方法：等比例取值、等间距取值、等指数比例取值等等，在程序中就是最简单的等距离取值。
（4）序列长度N，不能太少，一般要求N>=1500
转自：http://blog.sina.com.cn/s/blog_132b778440102vpd7.html
页:
[1]