24小时热门版块排行榜    

查看: 981  |  回复: 3

Love_Family

新虫 (初入文坛)

[求助] 数值解的误差曲面问题 已有2人参与

我用向前差分方法计算了一个抛物方程的初边值问题,结果要求画出误差曲面,但是我总是画不对;
我的计算步骤如下,请高人指点:
1、方程的经典解是u(x,t)=exp(x+t);  0=<x<=1, 0=<t<=1;
2、计算的数值解为 U(xi,tk) (总共计算出11*11个节点的数值)
3、matlab编程如下:
x=0:0.01:1;
y=x;
[X,Y]=meshgrid(x,y);
for i=1:11
   for j=1:11
   E(i,j)=abs(U(i,j)-u(i,j));   %数值解与经典解得误差
   end
end
surf(X,Y,E)
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

No Complain,Only work hard!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongcailiu

金虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
Love_Family: 金币+1, 有帮助 2014-09-29 21:43:53
没有计算数值解的过程啊,另外,变量x的步长和后面的也不符啊
2楼2014-09-29 08:28:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

至尊木虫 (职业作家)

科苑小木虫

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
Love_Family: 金币+1, 有帮助 2014-09-29 21:44:03
feixiaolin: 金币+20, 2014-4 优秀板油 2014-12-07 01:00:25
x=0:0.01:1和你的11个点不符。改了看看
板凳要做十年冷文章不发一个字
3楼2014-09-29 10:10:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Love_Family

新虫 (初入文坛)

步长是0.1,我写帖子里写错了;数值解的计算过程我没有写出来,数值解的计算没有问题,现在的问题是我画的曲面和课本不一样。
No Complain,Only work hard!
4楼2014-09-29 21:46:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Love_Family 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[教师之家] 跳槽后在研项目怎么办? +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见