克罗内克积

zz: http://zh.wikipedia.org/zh/%E5%85%8B%E7%BD%97%E5%86%85%E5%85%8B%E7%A7%AF

定义

如果A是一个 m × n 的矩阵,而B是一个 p × q 的矩阵,克罗内克积A /otimes B则是一个 mp × nq 的分块矩阵

 A /otimes B = /begin{bmatrix} a_{11} B & /cdots & a_{1n}B // /vdots & /ddots & /vdots // a_{m1} B & /cdots & a_{mn} B /end{bmatrix}.

更具体地可表示为

 A /otimes B = /begin{bmatrix}
   a_{11} b_{11} & a_{11} b_{12} & /cdots & a_{11} b_{1q} & 
                   /cdots & /cdots & a_{1n} b_{11} & a_{1n} b_{12} & /cdots & a_{1n} b_{1q} //
   a_{11} b_{21} & a_{11} b_{22} & /cdots & a_{11} b_{2q} & 
                   /cdots & /cdots & a_{1n} b_{21} & a_{1n} b_{22} & /cdots & a_{1n} b_{2q} //
   /vdots & /vdots & /ddots & /vdots & & & /vdots & /vdots & /ddots & /vdots //
   a_{11} b_{p1} & a_{11} b_{p2} & /cdots & a_{11} b_{pq} & 
                   /cdots & /cdots & a_{1n} b_{p1} & a_{1n} b_{p2} & /cdots & a_{1n} b_{pq} //
   /vdots & /vdots & & /vdots & /ddots & & /vdots & /vdots & & /vdots //
   /vdots & /vdots & & /vdots & & /ddots & /vdots & /vdots & & /vdots //
   a_{m1} b_{11} & a_{m1} b_{12} & /cdots & a_{m1} b_{1q} & 
                   /cdots & /cdots & a_{mn} b_{11} & a_{mn} b_{12} & /cdots & a_{mn} b_{1q} //
   a_{m1} b_{21} & a_{m1} b_{22} & /cdots & a_{m1} b_{2q} & 
                   /cdots & /cdots & a_{mn} b_{21} & a_{mn} b_{22} & /cdots & a_{mn} b_{2q} //
   /vdots & /vdots & /ddots & /vdots & & & /vdots & /vdots & /ddots & /vdots //
   a_{m1} b_{p1} & a_{m1} b_{p2} & /cdots & a_{m1} b_{pq} & 
                   /cdots & /cdots & a_{mn} b_{p1} & a_{mn} b_{p2} & /cdots & a_{mn} b_{pq} 
/end{bmatrix}.

[编辑]例子


  /begin{bmatrix} 
    1 & 2 // 
    3 & 1 // 
  /end{bmatrix}
/otimes
  /begin{bmatrix} 
    0 & 3 // 
    2 & 1 // 
  /end{bmatrix}
=
  /begin{bmatrix} 
    1/cdot 0 & 1/cdot 3 & 2/cdot 0 & 2/cdot 3 // 
    1/cdot 2 & 1/cdot 1 & 2/cdot 2 & 2/cdot 1 // 
    3/cdot 0 & 3/cdot 3 & 1/cdot 0 & 1/cdot 3 // 
    3/cdot 2 & 3/cdot 1 & 1/cdot 2 & 1/cdot 1 // 
  /end{bmatrix}
=
  /begin{bmatrix} 
    0 & 3 & 0 & 6 // 
    2 & 1 & 4 & 2 //
    0 & 9 & 0 & 3 //
    6 & 3 & 2 & 1
  /end{bmatrix}.

/begin{bmatrix}
a_{11} & a_{12} //
a_{21} & a_{22} //
a_{31} & a_{32}
/end{bmatrix}
/otimes
/begin{bmatrix}
b_{11} & b_{12} & b_{13} //
b_{21} & b_{22} & b_{23}
/end{bmatrix}
=
/begin{bmatrix}
a_{11} b_{11} & a_{11} b_{12} & a_{11} b_{13} & a_{12} b_{11} & a_{12} b_{12} & a_{12} b_{13} //
a_{11} b_{21} & a_{11} b_{22} & a_{11} b_{23} & a_{12} b_{21} & a_{12} b_{22} & a_{12} b_{23} //
a_{21} b_{11} & a_{21} b_{12} & a_{21} b_{13} & a_{22} b_{11} & a_{22} b_{12} & a_{22} b_{13} //
a_{21} b_{21} & a_{21} b_{22} & a_{21} b_{23} & a_{22} b_{21} & a_{22} b_{22} & a_{22} b_{23} //
a_{31} b_{11} & a_{31} b_{12} & a_{31} b_{13} & a_{32} b_{11} & a_{32} b_{12} & a_{32} b_{13} //
a_{31} b_{21} & a_{31} b_{22} & a_{31} b_{23} & a_{32} b_{21} & a_{32} b_{22} & a_{32} b_{23}
/end{bmatrix}.

[编辑]特性

[编辑]双线性和结合律

克罗内克积张量积的特殊形式,因此满足双线性结合律

 A /otimes (B+C) = A /otimes B + A /otimes C /qquad /mbox{(if } B /mbox{ and } C /mbox{ have the same size)},
 (A+B) /otimes C = A /otimes C + B /otimes C /qquad /mbox{(if } A /mbox{ and } B /mbox{ have the same size)},
 (kA) /otimes B = A /otimes (kB) = k(A /otimes B),
 (A /otimes B) /otimes C = A /otimes (B /otimes C),

其中,AB 和 C 是矩阵,而 k 是常量。

克罗内克积不符合交换律:通常,A /otimes B 不同于 B /otimes A

A /otimes BB /otimes A是排列等价的,也就是说,存在排列矩阵PQ,使得

 A /otimes B = P /, (B /otimes A) /, Q.

如果AB是方块矩阵,则A /otimes BB /otimes A甚至是排列相似的,也就是说,我们可以取P = QT

[编辑]混合乘积性质

如果ABCD是四个矩阵,且矩阵乘积ACBD存在,那么:

 (/mathbf{A} /otimes /mathbf{B})(/mathbf{C} /otimes /mathbf{D}) = /mathbf{AC} /otimes /mathbf{BD}.

这个性质称为“混合乘积性质”,因为它混合了通常的矩阵乘积和克罗内克积。于是可以推出,A /,/otimes/, B可逆当且仅当AB是可逆的,其逆矩阵为:

 (/mathbf{A} /otimes /mathbf{B})^{-1} = /mathbf{A}^{-1} /otimes /mathbf{B}^{-1}.

[编辑]克罗内克和

如果An × n矩阵,Bm × m矩阵,/mathbf{I}_k表示k × k单位矩阵,那么我们可以定义克罗内克和/oplus为:

 /mathbf{A} /oplus /mathbf{B} = /mathbf{A} /otimes /mathbf{I}_m + /mathbf{I}_n /otimes /mathbf{B}.


### 克罗内克积的定义与运算法则 克罗内克积是一种特殊的矩阵乘法操作,记作 \( \otimes \)。给定两个矩阵 \( A \in \mathbb{R}^{m \times n} \) 和 \( B \in \mathbb{R}^{p \times q} \),它们的克罗内克积是一个更大的矩阵,其维度为 \( (mp) \times (nq) \)[^3]。 具体来说,如果矩阵 \( A = [a_{ij}] \) 是一个 \( m \times n \) 的矩阵,则克罗内克积的结果可以写成如下形式: \[ A \otimes B = \begin{bmatrix} a_{11}B & a_{12}B & \cdots & a_{1n}B \\ a_{21}B & a_{22}B & \cdots & a_{2n}B \\ \vdots & \vdots & \ddots & \vdots \\ a_{m1}B & a_{m2}B & \cdots & a_{mn}B \end{bmatrix}. \] 这意味着对于矩阵 \( A \) 中的每一个元素 \( a_{ij} \),将其替换为 \( a_{ij}B \),即该元素与整个矩阵 \( B \) 的数乘结果[^3]。 #### 数学定义的例子 假设矩阵 \( A \) 和 \( B \) 如下所示: \[ A = \begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix}, \quad B = \begin{bmatrix} 5 & 6 \\ 7 & 8 \end{bmatrix}. \] 那么它们的克罗内克积 \( A \otimes B \) 可以通过以下步骤计算得出: \[ A \otimes B = \begin{bmatrix} 1 \cdot B & 2 \cdot B \\ 3 \cdot B & 4 \cdot B \end{bmatrix} = \begin{bmatrix} 1 \cdot \begin{bmatrix} 5 & 6 \\ 7 & 8 \end{bmatrix} & 2 \cdot \begin{bmatrix} 5 & 6 \\ 7 & 8 \end{bmatrix} \\ 3 \cdot \begin{bmatrix} 5 & 6 \\ 7 & 8 \end{bmatrix} & 4 \cdot \begin{bmatrix} 5 & 6 \\ 7 & 8 \end{bmatrix} \end{bmatrix}. \] 进一步展开得: \[ A \otimes B = \begin{bmatrix} 5 & 6 & 10 & 12 \\ 7 & 8 & 14 & 16 \\ 15 & 18 & 20 & 24 \\ 21 & 24 & 28 & 32 \end{bmatrix}. \][^3] #### 使用场景 克罗内克积广泛应用于线性代数、信号处理以及量子力学等领域。例如,在张量分解和图像处理中,它常被用来构建复合系统的状态空间模型;在控制理论中,用于描述多个子系统之间的耦合关系。 以下是 Python 实现的一个简单示例代码片段,展示如何利用 NumPy 库完成克罗内克积运算: ```python import numpy as np # 定义输入矩阵 A = np.array([[1, 2], [3, 4]]) B = np.array([[5, 6], [7, 8]]) # 计算克罗内克积 result = np.kron(A, B) print(result) ``` 运行此程序会输出上述手动计算所得的结果矩阵。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值