CC方法求时间延迟程序设计.doc

精品文档
精品文档
芃1
虿
精品文档
C-C方法求时间延迟程序设计
一、子程序设计：
本程序需调用四个子函数：
1）：用来求解Heaviside函数的值；
2）：用来进行精品文档
精品文档
芃1
虿
精品文档
C-C方法求时间延迟程序设计
一、子程序设计：
本程序需调用四个子函数：
1）：用来求解Heaviside函数的值；
2）：用来进行相空间重构；
（3）：用来将时间序列分拆为t个不有关的时间序列；
4）：用来计算时间序列的关系积分。
程序说明：
1）：
(x)
0,x
0
1,x
0
输入r和d值的，根据公式求得
heaviside函数的值sita并输出。
（2）
：根据给定的时间延迟
t和嵌入维数m重构m维的相空间
Y(ti)
[x(ti),x(ti
),x(ti2
),
x(ti
(m1))],i
1,2,,M
M
N(m1)
程序输出一个m
M的矩阵；
（3）
：将时间序列分红t个不相交的时间序列，长度
lN/t
{x1,xt1,x2t1,}
{x2,xt2,x2t2,}
{xt,x2t,x3t,}
程序输出为一个tl的矩阵
4）：
2
(r
dij
),r0
关系积分定义为：
C(m,N,r,t)
1)1
M(M
ijM
其中：dijXi
Xj
精品文档
精品文档
芃2
虿
精品文档
(x)
0,x
0
0
1,x
程序流程：
1、输入数据X：X为根据给定的时间延迟
t和嵌入维数
m进行相空间重构后获得
的mM的矩阵；
2、对相空间中所有点求距离，并保留在数组
d(i,j)中；
3、调用Heaviside函数计算所有的
Heaviside函数值并求和获得
sum_h
4、计算出关系积分的值，并输出
C_I。
二、主程序设计：
1、读入数据；
2、计算时间序列的标准差；
3、让r从sigma/2到2sigma变化，m从2到5变化，t从1到200
变化；
4
、调用disjoint
函数将时间序列分拆成t
个不相交的子列，并调用
correlation_integral函数计算C(1,N,r,t);
5
、调用reconstitution函数对子列进行相空间重构，并调用
correlation_integral
函