基于MATLAB科学计算—解线性方程组迭代法.docVIP

科学计算—理论、方法 及其基于MATLAB的程序实现与分析 三、 解线性方程组的迭代法（Iteration） 线性方程组的理论求解公式 　　　　　　　　　（１） 在应用于实际问题的计算时，通常面临两方面的问题 １、计算过程复杂， ２、不能保证算法的稳定性； 此外，当初始数据（可能）存在误差时，按公式（１）即使求出了“精确解”意义也不大，因此，对于存在初始数据误差、特别是大型的线性方程组求解，需要寻求能达到精度要求的、操作和计算过程相对简单的求解方法。下面将要介绍的迭代法就属于这类方法。 迭代法求解线性方程组的基本思想是 不追求“一下子”得到方程组的解，而是在逐步逼近方程组的精确解的迭代过程中获得满足精度要求的近似解，这一点与直接法不同； 通过对问题的转化，避免（困难的）矩阵求逆运算。 用迭代法求解线性方程组，首先要把线性方程组写成等价的形式 　　　　　　　（２） 式（２）的右端称为迭代格式，由迭代格式（２）确定如下的迭代算法： 　　　　　　（３） 对于给定的线性方程组，可以写成不同的（无穷多）迭代格式，有意义的（可用的）迭代格式应具有收敛性―生成的解向量序列收敛于方程组的解；而好的迭代法应具有较高的收敛速度。 关于迭代法收敛性的两个判别条件： a、充分必要条件是：矩阵的谱半径 b、充分条件是：矩阵的某个算子范数。 　　设是方程组（２）的解，是迭代法（３）生成的任一序列，因为 ，