c语言doolittle分解矩阵的方法,矩阵Doolittle分解的快速算法.pdf

第38卷第1期 东华理工大学学报(自然科学版) V01．38No．1

EAST

2015年3月 JOURNALOF CHINAINSTITUTEOFTECHNOLOGYMar．2015

doi：10．3969／j．issn．1674-3504．2015．01．018

矩阵Doolittle分解的快速算法

吴光文， 黄乡生， 胡文龙

(东华理工大学机械与电子工程学院，江西抚州344000)

摘要：选列主元的Doohttle分解是一种稳定的矩阵三角分解算法。在使用计算机运算时，能将运算过程中产生的机器数

截断误差缩小。但是在进行选列主元的过程中，包含大量的数据交换操作，降低了运算的效率。本文提出一种使用交换

指针代替交换整行元素的方法来提高运算的效率，可以将交换数据的工作量减少到原来的1／n(n为矩阵的阶数)。通过

实验程序运行的结果证明，运算结果正确，算法执行过程稳定、可靠。

关键词：Doolittle分解；三角分解；截断误差

中图分类号：TP311．1文献标识码：A 文章编号：1674-3504(2015)01-0116-04

吴光文，黄乡生，胡文龙．2015．矩阵Doolittle分解的快速算法[J]．东华理工大学学报：自然科学版，38(1)：116．119．

Wu novelfast formatrixdoolittle of

Guang—wen，HungXiang-sheng，Huwen—long．2015．Aalgori出m decomposition[J]．Journal

EastChinaInstituteofTechnology(NaturalScience)，38(1)：116-119．

三角矩阵的特殊形状决定了它们有较好的运 者导致计算机进入除法溢出中断程序，均不能得到

算性质，这些性质使得它们在矩阵的行列式求解、 正确的结果。为了改善这种情况，实际操作中采用

矩阵逆的运算和求线性方程组的解等问题中发挥 选列主元或者全选主元的方法，通过变换主元所在

重要的作用。使得矩阵在阶数比较大时复杂的求 的行和列后再进行矩阵的分解运算。当矩阵的阶

解过程变得可以实现，且能提高效率。很多学者针 数比较大时，大量的行列数据交换运算会大大降低

对自己的应用领域对矩阵的分解做了大量的研究， 算法效率。

陈建平等(2001)对通用的基本线性代数子程序库 本文提出一种新的列选主元的Doolittle分解

linear

(basic 算法，达到不必移动数据即可进行主元行列交换的

algebrasubprograms，BIAS)的发展过

程进行了介绍。Chen(2003)对m×n矩阵的三角

分解进行了研究。Liu(2003)给出了Crout方法在阵的阶数。

线性方程组中的应用。智东杰等(2007)给出了

1 Doolittle分解的原理及实现方法

Doolittle递归分解算法的程序设计方法。鲁铁定等

(2009)探讨了采用QR分解进行误差方式式进行定理 设矩阵A∈C““，若其中的rt一1阶顺

直接解算的算法。李辉等(2013)探讨了矩阵的奇

序子式均不为0，则A存在三角分解A=LU，其中L

异值分解的核磁共振T2谱反演问题。 为下三角矩阵，U为上三角矩阵。此外，如果A非奇

理论上可以通过Doolittle分解将C“4矩阵分异，则分解方法唯一(朱元国等，2010)。A为待分解

解为对角元素为1的下三角矩阵￡和上三角矩阵 的n×n阶矩

的u的乘积，但是在实际计算过程中，计算机会将 all