混沌理论学习的总结_第1页
混沌理论学习的总结_第2页
混沌理论学习的总结_第3页
混沌理论学习的总结_第4页
混沌理论学习的总结_第5页
已阅读5页,还剩1页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、混沌理论学习的总结刘年东混沌的概念(名词解释):很难给出一个精确的定义,一般认为混沌就是指在确定系统中出现的一种貌似无规则的,类似随机的现象。混沌不是简单无序而是没有明显的周期和对称,但却是具有丰富的内部层次的有序结构,是非线性系统中一种新的存在形式。目前学者们一致的观点是:混沌是物质科学、数学科学和系统科学等多学科交义的边缘科学。混沌之所以有如此大的吸引力,因为它是一种关于过程的科学,而不是关于状态的科学;是关于演化的科学,而不是关于存在的科学,它使人们看到了运动演化中的生机和动力。混沌的主要特征:敏感依赖于初始条件伸长与折叠(从logistic映射形成混沌的过程可以看到,混沌具有伸长和折叠

2、的特性,这是形成敏感依赖于初始条件的主要机制。伸长是指系统内部局部不稳定引起的点之间的距离的扩大;折叠是指系统整体稳定所形成的点之间距离的稳定的限制经过多次的伸长和折叠,轨道被扰乱了,形成混沌)具有丰富的层次和门相似的结构(混沌区内有窗口(稳定的周期解),窗口里面还有混沌,这种结构无穷多次重复着,并具有各态历程和层次分明的特征。伸长和折叠使混沌运动具有大大小小的,无特征的各种尺度,这些都被称为自相似的结构非线性的耗散系统中存在混沌吸引子。(这是总体稳定和局部不稳定相结合的产物。通常的混沌吸引子都有负的(且无正的)Lyapunov指数,唯独混沌吸引子具有正的Lyapunov指数,而且混沌吸引子只

3、能用分数维來表征。)混沌行为最本质的特点是非线性系统对于初始-条件的极端敏感性。吸引子在一个确定耗散系统里,不属于其他任何更大的极限集,而且无轨道。在耗散系统里包含了拟周期吸引子、周期吸引子、混纯吸引子以及定常吸引子。几种典型的混沌吸引子图:iw*:*wLogistic混沌吸引子图Lorenz吸引子图Chen吸引子图混沌理论与其它预测技术的不同:传统的预测技术方法所基于的原理可归纳为惯性原理、相关原理、类推原理和概率推断原理等,这些理论对非线性问题难以给出很好的解释,乂由于传统预测模型多基于全局建模的角度,缺失局部特征的代表性,例如回归分析法、指数平滑法、趋势外推法等,通常预测精度较低,且存在

4、模型参数确定困难及普适性较差等缺点。虽然传统预测技术方法的应用在社会、经济和技术的发展中起了非常重要的作用,取得了积极的成果,但是对众多纷繁复杂的技术经济现象,现有预测理论与方法还不能给予合理的解释和有效的预测,预测结果常常与现实情况存在较大偏差。而混沌理论是非线性理论的重要组成部分,能够很好地描述非线性系统运动的变化规律,从而为预测模型的研究开辟了新思路。其对预测的影响不是对“原有理论与方法的修修补补”,而是要求根本性的变革,即进行预测范式的转换。判断系统是否具有混沌特征(需先求出T和m)常用表征系统是否具有混沌特征一般有两类方法:定性方法(功率谱)和定量方法(最大Lyapimov指数)。利

5、用Fourier分析法求出时间序列的功率谱,从而可以识别该时间序列表征的动力系统的规则性态与不规则性态。若时间序列具有混沌特征,则其功率谱具有连续性、噪声背景和宽峰特征等图形特征;若时间序列是确定性的周期系统,则其功率谱是仅包含有基频和其谐波或分频的离散波形;若时间序列是确定性的准周期系统,则其功率谱是包括不同层次频率的离散波形,但谱线并不像周期运动那样以某间隔的频率分离。最大Lyapunov指数是评判和表征非线性时间序列混沌特性的重要参数,是一个非常关键的混沌不变量。Lyapunov指数是用來描述混沌系统内部相邻相点间辐散的平均速率(其中正Lyapunov指数值(Lyapimov指数0)评判

6、两个相邻轨道的平均指数分离程度,负Lyapunov指数值(Lyapimov指数V0)评判两个相邻轨道的平均指数靠拢程度)。如果一个非线性系统是离散的,那么正Lyapimov扌旨数则是衡量系统是否混沌的一个重要指标。(PS:Lyapunov指数的倒数就是有效预测步数!)从时间序列的角度來研究混沌,我们知道对于决定系统长期演化的任一变量均包含了系统所有变量长期演化的所有信息。因此,我们可以通过决定系统的长期演化的任一单变量时间序列來研究系统的混沌行为,于是帕卡徳(Packard)等人提出的重构相空间理论。相空间的重构:混沌动力学系统分析的第一步是相空间重构。由Takens定理可知系统中任意一个分量

7、的演化均是由与它相互作用的其它分量所决定的。所以,这些相关分量的数据信息隐含在任意一分量的变化过程中,系统相空间的重构只需要考察其中的一个分量,再通过某些延时点上的观测数值找到如Y=x(i),x(iHr),.,x(i+(ml)i:,x(i+(m-l)i)所示的m维向量,就能重构出一个等价的相空间用于恢复原有的动力学系统。从Y中可以看出,未知参数只有m和兀,所以,如何选择适当的嵌入维数m和延迟时间t是相空间重构的主要研究内容。举例:设时间序列为X=1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,1920,假设算出此时间序列的兀=3,m=5o则相空间重构有

8、:M=N(m-1)T.1591317261014183711151948121620嵌入空间维度m的选择为了保证准确计算各种混沌不变量,乂尽量降低运算量和噪声的影响,需要选择合适的嵌入维数m。如果m取值过小,吸引子会发生折叠共至在某些地方会出现自相交,重构吸引子的几何形状和原始状态吸引子可能完全不同;如果m取值过大,吸引子的儿何结构会完全打开,混沌不变量计算精度提高,这在理论上是可行的,但会增加计算量(关联维数和Lyapunov指数),也放大了噪声的影响。目前,嵌入维数的确定方法很多,如饱和关联维数法、邻近点维数法、虚假邻近点法、奇异值分解法、累积局部变形法、轨道扩张法等。这些方法确定最佳嵌入

9、维的原理大致相同,都是依据随嵌入维数m升高而逐步收敛的情况,或利用预测效果來确定最优嵌入维数。本文选用上述方法中最为常用的饱和关联维数法和虚假邻近点法对供水系统的最优嵌入维数m进行研究。时间延迟兀的选择如果延迟时间T取值太小,那么在相空间中各向量的分量间儿乎不包括新信息,从而会低估关联维数(会增加运算量,延迟运算时间)。相反,如果延迟时间T取值太大,相空间重构的相关有效信息会遗漏,这样会高估关联维数。只有在最小的嵌入相空间内能够对相邻轨道实现最优分离的延迟时间T才是重构相空间的最佳延迟时间。目前,确定延迟时间的方法很多,其中由于计算简单,自相关函数法使用最为广泛,而FraserWSwiimey

10、认为自相关函数法对非线性系统的分析可能并不适合,因为自相关法测量的是连续的时间序列之间的线性依托关系,而互信息法则能提取时间序列的非线性特征,通过重构相空间來定量地分析吸引子动力学特征,所以建议选用互信息法來确定延迟时间。当然有些学者认为,m和i的取值是相互关联的,CC算法就可同时求出m和口预测模型:1、最大Lyapunov指数的预测步骤为:对时间序列x(tj.i二1,2,N进行FFT变换,计算平均周期P。(PS:这个建议使用平均法)用CC方法同时计算出嵌入维数m和时间延迟T。根据时间延迟T和嵌入维数m重构相空间Y“i二1,2,M(M=N(m+l)z)找相空间中每个点Y,的最近邻点Yi,并限制

11、短暂分离。(保证了邻近点的变化趋势与中心点、的变化趋势有较高的一致性)|YfYj|,jjPI对相空间中每个点Y”计算出该邻点对的i个离散时间的步后的距离dj(i)o込(=|+门+訂z=1,2rmin(Mj,Mj对每个i,求出所有j的lndj(i)平均y(i)。曲)=丘伞引)At为样本周期,q是非零d,i)的数目,并用最小二乘法作出回归直线,该直线的斜率就是最大Lyapunov指数入1。2、BP网络的混沌时间序列预测的具体步骤为:(1)建立网络,根据混沌时间序列计算出嵌入维数in,用m作为网络的输入个数学习阶段,将N-1个数值序列顺序输人,通过前向过程得到一个输出结果,将该结果与目标模型比较,如

12、果存在误差,立即进行反向传播过程,并修正网络权值,以减小误差,正向输出计算和反向权值修改交替进行,直到误差控制在允许范围内。(这N-1个数值不是单单的一个数字,也可能是一个单元)预测阶段,将第N个数值输人到神经网络中,此时的输出就是该序列的实际预测值。3、加权一阶局域法重构相空间;其基本思路是将一维混沌时间序列映射到高维的空间,目的是恢复有规律的吸引子,从而使蕴藏在时间序列中的信息显露。设混沌序列X】、X2、则重构相空间为:Y=x(i),x(i+r),.x(i+(ml)i:x(i+(m-l)i)(1)(i=l,2,-N-(m-l)r)式中t为时间延迟,m为嵌入维数。选取邻近点。设中心点XM的K

13、个邻近相点为(XMi二1,2,K),到中心点XM的欧式距离为4,设是中最小值,定义乂血的权值为P,则4=时沁(2)e(F(y)E=式中,a为权重调节系数,一般可取(3)进行计算预测。加权一阶局域法线性拟合为:XM1+1eXM1XM2+1=eXM2H.XMk+l.eXMk.其中e=l11:Xmi+】是Xmi的一步演化后的相点,a和b为待求的拟合参数。鉴于大多数文献中a、b计算式冗长、不易理解,本文给出a、b的矩阵表达式,使其便于理解且易于XM1+1eXM1其中Y=XM2+1,A=eXM2:XMk+l-编程。设线性方程组(4)的矩阵表达式为:Y=Ax(5)若定义目标函数为:f(x)=(Ax-Y)tP(Ax-Y)(6)其中P为权重系数矩阵;依据多元函数的极值理论,目标函数取得最小解的充要条件是:=2ArP(Ax-Y)=0比nArPAx=Xpy则可得参数a、b的矩阵表达式::(8)(4)预测计算;将求解的参数a、b带入式(4),则能计算下一步预测点

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论