线性方程组直接法例题.ppt

* 例 用列主元高斯消去法求解方程组(用三位有效数字计算) 解 选主元 选主元 消元过程完成，得到上三角形方程组 再作回代可求得 例 设 讨论a 取何值时，矩阵A可作LU分解。 解 当A的顺序主子式不为零，矩阵A有LU分解。 所以，有 例 紧凑格式 例题 试判断迭代公式是否收敛。 解 迭代矩阵 , 故迭代发散。 雅可比迭代收敛。 G-S迭代收敛。 例4.7 设 ,证明, 求解方程组 的Jacobi迭代与G-S迭代同时收敛或发散 证:雅可比迭代矩阵 其谱半径 证: G-S迭代矩阵 其谱半径 显然, 和 同时小于、等于或大于1,因而Jacobi迭代法与G-S迭代法具有相同的收敛性 当时?a?1时,Jacobi矩阵??GJ??∞1,对初值x(0)均收敛 例4.11 设 方程组 写出解方程组的Jacobi迭代公式和迭代矩阵 并讨论迭代收敛的条件。 写出解方程组的Gauss-Seidel迭代矩阵,并讨 论迭代收敛的条件。 解 ① Jacobi迭代公式和Jacobi矩阵分别为 例4.11设 方程组 写出解方程组的Gauss-Seidel迭代矩阵，并讨论 迭代收敛的条件。 解 ② Gauss-Seidel矩阵为 当时?a?1时, Gauss-Seidel矩阵 ??Gs??∞1, 所以对任意初值x(0)均收敛。 也可用矩阵的谱半径p(GS)1来讨论 求解AX=b,当?取何值时迭代收敛？ 解:所给迭代公式的迭代矩阵为 例4.13 给定线性方程组 AX= b 用迭代公式 X(K+1)=X(K)+?(b-AX(K)) (k=0,1,…） 即 ?2-(2-5 ?)?+1- 5 ?+4 ?2=0 ?2-(2-5 ?)?+(1- ? )(1-4?)=0 [?-(1-?)][?- (1-4?)]=0 ?1=1-? ?2=1-4? ?(B)=max{|1- ?|, |1-4?|}1 取0 ?1/2迭代收敛