资源描述：
《科学与计算工程办法第3章.pdf》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、第三章线性方程组的解法§3.0引言§3.1雅可比(Jacobi)迭代法§3.2高斯-塞德尔(Gauss-Seidel)迭代法§3.3超松驰迭代法§3.4迭代法的收敛性§3.5高斯消去法§3.6高斯主元素消去法§3作业讲评3§3.7三角分解法§3.8追赶法§3.9其它应用§3.10误差分析§3.11总结§3.0引言重要性：解线性代数方程组的有效方法在计算数学和科学计算中具有特殊的地位和作用.如弹性力学、电路分析、热传导和振动、以及社会科学及定量分析商业经济中的各种问题.分类：线性方程组的解法可分为直接法和迭代法两种方法.(a)直接法：对于给定的方程组
2、，在没有舍入误差的假设下，能在预定的运算次数内求得精确解.最基本的直接法是Gauss消去法，重要的直接法全都受到Gauss消去法的启发.计算代价高.(b)迭代法：基于一定的递推格式，产生逼近方程组精确解的近似序列.收敛性是其为迭代法的前提，此外，存在收敛速度与误差估计问题.简单实用,诱人.§3.1雅可比Jacobi迭代法(AX=b)1基本思想:与解f(x)=0的不动点迭代相类似,将AX=b改写为X=BX+f的形式,建立雅k+1(k)可比方法的迭代格式：X=BX+f,其中，B称为迭代矩阵.其计算精度可控,特别适用于求解系数为大型稀疏矩阵(sparse
3、matrices)的方程组.2问题:(a)如何建立迭代格式？k(b)向量序列{X}是否收敛以及收敛条件?3例题分析:⎧10x1−x2−2x3=2.7⎪考虑解方程组⎨−x1+10x2−2x3=3.8(1)⎪⎩−x1−x2+5x3=2.4*其准确解为X={1,1.2,1.3}.建立与式(1)相等价的形式：⎧x1=1.0x2+2.0x3+.072⎪⎨x2=1.0x1+2.0x3+.083(2)⎪⎩x3=1.0x1+2.0x2+.084据此建立迭代公式：(k+)1(k)k⎧x=1.0x+2.0x+.072123⎪(k+)1(k)(k)⎨x2=1.0x1+2
4、.0x3+.083(3)⎪(k+)1(k)(k)x=1.0x+2.0x+.084⎩312)0()0()0(取迭代初值x=x=x=0,迭代结果如下表.123JocabiMethodP31.cpp迭代次数x1x2x3000010.720.830.8420.9711.071.1531.0571.15711.248241.085351.185341.2828251.0950981.1950991.29413861.0983381.1983371.29803971.0994421.1994421.29933581.0998111.1998111.299777
5、91.0999361.1999361.299924101.0999791.1999791.299975111.0999931.1999931.299991121.0999981.1999981.299997131.0999991.1999991.299999141.11.21.3151.11.21.34Jocobi迭代公式:设方程组AX=b,通过分离变量的过程建立Jocobi迭代公式，即n∑ax=b,a≠0(i=,2,1?,n)ijjiiii=11nx=(b−∑ax()i=,2,1?,n)iiijjaj=1iij≠i由此我们可以得到Jacobi迭代
6、公式：1n(k+)1kx=(b−∑ax()i=,2,1?,n)iiijiaj=1iij≠i[Jacobi迭代公式的算法]1:初始化.n,(aij),(bj),(x1),M.2:执行k=1直到M为止.①执行i=1直到n为止.nui←(bi−∑aijxj/)aii;j=1j≠i②执行i=1直到n为止.x←u;ii③输出k,(xi).另外，我们也可以建立Jacobi迭代公式的矩阵形式.设方程组AX=b,其中，A=(aij)n为非奇异阵，TX=(x1,x2,…,xn),Tb=(b1,b2,…,bn)将系数阵A分解为:A=U+D+L,U为上三角矩阵，D为对角
7、矩阵，L为下三角矩阵.于是AX=b可改写为-1-1(U+D+L)X=b⇔X=Db-D(U+L)X由此可得矩阵形式的Jocobi迭代公式：k+1(k)X=BX+f□§2高斯-塞德尔Gauss-Seidel迭代法(k+)1注意到利用Jocobi迭代公式计算x时，已经计算好i(k)(k)(k)x,x,?,x的值，而Jocobi迭代公式并不利用这些最新的近似值计12i−1(k)(k)(k)算，仍用x,x,?,x.这启发我们可以对其加以改进,即在每个分量的12i−1计算中尽量利用最新的迭代值，得到1i−1n(k+)1(k+)1kx=(b−∑ax−∑ax()i
8、=,2,1?,n)iiijjijjaj=1j=i+1ii上式称为Gauss-Seidel迭代法.其矩阵形式是-1-1X=-
