四元数矩阵方程的复转化及保结构算法.pdfVIP

维普资讯 2008年6月 纯粹数学与应用数学 Jun．2008 第24卷 第2期 PureandApplied Mathematics V0l1．24No．2 四元数矩阵方程的复转化及保结构算法 黄敬频 ，于 艳 (广西民族大学 数学与计算机科学学院 南宁530006) 摘 要：给出四元数矩阵复表示运算定义及其相关性质，并运用复表示运算的保结构 特性，讨论了四元数矩阵Moore-Penrose逆计算以及两类四元数矩阵方程 AXB=C 和A —XB=C 的数值求解方法．数值算例检验 了所给算法的可行性． 关 键 词：四元数矩阵；矩阵方程；复表示；保结构算法；Moore-Penrose逆 中图分类号：0151．21 文献标识码：A 文章编号：1008—5513(2008)02—0321—06 1 引 言 20世纪 80年代 以来，随着刚体动力学的发展，人们发现利用 四元数和四元数矩阵可以较 好地处理刚体运动 问题，从而使 四元数方法在理论力学中开始获得应用．现在人们发现四元数 和 四元数矩阵的应用范围越来越广泛．比如，在计算机图形学中，对图形旋转和球面线性插值可 用 四元数来构造强制变换；在航天领域，可用 四元数来描述并处理航天器的姿态控制 问题；在雷 达和无线通信中，可用信道传输的四元数模型处理用信息传输时延和方向角进行定位 问题 [1-3]． 关于 四元数矩 阵的理论研究方面，近三十年来，取得 了许多重要成果 4【-6J．然而，目前关于 四元数矩阵数值计算方面的成果并不多，主要原因在于四元数体上元素乘法的非可换性，给相 关问题的计算机处理带来困难．但随着 四元数矩阵应用范围的不断扩大，对于它的数值计算 问 题是不可回避的．因此，研究四元数矩阵问题的数值算法具有较大的实际意义． 本文主要利用 四元数矩阵复表示运算及其性质，讨论 四元数矩阵 M—P逆 的计算、四元数 矩阵方程AXB=C和 A — B= 的复转化、可解性判别与保结构数值求解方法． 2 四元数矩阵复表示运算及性质 我们分别用 ，C，Q表示实数域、复数域和 四元数体； ，Cm ，Qm×n分别表示全 体 仇 Xn实矩阵、复矩阵和四元数矩阵的集合； ，AH分别表示矩阵 的共轭和共轭转置矩 阵；p(A)表示复方阵A的谱半径；lIAIIr表示 仇Xn四元数矩阵 的Frobenius范数．四元数 矩阵的复表示运算定义如下． 定义 1 设A=A1+A2j∈Q (A1，A2∈C )是 在复数域 C上的分解式，则称 r LA()=J【 z]∈c2mx2n 2 1 j 为四元数矩阵 的复表示矩阵．这里，我们称 为四元数矩阵的复表示运算 下面给出四元数矩阵复表示运算 的一些性质． 收稿日期：2006-11—08． 基金项 目：广西自然科学基金项目(2008J32)，广西教育厅科研基金项 目(2006J26) 作者简介：黄敬频 (1964一)，教授．研究方向：矩阵计算及其应用． 维普资讯 322 纯粹数学与应用数学 第 24卷 性质 14{71设 Q上的矩阵 ， 对于下列运算有意义，则 (a)L(A)=L(B)甘A=B，L(A-4-B)=L(A)-4- ()，L(AB)= () ()； (b)A∈Qn×非奇异甘 ()非奇异，JtL(A )=(()) ； (C)，rank(L(A))=2rank(A)，即L(A)的秩为偶数； (d)L(A爿)= (())爿；A是四元数酉阵 ()是复酉阵；