矩阵Doolittle分解的快速算法.pdfVIP

第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