关于最大Lyapunov指数随参数变化的曲线绘制问题！

octopussheng

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

octopussheng

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

咕噜噜

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

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

flybaly

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


%    The alogrithm employed in this toolbox for determining Lyapunov
%    exponents is according to the algorithms proposed in
%
%    [1] 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.
%
%    [2] 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 [1] is used for first-order systems while
%    the QR-based algorithm proposed in [2] is applied for higher order
%    systems.

octopussheng

:@L :@L :@L

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

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

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

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

octopussheng

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

                               
登录/注册后可看大图

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

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

flybaly

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

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

咕噜噜

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

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

octopussheng

原帖由 <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 编辑 ]

咕噜噜

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

octopussheng

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

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

咕噜噜

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

octopussheng

有这个可能性啊！

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

咕噜噜

好的，记得一起探讨分享

hohoo

在let程序中的Q的初值是一个单位矩阵 根据你方程的自由度数而定 你的方程是一个两自由度的微分方程 因此Q的初值是一个4X4的单位矩阵 在let程序的界面上有一个选项No of linearlized ODEs  应该输入16即可