四元数矩阵方程axah+基本知识的最小二乘问题.docxVIP

四元数矩阵方程axah+基本知识的最小二乘问题 0 元数矩阵的解 哈密顿的四位数是1843年发现的数学概念。四位数的乘法不符合交换法，因此它们似乎破坏了科学知识的最基本原则，并且长期未被纳入。最近，人们将恢复平面推广到四维空间，并在多个领域应用了四元数理论。 四元数乘法的非可交换性,造成在已有的数域上矩阵方程的一些解法难以推广到四元数体上.文献中用四元数矩阵分解的方法给出了矩阵方程有解和约束解的条件及解的表达式.在文献中,借助于四元数矩阵的实表示,解决了四元数矩阵方程AXB-CYD=E的求解,并得到了一些四元数体上的结论.在文献中已经得到实数域上的矩阵方程AXAT+BYBT=C的对称与反对称最小范数最小二乘解. 本文借助于四元数矩阵的实表示,将四元数矩阵方程等价地表示为实数矩阵方程,并利用实数域上已有的结论来解决下面问题: 给定A∈Qm×n,B∈Qm×l,C∈Qm×m,求～X∈SCn(Q)?～Y∈SCn(Q)X～∈SCn(Q)?Y～∈SCn(Q)使得 ∥A～XAΗ+B～YBΗ-C∥F=minX∈SCn(Q)Y∈SCL(Q)∥AXAΗ+BYBΗ-C∥F.∥AX～AH+BY～BH?C∥F=minX∈SCn(Q)Y∈SCL(Q)∥AXAH+BYBH?C∥F. SCn(Q)为Q上的全体n×n自共轭矩阵. 1 元数矩阵的定义 定义1 设q=a+bi+cj+dk,a、b、c、d∈R,其中i、j、k满足i2=j2=k2=-1,ij=-ji=k,jk=-kj=i,ki=-ik=j,则称形式如q=a+bi+cj+dk的数为四元数. 称a为q的实部,记为Re(q)=a,称bi+cj+dk为q的虚部,记为Im(q)=bi+cj+dk. 四元数的全体记为Q,即Q={a+bi+cj+dk|a,b,c,d∈R}.容易验证四元数的加法满足结合律和交换律,乘法只满足结合律,无交换律. 定义2 设矩阵A=(aij)m×n,aij∈Q,则称A为四元数矩阵,m×n矩阵的全体记为Qm×n.单位矩阵In=diag(1,…,1).四元数矩阵的基本运算(加法,乘法,数乘)的定义与常规矩阵一样. 定义3 对任意四元数矩阵A=A1+A2i+A3j+A4k∈Qm×n,其中Ai∈Rm×n,i=1,2,3,4.有 Aσ=(A1-A2-A3-A4A2A1-A4A3A3A4A1-A2A4-A3A2A1)∈R4m×4n, 实矩阵Aσ称为四元数矩阵A的实表示. 2 q-现代类,有3个主要的性质 性质1 设A,B∈Qm×n,C∈Qn×s,α∈R,则 (A+B)σ=Aσ+Bσ, (1) (αΑ)σ=αAσ, (2) (AC)σ=AσCσ, (3) (AH)σ=(Aσ)T, (4) ‖Aσ‖F=2‖A‖F. (5) 性质2 设 Ρt=(Ιt-ΙtΙt-Ιt)?Qt=(-ΙtΙtΙt-Ιt)?Rt=(-Ιt-ΙtΙtΙt),St=(-ΙtΙt-ΙtΙt), 则Q-1mAσQm=Aσ,R-1mAσRm=Aσ,S-1mAσSm=Aσ. 3 rsa-r-sw13w3waxa-1823w33 对于给定A∈Rm×n,B∈Rm×l,C∈Rm×m,先作如下假设: 假设1CT=C(对于一般的矩阵C可以作变换C′=12(CΤ+C),则矩阵C′是对称矩阵,代入即可). 假设2na=rank(A)和nb=rank(B),且设矩阵A和B的简化奇异值分解分别为 A=UADAVTA,B=UBDBVTB, 其中:DA和DB分别为na和nb阶满秩对角方阵;UA和VA均为列正交矩阵;UB和VB也均为列正交矩阵. 假设3 (UTAUB)的奇异值分解为 (UTAUB)=UABDABVTAB, 其中UAB与VAB为正交矩阵;DAB=diag(Ir,Σ,0)∈Rna×nb为块对角矩阵,Ir为r阶单位矩阵,Σ=diag(σ1,…,σs);且正整数r和s满足r+s=rank(UTAUB). 假设4 UΤABUΤACUAUAB=(W11W12W13WΤ12W22W23WΤ13WΤ23W33)rsna-r-s,UΤABUΤBCUBVAB=(Ζ11Ζ12Ζ13ΖΤ12Ζ22Ζ23ΖΤ13ΖΤ23Ζ33)rsnb-r-s. 假设5 ～X=(R?×(W12-Ζ12∑)W13(?×(W12-Ζ12∑))Τψ×(W22-∑Ζ22∑)W23WΤ13WΤ23W33),～Y=(W11-R?×(Ζ12-W12∑)Ζ13(?×(Ζ12-W12∑))Τψ×(Ζ22-∑W22∑)Ζ23ΖΤ13ΖΤ23Ζ33). 其中R∈Sr(R)为任意矩阵, ?(φij)r×s=(11-σ2j)r×s?ψ=(ψij)s×s=(11-σiσ2j)s×s. 引理1给定A∈Rm×n,B∈Rm×l,C∈Rm×m,则存在～X∈Sn(R)?～Y