求助大家一个龙格库塔的简单问题（在线等）

radicle

我在用三阶龙格库塔解迁移方程时出现了耗散现象，时间推进越长，幅值越小，这是为什么啊？

zinsser

这是正常的，离散迁移方程的格式有耗散，幅值当然会降低。

lixl-imech

应当是数值耗散引起的。 三阶RK 方法的数值耗散相对较小，估计可能是空间离散引起的耗散。你可以采用低耗散或无耗散的差分格式试一下。 例如把空间导数的差分由原先的一阶迎风换成2阶中心或高阶迎风格式。  
    如果不改变数值方法，单纯减小时间步长及空间步长也会改善结果。