- 1、本文档共32页,可阅读全部内容。
- 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 二、实验内容‥ ‥ ‥ ‥ ‥ ‥ ‥ 4 三、实验要求‥ ‥ ‥ ‥ ‥ ‥ ‥ 5 四、实验过程‥ ‥ ‥ ‥ ‥ ‥ ‥ 6 五、结果分析‥ ‥ ‥ ‥ ‥ ‥ ‥ 19 六、程序代码‥ ‥ ‥ ‥ ‥ ‥ ‥ 21 七、课程信息‥ ‥ ‥ ‥ ‥ ‥ ‥ 30 一、实验分析 关于线性方程组的数值解法一般有两类:直接法和迭代法。而对于系数矩阵较特殊的一类线性方程组,解法也相对有所不同。在实际问题中,例如解常微分方程边值问题,解热传导问题等,都会要求解系数矩阵为对角占优的三对角方程组。 如何求解三对角方程组就成为关键。三对角线性方程组是系数矩阵为三对角矩阵的一类方程组,其数值解法有:追赶法,Jacobi迭代法,Gauss迭代法和SOR法。 二、实验内容 考虑线性方程组: 其中A为三对角矩阵,编制程序求解该方程组 三、实验要求 (1)考虑不同的数值解法,对方程组进行编程,编写其通用代码; (2)对不同的程序代码用实例进行验证,取矩阵: 方程初值 : 。 (3)比较不同方法对求解三对角线性方程组的优异性。 四、实验过程 (一)追赶法 编程思想: 矩阵的三对角线性方程组是由矩阵的直接三角分解法来推到的来的,求解等价于解两个三角形方程组: (1) ,求y;(2) ,求x 从而得到解三对角线方程组的追赶法公式。 1. 计算的递推公式: (i=2,3,……n-1); 2. 解 (i=2,3,……n); 3. 解 (i=n-1,n-2,……,2,1)。 运行结果:(点击查看源程序) a=-1*ones(1,4); c=a; b=2*ones(1,5); f=[1 0 1 0 0]; [x]=ZhuiGan(a,b,c,f)‘ n =5 x = 1.33333333333333 1.66666666666667 2.00000000000000 1.33333333333333 0.66666666666667 注:其中a为对角下向量,b为对角向量,c为对角上向量,f为方程常 系数。x为方程组的解。 (二)迭代法 对于 ,A为非奇异矩阵,除了可用追赶法计算之外,也可建立迭代法求解。迭代法的一般格式为: ,其中B为迭代矩阵,对由迭代格式产生迭代序列 ,若 ,则 即 为 方程 的解。 迭代法有以下几种方法:Jacobi迭代法,Gauss迭代法和SOR法。 Jacobi迭代法 1.迭代过程: 对于 ,设A可逆,且对于 令D为A的对角元素部分,则 ,继而对于方程组有 ,所以有解 ,令 为迭代矩阵, 。则Jacobi迭代法的迭代格式为: 。 2.算法设计: (1)将方程组系数矩阵用 表示,常数项用 表示。 (2)用 存储初值; (3)用库函数diag,tril,triu求取对角矩阵D,下三角矩阵L,上三角矩阵U; (4)利用得到的D,L,U,求取迭代矩阵J,和f; (5)循环判断, 存取循环值,用二范式迭代循环,当 时,终止循环。 3.运行结果:(点击查看源程序)
您可能关注的文档
- 《雪孩子》图片字体改良版.ppt
- 【公司金融】第讲 资本结构基本概念.ppt
- 【化学】《 化学平衡的移动》课件苏教版选修.ppt
- 【中医外科学】瘤岩.ppt
- 【人教版】必修一地理《地球的圈层结构》课件.ppt
- 「先期品质管理AQM基础」.ppt
- 【】 【】年美菱冰箱最新功能卖点汇总.ppt
- 【地理】湘教版选修 第二章 第二节 我国的干旱洪涝寒潮与台风课件.ppt
- 【地理】高考复习专题训练课件专题时空分析第节季节判断.ppt
- 【学渣君高考历史复习系列】世纪现代化模式的探索—年.doc
- 剧本杀市场新变革:2025年商业模式创新与消费者行为洞察.docx
- 2025数字藏品交易监管政策及市场分析报告.docx
- 2025年航空电子设备国产化替代产业链产业政策环境分析报告.docx
- 2025年新能源汽车电池回收利用政策法规与市场前景分析报告.docx
- 2025年航天军工企业民品市场拓展策略报告.docx
- 数字藏品市场准入机制与行业创新能力研究报告.docx
- 2025年短视频平台算法推荐机制对青少年认知发展模式的优化策略报告.docx
- 2025剧本杀行业人才流失原因及行业人才流失风险控制研究报告.docx
- 2025年智慧社区医疗平台健康管理服务与医疗资源整合报告.docx
- 新能源电动巡逻车在智慧城市建设中的AI算法应用前景报告.docx
文档评论(0)