第三章编程作业.docVIP

数值实验题三 3.用有限差分方法（五点差分格式）求解正方形域上的Poisson方程边值问题 用MATLAB语言编写求解线性方程组的算法程序，采用下列方法计算，比较计算结果和算法性能，对计算结果给出结论。 （1）用Jacobi迭代法求解线性方程组。 （2）用块Jacobi迭代法求解线性方程组。 （3）用（预条件）共轭斜向量法求解线性方程组。 解：由差分格式可得: 写成矩阵形式：Au=f 其中： 其中： （1）Jacobi迭代法： function [U,k,er,t]=Jcb(n) % Jacobi迭代法 % U 表示方程组的解；h 表示步长； A 表示迭代矩阵；k 表示迭代次数；n 表示非边界点数 % f 表示线性方程组A*U=f的右端矩阵f ；e 表示允许误差界；er 表示迭代误差 % t 表示计算时间 h=1/(n+1); f(2:n+1,2:n+1)= 2*h^2; %初始化f A=zeros(n+2); %初始化A为n+2阶零矩阵 U=zeros(n+2); %初始化U为n+2阶零矩阵 e=0.000000001; %设置误差界 tic; for k=1:10000 %迭代求解 er=0; for i=2:n+1 for j=2:n+1