作用在兩個相同形狀的矩陣上的阿達瑪乘積,結果是第三個相同形狀的矩陣。
「
逐項積 」重新導向至此。關於函數的逐點乘積,請見「
逐點乘積 」。
在數學 中,阿達瑪乘積 (英語:Hadamard product ,又譯哈達瑪乘積 ),又名舒爾乘積 (Schur product )[ 1] 或逐項乘積 (entrywise product )[ 2] :ch. 5 ,是一個二元運算 ,其輸入為兩個相同形狀的矩陣 ,輸出是具有同樣形狀的、各個位置的元素等於兩個輸入矩陣相同位置元素的乘積的矩陣。此乘積歸功於法國數學家雅克·阿達馬 或德國數學家伊賽·舒爾 ,並以其命名。
定義
若兩個矩陣
A
{\displaystyle A}
和
B
{\displaystyle B}
具有相同的維度
m
×
n
{\displaystyle m\times n}
,則它們的阿達瑪乘積
A
∘
B
{\displaystyle A\circ B}
是一個具有相同維度的矩陣,其元素值為:
(
A
∘
B
)
i
j
=
(
A
)
i
j
(
B
)
i
j
.
{\displaystyle (A\circ B)_{ij}=(A)_{ij}(B)_{ij}.}
對於維度不相等的矩陣(m × n 矩陣和 p × q 矩陣,其中m ≠ p 或n ≠ q ),阿達瑪乘積沒有定義。
樣例
3
×
3
{\displaystyle 3\times 3}
矩陣A 與
3
×
3
{\displaystyle 3\times 3}
矩陣B 的阿達瑪乘積為:
[
a
11
a
12
a
13
a
21
a
22
a
23
a
31
a
32
a
33
]
∘
[
b
11
b
12
b
13
b
21
b
22
b
23
b
31
b
32
b
33
]
=
[
a
11
b
11
a
12
b
12
a
13
b
13
a
21
b
21
a
22
b
22
a
23
b
23
a
31
b
31
a
32
b
32
a
33
b
33
]
.
{\displaystyle {\begin{bmatrix}a_{11}&a_{12}&a_{13}\\a_{21}&a_{22}&a_{23}\\a_{31}&a_{32}&a_{33}\end{bmatrix}}\circ {\begin{bmatrix}b_{11}&b_{12}&b_{13}\\b_{21}&b_{22}&b_{23}\\b_{31}&b_{32}&b_{33}\end{bmatrix}}={\begin{bmatrix}a_{11}\,b_{11}&a_{12}\,b_{12}&a_{13}\,b_{13}\\a_{21}\,b_{21}&a_{22}\,b_{22}&a_{23}\,b_{23}\\a_{31}\,b_{31}&a_{32}\,b_{32}&a_{33}\,b_{33}\end{bmatrix}}.}
性質
阿達瑪乘積滿足交換律 (當其元素屬於交換環時), 結合律 和對加法的分配律 :
A
∘
B
=
B
∘
A
,
A
∘
(
B
∘
C
)
=
(
A
∘
B
)
∘
C
,
A
∘
(
B
+
C
)
=
A
∘
B
+
A
∘
C
.
{\displaystyle {\begin{aligned}&\mathbf {A} \circ \mathbf {B} =\mathbf {B} \circ \mathbf {A} ,\\&\mathbf {A} \circ (\mathbf {B} \circ \mathbf {C} )=(\mathbf {A} \circ \mathbf {B} )\circ \mathbf {C} ,\\&\mathbf {A} \circ (\mathbf {B} +\mathbf {C} )=\mathbf {A} \circ \mathbf {B} +\mathbf {A} \circ \mathbf {C} .\end{aligned}}}
在阿達瑪乘積意義下,m × n 矩陣的單位元 是全部元素均為1的m × n 矩陣 。這跟普通矩陣乘法的單位元 只有主對角線上的元素為1的單位矩陣 不同。此外,若且唯若沒有任何元素等於 0 時,矩陣的阿達瑪乘積有逆矩陣。[ 3]
對於向量x 和y ,以及以這些向量作為主對角線的對應對角矩陣D x 和D y ,以下恆等式成立:[ 2] :479
x
∗
(
A
∘
B
)
y
=
tr
(
D
x
∗
A
D
y
B
T
)
,
{\displaystyle \mathbf {x} ^{*}({A}\circ {B})\mathbf {y} =\operatorname {tr} \left({D}_{\mathbf {x} }^{*}{A}{D}_{\mathbf {y} }{B}^{\mathsf {T}}\right),}
其中x * 指的是x 的共軛轉置 。特別的,使用全1向量,可以發現阿達瑪乘積的所有元素求和是AB T (上標T表示矩陣轉置 )的跡 。對於方陣 A 和 B ,一個相似的結論是他們的阿達瑪乘積按行求和後得到的向量是AB T 的對角元素組成的向量:[ 4]
∑
i
(
A
∘
B
)
i
j
=
(
B
T
A
)
j
j
=
(
A
B
T
)
i
i
.
{\displaystyle \sum _{i}(A\circ B)_{ij}=\left(B^{\mathsf {T}}A\right)_{jj}=\left(AB^{\mathsf {T}}\right)_{ii}.}
相似地,
(
y
x
∗
)
∘
A
=
D
y
A
D
x
∗
{\displaystyle \left(\mathbf {y} \mathbf {x} ^{*}\right)\circ {A}={D}_{\mathbf {y} }{A}{D}_{\mathbf {x} }^{*}}
此外,阿達瑪乘積矩陣與向量的積可以表示為:
(
A
∘
B
)
y
=
diag
(
A
D
y
B
T
)
{\displaystyle (A\circ B)\mathbf {y} =\operatorname {diag} \left(AD_{\mathbf {y} }B^{\mathsf {T}}\right)}
其中
diag
(
M
)
{\displaystyle \operatorname {diag} (M)}
是M 的對角元素組成的向量。
阿達瑪乘積是克羅內克乘積 的主要子矩陣 。
阿達瑪乘積滿足秩不等式
rank
(
A
∘
B
)
≤
rank
(
A
)
rank
(
B
)
{\displaystyle \operatorname {rank} ({A}\circ {B})\leq \operatorname {rank} ({A})\operatorname {rank} ({B})}
如果A 和B 是正定矩陣 ,那麼下列不等式成立:[ 5]
∏
i
=
k
n
λ
i
(
A
∘
B
)
≥
∏
i
=
k
n
λ
i
(
A
B
)
,
k
=
1
,
…
,
n
,
{\displaystyle \prod _{i=k}^{n}\lambda _{i}({A}\circ {B})\geq \prod _{i=k}^{n}\lambda _{i}({A}{B}),\quad k=1,\ldots ,n,}
其中λi (A ) 是A 的第i 大的特徵值 。
如果D 與E 是對角矩陣 ,那麼[ 6]
D
(
A
∘
B
)
E
=
(
D
A
E
)
∘
B
=
(
D
A
)
∘
(
B
E
)
=
(
A
E
)
∘
(
D
B
)
=
A
∘
(
D
B
E
)
.
{\displaystyle {\begin{aligned}{D}({A}\circ {B}){E}&=({D}{A}{E})\circ {B}=({D}{A})\circ ({B}{E})\\&=({A}{E})\circ ({D}{B})={A}\circ ({D}{B}{E}).\end{aligned}}}
兩個向量
a
{\displaystyle \mathbf {a} }
和
b
{\displaystyle \mathbf {b} }
的阿達瑪乘積與一個向量和另一個向量對應的對角矩陣 做矩陣乘法得到的結果相同:
a
∘
b
=
D
a
b
=
D
b
a
{\displaystyle \mathbf {a} \circ \mathbf {b} =D_{\mathbf {a} }\mathbf {b} =D_{\mathbf {b} }\mathbf {a} }
將向量映射到對角矩陣的
diag
{\displaystyle \operatorname {diag} }
運算 可以用阿達瑪乘積寫為:
diag
(
a
)
=
(
a
1
T
)
∘
I
{\displaystyle \operatorname {diag} (\mathbf {a} )=(\mathbf {a} \mathbf {1} ^{T})\circ I}
其中
1
{\displaystyle \mathbf {1} }
是全
1
{\displaystyle 1}
向量 ,
I
{\displaystyle I}
是單位矩陣 。
參考資料
^ Davis, Chandler. The norm of the Schur product operation. Numerische Mathematik. 1962, 4 (1): 343–44. doi:10.1007/bf01386329 .
^ 2.0 2.1 Horn, Roger A.; Johnson, Charles R. Matrix analysis. Cambridge University Press. 2012.
^ Million, Elizabeth. The Hadamard Product (PDF) . [2 January 2012] . (原始內容 (PDF) 存檔於2013-06-12).
^ Styan, George P. H., Hadamard Products and Multivariate Statistical Analysis, Linear Algebra and Its Applications, 1973, 6 : 217–240, doi:10.1016/0024-3795(73)90023-2 , hdl:10338.dmlcz/102190
^ Hiai, Fumio; Lin, Minghua. On an eigenvalue inequality involving the Hadamard product. Linear Algebra and Its Applications. February 2017, 515 : 313–320. doi:10.1016/j.laa.2016.11.017 .
^ Project (PDF) . buzzard.ups.edu. 2007 [2019-12-18 ] . (原始內容 (PDF) 存檔於2013-06-12).