octopussheng 发表于 2007-7-9 13:23

啊?怎么这样呢?奇怪中!我再试试!咕噜噜是用let计算的吧?

octopussheng 发表于 2007-7-9 13:30

回复咕噜噜

咕噜噜你好,下面的图是我截的算lorenz系统的LE图,lorenz方程的建立函数就和上面一样的,不用定义Q矩阵!

咕噜噜 发表于 2007-7-9 14:52

^_^,不好意思,是我这里出了点问题
又试过了,那个lorenz方程例子没问题,可我也找不到你自己的问题在那里:@L再找找看啊我

[ 本帖最后由 咕噜噜 于 2007-7-9 15:31 编辑 ]

flybaly 发表于 2007-7-9 14:58

强烈建议调程序之前请看下面参考文献把理论弄清楚!!!
%    The algorithm given in is used for first-order systems while
%    the QR-based algorithm proposed in is applied for higher order
%    systems.


%    The alogrithm employed in this toolbox for determining Lyapunov
%    exponents is according to the algorithms proposed in
%
%    A. Wolf, J. B. Swift, H. L. Swinney, and J. A. Vastano,
%      "Determining Lyapunov Exponents from a Time Series," Physica D,
%      Vol. 16, pp. 285-317, 1985.
%
%    J. P. Eckmann and D. Ruelle, "Ergodic Theory of Chaos and Strange
%      Attractors," Rev. Mod. Phys., Vol. 57, pp. 617-656, 1985.
%
%    The algorithm given in is used for first-order systems while
%    the QR-based algorithm proposed in is applied for higher order
%    systems.

octopussheng 发表于 2007-7-9 15:41

回复flybaly

:@L :@L :@L

求连续系统LE的基本理论看了n遍了,就是一个公式啊!拜托不要老是拿这些什么“请看下面参考文献把理论弄清楚”这些话来压人!

如果你把我的程序存在的问题一一指出来,告诉我为什么错了,然后再跟我说“请看下面参考文献把理论弄清楚”,我肯定很虚心接受的,你就这么来一句,究竟是炫你的理论基础雄厚,还是看不起我们这些浅薄之士呢!

每个人在做课题、项目的时候都会遇到问题,论坛的作用就是通过会员间的交流、互相帮助才能有更大的发展的,如果新会员问问问题,就被你这句话压一下,会认为论坛有作用嘛!

建议各位兄弟姐妹们多向咕噜噜、无水、万老师学习!

octopussheng 发表于 2007-7-9 15:42

原帖由 咕噜噜 于 2007-7-9 14:52 发表 http://www.chinavib.com/forum/images/common/back.gif
^_^,不好意思,是我这里出了点问题
又试过了,那个lorenz方程例子没问题,可我也找不到你自己的问题在那里:@L再找找看啊我


感谢咕噜噜的帮忙啊!
实在不行我就用小数据量法挨个算了,呵呵,我自己也再检查检查!

flybaly 发表于 2007-7-9 16:43

呵呵,既然你理论都知道了,而且把他归结为一个公式,我也没有什么好说的。从语气上看,我是非常的不对了。
讨论问题从来不应该心存芥蒂的。

[ 本帖最后由 flybaly 于 2007-7-9 16:53 编辑 ]

咕噜噜 发表于 2007-7-9 16:47

回复 #37 flybaly 的帖子

fly,嘿嘿,你小子啊,他这个东西我找不出问题来,你是高手帮个忙,我对这方面东西还真是欠缺一些的,学习中

[ 本帖最后由 咕噜噜 于 2007-7-9 20:17 编辑 ]

octopussheng 发表于 2007-7-9 19:34

原帖由 <i>flybaly</i> 于 2007-7-9 16:43 发表 <a href="http://www.chinavib.com/forum/redirect.php?goto=findpost&pid=254396&ptid=47371" target="_blank"><img src="http://www.chinavib.com/forum/images/common/back.gif" border="0"   alt="" /></a><br />
呵呵,既然你理论都知道了,而且把他归结为一个公式,我也没有什么好说的。从语气上看,我是非常的不对了。<br />
讨论问题从来不应该心存芥蒂的。


首先向flybaly致歉,今天下午确实因为问题没有解决而很郁闷,说的话太过激了,请你不要放在心上!

算这个东西也有段时间了,一直出不来好的结果,心里也很急啊!如果这条路走不通,只有换别的方法试试了!

如果flybaly不介意的话,也请你一起帮忙探讨探讨!谢谢

[ 本帖最后由 无水1324 于 2007-7-11 12:13 编辑 ]

咕噜噜 发表于 2007-7-9 20:17

回复 #39 octopussheng 的帖子

flybaly是不介意的,他大度的很,是个专家

octopussheng 发表于 2007-7-9 21:23

回复flybaly、咕噜噜

今天确实比较郁闷,有点口不择言,希望flybaly大人大量啦!再次道歉

不过我这个问题确实很奇怪的,还没有遇到分数维等于0的情况,不知道咕噜噜以前有没有碰到!

咕噜噜 发表于 2007-7-9 21:44

回复 #41 octopussheng 的帖子

没有遇到过,会不会是你的参数选取的问题

octopussheng 发表于 2007-7-10 08:00

有这个可能性啊!

我再试试别的方法吧,谢谢咕噜噜啦!有结果了一定马上告知!

咕噜噜 发表于 2007-7-10 08:10

回复 #43 octopussheng 的帖子

好的,记得一起探讨分享

hohoo 发表于 2007-7-10 19:03

在let程序中的Q的初值是一个单位矩阵 根据你方程的自由度数而定 你的方程是一个两自由度的微分方程 因此Q的初值是一个4X4的单位矩阵 在let程序的界面上有一个选项No of linearlized ODEs应该输入16即可
页: 1 2 [3] 4 5
查看完整版本: 关于最大Lyapunov指数随参数变化的曲线绘制问题!