版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
雅可比迭代法求解线性方程组的实验报告一、实验题目分别利用雅可比迭代法和高斯-塞德尔迭代法求解以下线性方程组:使得误差不超过0.00001。二、实验引言1.实验目的=1\*GB3①掌握用迭代法求解线性方程组的基本思想和步骤,熟悉计算机fortran语言;=2\*GB3②了解雅可比迭代法在求解方程组过程中的优缺点。2.实验意义雅克比迭代法就是众多迭代法中比较早且较简单的一种,求解方便实用。三、算法设计1.雅可比迭代法原理:
设有线性方程组Ax=b满足,将方程组变形为:x=Bx+f,则雅可比(Jacobi)迭代法是指,即由初始解逐步迭代即可得到方程组的解。算法步骤如下:步骤1.给定初始值,精度e,最大容许迭代次数M,令k=1。步骤2.对i=1,2,…,n依次计算步骤3.求出,若,则输出结果,停止计算。否则执行步骤4.步骤4.若转步骤2继续迭代。若表明迭代失败,停止计算。2.算法流程图四、程序设计programjacobiimplicitnoneinteger::i,jinteger::ksavekreal,parameter::e=0.001integer,parameter::n=3real::x(n),y(n),b(n)datab/7.2,8.3,4.2/real::Dreal::a(n,n)五、结果及讨论1.实验结果**********矩阵A的形式为**********10.00-1.00-1.00-1.0010.00-1.00-2.00-2.005.00迭代次数为:1迭代次数为:2迭代次数为:3迭代次数为:4迭代次数为:5迭代次数为:6迭代次数为:7****************************************用jacobi方法解得的结果X[t]为:1.101.201.302.讨论分析(1)误差从上述输出结果中可以看出,当迭代次数k增大时,迭代值x1,y1,z1会越来越逼近方程组的精确解x=1.0,y=1.2,z=1.3。(2)收敛性在本题目中,用雅可比迭代法和高斯-塞德尔迭代法分别求解该线性方程组,得到的近似根是收敛的六、算法评价优点:迭代法算法简单,编制程序比较容易。缺点:迭代法要求方程组的系数矩阵有某种特殊性质(譬如是所谓对角占优阵)以保证过程的收敛性。高斯—塞德尔迭代法比雅可比迭代法收敛快(达到同样的精度所需迭代次数少),但这个结论,在一定条件下才是对的,甚至有这样的方程组,雅可比方法收敛,而高斯—塞德尔迭代法却是发散的。在雅可比迭代法求解线性方程组时,只要误差截断设计的合理,原则上可以得到很正确的解。而通常我们选取设计误差限或设计最大迭代次数的方法来控制。由于它的准确性,故在实际应用中比较常见,对于解一般线性方程组非常有效准确。通过该算法以及编程对求解的过程,我们不难发现,\o"编辑本段"雅克比迭代法的优点明显,计算公式简单,每迭代一次只需计算一次矩阵和向量的乘法,且计算过程中原始矩阵A始终不变,比较容易并行计算。然而这种迭代方式收敛速度较慢,而且占据的存储空间较大,所以工程中一般不直接用雅克比迭代法,而用其改进方法。附:高斯—赛德尔程序programG-Simplicitnoneinteger::i,jinteger::ksavekreal,parameter::e=0.001integer,parameter::n=3real::x(n),y(n),b(n)datab/7.2,8.3,4.2/real::Dreal::a(n,n)open(unit=10,file='1.txt')dataa/10,-1,-1,-1,10,-1,-2,-2,5/write(10,*)"**********矩阵A的形式为**********"write(10,"(1x,3f6.2,/)")aforall(i=1:n)x(i)=0endforallk=0100D=0doi=1,n y(i)=b(i) doj=1,n if(i<j)y(i)=y(i)-a(i,j)*x(j) if(i>j)y(i)=y(i)-a(i,j)*y(j) enddo y(i)=y(i)/a(i,i) enddo doj=1,n D=abs(x(j)-y(j)) enddoforall(i=1:n)x(i)=y(i)endforall if(D>=e)then k=k+1 write(10,*)"迭代次数为:",k goto100 else goto200endif200write(10,*)"****************************************"write(10,*)"用Gauss-seidel方法解得的结果X[t]为:" write(10,"(1x,3f6.2,/)")x(:) stop endprogram**********矩阵A的形式为**********10.00-1.00-1.00-1.0010.00-1.00-2.00-2.005.00迭代次数为:1迭代次数为:2
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 会计业务知识测试题库含答案
- 中小学生预防校园欺凌和处理制度措施
- 陕西省西安市长安第一中学2025年高考复习全程精练模拟卷(全国卷)物理试题含解析
- 煤矿自动化控制技术在皮带运输中的应用
- 陕西韩城象山中学2025年高三下学期月考(一)物理试题试卷含解析
- 山西省运城市永济市重点名校2024-2025学年初三下学期半期联合考试数学试题含解析
- 河南省商丘市民权县2023-2024学年八年级上学期期中英语试卷
- 企业去杠杆、投资与创新:一个投资者保护的视角
- 《ISO56001-2024创新管理体系 - 要求》之16:“7支持-7.1资源”解读和应用指导材料(雷泽佳编制-2024)
- 工艺美术专业工艺品设计方案课程标准
- 南京市2021-2022学年高一上学期地理期末统考试卷(原卷版)
- 科研伦理与学术规范-课后作业答案
- 钣金件设计工艺设计规范
- 湖南文艺出版社小学六年级上册(简谱)音乐全册教学课件
- (新版)非阿片类镇痛药治疗慢性疼痛病中国指南
- 不变与万变:葛剑雄说国史
- 新版入团志愿书表格(含申请书范本)
- 《行路难》-初中语文
- 病理生理学教案大全
- (2023)垂直轴微风发电机制造生产建设项目可行性研究报告写作模板-立项备案(一)
- 金属工艺学第一篇 金属材料的基本知识
评论
0/150
提交评论