我正在尝试使用有限差分法求解微分方程组。
的项很少,和 r 在边界处都归零。从分析上讲,该术语定义明确,但从数值上讲,这些术语会导致错误,从而很快破坏模拟。
如果这很重要,我正在使用跨越式方法。我在某处读到过蛙式方法具有内置的耗散,因此它应该是稳定的,但显然情况并非如此。
关于如何进行的任何建议?我应该改变我的算法吗?还是使用高阶方案?
提前致谢
我正在尝试使用有限差分法求解微分方程组。
的项很少,和 r 在边界处都归零。从分析上讲,该术语定义明确,但从数值上讲,这些术语会导致错误,从而很快破坏模拟。
如果这很重要,我正在使用跨越式方法。我在某处读到过蛙式方法具有内置的耗散,因此它应该是稳定的,但显然情况并非如此。
关于如何进行的任何建议?我应该改变我的算法吗?还是使用高阶方案?
提前致谢