当前位置:   article > 正文

【矩阵论】4. 矩阵运算——张量积

张量积

矩阵论的所有文章,主要内容参考北航赵迪老师的课件

[注]由于矩阵论对计算机比较重要,所以选修了这门课,但不是专业搞数学的,所以存在很多口语化描述,而且对很多东西理解不是很正确与透彻,欢迎大家指正。我可能间歇性忙,但有空一定会回复修改的。

矩阵论
1. 准备知识——复数域上矩阵,Hermite变换
1.准备知识——复数域上的内积域正交阵
1.准备知识——Hermite阵,二次型,矩阵合同,正定阵,幂0阵,幂等阵,矩阵的秩
2. 矩阵分解——SVD准备知识——奇异值
2. 矩阵分解——SVD
2. 矩阵分解——QR分解
2. 矩阵分解——正定阵分解
2. 矩阵分解——单阵谱分解
2. 矩阵分解——正规分解——正规阵
2. 矩阵分解——正规谱分解
2. 矩阵分解——高低分解
3. 矩阵函数——常见解析函数
3. 矩阵函数——谱公式,幂0与泰勒计算矩阵函数
3. 矩阵函数——矩阵函数求导
4. 矩阵运算——观察法求矩阵特征值特征向量
4. 矩阵运算——张量积
4. 矩阵运算——矩阵拉直
4.矩阵运算——广义逆——加号逆定义性质与特殊矩阵的加号逆
4. 矩阵运算——广义逆——加号逆的计算
4. 矩阵运算——广义逆——加号逆应用
4. 矩阵运算——广义逆——减号逆
5. 线性空间与线性变换——线性空间
5. 线性空间与线性变换——生成子空间
5. 线性空间与线性变换——线性映射与自然基分解,线性变换
6. 正规方程与矩阵方程求解
7. 范数理论——基本概念——向量范数与矩阵范数
7.范数理论——基本概念——矩阵范数生成向量范数&谱范不等式
7. 矩阵理论——算子范数
7.范数理论——范数估计——许尔估计&谱估计
7. 范数理论——非负/正矩阵
8. 常用矩阵总结——秩1矩阵,优阵(单位正交阵),Hermite阵
8. 常用矩阵总结——镜面阵,正定阵
8. 常用矩阵总结——单阵,正规阵,幂0阵,幂等阵,循环阵


在这里插入图片描述

4.2 张量积

4.2.1 定义

A = ( a i j ) m × n A=(a_{ij})_{m\times n} A=(aij)m×n B = ( b i j ) p × q B=(b_{ij})_{p\times q} B=(bij)p×q ,则称分块矩阵 ( a 11 B a 12 B ⋯ a 1 n B a 21 B a 22 B ⋯ a 2 n B ⋮ ⋮ ⋱ ⋮ a n 1 B a n 2 ⋯ a n n ) \left(

a11Ba12Ba1nBa21Ba22Ba2nBan1Ban2ann
\right) a11Ba21Ban1Ba12Ba22Ban2a1nBa2nBann A A A B B B 的张量积,记作 A ⊗ B = ( a i j B ) m p × n q A\otimes B=(a_{ij}B)_{mp\times nq} AB=(aijB)mp×nq

eg
A = ( a b c d ) , B = ( 2 3 ) A=\left(

abcd
\right),B=\left(
23
\right) A=(acbd),B=(23)

A ⊗ B = ( a B b B c B d B ) = ( 2 a 2 b 3 a 3 b 2 c 2 d 3 c 3 d ) , B ⊗ A = ( 2 A 3 B ) = ( 2 a 2 b 2 c 2 d 3 a 3 b 3 c 3 d )

AB=(aBbBcBdB)=(2a2b3a3b2c2d3c3d),BA=(2A3B)=(2a2b2c2d3a3b3c3d)
AB=(aBcBbBdB)= 2a3a2c3c2b3b2d3d ,BA=(2A3B)= 2a2c3a3c2b2d3b3d

  • 张量积不满足交换律 ,即 A ⊗ B ≠ B ⊗ A A\otimes B\neq B\otimes A AB=BA
定理
  1. 两个上三角的张量积也是上三角

  2. 两个对角阵的张量积是对角阵
    P − 1 A P = A 1 = ( λ 1 0 ⋱ 0 λ m ) , Q − 1 B Q = ( t 1 0 ⋱ 0 t n ) ( P ⊗ Q ) − 1 ( A ⊗ B ) ( P ⊗ Q ) = 吸收公式 ( P − 1 A P ) ⊗ ( Q − 1 B Q ) = A 1 ⊗ B 1 = ( λ 1 B 1 0 ⋱ 0 λ m B n )

    P1AP=A1=(λ100λm),Q1BQ=(t100tn)(PQ)1(AB)(PQ)\xlongequal(P1AP)(Q1BQ)=A1B1=(λ1B100λmBn)
    P1AP=A1= λ100λm ,Q1BQ= t100tn (PQ)1(AB)(PQ)吸收公式 (P1AP)(Q1BQ)=A1B1= λ1B100λmBn

  3. I n ⊗ I m = I m ⊗ I n = I m × n I_n\otimes I_m=I_{m}\otimes I_n=I_{m\times n} InIm=ImIn=Im×n

4.2.2 计算

a. 分块法

右进右出
( A B C D ) ⊗ F = ( A ⊗ F B ⊗ F C ⊗ F D ⊗ F ) ( A B ) ⊗ F = ( A ⊗ F B ⊗ F ) \left(

ABCD
\right)\otimes F=\left(
AFBFCFDF
\right)\\ (A\quad B)\otimes F=(A\otimes F\quad B\otimes F) (ACBD)F=(AFCFBFDF)(AB)F=(AFBF)

一般情况下: F ⊗ ( A B ) ≠ ( F ⊗ A F ⊗ B ) F \otimes (A\quad B)\neq (F\otimes A\quad F\otimes B) F(AB)=(FAFB)

在这里插入图片描述

b. 向量与向量张量积

列向量 α = ( a 1 a 2 ⋮ a n ) β = ( b 1 b 2 ⋮ b q ) , 则 α ⊗ β = ( a 1 ⊗ β a 2 ⊗ β ⋮ a n ⊗ β ) n q × 1 = ( a 1 b 1 a 1 b 2 ⋮ a 1 b q ⋮ a n b 1 a n b 2 ⋮ a n b q ) n q × 1 列向量\alpha=\left(

a1a2an
\right)\beta=\left(
b1b2bq
\right),则\alpha \otimes \beta=\left(
a1βa2βanβ
\right)_{nq\times 1}=\left(
a1b1a1b2a1bqanb1anb2anbq
\right)_{nq\times 1} 列向量α= a1a2an β= b1b2bq ,αβ= a1βa2βanβ nq×1= a1b1a1b2a1bqanb1anb2anbq nq×1

c. 向量与矩阵张量积

列向量 α = ( a 1 a 2 ⋮ a m ) , B n × q = ( β 1 , β 2 , ⋯   , β q ) , α ⊗ B = ( a 1 B a 2 B ⋮ a m B ) = ( a 1 b 11 a 1 b 12 ⋯ a 1 b 1 q a 1 b 21 a 1 b 22 ⋯ a 1 b 2 q ⋮ ⋮ ⋱ ⋮ a 1 b p 1 a 1 b p 2 ⋯ a 1 b p q a 2 b 11 a 2 b 12 ⋯ a 2 b 1 q a 2 b 21 a 2 b 22 ⋯ a 2 b 2 q ⋮ ⋮ ⋱ ⋮ a 2 b p 1 a 2 b p 2 ⋯ a 2 b p q ⋮ ⋮ ⋮ ⋮ a m b 11 a m b 12 ⋯ a m b 1 q a m b 21 a m b 22 ⋯ a m b 2 q ⋮ ⋮ ⋱ ⋮ a m b p 1 a m b p 2 ⋯ a m b p q ) m p × q 列向量\alpha=\left(

a1a2am
\right),B_{n\times q}=(\beta_1,\beta_2,\cdots,\beta_q),\alpha\otimes B=\left(
a1Ba2BamB
\right)=\left(
a1b11a1b12a1b1qa1b21a1b22a1b2qa1bp1a1bp2a1bpqa2b11a2b12a2b1qa2b21a2b22a2b2qa2bp1a2bp2a2bpqamb11amb12amb1qamb21amb22amb2qambp1ambp2ambpq
\right)_{mp\times q} 列向量α= a1a2am ,Bn×q=(β1,β2,,βq)αB= a1Ba2BamB = a1b11a1b21a1bp1a2b11a2b21a2bp1amb11amb21ambp1a1b12a1b22a1bp2a2b12a2b22a2bp2amb12amb22ambp2a1b1qa1b2qa1bpqa2b1qa2b2qa2bpqamb1qamb2qambpq mp×q

4.2.3 运算律

  1. 数乘: k ( A ⊗ B ) = ( k A ) ⊗ B = A ⊗ ( k B ) k(A\otimes B)=(kA)\otimes B=A\otimes (kB) k(AB)=(kA)B=A(kB)

  2. 分配律(右进右出): ( A + B ) ⊗ C = A ⊗ C + B ⊗ C (A+B)\otimes C=A\otimes C+B\otimes C (A+B)C=AC+BC C ⊗ ( A + B ) = C ⊗ A + C ⊗ B C\otimes(A+B)=C\otimes A+C\otimes B C(A+B)=CA+CB

  3. 结合律: ( A ⊗ B ) ⊗ C = A ⊗ ( B ⊗ C ) (A\otimes B)\otimes C=A\otimes (B\otimes C) (AB)C=A(BC)

    在这里插入图片描述

  4. 吸收律: ( A ⊗ B ) ( C ⊗ D ) = ( A C ) ⊗ ( B D ) (A\otimes B)(C\otimes D)=(AC)\otimes (BD) (AB)(CD)=(AC)(BD)

    张乘 ( 左边括号间乘法个数 ) → 乘张 ( 右边括号内乘法个数相等 ) 张乘(左边括号间乘法个数)\rightarrow 乘张(右边括号内乘法个数相等) 张乘(左边括号间乘法个数)乘张(右边括号内乘法个数相等)

    在这里插入图片描述

    推论:
    若 A = A m × m 为 m 阶方阵 , B = B n × n 为 n 阶方阵,则 ( A ⊗ B ) k = A k ⊗ B k ( A ⊗ I n ) ( I m ⊗ B ) = A ⊗ B ( A 1 ⊗ B 1 ) ( A 2 ⊗ B 2 ) ⋯ ( A k ⊗ B k ) = ( A 1 A 2 ⋯ A k ) ⊗ ( B 1 ⊗ B 2 ⋯ B k ) ( A 1 ⊗ A 2 ⊗ ⋯ ⊗ A k ) ( B 1 ⊗ B 2 ⊗ ⋯ ⊗ B k ) = ( A 1 B 1 ) ⊗ ( A 2 ⊗ ⋯ ⊗ A k ) ( B 2 ⊗ ⋯ ⊗ B k ) = ( A 1 B 1 ) ⊗ ( A 2 B 2 ) ⊗ ⋯ ⊗ ( A k B k )

    A=Am×mm,B=Bn×nn(AB)k=AkBk(AIn)(ImB)=AB(A1B1)(A2B2)(AkBk)=(A1A2Ak)(B1B2Bk)(A1A2Ak)(B1B2Bk)=(A1B1)(A2Ak)(B2Bk)=(A1B1)(A2B2)(AkBk)
    A=Am×mm阶方阵,B=Bn×nn阶方阵,则(AB)k=AkBk(AIn)(ImB)=AB(A1B1)(A2B2)(AkBk)=(A1A2Ak)(B1B2Bk)(A1A2Ak)(B1B2Bk)=(A1B1)(A2Ak)(B2Bk)=(A1B1)(A2B2)(AkBk)

    eg:
    A = A m × m ,证明 e A ⊗ I n = e A ⊗ I n A=A_{m\times m},证明e^{A\otimes I_n}=e^A\otimes I_n A=Am×m,证明eAIn=eAIn

e A ⊗ I n = ∑ k = 1 ∞ 1 k ! ( A ⊗ I ) k = ∑ k = 1 ∞ 1 k ! ( A k ⊗ I k ) = ( ∑ k = 1 ∞ 1 k ! A k ) ⊗ I k = e A ⊗ I

eAIn=k=11k!(AI)k=k=11k!(AkIk)=(k=11k!Ak)Ik=eAI
eAIn=k=1k!1(AI)k=k=1k!1(AkIk)=(k=1k!1Ak)Ik=eAI

  1. 转置与求逆公式:

( A ⊗ B ) H = A H ⊗ B H ( A ⊗ B ) − 1 = A − 1 ⊗ B − 1

(AB)H=AHBH(AB)1=A1B1
(AB)H=AHBH(AB)1=A1B1

  1. 若A与B都是U阵,则 A ⊗ B A\otimes B AB 为U阵

  2. 秩公式: r ( A ⊗ B ) = r ( A ) r ( B ) r(A\otimes B)=r(A)r(B) r(AB)=r(A)r(B)

    推论:
    若 X 1 、 ⋯ 、 X p 为 C m 中 p 个线性无关的列向量, Y 1 、 ⋯ 、 Y q 为 C n 中 q 个线性无关列向量,则 则 p q 个列向量 { X i ⊗ Y j } 线性无关

    X1XpCmp线Y1YqCnq线pq{XiYj}线
    X1XpCmp个线性无关的列向量,Y1YqCnq个线性无关列向量,则pq个列向量{XiYj}线性无关

    由于 r ( { X ⊗ Y } ) = r ( { X } ) r ( { Y } ) r(\{X\otimes Y\})=r(\{X\})r(\{Y\}) r({XY})=r({X})r({Y}) ,张量积的秩等于两个向量组的秩的乘积,所以张量积线性无关

4.2.4 张量积行列式

A = ( a i j ) ∈ C m × m , B = ( b i j ) ∈ C n × n A=(a_{ij}) \in C^{m\times m} ,B=(b_{ij})\in C^{n\times n} A=(aij)Cm×mB=(bij)Cn×n ,则
t r ( A ⊗ B ) = t r ( A ) ⋅ t r ( B ) ∣ A ⊗ B ∣ = ∣ A ∣ n ∣ B ∣ m

tr(AB)=tr(A)tr(B)|AB|=|A|n|B|m
tr(AB)=tr(A)tr(B)AB=AnBm

证明:

在这里插入图片描述

由许尔公式,存在可逆阵 P 使 P − 1 A P = ( λ 1 ∗ ⋱ 0 λ m ) 为上三角,且 P ⊗ I 为可逆阵,构造一个新公式 ( P ⊗ I ) − 1 A ⊗ B ( P ⊗ I ) = ( P − 1 A P ) ⊗ B = A 1 ⊗ B = ( λ 1 B ( ∗ B ) λ 2 B ⋱ 0 λ m B ) 故 ∣ A ⊗ B ∣ = ∣ λ 1 n B ∣ ∣ λ 2 n B ∣ ⋯ ∣ λ m n B ∣ m = ∣ λ 1 n λ 2 n ⋯ λ m n ∣ ∣ B ∣ m = ∣ A ∣ n ∣ B ∣ m

P使P1AP=(λ10λm)PI(PI)1AB(PI)=(P1AP)B=A1B=(λ1B(B)λ2B0λmB)|AB|=|λ1nB||λ2nB||λmnB|m=|λ1nλ2nλmn||B|m=|A|n|B|m
由许尔公式,存在可逆阵P使P1AP= λ10λm 为上三角,且PI为可逆阵,构造一个新公式(PI)1AB(PI)=(P1AP)B=A1B= λ1B0λ2B(B)λmB AB=λ1nB∣∣λ2nBλmnBm=λ1nλ2nλmn∣∣Bm=AnBm

4.2.5 张量积特征值特征向量

a. 特征值

A = A m × m A=A_{m\times m} A=Am×m 的特征根为 λ 1 , λ 2 , ⋯   , λ m , B = B n × n \lambda_1,\lambda_2,\cdots,\lambda_m,B=B_{n\times n} λ1,λ2,,λmB=Bn×n 的特征根为 t 1 , t 2 , ⋯   , t n t_1,t_2,\cdots,t_n t1,t2,,tn

  • A ⊗ B A\otimes B AB 的全体特征根为 m n mn mn 个数 { λ i t j } \{\lambda_it_j\} {λitj} ,有重根

  • A ⊗ I n ± I m ⊗ B A\otimes I_n\pm I_m\otimes B AIn±ImB 是的全体特根为 m n mn mn 个数 { λ i ± t j } \{\lambda_i\pm t_j\} {λi±tj}

b. 特征向量

{ X 1 , ⋯   , X p } \{X_1,\cdots,X_p\} {X1,,Xp} A ∈ C m × m A\in C^{m\times m} ACm×m 关于 λ \lambda λ 的线性无关的特征向量, { Y 1 , ⋯   , Y q } \{Y_1,\cdots,Y_q\} {Y1,,Yq} B ∈ C n × n B\in C^{n\times n} BCn×n 关于t的线性无关的特征向量

  • p q pq pq 个向量 { X i ⊗ Y j } \{X_i\otimes Y_j\} {XiYj} A ⊗ B A\otimes B AB 关于 λ t \lambda t λt 的特征向量
  • { X ⊗ Y } \{X\otimes Y\} {XY} A ⊗ I n ± I m ⊗ B A\otimes I_n\pm I_m\otimes B AIn±ImB 关于 λ + t \lambda+t λ+t 的一个特征向量

在这里插入图片描述

eg

在这里插入图片描述

A = ( 2 2 1 3 ) ⊗ ( 1 − 1 0 1 ) = B ⊗ D , 且 B 为行和等矩阵, D 为上三角阵,则 λ ( A ) = { 4 , t r ( A ) − 4 } = { 4 , 1 } , λ ( B ) = { 1 , 1 } ∴ A ⊗ B = ∏ λ A λ B = { 4 , 4 , 1 , 1 } , 特式 ∣ λ I − A ∣ = ( λ − 4 ) 2 ( λ − 1 ) 2 可知, ( 1 1 ) ( − 2 1 ) 是 B 的特征向量, ( 1 0 ) 是 D 的特征向朗 , 故 A ⊗ B 的特征向量为 ( 1 0 1 0 ) , ( − 2 0 1 0 )

A=(2213)(1101)=BD,BDλ(A)={4,tr(A)4}={4,1},λ(B)={1,1}AB=λAλB={4,4,1,1},|λIA|=(λ4)2(λ1)2(11)(21)B(10)D,AB(1010),(2010)
A=(2123)(1011)=BD,B为行和等矩阵,D为上三角阵,则λ(A)={4,tr(A)4}={4,1},λ(B)={1,1}AB=λAλB={4,4,1,1},特式λIA=(λ4)2(λ1)2可知,(11)(21)B的特征向量,(10)D的特征向朗,AB的特征向量为 1010 , 2010

c. 矩阵函数的特根特向

设 A 和 B 分别是 m 阶与 n 阶方阵, f ( x , y ) f(x,y) f(x,y) 是二元多项式, f ( x , y ) = ∑ i , j = 1 p c i j x i y j f(x,y)=\sum_{i,j=1}\limits^pc_{ij}x^iy^j f(x,y)=i,j=1pcijxiyj

定义 mn 阶矩阵 f ( A , B ) f(A,B) f(A,B) 如下:
f ( A , B ) = ∑ i , j = 1 p c i j A i ⊗ B j , ( A 0 = I m , B 0 = I n ) f(A,B)=\sum_{i,j=1}\limits^pc_{ij}A^i\otimes B^j,(A^0=I_m,B^0=I_n) f(A,B)=i,j=1pcijAiBj,(A0=Im,B0=In)

在这里插入图片描述

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/你好赵伟/article/detail/798373
推荐阅读
相关标签
  

闽ICP备14008679号