_解线性方程组的直接法解说.ppt

第二章 求解线性方程组的数值解法 在第k步消元前，在系数矩阵第k列的对角线以下 的元素中找出绝对值最大的元。 回顾：对称正定阵A的几个重要性质 （1）A?1 亦对称正定，且 aii 0 （2）A 的顺序主子阵 Ak 亦对称正定 （3）A 的特征值 ?i 0 （4）A 的全部顺序主子式 det ( Ak ) 0 §2.1.3 有关定理 初等变换矩阵与Gauss变换矩阵 2. 消去过程系数矩阵的Gauss型变换 每一步消去过程相当于左乘Gauss变换矩阵 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 A 的 LU 分解 ( LU factorization ) -1 -1 -1 -1 -1 -1 -1 定理2.1.2： Gauss 消去法可顺利进行 A的顺序主子式 . Cor1 Cor2 3. 列主元Gauss消去法与矩阵的三角分解 定理2.1.3： 4. 矩阵三角分解的存在性与唯一性 LDU 分解 定理2.1.4： 非奇异阵A 存在唯一的LDU 分解 A的所有顺序主子式 均不为0. 注：此结果对|A|=0亦成立。 此时|Ak|=0(k=1,…n-1) §2.1.4 求解正定方程组的Cholesky方法/平方根法 定理2.1.5 设矩阵A对称正定，则存在唯一的单位下 三角阵L和对角元全为正的对角阵D 使得 A＝LDLT. 证明： 定理2.1.6 设矩阵A对称正定，则存在唯一的对角 元全为正的下三角阵G 使得 A＝GGT 证明： * * 用克莱姆法则（Cramer rule）求解n阶线性方程组 每个行列式由n!项相加，而每项包含了n个因子相乘，乘法运算次数为(n-1)n !次. 仅考虑乘(除)法运算,计算解向量包括计算n+1个行列式和n次除法运算,乘(除)法运算次数N=(n+1)(n-1)n!+n. 解线性方程组的两类数值方法: 直接法: 经过有限次运算后可求得方程组精确解的方法(不计舍入误差!) 迭代法：从解的某个近似值出发，通过构造一个无穷序列去逼近精确解的方法(一般有限步内得不到精确解) §2.1 解线性方程组的直接法 一、高斯消去法（简单高斯消去法/顺序高斯消去法） 消元法的基本思想就是通过对方程组作初等变换，把一般形式的线性方程组化为等价的易于求解的三角方程组。 首先将方程组Ax=b 化为上三角方程组，此过程称为消去过程，再求解上三角方程组，此过程称为回代过程. §2.1.1 高斯消去法和选主元高斯消去法 将增广矩阵的第 i 行 + li1 ? 第1行，得到： 消去过程： 第一步：设 ，计算因子 其中 k-1 步消元后 第k步：设 ，计算因子 且计算 共进行 n ? 1步，得到 回代过程： 运算次数： 乘法： 加法： 除法： 二、 选主元消去法 为避免这种情况的发生, 可通过交换方程的次序，选取 绝对值大的元素作主元. 基于这种思想导出了主元素法 在高斯消去法消去过程中可能出现 的情况，这时 高斯消去法将无法进行；即使主因素 但很小， 其作除数 ，也会导致其它元素数量级的严重增长和舍 误差的扩散 例 若p≠k，交换第k个与第p个方程后，再继续消去计算. 这种方法称为列主元Gauss消去法。 列主元Gauss消去法保证了｜lik｜≤1 (i=k+1,k+2,…，n)，提高了数值解的精确度。 列主元消去法 全主元消去法 在第k步消去前， 在系数矩阵右下角的n-k+1阶主子阵中，选绝对值最大的元素作为主元素。 (1) If p ? k then 交换第 k 行与第p行; If q ? k then 交换第 k 列与第 q 列; (2) 消元 注：列交换改变了 xi 的顺序，须记录交换次序， 解完后再换回来。 保证稳定，但费时 算法描述 顺序高斯消去法 function mygauss1(n,a,b,eps) for k=1:n-1 if abs(a(k,k))eps fprintf(a(%d,%d)eps, 算法失败\n,k,k); return; else for i=k+1:n l(i,k)=-a(i,k)/a(k,k); b(i)=b(i)+l(i,k)*b(k); for j=k+1:n a(i,j)=a(i,j