下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
周期脉冲作用下logisic映射的复杂动力学行为
1周期脉冲作用下离散映射的复杂动力学行为及其分岔机理从数学的形式来看,logistic映射是一个非常简单的二次映射。然而,该系统的动态行为极其复杂。20世纪50年代,生态科学家开始使用这一简单的差分方程来描述群体的变化。目前,对logistic映射的研究取得了许多成果,并在人口预测、生态环境、保密通信等领域得到了应用。脉冲控制具有控制量小、实现方便等优点因此被广泛应用于大型航天器的减震装置、卫星轨道的转换、混沌系统及复杂网络的控制与同步等本文以典型的Logistic映射为例,研究了周期脉冲作用下离散映射的复杂动力学行为及其分岔机理.首先分别研究了两种周期脉冲作用下离散映射的复杂动力学行为,即常量与比例两种形式;然后通过构造Poincar´e映射,对周期脉冲作用下Logistic映射进行了分岔分析;最后基于Floquet理论揭示了该系统周期解的分岔机理.Logistic映射的方程为其中x(k)∈N,k∈N为系统状态,N为自然数集,r为系统参数.设Logistic映射受到周期τ的脉冲作用,即其中m,τ∈N,τ2,g(x(k))为函数.未加脉冲作用时Logistic映射的分岔图与Lyapunov指数如图1所示.从图1可知,当0r<1时,系统(1)有稳定的平衡点,即x=0;当r=1时,系统经鞍结分岔变为稳定的周期一解;当r=3时,系统经倍周期分岔,由周期一解变为稳定的周期二解;当r=3.449时,系统经倍周期分岔,由周期二解变为稳定的周期四解.随着r的逐渐增加,系统经倍周期分岔到达混沌.在混沌区域中存在周期窗口.2系统对的影响本文考虑周期脉冲作用的两种典型情形:(I)g(x)=α;(II)g(x)=αx,即常量和比例两种脉冲作用.本文首先固定周期τ,α,取r为分岔参数,利用频闪映射研究参数r对系统的影响.针对r的某个变化范围中的每一个值,求出系统的数值解,然后取出系统状态的频闪映射的点(即周期脉冲作用后的点)来得到r对系统的影响.对每一个r的值,求出长度为250个脉冲周期的数值解;然后对每一个r的值,给出该系统的最后50个频闪点,从而得到该系统的分岔图.根据频闪映射与原系统的关系得知:频闪映射的不动点对应于原系统存在一个以脉冲周期为周期的周期解,即1T周期解;频闪映射的k周期点对应于原系统存在一个以k倍脉冲周期为周期的周期解,即kT周期解;频闪映射的不变环对应于原系统的拟周期解或圆环面2.1空心点与空心点的分岔机理固定周期τ=3,α=0.1,系统(3)随r变化的分岔图与Lyapunov指数图如图2所示.分岔图上的点为周期脉冲作用后Logistic映射的状态,而文献[20]中的分岔图上的点为Logistic映射的所有状态,即本文采用频闪映射进行分析,这样有利于揭示分岔机理.从图2知,当0r<3.018时系统为稳定的1T周期解,如图3所示当r=3.018时系统经倍周期分岔由1T周期解变为2T周期解,如图4所示.随着r的逐渐增加,系统经倍周期分岔由4T周期解,8T周期解,一直到达混沌状态,如图5和图6所示.当r=3.405时,系统产生混沌危机,直接由混沌状态到达1T周期解.最后经过倍周期分岔再次到达混沌.图3—6中空心圆圈代表周期脉冲作用点,实心点代表Logistic映射作用点.2.2周期脉冲作用的影响固定周期τ=4,α=0.125,系统(3)随r变化的分岔图与Lyapunov指数图如图7所示.从图7知,当0r<0.961时,系统有稳定的平衡点,即x=0.当r=0.961时,系统经鞍结分岔,由平衡解变为稳定的1T周期解,如图8所示.当r=2.927时,系统经倍周期分岔,由1T周期解变为稳定的2T周期解,如图9所示.当r=3.094时,系统经倍周期分岔,由2T周期解变为稳定的4T周期解,如图10所示.随着r的逐渐增加,系统经倍周期分岔到达混沌,如图11所示.下面固定r=2及周期τ,取g(x)=α,并以α为分岔参数,研究脉冲作用对系统的影响.系统(3)随α变化的分岔图如图12所示.从图12可知,在脉冲作用周期固定的情况下,当α=0时,未受脉冲作用的Logistic映射为周期一解;随着α的增加,Logistic映射受到常量脉冲作用,且经历周期解、混沌等复杂动力学现象.从图12(a)—(d)可知,随着脉冲作用周期τ的增加,系统(3)的倍周期分岔过程逐渐被压缩.所以脉冲作用强度及脉冲作用周期对系统动力学具有显著影响.总之周期脉冲作用下的Logistic映射走向混沌的演化过程为经倍周期分岔到达混沌,中间伴随阵发性混沌出现;在混沌区域中存在周期窗口.3非润滑连续系统分岔分析由于周期脉冲作用下的Logistic映射为非光滑离散系统,故不能用一般离散系统的分岔分析方法进行分析.下面将文献[11,12,22]所提出的非光滑连续系统的分岔分析方法推广到非光滑离散系统并对周期脉冲作用下的Logistic映射的动力学行为进行分岔分析.令f(x)=rx(1-x).定义映射P下面基于Floquet理论揭示周期脉冲作用下的Logistic映射的分岔机理.3.1系统特征乘子的floquet特征乘子分析考虑情形I,当3.4047<r<3.49时,系统的Floquet特征乘子在单位圆内;而在r=3.4047处系统的Floquet特征乘子通过1穿出单位圆,如图13所示.根据Floquet理论可判断当3.4047<r<3.49时,系统为稳定的1T周期解,而在r=3.4047处系统发生鞍结分岔,系统由1T周期解变为混沌,此结果与图2(a)的分岔图是一致的.考虑情形II,当0<r<0.9620时,系统的Floquet特征乘子在单位圆内;而在r=0.9620处系统的Floquet特征乘子通过1穿出单位圆,如图14所示.根据Floquet理论可判断当0<r<0.9620时系统为稳定的平衡解,而在r=0.9620处系统发生鞍结分岔,系统由平衡解变为稳定的1T周期解,此结果与图7(a)的分岔图是一致的.总之,脉冲作用下的Logistic映射经鞍结分岔可能由平衡解变为周期解,也可能由周期解直接到达混沌.这两种情形均是发生了鞍结分岔,即系统的Floquet特征乘子通过1穿出单位圆,然而这两种情形的动力学机理并不完全相同.从Lyapunov指数来看,情形II中由平衡解变为周期解时系统的Lyapunov指数随着参数r的增加达到零后就逐渐减少,而情形I中由周期解变为混沌时系统的Lyapunov指数随着参数r的减少达到零后仍逐渐增加.特别地,情形I中随着参数r的增加,系统本质上是发生了混沌危机,即混沌吸引子突变为周期解.所谓混沌危机是指系统混沌吸引子不连续地改变其大小、突然出现或突然消失.3.2t周期解的floquet特征乘子分析考虑情形I,当0<r<3.0180时,系统的Floquet特征乘子在单位圆内;而在r=3.0180处系统的Floquet特征乘子通过-1穿出单位圆,如图15所示.根据Floquet理论可判断当0<r<3.0180,系统做稳定的1T周期解运动,而在r=3.0180处系统发生倍周期分岔,系统由1T周期解变为稳定的2T周期解,此结果与图2(c)的分岔图是一致的.考虑情形II,当2.927<r<3.094时,系统的Floquet特征乘子在单位圆内;而在r=3.094处系统的Floquet特征乘子通过-1穿出单位圆,如图16所示.根据Floquet理论可判断当2.927<r<3.094时,系统做稳定的2T周期解运动,而在r=3.094处系统发生倍周期分岔,系统由2T周期解变为稳定的4T周期解,此结果与图7(c)的分岔图是一致的.总之,周期脉冲作用下Logistic映射可能会发生倍周期分岔,即系统的某一周期解变为另一两倍周期的周期解,并会发生级联倍周期分岔,即系统的周期解经倍周期分岔到达混沌.4floquet理论分岔行为的理论分析本文研究了周期脉冲作用下Logistic映射
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 医院护理后备人才竞聘
- 2024自卸车买卖合同模板
- 2024至2030年中国花岗石Ⅴ形块行业投资前景及策略咨询研究报告
- 2024至2030年中国电感识别系统数据监测研究报告
- 2024至2030年中国混合环己酮数据监测研究报告
- 2024至2030年中国文字校对系统数据监测研究报告
- 2024至2030年中国小花盆行业投资前景及策略咨询研究报告
- 2024至2030年中国喷油器清洗流量测试台行业投资前景及策略咨询研究报告
- 天津市人教版初中物理八年级上册第一章综合经典复习题(含解析)
- 医疗行业年会
- 建筑设计防火规范
- 4D厨房设备设施管理责任卡
- GB/T 3655-2022用爱泼斯坦方圈测量电工钢带(片)磁性能的方法
- GB/T 25420-2021驱动耙
- 特应性皮炎积分指数AD 病情严重程度积分法(SCORAD)
- GB/T 19520.1-2007电子设备机械结构482.6mm(19in)系列机械结构尺寸第1部分:面板和机架
- GB/T 16762-2020一般用途钢丝绳吊索特性和技术条件
- 2023年北京市昌平区广播电视台(融媒体中心)招聘笔试题库及答案解析
- 主要耗能设备管理台账
- 2018年木地板公司组织架构及部门职能
- 露天矿山开采课件
评论
0/150
提交评论