tony 发表于 2024-7-13 11:38:29

Matlab混沌四阶电力系统

非线性微分方程如下，我觉得方程没有问题，但是在Q=10以及以上情况下就会出现这个问题：警告: 在 t=3.262927e-02 处失败。在时间 t 处，步长必须降至所允许的最小值(1.110223e-16)以下，才能达到积分容差要求。 > 位置：ode45 (第 352 行)位置: solve_and_plot_chaotic_system (第 10 行) ，有解决方法吗

花月mmc 发表于 2024-7-13 14:15:58

ode45换个求解的方法，鲁棒性强的？

tony 发表于 2024-7-13 19:21:30

花月mmc 发表于 2024-7-13 14:15
ode45换个求解的方法，鲁棒性强的？

这个我试了

花月mmc 发表于 2024-7-13 19:28:27

tony 发表于 2024-7-13 19:21
这个我试了

那看起来就是数值计算的精度要求了，要不就改小步长

tony 发表于 2024-7-13 20:05:31

花月mmc 发表于 2024-7-13 19:28
那看起来就是数值计算的精度要求了，要不就改小步长

这个减小步长怎么进行调整啊大佬，我方程应该是没有问题的，我和开源的类似的论文也对照了但就是跑不出来Q1=10以上的图像，应该是超过4就不行了

花月mmc 发表于 2024-7-13 20:28:21

tony 发表于 2024-7-13 20:05
这个减小步长怎么进行调整啊大佬，我方程应该是没有问题的，我和开源的类似的论文也对照了但就是跑不出来 ...

答案不是就在问题里？
在时间 t 处，步长必须降至所允许的最小值(1.110223e-16)以下，才能达到积分容差要求。
odeset这个选项里面不能设置嘛？

tony 发表于 2024-7-14 11:13:35

花月mmc 发表于 2024-7-13 20:28
答案不是就在问题里？
在时间 t 处，步长必须降至所允许的最小值(1.110223e-16)以下，才能达到积分容差要 ...

问题解决了 我把初值顺序搞错了
页: [1]
查看完整版本: Matlab混沌四阶电力系统