正余弦矩阵函数的精细积分算法.pdfVIP

第17卷第4期 广西工学院学报 V01．17No．4 0FGUANGXIUNIVERSITYOF Dec．2006 2006年12月 JOURNAL TECHNOLOGY 文章编号l004—64lo(2006)04一0089一03 正余弦矩阵函数的精细积分算法 陆 静1’2，韦笑梅1，齐辉2 (1．广西工学院汽车工程系，广西柳州545006；2．哈尔滨工程大学建筑工程学院，黑龙江哈尔滨150001) 摘要：利用正、余弦函数的倍角公式，提出了一种原理简单、实施容易的矩阵正、余弦函数的精细积分算法。该方法 不需要求矩阵的特征值和特征向量，避免了级数展开解法计算精度不高，效率较低的缺点，具有较高的计算精度和 效率，并用数值算例表明了该方法的有效性。 关键词：精细积分i正弦矩阵函数；余弦矩阵函数 中图分类号：0151．21 文献标识码：A O 引言 正、余弦矩阵函数是一种常见的矩阵函数，通常采用Jordan标准形以及级数展开等方法来求解。采用 值及其相应的特征向量，即使矩阵形式比较简单，这一工作也是相当费时的。求大型、复杂的矩阵的特征值和 特征向量具有一定的困难，只能借助数值计算的方法来求出它的近似解。所以，Jordan标准形法不适合解大 型矩阵的矩阵函数。级数展开的方法虽然不需要求特征值和特征向量，但存在截断误差，计算的精度不容易 保证，有时还可能不收敛，计算效率也不高。因此，探索一种用于求解大型的正、余弦矩阵函数的高精度、高效 率的算法就显得尤为重要。钟万勰先生在上世纪90年代提出的精细积分算法有效地解决了指数矩阵函数的 计算问题，得到了广泛的应用[1~3]。精细积分的原理非常简单，当矩阵和积分步长满足一定条件时，其计算结 果可以达到计算机所能表示的满精度，计算效率也很高。本文利用精细积分算法实现了正、余弦函数的数值 求解，数值算例表明该算法得出的结果具有很高的效率和精度。 1正余弦矩阵函数的精细积分算法 设A为挖阶方阵，取m=2Ⅳ，则正、余弦矩阵函数可化为 sinA—Sin(2Ⅳ曰) (1) cosA=cos(2Ⅳ8) 式中曰=A／川．利用正弦函数和余弦函数的倍角公式，可以得到： sin21艿一2sin口cos曰；cos21曰==2cos2曰一j sin(2‘B)=sin(2×2‘一1B)=2sin(2‘一1B)coS(2‘一1B) (2) cos(2‘召)=cos(2×2i一1B)一2cos2(2卜1口)一， cosA—cos(2ⅣB)=2cos2(2。Ⅳ一1曰)一， 收稿日期：2006—10一11 作者简介：陆静(1973一)，女，广西武鸣人，广西工学院汽车工程系讲师，哈尔滨工程大学在读研究生。 万方数据 90 广西工学院学报 第17卷 3项已可以保证精度了。以B矩阵为自变量的正弦、余弦函数的泰勒展开式为 s瑚：曰一等+等：s㈨ J： 。5 (3) c。s8：，一等+筹：J+c∞) 在计算机编程计算时，若将它们直接与单位矩阵相加，就会成为尾数，在计算机的舍