2楼: Originally posted by feixiaolin at 2013-10-23 19:50:13
可以调整步长，
或者引入递归滤波试一试: New_X(k+1)=0.9*X(k+1)+0.1*X(k);

4楼: Originally posted by 绿野渊 at 2013-10-25 14:55:26
楼主做的是瞬态计算，应该确保每一个时间步长内解都不发散。也就是说，每一个时间步长内每个物理量对应的残差曲线都要在（这个时间步长的）迭代结束时变成水平方向，即残差不再随着迭代的进行显著减小。
建议楼主把 ...