第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