




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、第4章 解线性方程组的迭代法 直接法得到的解是理论上准确的,但是我们可以看得出,它们的计算量都是n3数量级,存储量为n2量级,这在n比较小的时候还比较合适(n400),但是对于现在的很多实际问题,往往要我们求解很大的n的矩阵,而且这些矩阵往往是系数矩阵就是这些矩阵含有大量的0元素。对于这类的矩阵,在用直接法时就会耗费大量的时间和存储单元。因此我们有必要引入一类新的方法:迭代法。 迭代法具有的特点是速度快。与非线性方程的迭代方法一样,需要我们构造一个等价的方程,从而构造一个收敛序列,序列的极限值就是方程组的根对方程组bAx 做等价变换gGxxbMNxMxNxbMxbxNMbAx11)(如:令NM
2、A,则则,我们可以构造序列gxGxkk)()1( 若*)(xxkbAxgxGx* *同时:*)(*)()()1(xxGGxGxxxkkk*)()0(1xxGk0kG所以,序列收敛与初值的选取无关与初值的选取无关定义6.1:(收敛矩阵)0kG定理:矩阵G为收敛矩阵,当且仅当G的谱半径eps) x1=x2; for(i=0;in;i+) x2i=0; for(j=0;ji;j+) x2i += Aij*x1j for(j=i+1;jn;j+) x2i += Aij*x1j x2i=-(x2i-bi)/Aii 4、输出解x2 迭代矩阵迭代矩阵ULDA记nnaaD0011000001121nnnaaa
3、L000001112nnnaaaU易知,Jacobi迭代有bxULD)(bxULDx)(bDxULDx11)(bDgADIULDG111 , )( 收敛条件收敛条件 迭代格式收敛的充要条件是G的谱半径eps) for(i=0;in;i+) for(j=0;ji;j+) x2i += Aij*x2j for(j=i+1;jn;j+) x2i += Aij*x2j x2i=-(x2i-bi)/Aii 4、输出解x2 迭代矩阵迭代矩阵)()()1(1)1(bUxLxDxkkkbDUxDxLDIkk1)(1)1(1)(bDLDIUxDLDIxkk111)(111)1()()(bLDUxLDxkk1)(
4、1)1()()(bLDgULDG11)( , )( 是否是原来的方程的解?A=(D-L)-U 收敛条件收敛条件 迭代格式收敛的充要条件是G的谱半径eps) for(i=0;in;i+) temp-0 for(j=0;ji;j+) temp += Aij*x2j for(j=i+1;jn;j+) temp += Aij*x2j temp = -(x2i-bi)/Aii x2i = (1-omega)*x2i+omega*temp 4、输出解x2)()1 (1)(1)1(1)()1(bDUxDLxDxxkkkk 迭代矩阵迭代矩阵bDxUDIxLDIkk1)(1)1(1)1()(bDLDIxUDIL
5、DIxkk111)(111)1()()1()(bDLDIgUDILDIG111111)( )1()( 定理:松弛迭代收敛20定理:A对称正定,则松弛迭代收敛20是否是原来的方程的解? SORSOR方法收敛的快慢与松弛因子的选择有密切关系.但是如何选取最佳松弛因子,即选取=*,使( )达到最小,是一个尚未很好解决的问题.实际上可采用试算的方法来确定较好的松弛因子.经验上可取1.41.6. 定理定理 若SORSOR方法收敛, 则02. 证证 设SORSOR方法收敛, 则( )1,所以 |det( )| =|12 n|1而 det( ) =det(D-L)-1 (1-)D+U) =det(E-D-1L)-1 det(1-)E+D-1U) =(1-)n于是 |1-|1, 或 02 定理定理 设A是对称正定矩阵, 则解方程组Ax=b的SORSOR方法,当00 (Uy,y)=(y,Ly)=(Ly,y) =-i 0(Ay,y)=(Dy,y)-(Ly,y)-(Uy,y) =-2所以)()()1 (ii2222222)()(|当02时,有 (-+)2-(-)2= (2-)(2-) = (2-)(2-)0所以|21, 因此( )1,即S0R方法收敛.可得 =2/ 设是B的任一特征值, y是对应的特征向量, 则 (L+U)y=Dy于是 (Ly,y)+(Uy,y)=
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 垃圾分类试点协议书
- 企业清洁承包协议书
- 婚俗改革巡演协议书
- 申请暂缓就业协议书
- 运营车位转让协议书
- 民间协议书时效多久
- 客房转包协议书范本
- 自愿补偿修车协议书
- 装修验收申请协议书
- 老公外出喝酒协议书
- 土方回填施工记录表
- 旋挖钻机基坑支护工程施工隐患排查治理清单
- 空调维保质量保障体系及措施方案
- 平面向量在三角函数中的应用(学案)
- 中药的道地药材课件
- 幼儿园《3-6岁儿童学习与发展指南》健康领域知识试题及答案
- 国家职业技能标准 (2021年版) 婴幼儿发展引导员
- 幼儿园小班科学:《小鸡和小鸭》 PPT课件
- 伯努利方程-ppt课件
- 年产20吨阿齐沙坦原料药生产车间的设计和实现材料学专业
- 电子公章模板
评论
0/150
提交评论