24小时热门版块排行榜    

CyRhmU.jpeg
查看: 906  |  回复: 9
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

大葱1890

铜虫 (小有名气)

[交流] 【求助】有限差分法误差都是多少啊

用有限差分法计算二阶导数时的误差是多少啊?比如(d^2 x)/(dt^2 )假设t的步长取为h 那么误差是h的几次方项?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yzcluster

金虫 (著名写手)

小木虫扫盲人


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by mozhui at 2009-6-7 16:54:
还要除以h^2…

是的
9楼2009-06-07 16:56:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

mozhui


小木虫(金币+0.5):给个红包,谢谢回帖交流
这个和具体近似方法有关吧,比如,对于一阶导数的近似,采用向前和向后差商带来的截断误差都是h的一次方,而采用中心差商带来的截断误差则为h的平方。楼主是如何计算二阶导数的?
2楼2009-06-06 13:17:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大葱1890

铜虫 (小有名气)

(d^2 x)/(dt^2 )=(x_(i+1)-2x_i+x_(i-1))/h^2  其中h是时间步长,对于一阶的问题我算过了,可以用x_n+1在nh上展开级数得出,一阶向前差分的相当于把x_n+1项的h的平方项得省略,可是我用相同方法算不出来二阶的,很郁闷!
3楼2009-06-07 15:52:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yzcluster

金虫 (著名写手)

小木虫扫盲人


小木虫(金币+0.5):给个红包,谢谢回帖交流
如果我印象没错误差应该是O(h^{3})。其实很简单,你把f(x+h)和f(x-h)分别泰勒展开,两者做和(把步长h当作小量),你会发现所有h的单数项都被消去了,只剩下2,4,6,。。。等项。自然,误差就是h^{3}的高阶小。
4楼2009-06-07 16:02:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见