- 1、本文档共37页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
清华大学数值分析实验报告
数值分析实验报告 实验3.1题目:考虑线性方程组,,,编制一个能自动选取主元,又能手动选取主元的求解线性代数方程组的Gauss消去过程。(1)取矩阵,,则方程有解。取计算矩阵的条件数。分别用顺序Gauss消元、列主元Gauss消元和完全选主元Gauss消元方法求解,结果如何?(2)现选择程序中手动选取主元的功能,每步消去过程都选取模最小或按模尽可能小的元素作为主元进行消元,观察并记录计算结果,若每步消去过程总选取按模最大的元素作为主元,结果又如何?分析实验的结果。(3)取矩阵阶数n=20或者更大,重复上述实验过程,观察记录并分析不同的问题及消去过程中选择不同的主元时计算结果的差异,说明主元素的选取在消去过程中的作用。(4)选取其他你感兴趣的问题或者随机生成的矩阵,计算其条件数,重复上述实验,观察记录并分析实验的结果。1. 算法介绍首先,分析各种算法消去过程的计算公式,顺序高斯消去法:第k步消去中,设增广矩阵中的元素(若等于零则可以判定系数矩阵为奇异矩阵,停止计算),则对k行以下各行计算,分别用乘以增广矩阵的第行并加到第行,则可将增广矩阵中第列中以下的元素消为零;重复此方法,从第1步进行到第n-1步,则可以得到最终的增广矩阵,即;列主元高斯消去法:第k步消去中,在增广矩阵中的子方阵中,选取使得,当时,对中第行与第行交换,然后按照和顺序消去法相同的步骤进行。重复此方法,从第1步进行第n-1步,就可以得到最终的增广矩阵,即;完全主元高斯消去法:第k步消去中,在增广矩阵中对应的子方阵中,选取使得,若或,则对中第行与第行、第列与第列交换,然后按照和顺序消去法相同的步骤进行即可。重复此方法,从第1步进行到第n-1步,就可以得到最终的增广矩阵,即;接下来,分析回代过程求解的公式,容易看出,对上述任一种消元法,均有以下计算公式:2. 实验程序的设计一、输入实验要求及初始条件;二、计算系数矩阵A的条件数及方程组的理论解;三、对各不同方法编程计算,并输出最终计算结果。3. 计算结果及分析(1)先计算系数矩阵的条件数,结果如下,可知系数矩阵的条件数较大,故此问题属于病态问题, b或A的扰动都可能引起解的较大误差;采用顺序高斯消去法,计算结果为:最终解为x=(1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000001, 0.999999999999998, 1.000000000000004, 0.999999999999993, 1.000000000000012, 0.999999999999979, 1.000000000000028)T使用无穷范数衡量误差,得到=2.842170943040401e-14,可以发现,采用顺序高斯消元法求得的解与精确解之间误差较小。通过进一步观察,可以发现,按照顺序高斯消去法计算时,其选取的主元值和矩阵中其他元素大小相近,因此顺序高斯消去法方式并没有对结果造成特别大的影响。 若采用列主元高斯消元法,则结果为:最终解为x=(1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000)T同样使用无穷范数衡量误差,有=0; 若使用完全主元高斯消元法,则结果为最终解x=(1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000, 1.000000000000000)T同样使用无穷范数衡量误差,有=0;(2)若每步都选取模最小或尽可能小的元素为主元,则计算结果为最终解x=(1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000001 0.999999999999998 1.000000000000004 0.999999999999993 1.000000000000012 0.999999999999979 1.000000000000028)T使用无穷范数衡量误差,有为2.842170943040401e-14;而完
您可能关注的文档
- 法学英语(二)书上翻译题答案.doc
- 河流地貌的形成课时练习(含答案).doc
- 洋里学区洋里中心小学课堂教学设计综合实践.doc
- 流通概论形成性考核册答案(全).doc
- 浅析中国古代宰相制度的演变史.doc
- 浅xi谢灵运山水诗的艺术特色.doc
- 浅析互助激励在美术课堂上的效应.doc
- 流泗镇山洪灾害防御预案.doc
- 浅析商业地产的运营管理模式.docx
- 流星灯产品规格书.docx
- (高清版)DB41∕T 2772-2024省级公益林建设技术规程.pdf
- (高清版)DB37∕T 5305-2024 回收水泥浆在混凝土中应用技术规程.pdf
- (高清版)DB37∕T 4817—2025 沥青路面就地冷再生技术规范.pdf
- (高清版)DB37∕T 4601.1—2024 服务业用水定额 第1部分:批发零售、运输仓储、餐饮、居民服务、洗车及体育.pdf
- (高清版)DB37∕T 4715—2024 水利工程泵站信息化项目建设指南.pdf
- (高清版)DB37∕T 4722—2024 城市轨道交通工程建设BIM应用指南.pdf
- (高清版)DB37∕T 4744—2024 海洋自然资源资产核算技术指南.pdf
- (高清版)DB37∕T 4755.1—2024 公共数据共享+第1部分:基本要求.pdf
- (高清版)DB37∕T 4772—2024 公路工程废旧混凝土再生集料混凝土应用技术规程.pdf
- (高清版)DB37∕T 4784—2024 3岁以下婴幼儿生活照护托育服务规范.pdf
最近下载
- 长期卧床患者的体位护理技巧.pptx VIP
- JBT14159.3-2022 空气分离设备能效限额第3部分:液化设备.pdf VIP
- 《民用建筑设计通则》图示.docx VIP
- 灭火应急疏散预案.doc VIP
- 铁路内燃、电力司机考试题库.doc VIP
- 新概念英语第一册 Lesson 55-56 课件.ppt VIP
- 2025AI+BI峰会:DeepInsight Copilot 演进史以及未来探索.pdf VIP
- 《临床药学文献阅读汇报抗感染粘菌素在治疗鲍曼不动杆菌感染vap中的应用》课件模板.ppt VIP
- 2025年福建公需课答案.pdf VIP
- 2020四川园林绿化工程说明及计算规则.pdf VIP
文档评论(0)