工程数学(09用矩阵分解法解线性代数方程组.pptVIP

* 工程数学 工程数学 * 第三节 用矩阵分解法求解线性方程组 七、 三对角方程组的解法 lupqdsv.m %功能：调用全主元三角分解函数[LU,p,q]=lupqd(A) % 求解线性方程组Ax=b。 %解法：PAQ-1=LU, Ax=b←→(PAQ-1)(Qx)=Pb % LU(Qx)=Pb, z=Qx, y=Uz % Ly=f=Pb, f(i)=b(p(i)) % Uz=y, z=Qx , x(q(i))=z(i). %输入：方阵A，右端项b（行或列向量均可） %输出：解x（行向量） 定义1 若n 阶矩阵A=(aij)的元素满足:对于1p,qn的正整数p、q,有j≥i+p及i≥j+q时，aij=0，则A称为带状矩阵. 带宽为w=p+q-1。 A称为三对 角矩阵。 较常见带状矩阵为带宽为3（p=q=2,w=3)的矩阵。 系数矩阵为三对角矩阵的线性方程组称为三对角方程组。 七、 三对角方程组的解法 三对角线性方程组 应用追赶法求解三对角线性方程组。追赶法仍然 保持LU分解特性,它是一种特殊的LU分解。充分利用 了系数矩阵的特点，而且使之分解更简单,得到对三对 角线性方程组的快速解法。 定理 如果带宽为 w=p+q-1 的n阶带状矩阵A有LU 分解:A=LU，则L是带宽为p