词条 | 松弛法 |
释义 | songchifa 松弛法(卷名:数学) relaxation method 数值计算中解线性代数方程组的一类迭代法。当方程组的未知量个数甚多而又有大量的零系数时,常用这类方法求解。 基本迭代格式 设线性代数方程组为 ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() 若在迭代格式(2)中,将第i个方程迭代出的 ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() 收敛性 迭代格式(7)和(8)可统一地写为 ![]() ρ(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ω)之间有关系 ![]() ![]() ![]() ![]() ![]() 当取ω=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.) |
随便看 |
百科全书收录78206条中英文百科知识,基本涵盖了大多数领域的百科知识,是一部内容开放、自由的电子版百科全书。