1.6 Kronecker 积
定义1.19 设 A∈Rm×n,B∈Rp×q , 则 A 与 B 的 Kronecker 积定义为
A⊗B=a11Ba21B⋮am1Ba12Ba22B⋮am2B……⋱…a1nBa2nB⋮amnB∈Rmp×nq. Kronecker积也称为直积(direct product),或张量积(tensor product).
任意两个矩阵都存在 Kronecker 积, 且 A⊗B 和 B⊗A 是同阶矩阵. 通常 A⊗B=B⊗A , 但它们之间存在下面的关系式.
定理1.61 设 A∈Cm×m , B∈Cn×n , 则存在 mn 阶置换矩阵 P 使得
PT(A⊗B)P=B⊗A. 定理1.62 矩阵的Kronecker积有以下性质:
=(A1B1)⊗(A2B2)⊗⋯⊗(AkBk); (1) (αA)⊗B=A⊗(αB)=α(A⊗B),∀α∈R;
(2) (A⊗B)T=AT⊗BT, (A⊗B)∗=A∗⊗B∗;
(3) (A⊗B)⊗C=A⊗(B⊗C);
(4) (A+B)⊗C=A⊗C+B⊗C;
(5) A⊗(B+C)=A⊗B+A⊗C;
(6) 混合积: (A⊗B)(C⊗D)=(AC)⊗(BD)
(7) (A1⊗A2⊗⋯⊗Ak)(B1⊗B2⊗⋯⊗Bk)
(8) (A1⊗B1)(A2⊗B2)…(Ak⊗Bk)=(A1A2…Ak)⊗(B1B2…Bk)
(9) rank(A⊗B)=rank(A)rank(B) ;
推论1.63 设 A=Q1Λ1Q1−1∈Cm×m,B=Q2Λ2Q2−1∈Cn×n ,则
A⊗B=(Q1⊗Q2)(Λ1⊗Λ2)(Q1⊗Q2)−1. Kronecker积的特征值
定理1.64 设 A∈Rm×m,B∈Rn×n , 并设 (λ,x) 和 (μ,y) 分别是 A 和 B 的一个特征对, 则 (λμ,x⊗y) 是 A⊗B 的一个特征对. 由此可知, B⊗A 与 A⊗B 具有相同的特征值.
推论1.65设 A∈Cm×m,B∈Cn×n ,则 A⊗In+Im⊗B 的特征值为 λi+μj ,其中 λi 和 μj 分别为 A 和 B 的特征值.
定理1.66 设 A∈Rm×m , B∈Rn×n , 则
(1) tr(A⊗B)=tr(A)tr(B) ;
(2) det(A⊗B)=det(A)ndet(B)m ;
(3) 若 A 和 B 都非奇异, 则 (A⊗B)−1=A−1⊗B−1 .
Kronecker积与向量的乘积
定理1.67 设矩阵 X=[x1,x2,…,xn]∈Rm×n , 记 vec(X) 为 X 按列拉成的 mn 维列向量, 即
vec(X)=[x1⊤,x2⊤,…,xn⊤]⊤, 则有
vec(AX)=(I⊗A)vec(X),vec(XB)=(B⊤⊗I)vec(X), 以及
(A⊗B)vec(X)=vec(BXA⊤). 我们称 vec(X) 为向量化算子或拉直算子. 该结论对节省运算量和存储量都有好处.
Kronecker积与矩阵方程
Kronecker 积一个重要应用是可以将某些矩阵方程转化成一般的代数方程.
定理1.68 矩阵方程
AX+XB=D 等价于代数方程
(I⊗A+BT⊗I)vec(X)=vec(D).