计算方法 线性方程组的迭代解法优秀PPT.ppt

上传人:石*** 文档编号:65760495 上传时间:2022-12-08 格式:PPT 页数:50 大小:2.75MB
返回 下载 相关 举报
计算方法 线性方程组的迭代解法优秀PPT.ppt_第1页
第1页 / 共50页
计算方法 线性方程组的迭代解法优秀PPT.ppt_第2页
第2页 / 共50页
点击查看更多>>
资源描述

《计算方法 线性方程组的迭代解法优秀PPT.ppt》由会员分享,可在线阅读,更多相关《计算方法 线性方程组的迭代解法优秀PPT.ppt(50页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、计算方法 线性方程组的迭代解法现在学习的是第1页,共50页三、小结三、小结 二、线性方程组的解法二、线性方程组的解法一、线性方程组有解的判定条件一、线性方程组有解的判定条件回顾回顾线性代数线性代数中线性方程组的解法中线性方程组的解法现在学习的是第2页,共50页一、线性方程组有解的判定条件一、线性方程组有解的判定条件一、线性方程组有解的判定条件定理定理4 n 元线性方程组元线性方程组Ax=b(i)无解的充分必要条件是无解的充分必要条件是R(A)R(A,b);(ii)有唯一解的充分必要条件是有唯一解的充分必要条件是R(A)=R(A,b)=n;(iii)有无限多解的充分必要条件是有无限多解的充分必要

2、条件是R(A)=R(A,b)n.当R(A)=R(B)=rn时,n元线性方程组可由含有n-r个参数的解来表示,这是线性方程组的通解。现在学习的是第3页,共50页定理定理定理5 线性方程组线性方程组Ax=b有解的充分必要条件是有解的充分必要条件是R(A)=R(A,b).定理定理6 n元齐次线性方程组元齐次线性方程组Ax=0有非零解的充分必要条件是有非零解的充分必要条件是R(A)n.定理定理7 矩阵方程矩阵方程AX=B有解的充分必要条件是有解的充分必要条件是R(A)=R(A,B).定理定理8 设设AB=C,则,则R(C)minR(A),R(B).定理定理9 矩阵方程矩阵方程AX=0只有零解的充分必要

3、条件是只有零解的充分必要条件是R(A)=n现在学习的是第4页,共50页小结小结有唯一解有唯一解bAx=()()nBRAR=()()nBRAR=有无穷多解有无穷多解.bAx=齐次线性方程组:系数矩阵化成行最简形矩阵,便可写出齐次线性方程组:系数矩阵化成行最简形矩阵,便可写出其通解;其通解;非齐次线性方程组:增广矩阵化成行阶梯形矩阵,便可非齐次线性方程组:增广矩阵化成行阶梯形矩阵,便可判断其是否有解若有解,化成行最简形矩阵,便可写判断其是否有解若有解,化成行最简形矩阵,便可写出其通解;出其通解;当当R(A)=R(B)=rn时,由于含有时,由于含有n-r个参数的解可表示线性方程个参数的解可表示线性方

4、程组的任一解,因此称为线性方程组的组的任一解,因此称为线性方程组的通解通解。现在学习的是第5页,共50页二、线性方程组的解法例例1 1 求解齐次线性方程组求解齐次线性方程组解解二、线性方程组的解法二、线性方程组的解法现在学习的是第6页,共50页即得与原方程组同解的方程组即得与原方程组同解的方程组现在学习的是第7页,共50页由此即得由此即得方程组的通解是:现在学习的是第8页,共50页例例 求解非齐次线性方程组求解非齐次线性方程组解解对增广矩阵对增广矩阵B进行初等变换,进行初等变换,,3)(,2)(=BRAR由于,由于,故方程组无解故方程组无解现在学习的是第9页,共50页例例 求解非齐次方程组的通

5、解求解非齐次方程组的通解解解 对增广矩阵对增广矩阵B进行初等变换进行初等变换现在学习的是第10页,共50页故方程组有解,且有故方程组有解,且有现在学习的是第11页,共50页所以方程组的通解为所以方程组的通解为k1k2现在学习的是第12页,共50页例例 解证解证对增广矩阵对增广矩阵B进行初等变换,进行初等变换,方程组的增广矩阵为方程组的增广矩阵为现在学习的是第13页,共50页现在学习的是第14页,共50页由于原方程组等价于方程组由于原方程组等价于方程组由此得通解:由此得通解:现在学习的是第15页,共50页例例 设有线性方程组设有线性方程组解解现在学习的是第16页,共50页现在学习的是第17页,共

6、50页其通解为其通解为现在学习的是第18页,共50页这时又分两种情形:这时又分两种情形:现在学习的是第19页,共50页现在学习的是第20页,共50页此题也可以先求系数行列式。此题也可以先求系数行列式。现在学习的是第21页,共50页三、小结()()nBRAR=()()nBRAR|aij|i=1,2,n,j=1,ji 则称方阵则称方阵A是严格是严格(行行)对角占优的对角占优的.a11 a12 a13 a1n a21 a22 a23 a2n A=L+D+U an1 an3 an4 ann -4 2 1例 矩阵 A=1 -9 7 2 -6 10ULD现在学习的是第31页,共50页Jacobi 迭代一:

7、设有方程组 a11x1+a12x2+a1nxn=b1 a21x1+a22x2+a2nxn=b2 .an1x1+an2x2+annxn=bn用矩阵表示:Ax=b (A 为系数矩阵,非奇异;b为右端,x为解向量)现在学习的是第32页,共50页假设 aii0 令 cij=-aij/aii (ij)gi=bi/aij,i=1,2,3,n 则 x1(k+1)=c12x2(k)+c13x3(k)+c1nxn(k)+g1 x2(k+1)=c21x1(k)+c23x3(k)+c2nxn(k)+g2 。xn(k+1)=cn1x1(k)+cn2x2(k)+cn(n-1)xn-1(k)+gn Jacobi迭代格式迭

8、代格式 若令若令 0 c12 c13 c1n x1 c21 0 c23 c2n x2BJ=x=.cn1 cn3 cn4 0 xn a11 g1 a22 g=g2 易看出:BJ=D-1(L+U)P148,D=.ann gn 把方程组写成容易迭代的形式:现在学习的是第33页,共50页Jacobi迭代公式迭代公式现在学习的是第34页,共50页Gauss-Seidel迭代法迭代法 为了加快收敛速度,同时为了节省计算机的内存,我们作如下的改进:每算出一个分量的近似值,立即用到下一个分量的计算中去,即用迭代格式:这样所得的迭代法就称为Gauss-Seidel迭代法,也称为“异步迭代法”,简称为GS迭代法利

9、用Ax=b 及A=L+D+U,其中D为对角矩阵,L,U分别为严格下,上三角矩阵则有,GS迭代法的矩阵形式为:现在学习的是第35页,共50页Seidel迭代法的具体形式迭代法的具体形式Seidel迭代格式迭代格式 x1(k+1)=c12x2(k)+c13x3(k)+c1nxn(k)+g1 x2(k+1)=c21x1(k+1)+c23x3(k)+c2nxn(k)+g2 。xn(k+1)=cn1x1(k+1)+cn2x2(k+1)+cn(n-1)xn-1(k+1)+gn 假设 aii0 令 cij=-aij/aii (ij)gi=bi/aij,i=1,2,3,n 现在学习的是第36页,共50页收敛性

10、及误差估计Jacobi迭代和GS迭代格式可表述为统一形式:对于其收敛性,我们有如下定理:定理:对任意初始向量x(0)及任意右段向量 g,由此产生的迭代向量序列x(k)收敛的充要条件是证明:必要性:设x(k)收敛,其极限为 x*,则现在学习的是第37页,共50页因上式对任意均成立,故Bk 0 (k ),(B)1 充分性:设(B)1,则IB为非奇异阵,且 =0,所以 有唯一解,记为 则 定理:若|B|1,则迭代法收敛.推论若 满足下列条件之一:()现在学习的是第38页,共50页推论2:如果A为严格对角占优阵,则其 Jacobi迭代和Seidel迭代对任何初始向量x(0)都收敛。推论3:如果A为对称

11、正定阵,则其 Seidel迭代对任何初始向量x(0)都收敛。现在学习的是第39页,共50页三例题及求解三例题及求解例:用迭代法解方程组AX=b,其中已知该方程的解为 解:本题分别用简单迭代法(Jacobi迭代法)和GS迭代法来解(1)简单迭代法现在学习的是第40页,共50页现在学习的是第41页,共50页表现在学习的是第42页,共50页现在学习的是第43页,共50页表2现在学习的是第44页,共50页四.相关程序设计原始数据(A,b)可用一个二维数组存储,也可将A用一个二维数组,b用一个一维数组分别存储,存储需要一个一维数组程序中应方便地对迭代方法和终止条件的选择以及对初始向量和值的设置在迭代过程

12、中,为反映迭代情况,可设置一些中间数据的输出,如迭带次数,迭代向量,迭代残向量等当然不需要每迭代一次都作输出这可作为收敛情况或不收敛情况的分析作为不收敛的判定,可设置一个大的整数,当迭带次数超过该数时作为不收敛处理GS 迭代法的计算公式为:现在学习的是第45页,共50页开始TFTFT现在学习的是第46页,共50页SOR法介绍法介绍利用利用Seidel Seidel 迭代法迭代法 ,考虑如下,考虑如下SOR SOR 迭代格式迭代格式一步一步SeidelSeidel迭代迭代松弛因子松弛因子SORSOR迭代格式也是线性迭代形式迭代格式也是线性迭代形式现在学习的是第47页,共50页现在学习的是第48页,共50页Seidel迭代格式为 从式中解出 得 故可得Seidel迭代矩阵为现在学习的是第49页,共50页迭代法的特点优点:方法简单,每次迭代都是简单的重复运算,易于编制程序;与求解线性方程的精确法相比,简单迭代法对于字长位数较少的计算机更为适用,它可以用增加迭代次数来弥补字长位数少的不足.初值可以任取,因而中间结果偶然错误不影响最后结果的获得。缺点:迭代速度较慢。用计算机计算时,。就其收敛性而言,某些用Seidel迭代法不能收敛而无法得出结果的线性代数方程组,用Jacoai迭代法却能进行收敛计算,反之亦然.现在学习的是第50页,共50页

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 生活休闲 > 资格考试

本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

工信部备案号:黑ICP备15003705号© 2020-2023 www.taowenge.com 淘文阁