chp3 线性方程组的求解.ppt

3.7带状方程组的三角分解法 基本原理 在很多实际问题中,例如:在精馏塔计算、传热微分方程数值解等问题中都会遇到求解大型带状线性方程组,其系数矩阵的特点是含有大量的零元素,且其非零元素均分布在对角线元素的两侧成带状,对于这类特殊线性代数方程组,可采用改进三角分解法，节约内存，从而有效地实现其求解运算。 设 , 对于小于n的正整数p和q,当 或 时, ,则称A是具有上带宽 q 和下带宽 p 的带状矩阵,以带状矩阵为系数矩阵的方程组AX=t,称为带状方程组.特别,当p=q=s时,称s为A的半带宽, 2s+1称为A的总带宽. 基本原理 一般,n较大,且p,qn,带状矩阵的LU分解只需计算非零元素,且可证明L和U仍保持A的带状结构,即有 A= = =LU A的LU分解可以像3.6节的公式一样计算，但不需要计算L和U中的零元,即对L中的第k列,只需计算第k+1至k+p行的元，对U的第k行只需计算