关于对偶复矩阵的笛卡尔分解的酉不变范数不等式
Unitary Invariant Norm Inequalities for Cartesian Decompositions of Dual Complex Matrices
摘要: 研究了对偶复矩阵的笛卡尔分解的酉不变范数不等式。给出了对偶复矩阵的笛卡尔分解的定义。利用对偶向量间的优超关系以及对偶复矩阵的Mirsky定理证明了对偶复矩阵的笛卡尔分解的一个酉不变范数不等式。该不等式揭示了一个对偶复矩阵与其笛卡尔分解中两个对偶Hermite复矩阵的特征值的酉不变范数的数量关系。
Abstract: Unitarily norm inequalities for the Cartesian decomposition of dual complex matrices are studied. The definition of Cartesian decomposition of a dual complex matrix is defined. By using the majorization relation between dual vectors and Mirsky theorem of dual complex matrices, a unitarily norm inequality for the Cartesian decomposition of dual complex matrices is proved, which reveals the quantity relation between the unitarily norm of a dual complex matrix and the two dual Hermitian complex matrices in its Cartesian decomposition.
文章引用:刘霞. 关于对偶复矩阵的笛卡尔分解的酉不变范数不等式[J]. 理论数学, 2026, 16(2): 154-165. https://doi.org/10.12677/pm.2026.162044

1. 引言

1873年,Clifford [1]给出了对偶四元数的概念。如今,对偶四元数代数已成为Clifford代数的一个非常活跃的子代数。作为刚体运动中一个有力且简洁的数学工具,四元数代数在许多应用领域有重要应用,如机器人[2] [3]、空间运动学[4]-[6]、航天学[7]-[9]等。近年来,关于对偶四元数代数在空间运动学和机械系统等领域中的应用得到大量研究,如Yang和Freudenstein [10]介绍了对偶四元数在空间机构中的应用,且Yang [11]将对偶四元数应用于队列机制;Yacob和Semere [12]利用对偶四元数来补偿机械加工工艺变化;文[13] [14]利用对偶四元数标定空间运动学同时在串联机器人和平行机器人中实现;文[15]利用对偶四元数提供了由天线护送的地面车辆组成的编队控制器。

一个对偶四元数矩阵是元素全为对偶四元数的矩阵。目前,对偶矩阵的理论和计算已成功应用于多个领域。由于对偶四元数代数的应用要求建立对偶四元数代数的理论基础,且相比于实矩阵和复矩阵,对偶四元数矩阵的理论基础尚不成熟,所以迫切需要建立对偶四元数矩阵的基础理论和计算方法,如特征值、奇异值、广义逆、范数、方程等。Luo和Qi [16]给出了对偶Hermite四元数矩阵的谱分解和一般对偶四元数矩阵的奇异值分解;Ling,Qi和Yan [17]建立了对偶Hermite四元数矩阵的特征值的极大极小定理;Ding等[18]给出了计算。

对偶四元数矩阵的奇异值分解的一种可行方法;Ding,Li和Wei [19]给出了对偶Hermite四元数矩阵的特征值分解算法;Cui和Qi [20]给出了计算对偶Hermite四元数矩阵的主特征值的幂方法,并将此应用于同时定位和映射问题;Ling,Pan和Qi [21]介绍了对偶四元数矩阵的一种新的度量函数,并引入了两种可实现的近端点算法求解对偶四元数超定方程的逼近解;Ling,He和Qi [22]研究了对偶四元数的一些基本性质,包括极分解定理、奇异值的极小极大定理和Weyl型单调不等式、谱范数和毕达哥拉斯定理,以及最佳低秩逼近;Ling等[23]建立了一般对偶四元数矩阵的Von Neumann型迹不等式和Hoffman-Wielandt型不等式;此外,关于对偶四元数矩阵方程 AXB=C 的解的研究可见文[24]-[26]

笛卡尔分解是实矩阵和复矩阵中一种常见且重要的矩阵分解,但在对偶矩阵和对偶四元数矩阵中较少提及。文[4]介绍了对偶矩阵的笛卡尔分解在空间运动学中的应用。例如,令 A= A s + A ss A 的笛卡尔分解,则在刚体有限位移中, A 的反Hermite部分常用于构造旋转矩阵 Q=1+sinϕ A ss +( 1cosϕ ) A ss 2 = e ϕ A ss 。作为定义对偶复矩阵和对偶复向量的度量性质的重要工具,酉不变范数在最小二乘问题和最佳低秩逼近问题中发挥着重要作用。此外,酉不变范数也常用于分析空间运动学中的刚体运动,如最小范数位移。近年来,对偶四元数矩阵的酉不变范数得到关注,如Cheng和Hu [27]研究了对偶四元数矩阵的酉不变范数,并引入了对偶四元数矩阵的对称规度函数,同时研究了一些特殊的酉不变范数,如Schatten-p范数和Fan-k范数。本文利用对偶向量间的优超关系和Mirsky定理研究涉及对偶复矩阵的笛卡尔分解的酉不变范数不等式,推广了复矩阵的情形(见文[28])。

2. 预备知识

为四元数组成的集合。一个四元数具有形式 q= q 0 + q 1 i+ q 2 j+ q 3 k ,其中 q 0 , q 1 , q 2 , q 3 是实数, i,j,k 是四元数的三个虚数单位,满足 i 2 = j 2 = k 2 =1 ij=ji=k jk=kj=i ki=ik=j 。一个对偶四元数具有形式 q ^ = q 1 + q 2 ϵ ,其中四元数 q 1 q 2 分别称为 q ^ 的可观部和无穷小部, ϵ 是对偶元,满足 ϵ0 ϵ 2 =0 。若 q 1 0 ,则称 q ^ 为可观的;若 q 1 =0 ,则 q ^ 为无穷小的。特别地,当 q 1 q 2 是复数时,称 q ^ 为对偶复数;当 q 1 q 2 是实数时,称 q ^ 为对偶数。注意到一个对偶数 a= a 1 + a 2 ϵ 可逆当且仅当 a 1 0 ,此时, a 1 = a 1 1 a 1 2 a 2 ϵ 。分别记对偶四元数、对偶复数和对偶数的集合为 D D D 。元素分别取自 D D D m×n 的矩阵集合记为 D m×n D m×n D m×n 。文[29]引入了两个对偶数间的全序关系,令 a= a st + a in ϵ,b= b st + b in ϵD ,其中 a st b st 分别是 a b 的可观部, a in b in 分别是 a b 的无穷小部。若 a st > b st ,或 a st = b st a in > b in ,则有 a>b a=b 当且仅当 a st = b st a in = b in 。根据对偶数的全序关系,可定义正对偶数和非负对偶数,即若 a>0 ,则称 a 为正对偶数;若 a0 ,则称 a 为非负对偶数。若 a 为正可观的对偶数,则 a= a st + a in ϵ 的平方根为

.

一个复方阵 A 总可以表示为 A= A+ A * 2 +i A A * 2i ,其中 A * A 的共轭转置矩阵,该分解称为 A 的笛卡尔分解。同样地,一个对偶复方阵 A= A st + A in ϵ 可以表示为

A= A st + A st * 2 +i A st A st * 2i +( A in + A in * 2 +i A in A in * 2i )ϵ =( A st + A st * 2 + A in + A in * 2 ϵ )+i( A st A st * 2i + A in A in * 2i ϵ ) G+iH,

该表示称为对偶复方阵 A 的笛卡尔分解,其中 G H 为对偶Hermite矩阵,即满足 G= G * H= H *

矩阵 UD m×n 称为酉矩阵,如果 U * U=I ,其中 I 是单位矩阵。 D m×n 上的范数 . 称为酉不变的,若 A = UAV 对任意的 AD m×n 和任意的酉矩阵 UD m×m VD n×n 成立。用 Rank( A ) 表示矩阵 A 的秩, σ k ( A ) 表示 A 的第 k 大的奇异值, σ( A )= ( σ 1 ( A ), σ 2 ( A ),, σ n ( A ) ) T 表示 AD n×n 的奇异值向量。 λ k ( A ) 表示对偶Hermite复矩阵 A 的第 k 大的特征值, λ( A )= ( λ 1 ( A ), λ 2 ( A ),, λ n ( A ) ) T 表示 AD n×n 的特征值向量。本文总是将 AD m×n 的奇异值按降序排列为 σ 1 ( A ) σ 2 ( A ) σ min{ m,n } ( A ) diag( d 1 , d 2 ,, d n ) 表示以 d 1 , d 2 ,, d n 为对角元的对角矩阵, tr( A )= i=1 n a ii 表示 A 的迹。

Luo和Qi [16]给出了对偶Hermite四元数矩阵的谱分解,并说明对偶Hermite四元数矩阵的特征值全为对偶数。由此可知,作为特殊的对偶四元数矩阵,对偶复方阵的特征值也全为对偶数。在本文中,对于 n 阶复方阵 A=G+iH ,总是记Hermite矩阵 G H 的特征值为 α j β j 1jn 。下面介绍几种特殊的酉不变范数。对 1p AD m×n 的Schatten-p范数[27]定义为

A ={ ( j=1 min{ m,n } σ j p ( A ) ) 1 p , A st 0 ( j=1 min{ m,n } σ j p ( A in ) ) 1 p , A st =0

AD m×n 的Fan-k范数[27]定义为

A ( k ) = j=1 k σ j ( A ) , 1kmin{ m,n } .

特别地,当 k=1 时为 A 的谱范数,记为 A 2 ;当 k=min{ m,n } 时,为 A 的迹范数,记为 A 1

定义2.1 令 A D n×n 是非负对偶矩阵,若 A 的每行元素之和及每列之和都等于1,则称 A 为双随机对偶矩阵。

A= A st + A in ϵ D n×n 。若 A 是双随机对偶矩阵,则 A st 是双随机矩阵,且 A in 是每行元素及每列元素之和为0的实矩阵。不难看出双随机对偶矩阵关于矩阵乘法是封闭的,即两个同阶双随机对偶矩阵的乘积仍是双随机对偶矩阵。

定义2.2 给定向量 x,y D n 满足 x 1 x 2 x n y 1 y 2 y n 。若对任意的正整数 1kn ,有

i=1 k x i i=1 k y i ,

则称 x y 弱优超,记为 x w y 。若 x w y i=1 n x i = i=1 n y i ,则称 x y 优超,记为 xy

下面给出本文所需的一些引理。

引理2.1 [30] AD m×n p=min{ m,n } 。则对任意的正整数 1kp ,有

σ k ( A )=min{ AX 2 |Rank( X )k1,XD m×n } .

引理2.2 [30] A=diag( d 1 , d 2 ,, d n ) D n×n 。则 A 2 = max 1jn | d j |

引理2.3 [17] 对任意的 AD n×n 1in λ i ( A+ A * 2 ) σ i ( A )

引理2.4 [22] AD m×n ,则

Rank( A )=min{ l:A=BC,BD m×l ,CD l×n } .

引理2.5 [30] A,BD m×n ,则

diag( σ( A ) )diag( σ( B ) ) AB .

下面引理的证明和复矩阵类似,这里略去证明。

引理2.6 设 AD m×n 。则对于任意的正整数 1in ,有

A k =min{ X 1 +k Y 2 :A=X+Y,AD m×n } .

引理2.7 [27] A,BD m×n 。若 A ( k ) B ( k ) 对所有的 1kmin{ m,n } 成立,则 A = B 对所有 D m×n 上的酉不变范数成立。

引理2.8 设 x,y D n ,若 xy ,则存在 n 阶双随机对偶矩阵 A 使得 x=Ay

证: xy ,则对 n 用数学归纳法来证明存在双随机对偶矩阵 A 使得 x=Ay 。当 n=1 时结论平凡地成立。假设结论对 D n1 中的向量成立。令 x= ( x 1 , x 2 ,, x n ) T y= ( y 1 , y 2 ,, y n ) T ,且不失一般性,假设 x 1 x 2 x n y 1 y 2 y n 。首先由 xy 可知 y n x 1 y 1 ,因此存在某个正整数 2kn 使得 y k x 1 y k1 ,从而有 y k x 1 y 1 对某个正整数 2kn 成立。若 y 1 y k 是无穷小的,即 y 1 y k 都是无穷小的,或 y 1 y k 是有相同可观部的可观对偶数,则不难看出存在实数 0t1 使得 x 1 =t y 1 +( 1t ) y k 。若 y 1 y k 是可观的,令 t= y 1 x 1 y 1 y k ,则 tD ,且 x 1 =t y 1 +( 1t ) y k 。总言之,存在对偶数 0t1 使得 x 1 =t y 1 +( 1t ) y k

x ¯ = ( x 2 ,, x n ) T ,并记

y ¯ = ( y 2 ,, y k1 ,( 1t ) y 1 +t y k , y k+1 ,, y n ) T = ( y ¯ 2 ,, y ¯ n ) T .

下面验证 x ¯ y ¯ 。因为 y 1 y k1 x 1 x 2 x n ,所以对于 2mk1 ,有

j=2 m x j j=2 m y j .

对于 kmn

j=2 m x j = j=1 m x j x 1 j=1 m y j x 1 = j=1 m y j t y 1 ( 1t ) y k = j=2 k1 y j +[ ( 1t ) y 1 +t y k ]+ j=k+1 m y j = j=2 m y ¯ j

m=n 时, j=2 m x j = j=2 m y ¯ j 。又因为 xy ,所以有 x ¯ y ¯ 。由归纳假设,存在 n1 阶双随机对偶矩阵 A ¯ 使得 x ¯ = A ¯ y ¯ 。定义 n 阶矩阵 B=( b ij )

b ij ={ t, i=j=1k, 1t, i=1,j=ki=k,j=1, 1, i=j,i1ik, 0, .

则由于 t 是对偶数,不难验证 B 为双随机对偶矩阵。令 A=( 1 0 0 A ¯ )B ,则 A 为双随机对偶矩阵,且 x=Ay

3. 主要结果

本节利用对偶向量间的优超关系,以及Mirsky定理证明关于对偶复矩阵的笛卡尔分解的酉不变范数不等式。

定理3.1 A=diag( d 1 , d 2 ,, d n ) D n×n 满足 | d 1 || d 2 || d n | 。则 σ k ( A )=| d k | 对任意的正整数 1kn 成立。

证:由引理2.3,对任意的正整数 1kn λ k ( A )= d k σ k ( A ) 。同理, λ k ( A )= d k σ k ( A )= σ k ( A ) 。于是有 | d k | σ k ( A )

另一方面,令 D k1 =diag( d 1 , d 2 ,, d k1 ) ,则矩阵 A k1 ( D k1 0 0 0 ) 具有分解

A k1 =( D k1 0 )( I k1 0 ) .

从而由引理2.4可得 Rank( A k1 )k1

再由引理2.1和引理2.2可得

σ k ( A ) AX =| d k | .

综上所述, σ k ( A )=| d k | 对任意的正整数 1kn 成立。

推论3.1 对于任意的Hermite矩阵 AD n×n σ k ( A )=| λ k ( A ) | 对任意的正整数 1kn 成立。

N[ 0 ] 表示所有可观向量的集合。下面给出对偶Hermite复矩阵向量的一个优超关系。

定理3.2 设 A,BD n×n 为Hermite矩阵,则

λ( A )λ( B )λ( AB ) .

证:由引理2.5,Mirsky定理可以改写成

σ( diag( σ 1 ( A ) σ 1 ( B ), σ 2 ( A ) σ 2 ( B ),, σ n ( A ) σ n ( B ) ) ) w σ( AB ) . (1)

再由定理3.1可知(1)式等价于

( | σ 1 ( A ) σ 1 ( B ) |,| σ 2 ( A ) σ 2 ( B ) |,,| σ n ( A ) σ n ( B ) | ) w σ( AB ) . (2)

a max 1in | λ i ( A ) | bmax{ a+ max xN[ 0 ] x * ( AB )x x * x , max 1in | λ i ( B ) | } ,则 B+bI A+aI 均为半正定矩阵,且 A+aI( B+bI ) 也为半正定矩阵。

在(2)式中分别用 B+bI A+aI 代替 B A 可得

( | σ 1 ( A+aI ) σ 1 ( B+aI ) |,,| σ n ( A+aI ) σ n ( B+aI ) | ) w σ( AB+( ab )I ) .

对于半正定对偶复方阵,其奇异值等于特征值,所以由上式有

( | λ 1 ( A+aI ) λ 1 ( B+aI ) |,,| λ n ( A+aI ) λ n ( B+aI ) | ) w λ( AB+( ab )I ) ,

( | λ 1 ( A ) λ 1 ( B )+ab |,,| λ n ( A ) λ n ( B )+ab | ) w λ( AB )+ab .

又因为显然有

( λ 1 ( A ) λ 1 ( B )+ab,, λ n ( A ) λ n ( B )+ab ) w ( | λ 1 ( A ) λ 1 ( B )+ab |,,| λ n ( A ) λ n ( B )+ab | ),

综合可得

( λ 1 ( A ) λ 1 ( B )+ab,, λ n ( A ) λ n ( B )+ab ) w λ( AB )+ab ,

λ( A )λ( B ) w λ( AB ) 。又因为 tr( AB )=tr( A )tr( B ) ,从而有 λ( A )λ( B )λ( AB ) ,证毕。

推论3.2 设 G,HD n×n 为Hermite矩阵,则对任意的正整数 1jn

{ λ j ( G )+ λ nj+1 ( H ) }λ( G+H )λ( G )+λ( H ) . (3)

定理3.3 令 AD n×n A=G+iH A 的笛卡尔分解,则有

{ | α j +i β nj+1 | 2 }{ σ j 2 } ,

{ σ j 2 + σ nj+1 2 2 }{ | α j +i β j | } .

证:在(3)式中分别用 G 2 H 2 代替 A B 可得

{ λ j ( G 2 )+ λ nj+1 ( H 2 ) }λ( G 2 + H 2 )λ( G 2 )+λ( H 2 ) ,

{ α j 2 + β nj+1 2 }λ( G 2 + H 2 ){ α j 2 + β j 2 } . (4)

α j = q 1j + q 2j ϵ, β j = q 3j + q 4j ϵ ,则

| α j +i β j |=| q 1j + q 3j +( q 2j + q 4j )ϵ | ={ q 1j 2 + q 3j 2 + q 1j q 2j + q 3j q 4j q 1j 2 + q 3j 2 ϵ, q 1j +i q 3j 0, q 2j 2 + q 4j 2 ϵ, q 1j +i q 3j =0.

由此可知, | α j +i β nj+1 | 2 = α j 2 + β nj+1 2 | α j +i β j | 2 = α j 2 + β j 2 。于是(4)式可改写成

{ | α j +i β nj+1 | 2 }λ( G 2 + H 2 ){ | α j +i β j | 2 } . (5)

注意到 A * A 是完美Hermite矩阵[1],其特征值要么是正可观对偶数,要么是0,且A的正可观奇异值为 A * A 的正可观特征值的算数平方根,其余奇异值要么是正无穷小对偶数,要么是0,从而有 λ j ( A * A )= σ j 2 ( A ) 。又因为 σ j ( A )= σ j ( A * ) ,于是也有 λ j ( A A * )= σ j 2 ( A )

不难看出, G 2 + H 2 = A * A+A A * 2 。在(4)式中分别用 A * A 2 A A * 2 代替 G H 可得

{ λ j ( A * A 2 )+ λ nj+1 ( A A * 2 ) }λ( G 2 + H 2 ) λ j ( A * A 2 )+ λ j ( A A * 2 ) ,

{ σ j 2 + σ nj+1 2 2 }λ( G 2 + H 2 ){ σ j 2 } . (6)

于是从(5)和(6)可得

{ | α j +i β nj+1 | 2 }{ σ j 2 } .

{ σ j 2 + σ nj+1 2 2 }{ | α j +i β j | 2 } .

下面给出本文的主要结论。

定理3.4 设 A=G+iHD n×n G,H 为Hermite矩阵。设 α j β j 分别是 G H 的特征值,且满足 | α 1 || α 2 || α n | | β 1 || β 2 || β n | ,则

diag( α 1 +i β 1 ,, α n +i β n ) 2 A (7)

对任意的酉不变范数都成立。

证:首先证明不等式(7)对于谱范数成立。由引理2.2可知 G 2 =| α 1 | H 2 =| β 1 | 。此外,由

G= A+ A * 2 ,H= A A * 2i

可得 | α 1 |= G 2 A 2 | β 1 |= H 2 A 2

α 1 = q 1 + q 2 ϵ β 1 = q 3 + q 4 ϵ 。则

| α 1 |={ | q 11 |+sgn( q 1 ) q 2 ϵ, q 11 0, | q 2 |ϵ, q 11 =0. | β 1 |={ | q 3 |+sgn( q 3 ) q 4 ϵ, q 3 0, | q 4 |ϵ, q 4 =0.

其中 sgn( x ) 为符号函数,且

q 1 +i q 3 =0 ,则 | α 1 |=| q 2 |ϵ A 2 | β 1 |=| q 4 |ϵ A 2 。此时, | α 1 +i β 1 |= q 2 2 + q 4 2 ϵ 2 A 2

q 1 +i q 3 0 q 1 q 3 中有一个为0,不妨设 q 1 =0 ,则 | α 1 +i β 1 |=| β 1 | A 2 。若 q 1 q 3 均不为0,设 A 2 = p 1 + p 2 ϵ ,由 | α 1 | A 2 | β 1 | A 2 可得 p 1 max{ | q 1 |,| q 3 | } 。若 | q 1 || q 3 | ,则显然有 q 1 2 + q 3 2 2 p 1 ,从而有 | α 1 +i β 1 | 2 A 2 。若 | q 1 |=| q 3 |< p 1 ,同样显然有 | α 1 +i β 1 | 2 A 2 。若 | q 1 |=| q 3 |= p 1 ,则只需比 | α 1 +i β 1 | 2 A 2 的无穷小部,即

q 1 q 2 + q 3 q 4 q 1 2 + q 3 2 = q 1 q 2 + q 3 q 4 2 | q 1 | = 1 2 sgn( q 1 ) q 2 + 1 2 sgn( q 3 ) q 4 1 2 p 2 + 1 2 p 2 = 2 p 2 .

综上可得, | α 1 +i β 1 | 2 A 2 总是成立的,即不等式(7)对谱范数是成立的。

接下来证明不等式(7)对于迹范数成立。为此,我们先证明这样的事实:对于两个对偶向量 x= ( x 1 , x 2 ,, x n ) T y= ( y 1 , y 2 ,, y n ) T ,其中 x i y i 要么为正可观对偶数,要么为0。若 xy ,则

( y 1 , y 2 ,, y n ) T w ( x 1 , x 2 ,, x n ) T . (8)

事实上,若 xy ,则由定理2.8可知,存在双随机对偶矩阵 A 使得 x=Ay ,即对 1in ,有 x i = j=1 n a ij y j 。由此不难验证 j=1 n a ij y j x i ,即 x i j=1 n a ij y j 。于是有

( x 1 x n )A( y 1 y n )( y 1 y n ) .

所以, ( x 1 , x 2 ,, x n ) T w ( y 1 , y 2 ,, y n ) T ,这就证明了(7)。

由定理3.3有 { σ j 2 + σ nj+1 2 2 }{ | α j +i β j | 2 } 。将(8)应用于 { σ j 2 + σ nj+1 2 2 }{ | α j +i β j | 2 } 可得

{ | α j +i β j | } w { 1 2 σ j 2 + σ nj+1 2 } .

特别地,有 j=1 n | α j +i β j | 1 2 j=1 n σ j 2 + σ nj+1 2 = 2 j=1 n σ j 。这就证明了不等式(7)对迹范数是成立的。

由引理2.7可知,要证明不等式(7),只需证明不等式对Fan-k范数成立即可。由引理2.6,存在对偶复矩阵 X,Y ^ n×n 使得 A=X+Y 使得 A ( k ) = X 1 +k Y 。令

X=M+iN,Y=P+iQ

分别为 X Y 的笛卡尔分解。所以 A=M+P+i( N+Q ) A 的笛卡尔分解。于是

2 A ( k ) = 2 X 1 +k 2 Y diag( α 1 ( M )+i β 1 ( N ),, α n ( M )+i β n ( N ) ) 1 + diag( α 1 ( P )+i β 1 ( Q ),, α n ( P )+i β n ( Q ) ) 2 j=1 k | α j ( M )+i β j ( N ) |+k| α 1 ( P )+i β 1 ( Q ) | = j=1 k | σ j ( M )+i σ j ( N ) |+k| σ 1 ( P )+i σ 1 ( Q ) | = j=1 k { | σ j ( M )+i σ j ( N )|+| σ 1 ( P )+i σ 1 ( Q ) | } j=1 k [ σ j ( M )+ σ 1 ( P ) ]+i[ σ j ( N )+ σ 1 ( Q ) ]

由文[30]中引理8和引理9可知

max 1jn | σ j ( A ) σ j ( B ) | σ 1 ( AB ) .

对任意的 A,BD n×n 成立,即,对任意的 1jn σ j ( A+B ) σ j ( A )+ σ 1 ( B ) .

所以,由前面的验证可得

2 A ( k ) j=1 k | [ σ j ( M )+ σ 1 ( P ) ]+i[ σ j ( N )+ σ 1 ( Q ) ] | j=1 k | σ j ( M+P )+i σ j ( N+Q ) | = diag( α 1 +i β 1 ,, α n +i β n ) ( k ) .

这就完成了不等式(7)的证明。

下面给出例子说明不等式(7)对于两类特殊的酉不变范数成立。

例3.1 令

G=( 0 1+iϵ 12iϵ 1iϵ 0 1iϵ 1+2iϵ 1+iϵ 0 ),H=( 0 1+iϵ 1 1iϵ 0 1iϵ 1 1+iϵ 0 )

A=G+iH=( 0 1+i+( i1 )ϵ 1+i2iϵ 1+i+( 1i )ϵ 0 1+i+( 1i )ϵ 1+i+2iϵ 1+i+( i1 )ϵ 0 ) .

经计算可得 G 的特征值为

α 1 =2, α 2 =11.1547ϵ, α 3 =1+1.1547ϵ .

H 的特征值为

β 1 =2, β 2 =1, β 3 =1 .

于是

| α 1 +i β 1 |=2 2 ,| α 2 +i β 2 |= 2 + 1.1547 2 ϵ,| α 3 +i β 3 |= 2 1.1547 2 ϵ .

A 的奇异值的可观部分别为

( σ 1 ) st = 7.4641 =2.7321, ( σ 2 ) st = 2 =1.4140, ( σ 3 ) st = 0.5359 =0.7321 .

所以,

diag( α 1 +i β 1 , α 2 +i β 2 , α 3 +i β 3 ) 2 =| α 1 +i β 1 |=2 2 =2.8280< 2 ( σ 1 ) st =3.8632 ,

diag( α 1 +i β 1 , α 2 +i β 2 , α 3 +i β 3 ) 2 < 2 A 2

此外,

diag( α 1 +i β 1 , α 2 +i β 2 , α 3 +i β 3 ) 1 =| α 1 +i β 1 |+| α 2 +i β 2 |+| α 3 +i β 3 | =4 2 =5.6560 < 2 ( ( σ 1 ) st + ( σ 2 ) st + ( σ 3 ) st )=6.8978.

diag( α 1 +i β 1 , α 2 +i β 2 , α 3 +i β 3 ) 1 < 2 A 1

一个有趣的问题是不等式(7)中等式何时成立,因此,我们提出如下问题:

问题3.1 刻画不等式(7)中等式成立的条件。

4. 总结与展望

本文围绕对偶复矩阵的笛卡尔分解及酉不变范数不等式展开研究。首先,明确界定了对偶复矩阵的笛卡尔分解定义。其次,基于对偶向量间的优超关系、对偶复矩阵的Mirsky定理,结合谱分解、奇异值分解等已有理论成果,通过推导,证明了一系列关键结论:

(1) 对角对偶矩阵的奇异值等于其对角元的绝对值;

(2) 对偶Hermite复矩阵满足 λ( A )λ( B )λ( AB ) 的优超关系;

(3) 在笛卡尔分解下,存在 { | α j +i β nj+1 | 2 }{ σ j 2 } { σ j 2 + σ nj+1 2 2 }{ | α j +i β j | 2 } 的优超关系;

(4) 核心成果:证明了对任意酉不变范数,

diag( α 1 +i β 1 ,, α n +i β n ) 2 A

恒成立,该不等式清晰揭示了对偶复矩阵与其笛卡尔分解中两个对偶Hermite 复矩阵特征值的酉不变范数数量关联,成功推广了复矩阵的相关结论。

基于上述工作,未来可进一步拓展该不等式的应用场景,如将其应用于空间运动学中的刚体运动分析、最小二乘问题等实际领域;还可探索更高维度对偶矩阵或其他类型对偶代数结构下的笛卡尔分解及范数不等式,完善对偶矩阵理论体系。

参考文献

[1] Clifford, (1871) Preliminary Sketch of Biquaternions. Proceedings of the London Mathematical Society, 1, 381-395. [Google Scholar] [CrossRef
[2] Perez, A. and McCarthy, J.M. (2003) Dual Quaternion Synthesis of Constrained Robotic Systems. Journal of Mechanical Design, 126, 425-435. [Google Scholar] [CrossRef
[3] Valverde, A. and Tsiotras, P. (2018) Spacecraft Robot Kinematics Using Dual Quaternions. Robotics, 7, Article 64. [Google Scholar] [CrossRef
[4] Angeles, J. (1998) The Application of Dual Algebra to Kinematic Analysis. In: Angeles, J. and Zakhariev, E., Eds., Computational Methods in Mechanical Systems, Springer, 3-32. [Google Scholar] [CrossRef
[5] Condurache, D. (2019) A Davenport Dual Angles Approach for Minimal Parameterization of the Rigid Body Displacement and Motion. Mechanism and Machine Theory, 140, 104-122. [Google Scholar] [CrossRef
[6] Leclercq, G., Lefèvre, P. and Blohm, G. (2013) 3D Kinematics Using Dual Quaternions: Theory and Applications in Neuroscience. Frontiers in Behavioral Neuroscience, 7, Article 7. [Google Scholar] [CrossRef] [PubMed]
[7] Abaunza, H., Castillo, P., Lozano, R. and Victorino, A. (2016) Quadrotor Aerial Manipulator Based on Dual Quaternions. 2016 International Conference on Unmanned Aircraft Systems (ICUAS), Arlington, 7-10 June 2016, 152-161. [Google Scholar] [CrossRef
[8] Valverde, A. and Tsiotras, P. (2018) Dual Quaternion Framework for Modeling of Spacecraft-Mounted Multibody Robotic Systems. Frontiers in Robotics and AI, 5, Article 128. [Google Scholar] [CrossRef] [PubMed]
[9] Qiao, B., Tang, S., Ma, K. and Liu, Z. (2013) Relative Position and Attitude Estimation of Spacecrafts Based on Dual Quaternion for Rendezvous and Docking. Acta Astronautica, 91, 237-244. [Google Scholar] [CrossRef
[10] Yang, A.T. and Freudenstein, F. (1964) Application of Dual-Number Quaternion Algebra to the Analysis of Spatial Mechanisms. Journal of Applied Mechanics, 31, 300-308. [Google Scholar] [CrossRef
[11] Yang, A.T. (1967) Application of Dual Quaternions to the Study of Gyrodynamics. Journal of Engineering for Industry, 89, 137-143. [Google Scholar] [CrossRef
[12] Yacob, F. and Semere, D. (2020) Variation Compensation in Machining Processes Using Dual Quaternions. Procedia CIRP, 93, 879-884. [Google Scholar] [CrossRef
[13] Li, G., Zhang, F., Fu, Y. and Wang, S. (2019) Kinematic Calibration of Serial Robot Using Dual Quaternions. Industrial Robot: The International Journal of Robotics Research and Application, 46, 247-258. [Google Scholar] [CrossRef
[14] Luo, J., Chen, S., Zhang, C., Chen, C. and Yang, G. (2023) Efficient Kinematic Calibration for Articulated Robot Based on Unit Dual Quaternion. IEEE Transactions on Industrial Informatics, 19, 11898-11909. [Google Scholar] [CrossRef
[15] Luo, J., Chen, S., Jiang, D., Zheng, T., Li, H., Fang, Z., et al. (2024) Efficient Kinematic Calibration for Parallel Manipulators Based on Unit Dual Quaternion. IEEE Transactions on Industrial Informatics, 20, 6791-6801. [Google Scholar] [CrossRef
[16] Qi, L.Q. and Luo, Z.Y. (2023) Eigenvalues and Singular Values of Dual Quaternion Matrices. Pacific Journal of Optimization, 19, 257-272.
[17] Ling, C., Qi, L. and Yan, H. (2023) Minimax Principle for Eigenvalues of Dual Quaternion Hermitian Matrices and Generalized Inverses of Dual Quaternion Matrices. Numerical Functional Analysis and Optimization, 44, 1371-1394. [Google Scholar] [CrossRef
[18] Ding, W., Li, Y., Wang, T. and Wei, M. (2024) Dual Quaternion Singular Value Decomposition Based on Bidiagonalization to a Dual Number Matrix Using Dual Quaternion Householder Transformations. Applied Mathematics Letters, 152, Article ID: 109021. [Google Scholar] [CrossRef
[19] Ding, W., Li, Y. and Wei, M. (2025) A New Structure-Preserving Method for Dual Quaternion Hermitian Eigenvalue Problems. Applied Mathematics Letters, 163, Article ID: 109432. [Google Scholar] [CrossRef
[20] Cui, C. and Qi, L. (2024) A Power Method for Computing the Dominant Eigenvalue of a Dual Quaternion Hermitian Matrix. Journal of Scientific Computing, 100, Article No. 21. [Google Scholar] [CrossRef
[21] Ling, C., Pan, C. and Qi, L. (2025) A Metric Function for Dual Quaternion Matrices and Related Least-Squares Problems. Linear and Multilinear Algebra, 73, 3056-3079. [Google Scholar] [CrossRef
[22] Ling, C., He, H. and Qi, L. (2022) Singular Values of Dual Quaternion Matrices and Their Low-Rank Approximations. Numerical Functional Analysis and Optimization, 43, 1423-1458. [Google Scholar] [CrossRef
[23] Chen, L., Hongjin, H., Liqun, Q. and Tingting, F. (2024) Spectral Norm and Von Neumann Type Trace Inequality for Dual Quaternion Matrices. Pacific Journal of Optimization, 2024, 229-247. [Google Scholar] [CrossRef
[24] Zhu, L., Wang, Q. and Kou, Z. (2025) The Least-Norm Solution to a Matrix Equation over the Dual Quaterion Algebra. Symmetry, 17, Article 267. [Google Scholar] [CrossRef
[25] Shi, L., Wang, Q., Xie, L. and Zhang, X. (2024) Solving the Dual Generalized Commutative Quaternion Matrix Equation AXB = C. Symmetry, 16, Article 1359. [Google Scholar] [CrossRef
[26] Chen, Y., Wang, Q. and Xie, L. (2024) Dual Quaternion Matrix Equation AXB = C with Applications. Symmetry, 16, 287. [Google Scholar] [CrossRef
[27] Cheng, S. and Hu, H.F. (2024) Unitarily Invariant Norms on Dual Quaternion Matrices. Pacific Journal of Optimization, 2024, 371-387. [Google Scholar] [CrossRef
[28] Zhan, X. (1999) Norm Inequalities for Cartesian Decompositions. Linear Algebra and its Applications, 286, 297-301. [Google Scholar] [CrossRef
[29] Qi, L., Ling, C. and Yan, H. (2022) Dual Quaternions and Dual Quaternion Vectors. Communications on Applied Mathematics and Computation, 4, 1494-1508. [Google Scholar] [CrossRef
[30] Zhong, J. and Zhong, P. (2025) A Mirsky-Type Unitarily Invariant Norm Inequality for Dual Quaternion Matrices and Its Applications. Symmetry, 17, Article 1355. [Google Scholar] [CrossRef