24小时热门版块排行榜    

查看: 1080  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

erwty

金虫 (正式写手)

[求助] 关于拉普拉斯方程的求解问题!求助 已有1人参与

 CODE:
\"关于拉普拉斯方程的求解问题!求助\"

w178h2650753_1433914995_653.jpg|bcs#opennewwindow\" class=\"content_pic\" target=\"_blank\">CODE:
\"关于拉普拉斯方程的求解问题!求助\"

w178h2650753_1433914995_653.jpg|bcs\" border=\"0\" style=\"max-width:700px;\">
关于拉普拉斯方程的求解问题!求助-1
IMAG0227_proc.jpg
回复此楼

» 猜你喜欢

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

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

oinkmasta

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
erwty: 金币+5, 有帮助 2015-06-17 18:09:08
真抱歉,上面式子发错了,我读了一阶导数的系数了。正确的应该是:
x方向上的二阶偏导项,三个点用forward是     1*u(i,j) - 2*u(i+1,j) + 1*u(i+2,j)
                                           用backward是  1u(i-2,j) - 2*u(i-1,j) + 1*u(i,j)
y方向上的二阶偏导项,三个点用forward是     1*u(i,j) - 2*u(i,j+1) + 1*u(i,j+2)
                                           用backward是  1*u(i,j-2) - 2*u(i,j-1) + 1*u(i,j)

需要留意的是forward和backward在用点的数量相同的情况下比central低一级精确度;也就是说正比于h而不是h^2。
16届本科在北美
3楼2015-06-11 01:48:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

oinkmasta

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
feixiaolin: 金币+10, 2015-06应助之星 2015-07-02 12:01:52
题里给的是三个点的central difference,
x方向上的二阶偏导项是  1*u(i-1,j) - 2*u(i,j) + 1*u(i+1,j)
y方向上的二阶偏导项是  1*u(i,j-1) - 2*u(i,j) + 1*u(i,j+1)
合起来就是 1*u(i-1,j) + 1*u(i+1,j) + 1*u(i,j-1) + 1*u(i,j+1) - 4*u(i,j)

边值上缺一些点,不能用central;左边和下边用forward difference,右和上用backward。
x方向上的二阶偏导项,三个点用forward是  −3/2*u(i,j) + 2*u(i+1,j) - 1/2*u(i+2,j)
                                           用backward是  −1/2*u(i-2,j) + 2*u(i-1,j) - 3/2*u(i,j)
y方向上的二阶偏导项,三个点用forward是  −3/2*u(i,j) + 2*u(i,j+1) - 1/2*u(i,j+2)
                                           用backward是  −1/2*u(i,j-2) + 2*u(i,j-1) - 3/2*u(i,j)
用这些把原来的项替换掉就可以了。

另:这个边值问题可以直接用PDE的方法,分析上有解的。
16届本科在北美
2楼2015-06-11 01:43:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856材料与化工,270求调剂 +3 YXCT 2026-03-01 3/150 2026-03-01 19:51 by 无懈可击111
[考研] 高分子化学与物理调剂 +6 好好好1233 2026-02-28 12/600 2026-03-01 19:48 by 好好好1233
[考研] 0805总分292,求调剂 +5 幻想之殇 2026-03-01 5/250 2026-03-01 19:48 by 无懈可击111
[考研] 化工299分求调剂 一志愿985落榜 +5 嘻嘻(*^ω^*) 2026-03-01 5/250 2026-03-01 19:47 by 无际的草原
[考研] 291分工科求调剂 +9 science饿饿 2026-03-01 10/500 2026-03-01 18:55 by 18137688336
[考研] 290求调剂 +9 材料专硕调剂; 2026-02-28 11/550 2026-03-01 17:21 by sunny81
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 285求调剂 +8 满头大汗的学生 2026-02-28 8/400 2026-03-01 16:47 by caszguilin
[考研] 311求调剂 +6 亭亭亭01 2026-03-01 6/300 2026-03-01 15:41 by 324616
[考研] 307求调剂 +5 wyyyqx 2026-03-01 5/250 2026-03-01 15:21 by Fff-1
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[考研] 311求调剂 +9 南迦720 2026-02-28 10/500 2026-03-01 10:55 by sunny81
[考博] 博士自荐 +4 kkluvs 2026-02-28 4/200 2026-03-01 10:19 by 馥安馥安
[硕博家园] 2025届双非化工硕士毕业,申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 10/500 2026-03-01 10:02 by 科研狗111
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[考研] 272求调剂 +4 田智友 2026-02-28 4/200 2026-03-01 06:43 by 刘兵
[考研] 材料调剂 +4 爱擦汗的可乐冰 2026-02-28 4/200 2026-03-01 00:38 by 猫猫球alter
信息提示
请填处理意见