| 记得很早之前，在做某单位的活儿的时候，来了一个博士，在闲聊时告诉我，“瞬态动力学我们都用欧拉法就够了，什么龙格库塔什么Newmark方法都是浮云！”，令人印象深刻。时间一长，竟也被潜移默化，没有深究个中深意，总觉得欧拉法精度差，那就把时间步取密一点不就好了？ 最近需要自己处理一些瞬态动力学过程，遇到了这样那样的问题。不得不老老实实翻书，清华的《计算动力学》一开始就从中心差分法的稳定性开始分析，而且有言：“欧拉法一般只能用于起步或与其他方法配合使用”。这些不得不使我好奇被博士尊为"就够了"的欧拉法未什么这么不受待见，自己用书上的方法推导了一下欧拉法的稳定性，与其他方法相对比，才知道欧拉法精度之差，之难以稳定…… 总结了一下，写了个小程序做了三个算例，算是在这个问题上彻底给自己一个交代，下次遇到学术问题，不要全信他人，哪怕他是博士……呵呵 [attach]49510[/attach] |
GMT+8, 2024-5-20 14:29 , Processed in 0.030010 second(s), 15 queries , Gzip On.
Powered by Discuz! X3.4
Copyright © 2001-2021, Tencent Cloud.
