松弛因子对SOR法收敛速度的影响之C++实现

一、问题描述

用SOR法求解方程组Ax=b, 其中

要求程序中不存系数矩阵A,分别对不同的阶数取w=1.1, 1.2, ...,1.9进行迭代，记录近似解x^(k)达到||x^(k)-x^(k-1)||<10^-6时所用的迭代次数k，观察松弛因子对收敛速度的影响，并观察当w <= 0或 w>=2会有什么影响？

二、计算结果与分析

     (1) 在收敛标准||x^(k)-x^(k-1)||<e=10^-6且最大迭代次数200前提下，分别取阶数n=5、20、35、50、65、80，取w=1.1, 1.2, ...,1.9进行SOR法迭代求解，计算结果见图1：

图1不同阶数下松弛因子w对收敛速度的影响(w=1.1, 1.2, ...,1.9)

 从图1的计算结果可以看出，在w=1.1, 1.2, ...,1.9范围内，不同阶数时SOR法的收敛速度均随收敛松弛因子w的增加而减慢&#x