线性方程组的元解法 - 吉林大学.pptVIP

计算方法 第一章解线性方程组的直接法 误差 绝对误差 相对误差 计算方法 第一章解线性方程组的直接法 有效数字 如果 　 　则说x*近似表示x准确到第n位，并从第n位起直到最左边的非零数字之间的一切数字都称为有效数字，并把有效数字的位数称为有效位数。 如果 　 　 　则说近似值x*具有n位有效数字，这里n为正整数，m为整数。 计算方法 第一章解线性方程组的直接法 定理 若近似值x*具有n位有效数字，则其相对误差满足 反之，若x*的相对误差满足 则x*至少具有n位有效数字 * 第一次消元 * 计算方法 第一章解线性方程组的直接法 第一章解线性方程组的直接法 计算方法 第一章解线性方程组的直接法 求解Ax=b直接法是指在无舍入误差存在的情况下,经过有限步运算即可求得精确解的算法,因此又称精确法.因为舍入误差的存在,精确解也是不精确的.直接法的典型代表是Gauss消元法 计算方法 第一章解线性方程组的直接法 消元过程 回代过程 计算方法 第一章解线性方程组的直接法 按照矩阵变换的观点来描绘消元的过程 分别表示n维实和复向量空间，用R n×n 表示n×n阶实矩阵空间 考虑线性方程组 Ax=b 其中 计算方法 第一章解线性方程组的直接法 矩阵形式为 其增广矩阵为 计算方法 第一章解线性方程组的直接法 第一步消元过程相当于用矩阵G1左乘增广矩阵，即 它可以记为　　　　　 其中　　　　　是将A中的第一列中的元素a11换成0而得到的列向量， 是单位矩阵I的第一列所形成的列向量，1/2是a11的倒数。这种形式的矩阵称为Gauss矩阵。 计算方法 第一章解线性方程组的直接法 同样，若取Gauss矩阵为 则有 从上述讨论看出，消元过程等价于将方程组的增广矩阵依次左乘相应的Gauss矩阵，将其化为上三角形式 计算方法 第一章解线性方程组的直接法 下面，按照上述思想推导消元法的一般算式 对于一般的矩阵A=[aij]，设a11≠0，令 构造Gauss矩阵 用G1左乘a1得 从而G1(A,b)具有下列形式， 其中 计算方法 第一章解线性方程组的直接法 一般地，如果已经利用Gauss矩阵G1, …,Gk-1得到 则当 时,取 计算方法 第一章解线性方程组的直接法 计算方法 第一章解线性方程组的直接法 如此继续下去,最后,当 时得到 其中 上述过程结束后,方程组化成了具有上三角系数矩阵的方程组 当U非奇异，即uii≠0(i=1,2, …,n)时，方程组容易求解。 计算方法 第一章解线性方程组的直接法 Gauss消元法的计算量 回代过程所用的乘除法次数: 消元过程所用的乘除法次数: 计算A(k+1)和b(k+1)时,需要计算(n-k)个lik和(n-k)(n-k+1)个元素 , 而计算每一个元素只需要用一次乘(除)法,所以第k步消元过程所用乘(除)法的次数为(n-k)(n-k+2) 这样，n-1步消元过程及回代过程所需的乘除法次数的总和为 计算方法 第一章解线性方程组的直接法 Gauss消元法能够进行到底的条件 从上述消元过程可以看出:方程组Ax=b能够按照自然顺序将消元过程进行到底的充要条件是 为了进一步讨论这个条件与矩阵A的关系，将下式中的A(k+1)写成分块形式 其中Ak是A(k+1)的k阶顺序主子矩阵 计算方法 第一章解线性方程组的直接法 由于 可将Gauss矩阵写成 其中， 不难验证 因而 计算方法 第一章解线性方程组的直接法 如果将A也写成分块形式 其中Akk为k阶矩阵,则显然有 由于Lk是单位下三角阵，其行列式det(Lk)=1, Ak是上三角阵，所以 计算方法 第一章解线性方程组的直接法 引理1.1　按自然顺序消元过程得到的 均不为零的充分必要条件是顺序主子矩阵 非奇异，并且 由上述引理，当k=n-1时得到下面的结果 定理1.1　按自然顺序消元过程可以实现的充分必要 条件是A的顺序主子矩阵　　　　　均为非奇异矩阵. 计算方法 第一章解线性方程组的直接法 1.2　矩阵的LU分解 定理1.2（LU分解）设A的前n-1个顺序主子矩阵非奇异，则存在单位下三角阵L及上三角阵U，使 　　　　　　　　A＝LU 　而且这样的分解是唯一的。 计算方法 第一章解线性方程组的直接法 LU分解（Doolittle分解） 计算方法 第一章解线性方程组的直接法 Doolittle分解算法的程序语言表述 对于i=1,2,…,n (1