数值计算方法（南京大学）第6章解线性方程组的迭代法.ppt

21 第六章 解线性方程组的迭代法 一、基本迭代法的格式及收敛性 二、几种实用的基本迭代法 1、Jacobi迭代法 2、Gauss-Seidel迭代法 3、超松弛迭代法（SOR) 4、对称超松弛迭代法（SSOR) 5、块超松弛迭代法（BSOR法） 例:讨论用Gauss-Seidel迭代法求解方程组Ax=b时的收敛性，已知 解:(1)对A：不是严格对角占优的矩阵，无法用充分准则I， (2)考虑充分准则II，计算Jacobi迭代矩阵BJ=D-1（L+U）=I-D-1A 不满足充分准则II，故无法判断。 先求出Gauss-Seidel迭代矩阵 BG=(D-L)-1U (3)考虑用定理2的充分条件 不满足定理2的充分条件，故无法判断。 (4)再用定理1的充要条件 例:讨论用Jacobi迭代法和Gauss-Seidel迭代法求解方程组Ax=b时的收敛性，如果收敛，并比较哪种方法收敛较快，其中 1、Jacobi 迭代 Jacobi迭代矩阵 推导其分量形式 第i个方程除以aii(i =1,2,…,n),得 Jacobi迭代的分量形式 则 x(k+1)=BJx(k)+g ， 这里 BJ=D-1(L+U) , g=D-1b Jacobi迭代公式（分量形式） 给出初始向量 x(0), 即可得到向量序列： x(1),x(2),…,x(k),… 若 x(k) → x*, 则x*是解。 例1：设方程组为 解：Jacobi迭代格式为 试写出其Jacobi分量迭代格式以及相应的迭代矩阵，并求解。 故Jacobi迭代矩阵为 取 x(0)=(0,0,0)t, e=10-3,终止准则：‖x(k)-x(k-1)‖<e x(14)= -3.9997 2.9998 1.9998 例2：设方程组为 解： Gauss-Seidel迭代格式为 试写出Gauss-Seidel迭代格式. 2、Gauss-Seidel迭代法 Gauss-Seidel迭代的分量形式 推导Gauss-Seidel迭代法的矩阵形式 Gauss-Seidel迭代矩阵 Gauss-Seidel迭代公式 给出初始向量 x(0), 即可得到向量序列： x(1),x(2),…,x(k),… 若 x(k) → x*, 则x*是解。 Ab.m a(1,1)=1/2+1/4+1/3;a(1,2)=-1/4;a(1,3)=-1/3; a(2,1)=a(1,2);a(2,2)=1/4+1/3+1/5;a(2,3)=-1/5; a(3,1)=a(1,3);a(3,2)=a(2,3);a(3,3)=1/3+1/5+1/3; b(1)=20/2;b(2)=0;b(3)=5/3; function [x,k]=gs(A,b) [n n]=size(A); x=zeros(1,n); for k=1:1000 error=0; for i=1:n s=0;xb=x(i); for j=1:n if i~=j,s=s+A(i,j)*x(j);end end x(i)=(b(i)-s)/A(i,i); error=error+abs(x(i)-xb); end if error/n<0.0001,break;end end fprintf('k.no.=%3.0f,error=%7.2e\n',k,error) |A|=12+4-15=1, |2D-A|=12-4-15=-7 * 数值计算方法【第二版】电子教案 南京大学林成森 数值计算方法 【第二版】电子教案 科学出版社21世纪高等院校教材电子教案系列 南京大学数学系 林成森教授 授课教材 迭代法适用于求解大型稀疏的线性方程组，其 基本思想是通过构造迭代格式产生迭代序列，由迭代 序列来逼近原方程组的解，因此，要解决的基本问题 是：1. 如何构造迭代格式 2.迭代序列是否收敛 一 . 基本迭代法的格式及收敛性 二 . 几种实用的基本迭代法 设有线性代数方程组 A=M+N M的逆好求。Ax =b (M+N)x=b Mx=-Nx+b x=-M-1Nx+M-1b 用矩阵表示： Ax =b A 为系数矩阵，非奇异且设aii≠0；b为右端，x为解向量 注：分解A是一个重要问题 在Rn中，点列的收敛等价于每个分量的收敛。即