时间序列相空间重构的及其应用研究.docVIP

时间序列相空间重构及其应用研究 摘要 时间序列的重构分析是从产生该序列的系统特性的角度提取该时间序列的特征量,在这种分析方法的应用过程中,关联积分和关联维的正确、快速计算是重要的第一步.本文对混沌时间序列相空间重构中最佳延迟时间间隔和嵌入维数的选取方法作了综述, 基于时间序列分析的方法,提出了一种神经网络时间序列预测及建模方法.关键词 时间序列 ,相空间重构,延迟时间间隔, 关联维,神经网络1 引言混沌是一种低阶确定性的非线性动力系统所表现出来的非常复杂的行为,它对现代科学具有广泛而深远的影响,几乎覆盖了一切学科领域,尤其是在物理学、天体力学、数学、生物学、经济学等方面得到了广泛的应用.在对混沌时间序列的各种分析中,如混沌预测(prediction of chaos)。动力学不变量(dynamical invariants)的估计。混沌信号的诊断(detection of chaos)等,所要进行的第一步工作是要对混沌信号进行相空间重构.1981年Takens提出了相空间重构的延时坐标法,奠定了相空间重构技术的基础,这种方法用单一的标量时间序列来重构相空间,包括吸引子、动态特性和相空间的拓扑结构.现已成为最主要、最基本的相空间重构方法[1]. 分形维是用来描述混沌信号的一个重要参数,目前主要流行是基于GP算法的关联维提取算法。2 G.P算法的描述自从人们发现延迟时间 对重构相空间的重要之后,便开始了探索确定延迟时间 的方法,并取了显著的成效，相空间重构理论认为,要保证相空间重构的正确性,所选用的延迟时间 必须使重构相空间的各个分量保持相互独立,选择的延迟时间如果太大, 就混沌吸引子而言,由于蝴蝶效应的影响,时间序列的任意两个相邻延迟坐标点将毫不相关,不能反映整个系统的特性;而延迟时间选择过小的话,时间序列的任意两个相邻延迟坐标点又非常接近,不能相互独立,将会导致数据的冗余。.因此我们需要一种方法来选择恰当的 ,于是围绕这一条件便先后出现了用自相关函数和互信息来确定延迟时间 的方法[3]。自相关函数能够提供信号自身与它的时延之间由冗余到不相关比较这种的度量,一般取自相关函数值首次出现零点时的时延为所要确定的时间延迟。现描述如下:对于单变量时间序列x1, x2, x3,…, xn取延迟时间为 ,则其自相关函数为: （9）其中,n为时间序列点数, 为时间序列的平均值.延迟时间 的选取原则是让时间序列内元素之间的相关性减弱,同时又要保证时间序列包含的原系统的信息不会丢失.研究表明,当关联函数C的值第一次为0(或近似为0)对应的延迟时间 比较合适[4].4 关联维m的选取在关联维计算中，重构相空间维数的取值依赖于所分析和研究的关联维数，数据序列的关联维数不同，所选取的重构相空间维数亦不同。由单变量的时间序列重构相空间时,为了保证该相空间能包含原状态空间吸引子的特征, 关联维应该取得足够大.Takens在1980年证明了嵌入维数大小的嵌入定理: （10）其中:m---重构相空间维数;D2---原状态空间吸引子所处空间的关联维数该定理表明,嵌入空间的维数至少是吸引子维数的两倍,这时重构的相空间和原系统的状态空间拓扑等价.但是由于在一般情况下缺乏原动力系统的先验知识,选择m则具有随意性系统特征量饱和法，是使重构相空间维数由小到大变化,计算每一个重构相空间的系统特征量,如果特征量达到饱和,则饱和时的重构相空间维数就是所求的维数.一般情况下,确定性系统一般都会收敛到低维的相空间,关联维数将不再随嵌入维数的增加而发生明显的改变,因此可以取关联维数作为系统特征量,逐渐增大嵌入维数,观察关联维数是否达到饱和. 在实际的分形分析中，关联维数是所要求取的对象，并不知道其具体数值，因此，需首先估计出所求关联维的取值范围，从而得出重构相空间维数m的粗略估计值。在m粗略估计值的范围内对m取不同的值，然后分别求取系统的关联维数，当关联维数达到饱和时的m取值，即为重构相空间的实际维数。5 应用展望对于平稳的时间序列来讲,利用传统的AR、MA、ARMA等模型通常可获得较好的预报结果.而对混沌时间序列而言,即使模型对数据匹配的很好,有时也无法做出准确的预测,未来趋势会在性质上与原有时间序列趋势发生根本不同的变化.因此,对混沌时间序列的预测研究我们需另找出路. 时间序列的量往往伴随有大量的噪声，神经网络比其它方法更能容忍噪声。任何一个时间序列都可以看成是一个由非线性机制确定的输入输出系统，神经网络可以较好地揭示非线性时间序列在时延状态空间中的相关性,从而达到预测目的。具体说来，可以这样处理：(1)神经网络用于单变量时间序列预测设一个单变量时间序列x1, x2, x3,…, xn对它进行预测的前提是认为其未来值与其前面的m个值之间有着某种函数关系,可描述如下:xn +k=F