线性代数方程组与矩阵.ppt

一般带型方程组 ： A x = d 系数矩阵 A 是带状的，即 对 j i+q 有 aij =0, 则称 A =(aij) 具有上带宽 q； 对 i j+p 有 aij =0, 则称 A =(aij) 具有下带宽 p。 2. 带型方程组的解法 注：许多应用中只要方程排序就能使每个未知数 xi 只出现在 与第 i 个方程相邻的几个方程中，即构成带状方程组。 带型方程组的解法 当带型矩阵中绝大部分为零元素时，为节省存储空间，可以采用压缩存储的方法。 带型矩阵的压缩 其中 A 为对称带型矩阵，h 称为半带宽， 2h+1 称为带宽。 为了进行压缩存储，可用 n 行、2h+1 列的二维数组 B (1: n,1: 2h+1 ) 来存放 A 中带区内的非零元素。 其存放原则如下： （1）带型矩阵 A 中的行与二维数组 B 中的行一一对应； （2）带型矩阵 A 中每一行上带区内的非零元素以左边对齐 顺序存放在二维数组 B 中相应的行中，而对于前行与 最后行中最右边的空余部分均填入0； 特点： （1）A 中的前 h+1 行对角线元素在 B 中的对角线上； （2）A 中的 h+1 行后对角线元素在 B 中的第 h+1列上。 带型矩阵的压缩 带型矩阵的压缩 例2.3 带型矩阵的压缩 解：此矩阵的半带宽 h = 2，带宽 2h+1 = 5 ? 可用一个 9?5 的二维数组来表示 例2.3 带型矩阵的压缩 在压缩存储下可用列选主元的同步高斯消去法 求解带型方程组 Ax = d （ 算法2.6，A ? B ）： (1) for k = 1 : n-1，同步归一化和消元高斯变换 列选主元（行交换，始终在第一列中选） b(k,1) 系数矩阵归一化 常数向量归一化 系数矩阵消元 压缩存储下带型方程组的求解 最后一列置零 0? b(i, 2h+1) 常数向量消元 (2) 进行回代（求上三角方程组） 解出 xn 回代逐个解出 xn-1,?, x2, x1，即, 压缩存储下带型方程组的求解 算法2.6的C程序 及应用例子 带状高斯消去法： 设 A∈ Rnxn的上带宽q ，下带宽p ， 则当 A 存在LU分解时， U的上带宽为q，L的下带宽为P。 算法中 i j 的 A(i, j) 被 L(i, j) 覆盖， 否则被 U(i, j) 覆盖。 带型方程组的LU分解求法 下带宽 i ? j +p 上带宽 j ? i + q 含U(i,j) 带型方程组的LU分解求法 求解带状下三角方程组：（带向前消去法） 令 L∈ Rnxn是下带宽 p 的单位下三角阵， 算法中以 Lx＝b 的解覆盖 b ∈ Rn x(i) 带型方程组的LU分解求法 下带宽 i ? j +p 求解上三角带状方程组：（带向后消去法） 令 U∈ Rnxn是上带宽q的非奇异上三角阵， 算法中以 Ux＝b 之解覆盖 b∈ Rn x(i) 一般带型方程组的求解 上带宽 j ? i + q 设一个三对角对称正定阵 A∈Rnxn 和 b∈Rn，（即p=1,q=1) 则分解 A=LDLT ，其中 L 为 p=1 的单位下三角矩阵， 算法中计算 Ax＝b 的解覆盖 b。 设 A 的对角线元素存在 d(1:n) 中， 下次对角线元素存在 e(1:n-1) 中， 三对角对称正定阵的求解 x(i) L(ei) D(di) y(i) 算法2.10的C程序 及应用例子 Please give me a moment!