松弛法编辑本段回目录
正文编辑本段回目录
数值计算中解线性代数方程组的一类迭代法。当方程组的未知量个数甚多而又有大量的零系数时，常用这类方法求解。
基本迭代格式 设线性代数方程组为
,x
,…,x
，就可求出x
,x
,…,x
,由后者又可求出x
，x
,…,x
,如此继续直到满足精度要求为止,例如对所有的i,
或
（若
）或剩余 
｜小于某给定的小数ε时就终止迭代而将
取作方程组(1)的解。迭代格式(2)称为同时超松弛法或JOR法，其中的ω为一实参数，称为松弛因子，当ω＝1时就是通常所谓的简单迭代法，或称雅可比法。
若在迭代格式(2)中，将第i个方程迭代出的
的值代替其后各方程中的
，就得到迭代格式 
,
为n维列向量，则(1)可写为
① 若A为严格对角优势矩阵（见对角优势矩阵）或不可约弱对角优势矩阵，|D -1B|为矩阵D -1B的元素取绝对值所成的矩阵，则JOR法和SOR法在
0<ω<2/(1+ρ)|D -1B|范围内收敛，并且在所述条件下,上式右端是一个大于1的数；
② 若A为埃尔米特矩阵，且αii(i＝1,2,…,n)均为正，则JOR法收敛的充分必要条件为A和2ω -1D-A正定；
③ SOR法收敛的必要条件为0<ω<2；
④ 若A为埃尔米特矩阵，αii(i=1,2,…,n)均为正，则SOR法收敛的充分必要条件为A正定和0<ω<2。
松弛因子的选取 ω 的值选取得适当可使松弛法有较好的收敛性，然而如何选取最优的ω,还是一个困难的问题。通常用五点差分格式解二维二阶椭圆型方程得到的线性代数方程组的系数矩阵是块三对角矩阵，主对角块是三对角阵，非主对角块为对角阵。对这种方程组若记简单迭代和超松弛迭代的迭代矩阵为J和Lω，则它们的特征值λ(J)和λ(Lω)之间有关系
)分别为矩阵J和的谱半径。这些结果还可推广到所谓“相容次序”的矩阵。尽管如此，在实际计算中仍难以精确确定ωb的值，而往往是通过一些试算或用其他方法先估计ρ(J)，再来估计ωb。也可用若干ω值试迭代,从中选取使敛速较快者。从(12)式还可看出，ωb在1与2之间,且从大于ωb的一侧选取ω,然后逐步减小，较为有利。
当取ω=1时,由(12)有λ(L1)=(λ(J))2。即对相容次序的矩阵来说，赛德尔迭代比简单迭代敛速快一倍。在很多情况下,赛德尔迭代比简单迭代收敛快。实际上,若A =M-N，且A -1、M -1和N的元素均非负，则迭代过程
此外,若将A分成块来形成分块矩阵,而将(6)(7)(8)中的D、B、L和U分别取为A的主对角块、－A的非主对角块、下三角块和上三角块所形成的矩阵，则迭代(15)的收敛速度可能更快。这样得出的迭代格式称为块松弛迭代，相应地有块简单迭代，块超松弛迭代等。对照之下，前面的迭代就称为点松弛迭代，上面关于最佳松弛因子的讨论对块松弛迭代也是适用的。
上述的一些迭代法有时收敛都很慢，这就需要用一些辅助方法来加速收敛，例如半迭代加速、共轭梯度法加速等。
参考书目
冯康等编：《数值计算方法》，国防工业出版社，北京，1978。
D.M.Young,Iterative Solution of large Linear Systems,Academic Press, New York,1971.
R.S.瓦格著,蒋尔雄、游兆永、张玉德译:《矩阵迭代分析》，上海科学技术出版社,上海,1966。(R.S.Varga,Matrix Iterative Analysis , Prentice - Hall,Englewood Cliffs, New Jersey, 1962.)
基本迭代格式 设线性代数方程组为
(1)
(2)
,x,…,x，就可求出x,x,…,x,由后者又可求出x，x,…,x,如此继续直到满足精度要求为止,例如对所有的i,或（若）或剩余 
｜小于某给定的小数ε时就终止迭代而将取作方程组(1)的解。迭代格式(2)称为同时超松弛法或JOR法，其中的ω为一实参数，称为松弛因子，当ω＝1时就是通常所谓的简单迭代法，或称雅可比法。 若在迭代格式(2)中，将第i个方程迭代出的
的值代替其后各方程中的，就得到迭代格式 
(3)
(4)
,为n维列向量，则(1)可写为 
(5)
(6)
, (7)
。 (8)
(9)
ρ(G)<1。
(10)由此可得收敛的一个充分条件为这里‖G‖是G的任意一‖G‖<1 ， (11)
种范数。对上面松弛法的迭代格式可以得出一些有关收敛性的命题如下：① 若A为严格对角优势矩阵（见对角优势矩阵）或不可约弱对角优势矩阵，|D -1B|为矩阵D -1B的元素取绝对值所成的矩阵，则JOR法和SOR法在
0<ω<2/(1+ρ)|D -1B|范围内收敛，并且在所述条件下,上式右端是一个大于1的数；
② 若A为埃尔米特矩阵，且αii(i＝1,2,…,n)均为正，则JOR法收敛的充分必要条件为A和2ω -1D-A正定；
③ SOR法收敛的必要条件为0<ω<2；
④ 若A为埃尔米特矩阵，αii(i=1,2,…,n)均为正，则SOR法收敛的充分必要条件为A正定和0<ω<2。
松弛因子的选取 ω 的值选取得适当可使松弛法有较好的收敛性，然而如何选取最优的ω,还是一个困难的问题。通常用五点差分格式解二维二阶椭圆型方程得到的线性代数方程组的系数矩阵是块三对角矩阵，主对角块是三对角阵，非主对角块为对角阵。对这种方程组若记简单迭代和超松弛迭代的迭代矩阵为J和Lω，则它们的特征值λ(J)和λ(Lω)之间有关系
。 (12)
。 (13)
， (14)
)分别为矩阵J和的谱半径。这些结果还可推广到所谓“相容次序”的矩阵。尽管如此，在实际计算中仍难以精确确定ωb的值，而往往是通过一些试算或用其他方法先估计ρ(J)，再来估计ωb。也可用若干ω值试迭代,从中选取使敛速较快者。从(12)式还可看出，ωb在1与2之间,且从大于ωb的一侧选取ω,然后逐步减小，较为有利。 当取ω=1时,由(12)有λ(L1)=(λ(J))2。即对相容次序的矩阵来说，赛德尔迭代比简单迭代敛速快一倍。在很多情况下,赛德尔迭代比简单迭代收敛快。实际上,若A =M-N，且A -1、M -1和N的元素均非负，则迭代过程
(15)
此外,若将A分成块来形成分块矩阵,而将(6)(7)(8)中的D、B、L和U分别取为A的主对角块、－A的非主对角块、下三角块和上三角块所形成的矩阵，则迭代(15)的收敛速度可能更快。这样得出的迭代格式称为块松弛迭代，相应地有块简单迭代，块超松弛迭代等。对照之下，前面的迭代就称为点松弛迭代，上面关于最佳松弛因子的讨论对块松弛迭代也是适用的。
上述的一些迭代法有时收敛都很慢，这就需要用一些辅助方法来加速收敛，例如半迭代加速、共轭梯度法加速等。
参考书目
冯康等编：《数值计算方法》，国防工业出版社，北京，1978。
D.M.Young,Iterative Solution of large Linear Systems,Academic Press, New York,1971.
R.S.瓦格著,蒋尔雄、游兆永、张玉德译:《矩阵迭代分析》，上海科学技术出版社,上海,1966。(R.S.Varga,Matrix Iterative Analysis , Prentice - Hall,Englewood Cliffs, New Jersey, 1962.)
配图编辑本段回目录
