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

CC方法求时间延迟程序设计.docCC方法求时间延迟程序设计
CC方法求时间延迟程序设计
2 / 4
芃2
虿
CC方法求时间延迟程序设计
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, }
程序输出为一个 t l的矩阵
4）：
2
(r
dij
),r0
关系积分定义为：
C(m,N,r,t)
1)1
M(M
ijM
CC方法求时间延迟程序设计
CC方法求时间延迟程序设计
2 / 4
芃2
虿
CC方法求时间延迟程序设计
此中：dijXi
Xj
(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
函
数计算C(m,N,r,t);