- 1、本文档共28页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
* * * * * 线性代数方程组的数值解法(1) Gauss 消去法 (Demos in Matlab: airfoil in 2D) 线性代数方程组的数值解法 直接法:Gauss 消去法,SuperLU 迭代法:定常迭代(Jacobi, GS, SOR, SSOR) Krylov 子空间方法(CG, MINRES , GMRES, QMR, BiCGStab) The Landscape of Ax=b Solvers Direct A = LU Iterative y’ = Ay Non- symmetric Symmetric positive definite More Robust Less Storage More Robust More General 刘徽 (约220-280) Gaussian elimination, which first appeared in the text Nine Chapters on the Mathematical Art written in 200 BC, was used by Gauss in his work which studied the orbit of the asteroid Pallas. Using observations of Pallas taken between 1803 and 1809, Gauss obtained a system of six linear equations in six unknowns. Gauss gave a systematic method for solving such equations which is precisely Gaussian elimination on the coefficient matrix. (The MacTutor History of Mathematics, http://www-history.mcs.st-andrews.ac.uk/history/index.html) Gauss(1777-1855) 今有上禾三秉,中禾二秉,下禾一秉,实三十九斗;上禾二秉,中禾三秉,下禾一秉,实三十四斗;上禾一秉,中禾二秉,下禾三秉,实二十六斗。问上、中、下禾实一秉各几何?答曰:上禾一秉九斗四分斗之一。中禾一秉四斗四分斗之一。下禾一秉二斗四分斗之三。-------《九章算术》 一个两千年前的例子 Basic idea: Add multiples of each row to later rows to make A upper triangular 一个两千年前的例子(2) Solving linear equations is not trivial. Forsythe (1952) After k=1 After k=2 After k=3 After k=n-1 Gauss 消去过程图示 用矩阵变换表达消去过程 利用Gauss 变换矩阵的性质: 用矩阵变换表达消去过程(2) 单位下三角形 用Gauss消去法求解 A x=b --- LU分解 A = L U (cost = 2/3 n3 flops) --- 求解 L y = b (cost = n2 flops) --- 求解 U x = y (cost = n2 flops) 版本一 for k = 1 to n-1 for i = k+1 to n m = A(i,k)/A(k,k) for j = k to n A(i,j) = A(i,j) - m * A(k,j) 算法实现: Gauss Elimination Algorithm for k = 1 to n-1 … 对第k列,消去对角线以下元素 … (通过每行加上第k行的倍数) for i = k+1 to n … 对第k行以下的每一行i for j = k to n … 第k行的倍数加到第 i 行 A(i,j) = A(i,j) - (A(i,k)/A(k,k)) * A(k,j) 版本二: 在
您可能关注的文档
- 任务一网络设备登录方式-console口登录.doc
- 任意分布高速仿核信号发生器.pdf
- 三、施工进度计划和主要技术措施.doc
- 三个月340分到590分非常有效学习方法.doc
- 第1篇_行列式.pdf
- 第1篇Visual c++ 2010开发环境简介.ppt
- 三角函数最值问题几种常见类型.doc
- 第1篇数字逻辑基础.ppt
- 三角网数字地面模型和程序实现.doc
- 三态缓冲器74系列芯片的型号区别及功能略表.pdf
- (4篇)XX区抓党建促基层治理培训心得体会汇编12.docx
- 汇编1154期-在培训班上的讲话汇编(3篇).doc
- 汇编1173期-专题党课讲稿汇编(3篇)112.doc
- 汇编1076期-主题党课讲稿汇编(3篇).doc
- 汇编1177期-学习心得体会汇编(3篇)112.doc
- (6篇)党和国家机构改革心得体会汇编.docx
- 汇编1166期-坚定理想信念、全面从严治党、担当作为专题党课讲稿汇编(3篇)112.doc
- 汇编1174期-专题党课讲稿汇编(3篇)112.doc
- 教育13期-主题教育学习心得体会、研讨发言材料参考汇编(3篇).doc
- 汇编1458期-心得体会研讨发言提纲参考汇编(3篇)123.doc
最近下载
- 高桩码头标准化建设成果图集.pdf VIP
- 施工管理标准化图集最终版.pdf VIP
- 包装设计智慧树知到期末考试答案章节答案2024年浙江艺术职业学院.docx
- 2025年河南交通职业技术学院单招职业倾向性测试题库(考点梳理).docx VIP
- 胃肠外科“三基”考试卷含答案.doc
- 深入学习贯彻中央八项规定精神学习研讨发言材料.docx VIP
- 2022年汽车修理工(技师)证考试题库汇总(含真题和典型题).docx
- 52号令 水利工程质量管理规定 20230130.pdf VIP
- 供电营业规则复习试题附答案.docx
- 2025年中央八项规定精神学习教育党课讲稿+以学习教育为契机,深入贯彻中央八项规定精神.doc VIP
文档评论(0)