chapter6解线性代数方程组直接法和迭代法2018.pptxVIP

第6章 解线性代数方程组咱从初中就开始会啦！关于期末考试答疑;在科学研究和工程计算中经常会遇到求解线性代数方程组的问题。 许多有效的数值计算方法中，作为关键的一步就需要求解线性代数方程组，而且有可能要求解的未知数数目很大。因此线性代数方程组的数值解法是数值计算重要内容。;;在计算机上解线性代数方程组的方法分为二类，一类是直接法，一类是迭代法。 由于未知数数目n很大，计算的每一步都含有不可避免的舍入误差，因此有必要介绍多种??线性代数方程组的计算方法，去适应各种实际问题的需要。 本章先讨论解线性方程组的直接法。 ;对于中小型方程组（150？？），直接法是有效的；大型带状矩阵也可直接法。 对于阶数较高而又不具有带状结构的方程组，常用迭代法； 特别是高阶方程组中，系数矩阵为大型无规律稀疏矩阵，直接法难以克服存储问题，迭代法不需要存储系数矩阵的零元素。 非线性方程数值求解，也主要用迭代法。;高性能科学计算 与方程组求解;全球尺度模拟需要巨量的网格。;;D02华南： 网格距△x=9km， 时间步长约△t=200s， 网格数148X160=23680;教师 姓名;Color surfaces (red is positive, blue is negative) of the instantaneous streamwise vorticity in DNS of flow over a cylinder at Re=300 ;城市风环境模拟，对网格有较大要求;15;16;研究揉面，是有背景的。 多孔材料的其中一个重要应用是作为能量吸收器使用。在碰撞事故中，多孔材料变形时能够保持较低应力水平同时吸收大量能量，从而保护人或设备。 多孔材料单位体积所吸收能量，通常用名义应力对应变的积分计算，这就需要监测整个压缩过程中材料应力及应变的变化，所以传统上评估多孔材料吸能能力是一个较为繁琐的工作。。。;举例、高性能科学计算：推荐 西安交大 陶文铨院士 数值传热学;1、计算出来的结果不一定是对的，需要验证，需要按照标准计算；三种论文引用较多！2、现状：模式研究多，实验研究少。（ 2012年初来大气的体验，审稿中的问题）3、有时，不得不先能计算，再求准确性。（差一个量级-?差2-3倍-? 20%以内。。。）;;以香港为例子，控制通风与隔离和接种疫苗的比较;§6.1 高斯消元法（Gauss）;;;;;高斯消元法的简单流程;算法的时间复杂性分析:;三对角系数矩阵的高斯消元法;高斯消元法的时间复杂度(三对角方程组);列主元消元法;存在易求解的3种线性方程组（上三角、下三角、对角）;把一个n阶矩阵A分解成结构简单的三角形矩阵的乘积称为矩阵的三角分解。;;高斯消元法可求解的充要条件:;;;;比较等式两边第二行;;已知u的前k-1行和L的前k-1列，比较两边第k 行;;摘自《数值计算基础》 沈剑华， 同济大学出版社;;; 算法分析: 三角分解法求解方程组的时间复杂度也为O(n3)量级。 应当指出A为非奇异矩阵时，在用直接三角分解法解方程组Ax=b的过程中，有可能出现绝对值ukk很小的情况，使计算严重失真。;6.3 线性方程组的直接法 的误差分析;§5.3 向量和矩阵的范数、方程组的性态、条件数;一、向量和矩阵的范数;;有了向量范数，就可以用它来表示向量误差，刻画近似解的可靠程度。设x*为真正解，x为计算解，后者的误差表示为||x-x*||或||x-x*||/||x||，其中||x-x*||称为绝对误差，||x-x*||/||x||称为相对误差。;矩阵的范数 定义;;;在应用问题中，矩阵和向量的乘积经常出现，因而应让所用的矩阵范数和向量范数有某种关系。 定义3 对任意矩阵A?Rn?n和任意向量x?Rn，若满足条件||Ax||?||A||?||x||，则所定义的矩阵范数与向量范数是相容的。 常用范数的相容关系： ;;二 方程组的性态、条件数;例：系数矩阵A为希尔伯特(Hilbert)矩阵;当一个方程组由于初始数据的小误差，而使计算结果严重失真时，这样的问题就是一个“坏”问题。 因此讨论线性代数方程组性态的办法就是给方程组的系数阵或右端一个小误差，然后考虑方程组解的变化。;只考虑b的误差时;只考虑A的误差时;;同时考虑A和b的误差时 （教材p164）; 由于Cond(A)=||A||?||A-1||≥||AA-1||=||I||=1,可见条件数是一个放大倍数。 上面讨论中所给出的系数矩阵A或右端b的扰动，是在形成线性代数方程组时，原始数据的误差和运算过程中的舍入误差造成的。因此求解的对象总是有扰动的方程组，得到的解是近似方程组的解，这个解的好坏完全由条件数确定。因而系数矩阵的条件数刻画了线性代数方程组的性态。;教材 p166;对一个确