数值分析大作业 超松弛迭代法如何选取最佳松弛因子.docVIP

超松弛迭代法如何选取最佳松弛因子 船建学院 B1301095 wj 一、课题背景 逐次超松弛迭代法是Gauss-Seidel方法的一种加速方法，是解大型稀疏矩阵方程组的有效方法之一，它具有计算公式简单，程序设计容易，占用计算机内存较少等优点，但需要选择好的加速因子（即最佳松弛因子）。 最佳松弛因子ω的确定是数值代数中的一个理论难题，对于不同的矩阵，其最佳松弛因子往往相差很大，没有统一的计算公式来确定ω。由于对称正定矩阵sor方法收敛的充分必要条件为w在0到2之间，故利用对称正定矩阵一定收敛的性质，本文提供一种针对于系数矩阵为对称正定矩阵时，如何选取合适的最佳松弛因子的方法。 二、课题研究流程图 SOR迭代公式　　逐次超松弛(Successive Over Relaxation)迭代法，简称SOR迭代法，它是在GS法基础上为提高收敛速度，采用加权平均而得到的新算法，设解方程的GS法记为　　　　　　(1)再由与加权平均得　　　 这里ω＞0称为松弛参数，将（1）式代入则得　　　　　　　(2)称为SOR迭代法，[WTBX]ω＞0称为松弛因子，当ω=1时(2)式即为GS法，将(2)式写成矩阵形式，则得　　　 即　　于是得SOR迭代的矩阵表示　　　　　　　(3) Matlab程序%sor法确定对称正定矩阵的最佳松弛因子w% clc;clear; n=100;