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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 统计、非线性&软物质 » 【求助】有限差分法误差都是多少啊
101/1返回列表
查看: 977  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

大葱1890

铜虫 (小有名气)

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

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

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

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的回帖

mozhui

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
把x(i+1)和x(i-1)都在i处泰勒展开到四次项带到你的计算公式里可以得到截断误差是h^2，有什么问题吗？
赞一下(1人) 回复此楼
5楼2009-06-07 16:02:19
已阅   回复此楼   关注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的回帖

yzcluster

金虫 (著名写手)

小木虫扫盲人
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
"误差就是h^{3}的高阶小"，自然是指包括4阶在内以上的误差，不包括3阶项的误差。
赞一下(1人) 回复此楼
7楼2009-06-07 16:49:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mozhui

还要除以h^2…
回复此楼
8楼2009-06-07 16:54:16
已阅   回复此楼   关注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的回帖

mozhui

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
另外一般情况下约定小写o{}表示比括号里高阶的无穷小，而大写O{}表示与括号里同阶的无穷小~

[ Last edited by mozhui on 2009-6-7 at 17:08 ]
赞一下(1人) 回复此楼
10楼2009-06-07 16:58:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 大葱1890 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 293求调剂 +6 世界首富 2026-03-11 6/300 2026-03-17 17:04 by ruiyingmiao
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 10/500 2026-03-17 15:07 by 一鸭鸭哟
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[基金申请] 今年的国基金是打分制吗？ 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 求调剂（材料与化工327） +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 一志愿西南交大，材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 283求调剂，材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
信息提示
关闭
请填处理意见
关闭