nebula1983 发表于 2007-12-13 17:36

求助用误差控制龙格－库塔－费尔博格变步长法求解方程组的问题

到现在为止我看见的用龙格－库塔－费尔博格变步长法求解都是求解单个常系数线性方程，步长的变化控制使用这个公式:delta=0.84*(tol/E)^(1/4);其中tol是需要的精度，E是局部截断误差，由方程两个为当前步长时得到的近似值相减得到（数值分析书上有相关公式）。
   现在的问题是：当方程变为多个（组）时，每个方程都会有一个E，那么delta该如何确定呢？我觉得不可能分别求，那样的话每个方程的步长变化就会不一致了，该如何解决这个问题呢？请各位解惑，谢谢！

[ 本帖最后由 sigma665 于 2007-12-13 17:43 编辑 ]

sigma665 发表于 2007-12-13 17:59

我想问一下,步长不一致,会有什么影响.
我没用过..

nebula1983 发表于 2007-12-13 18:49

我们把范围（0,L）用步长分段，如果每个方程的步长不一样的话，我们到达L的‘时间’估计也就不一样了，好像不允许这样吧？

nebula1983 发表于 2007-12-14 23:03

没人知道吗？请各位帮帮忙啊

花如月 发表于 2007-12-15 11:31

回复 #4 nebula1983 的帖子

请不要催贴，版主们并非什么问题都懂。这里主要解决一些编程问题而非专业问题， xjzuo对解方程很有研究。你可以看看他的示例帖子，如果还是没办法就给点耐心、、

nebula1983 发表于 2007-12-15 16:34

哦，知道了。我只是怕帖子沉下去。而不是催贴
页: [1]
查看完整版本: 求助用误差控制龙格－库塔－费尔博格变步长法求解方程组的问题