第三章 椭圆型方程的差分格式
第1页/共36页
§ 正方形区域中的Laplace 方程 Dirichlet边值问题的差分模拟
设Ω是 xy 平面中的具有正方形边界 的
一个有界区域,考虑Laplace方程的第一边值
Dirichlet )问题
第2页/共36页
网格节点(l , m)
处的二阶中心差商代替
二阶微商
l , m+1
l–1,m l , m l+1 , m
l , m–1
第3页/共36页
Laplace方程的五点差分格式()为
截断误差为O(h2)。
-Ul, m+1
-Ul–1,m 4U l , m -U l+1 , m
-Ul, m–1
第4页/共36页
令
则Laplac方程的五点差分格式为()
即
第5页/共36页
例1 用五点差分格式求解 Laplace方程
在区域
内的近似解,边界值为:
取 。
第6页/共36页
解 网格点如图所示
U7
U8
U9
U4
U5
U6
U1
U2
U3
u(1,0)=20 u(2,0)=20 u(3,0)=20
u(1,4)=180 u(2,4)=180 u(3,4)=180
u(0,3)=80
u(0,2)=80
u(0,1)=80
u(4,3)=0
u(4,2)=0
u(4,1)=0
第7页/共36页
第8页/共36页
第9页/共36页
矩阵方程AU=K,K由边界条件所确定,解得
U = [U1 U2 U3 U4 U5 U6 U7 U8 U9]’=A-1K
= [
]T
第10页/共36页
五点差分格式 来自淘豆网www.taodocs.com转载请标明出处.