| 查看: 1084 | 回复: 2 | ||
[求助]
关于拉普拉斯方程的求解问题！求助 已有1人参与
|
回复此楼» 本主题相关价值贴推荐，对您同样有帮助:
Partial Differential Equations（三卷全）【Michael E. Taylor】【已搜索，无置重】
已经有291人回复
- 关于一个偏微分方程的转化问题
已经有5人回复
- 各位大神，这个方程怎么解啊？？
已经有14人回复
- 100金币求助拉普拉斯方程球坐标下的形式推导问题
已经有16人回复
- 关于菲克第二定律在边界条件下的求解问题
已经有12人回复
- 哪位高手帮我看看这个定理？关于反应扩散方程的
已经有4人回复
- 用matlab自带PDETOOL求解拉普拉斯方程，后处理中关于边界上二次偏导的问题，急！
已经有4人回复
- 傅立叶变换、拉普拉斯变换、Z变换的联系，求解惑
已经有5人回复
- 点电荷的周围的真空到底满足泊松方程求还是拉普拉斯方程
已经有5人回复
- 请问这个偏微分方程组有一般的解法么？
已经有12人回复
- 关于laplace变换和fourier变换的问题
已经有14人回复
- matlab可以求解差分方程系统的解吗？
已经有6人回复
- 微分方程求解（通解中含hypergeom函数）
已经有17人回复
- Mathematica问题求助
已经有7人回复
- 一个椭圆型偏微分方程的数值解
已经有19人回复
- 反拉普拉斯变换积分先后问题
已经有7人回复
- 求教这个是什么类型的偏微分方程，应该用什么方法来求解！
已经有9人回复
- 偏微分中 有点像反写的6 的符号 是在什么情况下使用，谢谢。
已经有5人回复
- 【求助】拉普拉斯方程求解
已经有6人回复
- 【求助】急请微分方程高手看下这个方程组！
已经有14人回复
- 【交流】关于贝塞尔函数
已经有7人回复
- 【求助】关于共形映射的问题，不甚感激！
已经有12人回复
- 【求助】离散拉普拉斯变换
已经有6人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
feixiaolin: 金币+10, 2015-06应助之星 2015-07-02 12:01:52
感谢参与，应助指数 +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的方法，分析上有解的。 |
2楼2015-06-11 01:43:13
【答案】应助回帖
★ ★ ★ ★ ★
erwty: 金币+5, ★有帮助 2015-06-17 18:09:08
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。 |
3楼2015-06-11 01:48:12
5