求解非线性方程组的迭代方法.pptVIP

前面介绍的解线性方程组的直接法是解低阶稠 密方程组的有效方法。但是，在工程技术中常产生 大型稀疏矩阵方程组，例如由某些偏微分方程数值 解所产生的线性方程组Ax=b，A的阶数很大 ， 但零元素较多，迭代法是能够充分利用系数矩阵稀 疏性特点的有效算法。;迭代法的构造;§4.1 解线性方程组的三种迭代法 4.1.1．雅克比（Jacobi）迭代法（以三阶方程组为例） 设有方程组：;假设;把X(1) 再代入到(4.1)中得方程组的二次近似.;对于n阶方程组;若用矩阵来记录雅可比矩阵,可作如下的推导:;则有AX=DX-LX-UX=b.即DX=b+(L+U)X;在某种条件下，按雅可比迭代所产生的向量序列的极限会存在，且等于原方程组的解。这种求解方法被称为雅可比迭代法，或简单迭代法。;例 用简单迭代法解下列方程组 ;取初始值x(0)= 0,按迭代公式 ;4.1.2 高斯——赛德尔迭代法;这种迭代的方法用公式写出来就是:;对给定的初值,用此迭代公式求线性方程组的方法被称为高斯—塞德尔迭代法。（G—S）;用矩阵表示此方法为：;例 用赛德尔迭代法解方程组 ;幅耙职郸楚曙队港殷仍疡耿逸污凶氰霉火洞泅膝息迢腕酵闽谰嚼听示酪挂求解非线性方程组的迭代方法求解非线性方程组的迭代方法;例1：分别用两种迭代法求下列线性方程组。初值均取(0，0，0)T;%用雅可比法解P91例1 a=[9,-1,-1;-1,8,0;-1,0,9]; D=-(a-triu(a)-tril(a)); L=-(tril(a)-b); U=-(triu(a)-b); xo=[0;0;0];bo=[7;7;8]; ep=0.0001;dx=1;k=0; while dxep k=k+1; x=D\(L+U)*xo+D\bo; dx=abs(norm(x)-norm(xo)); xo=x; end k,x; 在多数情况下用高斯—赛德尔迭代法比雅克比迭代法收敛快。但也有相反的情况，即高斯—赛德尔迭代法比雅克比迭代法收敛慢，甚至还有雅克比迭代法收敛，高斯—赛德尔迭代法发散的情形。;4.1.3 超松弛迭代法;首先，由高斯—赛德尔迭代法求出一个值，记;用此公式求解线性方程组的方法称为带有松弛因子ω的松弛迭代法. 当ω1时称为超松弛迭代法;(SOR法) 当ω1时称为低松弛迭代法; 当ω=1时就是G—S迭代法. 当某些方程组用高斯—赛德尔迭代法不收敛时，可以用低松弛方法获得收敛，;将上式写成矩阵的形式，得： 于是得SOR迭代的矩阵表示 ;例 用SOR法求解方程;%用SOR法解P96例2 a=[4,-2,-4;-2,17,10;-4,10,9]; D=-(a-triu(a)-tril(a)); L=-(tril(a)-D); U=-(triu(a)-D); xo=[0;0;0];bo=[10;3;-7];omiga=1.46; ep=0.000001;dx=1;k=0;;Matlab的关于三种迭代法的通用程序;2.高斯_塞德尔迭代法的通用程序;3.SOR法解线性方程组的通用程序;4.2 迭代法的收敛条件;定义 设给定Rn中的向量序列{ }，即;证:;定义4.2 设x*是方程组Ax=b的解,对于给定的初始向量x(0) ,若由某种迭代法产生的向量序列{x(m)}有;4.2.2 迭代法收敛的判定定理;证:先证 若;由于E-B是非奇异的，所以方程组(E-B)x=f 的解存在且唯一.设为x*,即x*=Bx*+f,进而有;所以迭代过程收敛.又;再由于;例 写出用雅可比迭代法和G-S迭代法解线性方程组收敛的迭代格式。;由此得;定义4.3 如果方阵A满足;定理4.2 若方程组Ax=b的系数矩阵按行(列)严格对角占优,则雅可比迭代法收敛,G-S迭代法也收敛.;以上两个定理都是收敛的充分条件.下面给出一个充分必要条件：;例：写出用雅可比迭代法求解方程组一定收敛的迭代格式。;由BJ的特征多项式;其系数矩阵是严格对角占优的所以用雅可比迭代法求解收敛.其迭代格式为:;定理4.4 设方程??Ax=b的系数矩阵A为实对称正定矩阵,且0?2, 则松驰迭代法 收敛.;例5:讨论例2中的方程组用SOR方法求解的收敛性. 解:例2中方程组的系数矩阵为:方程组;目前，只有少数特殊类型的矩阵，才有确定的最佳松弛因子的理论公式。例如，当A为对称正定的三对角矩