第二章解线性方程组的直接法解说.pptVIP

* 例9 用改进的平方根法求解方程组 分析 对称很显然，需验证正定，可通过验证各顺序主子式大 解 1. 分解 2.求解计算 到方程组的解： 于零. * 追赶法(Thomas算法) 对角占优矩阵: 补充 * 有一类方程组,在今后要学习的插值问题和边值问题中 有着重要的作用,即三对角线方程组,其形式为: 其中 --------(1) * * 以下以Doolittle分解导出三对角线方程组的解法 以Crout分解的三对角线方程组的解法请自己思考 设 用紧凑格式的 Doolittle分解 * 因此 二对角阵 ---(2) * 由 --------(3) --------(4) --------(5) --------(6) 可得L和U的元素的计算公式 * 得 --------(7) * 得 --------(8) 也称Thomas法 追赶法的计算流程 这个循环称之为追的 过程,相当于消元过程 这个循环称之为赶的 过程,相当于回代过程。 求解 Ax = f * 例1. 用追赶法解三对角线方程组 解: * * * 因此原线性方程组的解为 总结 事实上，追赶法的求解过程就是将系数矩阵分解两个简单的二对角线矩阵，从而归结为求解两个简单三角形方程组的过程。 追赶法也可由高斯消去法计算过程导出，但 考虑到方程组的特点，计算时会把大量零元 素撇开，从而大大节省计算量。 ■④ 2.3矩阵的条件数和方程组的性态 * 首先引入向量和矩阵的范数。 设 是n维向量空间， 分别定义x的非负实函数如下： 可以验证它们都是 上的范数。 * 定理1 设 是 的某种范数，则 是分量 的连续函数。 定理2 设 及 是 上的两种向 量范数，则有与x无关的常数 ，使 定理3 上所有范数是等价的。 * 接下来看矩阵的范数。 定义 矩阵 的非负实函数 若满足 则称 是矩阵A的范数。 * 设 ，常用矩阵范数有 * 如果对矩阵范数 和向量范数 ，有 则称矩阵范数 和向量范数 相容。 另外介绍一种由向量范数导出的矩阵范数 定义 设 为向量范数， 定义 矩阵A的非负函数 可验证它是一矩阵范数，称为从属范数。 * 可以验证，向量范数 所导出的矩阵范 数 与该向量范数是相容的。 定理 矩阵范数 分别是向量 范数 的从属范数。 下面来讨论线性方程组的解对系数矩阵和 右端常数项的敏感性问题。 * 的准确解是（1，1）T。若A及b作微小的变化，扰动后的方程组为 其准确解为（-2，10）T 。 先看一个例子，说明方程组的解对扰动的敏感性。 例 方程组  此例中， A和b的微小变化引起x很大的变化， x对A和b的扰动 是敏感的。 这种现象的出现完全是由方程组的性态决定的。 * 定义 如果方程组Ax= b中，矩阵A和右端b的微小变化，引起解向量x的很大变化，则称A为关于解方程组和矩阵求逆的病态矩阵，称相应的方程组为病态方程组。否则， 称A为良态矩阵，称相应的方程组为良态方程组。 　　我们需要一种能刻画矩阵和方程组病态程度的标准。 定理 设有n阶非齐次方程组Ax=b，其中A非奇异。设系数矩阵和右端项有小扰动 和 ，扰动后的方程组为 则当 时，有 * 将上式两端取范数，则有 证. 将代入扰动方程组 ，整理后有 经整理后得 由于 ，即得所证。 再利用 ，则有 * 如果 此时有 如果 ，此时有 * 为矩阵A的条件数。 　　如果矩阵范数取2范数，则记 　　定义 设A∈Rn×n为可逆矩阵，按算子范数，称 Cond（A）= 同样可以定义cond∞（A）和cond1（A）。