luohy1725 发表于 2012-3-27 20:54

最大lyapunov指数图与分岔图不符合的问题?

对一个动力系统进行求解,用龙格库塔法求得了分岔图,用wolf法求得了最大lyapunov指数图,但是分岔图和最大lyapunov指数图在一些却不相符(椭圆圈处最不符合),请同学们看下问题大概出在什么地方?

kezairenjian 发表于 2012-4-8 14:49

1、这个系统是不是有吸引子共存的情况,如果有,那么共存区域有两个不同的MLE。比如,周期吸引子和混沌吸引子共存,对应不同的初值会有不同的MLE,一个小于零的MLE对应为周其吸引子,一个大于零对应为混沌吸引子。因此分岔图的初值一定要和MLE的初值相等。
2、算法的精度问题。

luohy1725 发表于 2012-4-16 22:23

回复 2 # kezairenjian 的帖子

谢谢您的回复。计算的时候,MLE和分岔图的初值都是一样的。用的是定步长的四阶龙格库塔法法,后改为变步长的四阶龙格库塔法,结果还是一样。

luohy1725 发表于 2012-4-16 22:44

更正一下,最大lyapunov指数求解用的不是wolf方法,而是波内廷的方法,具体见陈立群,刘延柱的非线性振动一书。

公子汀 发表于 2012-11-30 13:41

学长,能把您的求最大lyapunov指数的程序发给我学习一下吗?
我最近也在学习这一块的东西.1245230415@qq.com

犟牛 发表于 2012-11-30 15:17

公子汀 发表于 2012-11-30 13:41 static/image/common/back.gif
学长,能把您的求最大lyapunov指数的程序发给我学习一下吗?
我最近也在学习这一块的东西

论坛有很多相关的程序,搜索一下

公子汀 发表于 2012-11-30 16:22

犟牛 发表于 2012-11-30 15:17 static/image/common/back.gif
论坛有很多相关的程序,搜索一下

我有搜索,但对于齿轮系统,不知道用哪种方法算比较好,我解动力学方程用的也是龙格库塔法。

gghhjj 发表于 2012-12-14 05:27

本帖最后由 gghhjj 于 2012-12-14 05:30 编辑

公子汀 发表于 2012-11-30 16:22 static/image/common/back.gif
我有搜索,但对于齿轮系统,不知道用哪种方法算比较好,我解动力学方程用的也是龙格库塔法。
李指数的计算方法大体上分为wolf法jocobian法两类
由于你的是仿真结果,不存在噪声,所以用wolf法就行
当然,如果只是判断系统是否会混沌,则计算最大李指数即可,可以用小数据量法

页: [1]
查看完整版本: 最大lyapunov指数图与分岔图不符合的问题?