24小时热门版块排行榜    

查看: 266  |  回复: 0

liujk189

木虫 (初入文坛)

[交流] 有哪位兄弟做过多次外推的, 请进来看一下.........

最近在写关于外推的程序, 先写了一个一维的Helmholtz方程的, 如下:
u''+ku=f x\in[-1,1]
以及相应的边界条件.......
用有限差分格式, 二阶导数用中心差分格式逼近,
取网格步长分别为h, h/2, h/4 计算出三组解, U_h, U_h/2, U_h/4,
然后用一次外推得到两个四阶的逼近解 U1和U2
U1=(4U_h/2-U_h)/3
U2=(4U_h/4-U_h/2)/3
用这两个四阶的逼近再进行外推得到
U3=(16U2-U1)/15
这个解应该是一个六阶的逼近, 但是我现在只能得到四阶精度, 请大家帮忙看看是哪里错误了??
如果有更好的建议, 请提出来, 本人不胜感激..........
QQ或者邮件都可以联系, 直接回贴也可以.

QQ: 357017677
Email: liujk189@163.com
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 liujk189 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见