无水，你给我的分岔图程序怎么只出来一列点啊？

yh0247

final-value 是多少，出来的图就在那一列

[ 本帖最后由 咕噜噜 于 2007-6-14 18:56 编辑 ]

无水1324

你仔细看一下那个程序的每一行。明白它的意思就知道了

yh0247

无水，我总是遇到这个警告：
Warning: Failure at t=3.647831e-001.  Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.295970e-015) at time t.
在一些时候，我怎么改都不成，有没有经验可以教教我阿

shenyongjun

这应该是方程出现奇异或者刚性的情况！

无水1324

对你得仔细检查你的方程的参数正确否，
再则有可能出现申老师说的奇异或者刚性的问题，那你就需要换其他的算法了

octopussheng

请问无水，如果方程出现刚性的话，改用什么算法较好呢？

无水1324

可能使你的方程有问题，适当调低点精度

无水1324

原帖由 octopussheng 于 2007-6-28 12:02 发表

                               
登录/注册后可看大图

请问无水，如果方程出现刚性的话，改用什么算法较好呢？

刚性方程，不适合用ode45求解，你可以采用ode15s

octopussheng

但是我现在用0de15s求解也会出问题啊！
Warning: Failure at t=4.369739e-001.  Unable to meet integration tolerances without reducing the step size below the smallest value allowed (8.881784e-016) at time t.
> In <a href="error:C:\MATLAB7\toolbox\matlab\funfun\ode45.m,355,1">ode45 at 355</a>
  In <a href="error:C:\MATLAB7\work\ode_equation.m,6,1">solve_ode_equation at 6</a>

无水1324

检查方程的正确性，参数的准确性。调低计算精度

octopussheng

看来引起这个问题的原因主要还是参数的选取啊！

谢谢无水

无水1324

原帖由 octopussheng 于 2007-6-28 13:30 发表

                               
登录/注册后可看大图

看来引起这个问题的原因主要还是参数的选取啊！

谢谢无水

  说实话，我原来遇到这个问题也是参数出了问题

octopussheng

不知道无水算Lyapunov指数的时候，有没有遇到过分数维等于0的情况？