[理学]线性方程组求解的直接法.pptVIP

第三章 线性方程组求解的数值方法 线性方程组的基本概念 线性方程组的基本概念 线性方程组的基本概念 克莱姆（Cramer）法则 克莱姆（Cramer）法则 高斯消去法（Gaussian elimination） 高斯消去法步骤： 高斯消去法步骤： 高斯消去法步骤： 高斯消去法： 高斯消去法 高斯消去法 例2：求解线性方程组 全/列主元消去法 如果高斯消去过程中对角线出现“0”元素怎么办？ 通过行变换和列变换将对角线元素换为非零元素。 列主元消去法 例3：求解线性方程组 列主元消去法 例4：求解线性方程组 全主元消去法 例5：求解线性方程组 高斯消去法 例6，第43页 3.2，解方程组 矩阵LU分解法求解线性方程组 矩阵LU分解法求解线性方程组 例：利用LU分解求解如下方程： LU分解计算方法 Matlab的符号计算 直接计算LU分解方法 例，第67页，习题1 正定对称矩阵 正定矩阵：一个矩阵 A 称为正定阵，如果 对任意非零向量 x 都成立。 对称正定矩阵Cholesky分解： 对称正定矩阵Cholesky分解 说明：上一页的证明过程是一个构造式的证明方法，不但得到了结论，也给出了LU分解和Cholesky矩阵分解的关系： 已知LU分解求cholesky分解： C = D-1U， 已知cholesky分解，求LU分解： U = DC， L = CTD-1 。 如何计算D？ U对角线元素开方，或C的对角线元素。 对称正定矩阵Cholesky分解 正定对称矩阵 为什么要单独讨论正定对称矩阵？ 很多问题都是正定对称矩阵：协方差矩阵，非负函数的代数矩， 正定对称矩阵 线性方程组的定义。 线性方程组解的情况：无解，有唯一解，有无穷多解 线性方程组通解的计算：“null”函数 将一般方程组转化为线性方程组 高斯消去法 全/列主元消去法 LU分解法求解线性方程组 LU分解的计算方法 Matlab的符号计算功能 对称正定矩阵的cholesky分解 已知LU分解，如何计算cholesky分解。 作业： P. 88： 思考题：2 习题：4， * 课程回顾 绝对误差、绝对误差限、相对误差、相对误差限、有效数字； 浮点数的概念； 算法的定义 算法研究的基本问题： 算法的构造 算法的误差 算法的稳定性 算法的运算量 第一节 线性方程组求解的直接法 线性方程组的解： 无解，（过定方程） 有唯一解， 有无穷多解，（欠定方程） 本章研究的线性方程组： 矩阵A为方阵 方程有唯一解 通解求解方法： Matlab中“null”函数可计算欠定方程 Ax=0 的通解。 将一般方程组转化为线性方程组： 如：取log，将乘法转化为加法，得到线性方程组。 高斯消去法思路： 上/下三角矩阵方程易于求解； 故将一般矩阵方程化为等价的三角矩阵方程，然后反向代入求解。 等价变换（行初等变换）： 倍加变换：把某行换成它本身与另一行的和； 对换变换：把两行互换； 倍乘变换：把某行所有元素乘以同一非零数； 注意：不能进行列初等变换； 性质1：若A的所有顺序主子式均不为0，则高斯消元无需换行即可进行到底，得到唯一解。 性质2：只要 A 非奇异，即 A?1 存在，则可通过逐次消元及行交换，将方程组化为三角形方程组，求出唯一解。 什么样的矩阵可以采用高斯消去法？ 例1：求解线性方程组 注意： 一定要等式两边的系数同时变。 高斯消去法第一步：写出增广矩阵 全主元消去法： 每一步选子矩阵中绝对值最大的元素为主元素 主元选择公式： 列主元消去法： 省去换列的步骤，每次仅选子矩阵中第k列最大的元（仅作行变换）。 主元选择公式： 注意：列交换改变了 xi 的顺序，须记录交换次序，解完后再换回来。 注意： 如果对角线元素为0，则需要交换元素。列主元消去法值选择该列 k~n 元素中最大值。 怎么换行？是不是高斯消去算法的问题？方程解到底是多少？ 全主元消去法选择子阵中最大元素。 注意：列交换改变了 xi 的顺序 即使主元素不为0，但是比较小，也应该换行 高斯消去法推广： 将矩阵A分解为A=LU,其中，L为下三角矩阵，U为上三角矩阵 则方程Ax=b可写为：LUx=b，令y=Ux，方程可转换为两个方程组： 依次求解y和x的得到原方程的解。 解：1、将矩阵A进行LU分解 解：2、求解方程组：Ly=b 解：3、求解方程组：Ux=y 按颜色顺序依次计算 给定矩阵A，如何计算其对应的L矩阵和U矩阵？ 思考题：67页习题2，写出4阶三对角矩阵的LU分解公式 LU分解矩阵 LU分解： 不是所有矩阵都可分解为A=LU！更一般形式： 思考：Matlab中lu为LU分解函数。 调用方式：“[L,U,P] =