Poisson方程离散时不同的系数矩阵方式

参考文献:科学计算中的偏微分方程有限差分法

B=[4 -2 0 0 0;-1 4 -1 0 0;0 -1 4 -1 0;0 0 -1 4 -1;0 0 0 -2 4];
I=eye(5);
Z=zeros(5);
A=[B -2*I Z Z Z;-I B -I Z Z;Z -I B -I Z;Z Z -I B -I;Z Z Z -2*I B];

为了保持形式上的对称,试着将A改写为

A=[1.0/2*B -I Z Z Z;-I B -I Z Z;Z -I B -I Z;Z Z -I B -I;Z Z Z -1*I 1.0/2*B];

 

 

 

 

posted @ 2013-01-10 22:24  liang_l  阅读(352)  评论(0编辑  收藏  举报