2-第二章 线性方程组的数值解法.ppt

§3 迭代法的收敛性 例：讨论用雅可比迭代法和高斯-赛德尔迭代法解方程组 的收敛性，其中 解： ，所以方程组有唯一解。 雅可比迭代法的迭代矩阵为 * §3 迭代法的收敛性 经过计算得到其三个特征值都为0，因此 ，所以 雅可比迭代法收敛。 高斯-赛德尔迭代法的迭代矩阵为 经过计算得到其三个特征值为0，2，2，因此，高斯-赛德尔 迭代法不收敛。 * 例：设有线性方程组 讨论用雅克比迭代法和高斯-赛德尔迭代法求解时的收敛性。 解：雅克比迭代矩阵为 因为 ，有 ，从而 ，所以雅克比迭代法不收敛。 §3 迭代法的收敛性 * 高斯-赛德尔迭代矩阵为 因为 ，有 ， 从而 ，所以高斯-赛德尔迭代 法收敛。 §3 迭代法的收敛性 * 练习题 1、分别用高斯消去法和列主元消去法求解下列方程组。（计算取4位小数） 2、证明下列方程组的雅克比迭代法和相应的高斯-赛德尔迭代法收敛，并写出迭代公式： * 3、设有方程组 （1）证明用雅克比迭代法与高斯-赛德尔迭代法解此方程组均收敛； （2）取初始向量 ，分别用雅克比迭代法与高斯-赛德尔迭代法求解，要求 时迭代终止。 练习题 * 4、设线性方程组 的系数矩阵为 证明用雅克比迭代法收敛，而用高斯-赛德尔迭代法不收敛。 5、设线性方程组 的系数矩阵为 证明用雅克比迭代法不收敛，而用高斯-赛德尔迭代法收敛。 练习题 * * §2 线性方程组的迭代解法 则原方程组用矩阵形式写为 求解后的方程组写为 令 上式可写为 * §2 线性方程组的迭代解法 迭代过程： 首先取 ，将其代入（2.2.2）式右端 得到 即 * §2 线性方程组的迭代解法 再将 代入（2.2.2）式右端，求得 ，如此 反复进行迭代，得到 即 由此得一个向量序列 如果 ，则 就是方程组 的解，这种求解 线性方程组的方法称为雅可比迭代法。 * §2 线性方程组的迭代解法 便于程序编写的形式 * §2 线性方程组的迭代解法 例：用雅可比迭代法求解方程组 当 时，迭代终止（方程组的精 确解为 ） * §2 线性方程组的迭代解法 解：取 ，雅可比公式为 * clc; %清屏 close all; %关闭窗口 clear all; %清除所有变量 n=4; %变量数 A=[5 1 -1 -2;2 8 1 3;1 -2 -4 -1;-1 3 2 7]; %系数矩阵 b=[-2;-6;6;12]; %常数向量 X0=zeros(n,1); %未知数向量 X=zeros(n,1); k=0; %迭代次数 §2 线性方程组的迭代解法 * §2 线性方程组的迭代解法 while 1 k=k+1; %计数 for i=1:n X(i,1)=X0(i,1)+(b(i,1)-A(i,:)*X0)/A(i,i); %雅可比计算公式 end Err(k)=max(abs(X-X0)); %误差 if Err(k)0.00001 break; %退出循环 end X0=X; end * §2 线性方程组的迭代解法 0 0 0 0 0 1 -0.4 -0.75 -1.5 1.7142857 1.7142857 2 0.1357143 -1.1053571 -1.65357