Matlab第次上机作业.docxVIP

第一次上机作业数值实验题三2.用有限差分法（五点分格式）求解正方形域上的Poission方程边值问题用MATLAB语言编写求解线性方程组Au=f的算法程序，采用下列方法计算结果和算法性能，对计算结果给出结论。（1）用Jacobi迭代法求解线性方程组Au=f；（2）用块Jacobi迭代法求解线性方程组Au=f；（3）用（预条件）共轭斜量法求解线性方程组Au=f。解：1、基本原理说明：首先剖分求解区域。用平行于坐标轴的直线将求解区域分为网格。然后用数值微分公式对微分算子进行离散。即有其中，uij表示u(xi, yi)的近似值。得到在每个点(xi, yi)上的有限差分方程为： （3-74）又称五点差分格式。其中fij=f(xi, yi)。在边界上，有（3-75）对非边界点编号，顺序为对直线y=yi从上往下，对固定的一条直线的网点从左往右依次编号，即有相应的解向量和右端向量分别为由3-74和3-75得差分方程组将差分方程组写成矩阵形式为，其中。差分方程等价于Av=b。这样，原问题就可用Jacobi迭代法或块Jacobi迭代法求解方程组。块Jacobi迭代法公式可化简为接下来就可以编写适当的数值算法求解线性方程组。2、计算方法：（1）Jacobi迭代法：算法程序如下：function [u,k]=xsj(n)% xsj:用Jacobi迭代法求解线性方程组A*u=f% u:方程组的解； k迭代