计算实习任务1汇总.docxVIP

中国石油大学(北京) 计算实习任务 中国石油大学(北京) 计算实习任务1 PAGE PAGE # 中国石油大学(北京) 计算实习任务 中国石油大学(北京) 计算实习任务1 PAGE PAGE # 用有限差分法求解正方形域上的 Poisson方程边值问题 +—宴 =f (x, y) =1,0 ex, y 』件 cy丿 u(0, y)二 u(1,y) = y(1 - y), 0 曲叮 u(x,0) u(x,1) x(1—x), 0exe1 解：用五点差分格式对椭圆型第一边值问题得到线性方程组为 4ui, j - ui J, j - ui 1,j - ui ,j」 ui, j 1 = h fj 1 i, j N 在边界上 U0,j =UN，1,^ j*h*(1- j*h) Ui,0 二Ui,N 1 =i*h*(1—i*h),0S,j 乞 N 1 写成矩阵形式Au=f，其中 —1u =V2■A^NN / —1 u = V2 ■ A ^NN / VNJ 『4 _1 b2 -1 4匕 Aii = + + a 11 4 - 4 - -1 bN丿 —1 4 J ,VN = (U1,N , U2, N , ,UN ,N V1 二(％1,U2,1, ,UN,1), V2 二(u1,2 , u2,2 , UN,2), b1 ~ h ( f1,1, f2,1, f3,1 , f