微分方程数值解法实验.docVIP

PAGE PAGE 1 一、实验概述： 【实验目的】 熟练掌握求解方程第一边值问题的五点差分格式。 【实验原理】 1.取定沿X轴和Y轴的步长h1和h2，作两族与坐标轴平行的直线：，两族直线的交点称为网点或节点，记为，假定为正则内点。沿x,y方向分别用二阶中心差商代替,,则得，由于差分方程中只出现u在（i,j）及其四个邻点上的值，故称为五点差分格式。 2.迭代法 设线性方程组 （1） 的系数矩阵可逆且主对角元素均不为零，令，并将分解成 （2） 从而（1）可以写成，令，其中 （3） 以为迭代矩阵的迭代法 （4） 称为迭代法，用向量的分量来表示，（4）为 （5） 其中为初始向量。 3.迭代法 把矩阵分解成 （6） 其中，分别为的主对角元除外的下三角和上三角部分，于是方程组（1）可以写成，即，其中 （7） 以为迭代矩阵构成的迭代法 （8） 称为迭代法，用向量表示的形式为 4.迭代 【实验环境】 MATLAB R2014a 二、实验内容： 【实验方案】 1.求解边值问题 取步长h=k=1/64,1/128,做五点差分格式。用jacobi迭代，Gauss-Seidel迭代和SOR迭代（取）求解差分方程，以前后两次重合到小数点后四位的迭代值作为解的近似。比较三种解法的迭代次数以及差分解与精确解的精度。 2.在