材料090班杨东麟0604090113(Page41).docVIP

材料090班杨东麟03(Page41) 材料090班杨东麟03(Page41) PAGE 材料090班杨东麟03(Page41) １、编程实现以下科学计算算法，并举一例应用之。 （参考书籍《精 通ＭＡＬＡＢ科学计算》 ，王正林等著，电子工业出版社，２００９ 年） “Gauss-Seidel 迭代法线性方程组求解” （1）迭代解法的基本思想： 根据给定方程组，设计出一个迭代公式，构造一数组的序列xi0,代入迭代公式，计算出xi1,在代入迭代公式，经过k次迭代运算后得到xik,若xik收敛于某一极限数组xi，则xi就是方程组的近似解。 迭代过程本质上就是计算极限的过程，一般不能得到精确解。但迭代的优点是程序简单，适合大型方程组求解，然而，缺点是要判断迭代是否收敛和收敛速度的问题。 (2)算法说明： Gauss-Seidel迭代法与简单迭代法类似，只是迭代公式有所改进。 简单迭代法：； Gauss-Seidel迭代法：； 设方程组Ax＝b，其中A和b中的元素都为常数，且A为非奇异，则A分可写成：A=D-L-U。 其中D上网意义同Jacobi迭代法，L为下三角矩阵，U为上三角矩阵，他的迭代公式为： 在MATLAB中编程实现的Gauss-Seidel迭代法函数为：gauseidel。 功能：用Gauss-Seidel迭代求线性方程组ax=b的解。 调用格式：[x,n]=gauseidel(A