应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 统计、非线性&软物质 » 【求助】有限差分法误差都是多少啊
51/1返回列表
查看: 1034  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

大葱1890

铜虫 (小有名气)

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

用有限差分法计算二阶导数时的误差是多少啊?比如(d^2 x)/(dt^2 ）假设t的步长取为h 那么误差是h的几次方项？
回复此楼

» 猜你喜欢
1楼 2009-06-06 12:53:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大葱1890

铜虫 (小有名气)
引用回帖:
Originally posted by yzcluster at 2009-6-7 16:02:
如果我印象没错误差应该是O(h^{3})。其实很简单，你把f(x+h)和f(x-h)分别泰勒展开，两者做和（把步长h当作小量），你会发现所有h的单数项都被消去了，只剩下2，4，6，。。。等项。自然，误差就是h^{3}的高阶小。

我算了一下是这么做，不过好像最后是得出h的四次方的误差，O(h^{3})代表的是四次方的还是三次的我也不怎么清楚。不过还是谢谢你！
赞一下 回复此楼
6楼2009-06-07 16:46:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

mozhui

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
这个和具体近似方法有关吧，比如，对于一阶导数的近似，采用向前和向后差商带来的截断误差都是h的一次方，而采用中心差商带来的截断误差则为h的平方。楼主是如何计算二阶导数的？
赞一下(1人) 回复此楼
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}的高阶小。
赞一下(1人) 回复此楼
4楼2009-06-07 16:02:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭