logo

二次形式と双線形形式の行列微分 📂多変数ベクトル解析

二次形式と双線形形式の行列微分

公式

  • 2つのベクトル aRm\mathbf{a} \in \mathbb{R}^{m}, bRn\mathbf{b} \in \mathbb{R}^{n}行列 XRm×n\mathbf{X} \in \mathbb{R}^{m \times n}について、双線形形式 aTXb\mathbf{a}^{\mathsf{T}} \mathbf{X} \mathbf{b}に対する勾配行列は次の通りだ。

    X(aTXb)=(aTXb)X=abT(1) \nabla_{\mathbf{X}} (\mathbf{a}^{\mathsf{T}} \mathbf{X} \mathbf{b}) = \dfrac{\partial (\mathbf{a}^{\mathsf{T}} \mathbf{X} \mathbf{b})}{\partial \mathbf{X}} = \mathbf{a}\mathbf{b}^{\mathsf{T}} \tag{1} 従って、二次形式 aTXa\mathbf{a}^{\mathsf{T}} \mathbf{X} \mathbf{a}について次のことが成り立つ。 (aTXa)X=aaT \dfrac{\partial (\mathbf{a}^{\mathsf{T}} \mathbf{X} \mathbf{a})}{\partial \mathbf{X}} = \mathbf{a}\mathbf{a}^{\mathsf{T}}

  • 2つのベクトル aRn\mathbf{a} \in \mathbb{R}^{n}, bRn\mathbf{b} \in \mathbb{R}^{n}行列 XRm×n\mathbf{X} \in \mathbb{R}^{m \times n}について、 aTXTXb\mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}}\mathbf{X} \mathbf{b}に対する勾配行列は次の通りだ。

    X(aTXTXb)=(aTXTXb)X=X(abT+baT)(2) \nabla_{\mathbf{X}} (\mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}} \mathbf{X} \mathbf{b}) = \dfrac{\partial (\mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}} \mathbf{X} \mathbf{b})}{\partial \mathbf{X}} = \mathbf{X}(\mathbf{a} \mathbf{b}^{\mathsf{T}} + \mathbf{b}\mathbf{a}^{\mathsf{T}} ) \tag{2}

    もし b=a\mathbf{b} = \mathbf{a}なら、

    (aTXTXa)X=2XaaT(2) \dfrac{\partial (\mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}} \mathbf{X} \mathbf{a})}{\partial \mathbf{X}} = 2\mathbf{X}\mathbf{a} \mathbf{a}^{\mathsf{T}} \tag{2}

説明

多項式の微分と類似した結果だ。

  • (1)(1): 本質的には一次関数の微分と同じなので、係数だけが残る形になることが分かる。少し混乱するかもしれない点は、計算結果が行列でなければならないためaTb\mathbf{a}^{\mathsf{T}}\mathbf{b}ではなくabT\mathbf{a}\mathbf{b}^{\mathsf{T}}だということだ。
  • (2)(2): 二次形式の観点から見るとaTXa\mathbf{a}^{\mathsf{T}} \mathbf{X} \mathbf{a}aTXTXa\mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}}\mathbf{X} \mathbf{a}もどれも単に二次形式だが、微分する変数が何度も掛けられているため、異なる結果を示す。本質的には二次関数の微分に対応する。

以下の証明では直接計算で示したが、 🔒(25/07/17)トレーストリックとして知られる方法を使用すると、もっと簡単に計算できる。X\mathbf{X}に対する任意の形や、X\mathbf{X}が多く含まれた式の場合、実際直接計算するのは非常に難しく、トレーストリックを使用する必要がある。

より多くの公式は、🔒(25/07/19)スカラ関数の行列微分表で確認できる。

証明

(1)(1)

双線形形式はaTXb=i=1mj=1naixijbj\mathbf{a}^{\mathsf{T}} \mathbf{X} \mathbf{b} = \sum\limits_{i=1}^{m}\sum\limits_{j=1}^{n} a_{i} x_{ij} b_{j}で表せる。よって(i=1mj=1naixijbj)xk=akb\dfrac{\partial \left( \sum\limits_{i=1}^{m}\sum\limits_{j=1}^{n} a_{i} x_{ij} b_{j} \right)}{\partial x_{k\ell}} = a_{k}b_{\ell}であり、

X(aTXb)=[(aTXb)x11(aTXb)x1n(aTXb)xm1(aTXb)xmn]=[a1b1a1b2a1bna2b1a2b2a2bnamb1amb1ambn]=[a1an][b1bn]=abT=ab \begin{align*} \nabla_{\mathbf{X}} (\mathbf{a}^{\mathsf{T}} \mathbf{X} \mathbf{b}) &= \begin{bmatrix} \dfrac{\partial (\mathbf{a}^{\mathsf{T}} \mathbf{X} \mathbf{b})}{\partial x_{11}} & \cdots & \dfrac{\partial (\mathbf{a}^{\mathsf{T}} \mathbf{X} \mathbf{b})}{\partial x_{1n}} \\ \vdots & \ddots & \vdots \\ \dfrac{\partial (\mathbf{a}^{\mathsf{T}} \mathbf{X} \mathbf{b})}{\partial x_{m1}} & \cdots & \dfrac{\partial (\mathbf{a}^{\mathsf{T}} \mathbf{X} \mathbf{b})}{\partial x_{mn}} \end{bmatrix} \\ &= \begin{bmatrix} a_{1}b_{1} & a_{1}b_{2} & \cdots & a_{1}b_{n} \\ a_{2}b_{1} & a_{2}b_{2} & \cdots & a_{2}b_{n} \\ \vdots & \vdots & \ddots & \vdots \\ a_{m}b_{1} & a_{m}b_{1} & \cdots & a_{m}b_{n} \end{bmatrix} \\ &= \begin{bmatrix} a_{1} \\ \vdots \\ a_{n} \end{bmatrix} \begin{bmatrix} b_{1} & \cdots & b_{n} \end{bmatrix} \\ &= \mathbf{a} \mathbf{b}^{\mathsf{T}} = \mathbf{a} \otimes \mathbf{b} \end{align*}

\otimesは2つのベクトルの外積、または2つの行列のクロネッカー積である。

(2)(2)

直接計算

2つの行列の積は[AB]ij=k=1naikbkj[AB]_{ij} = \sum\limits_{k=1}^{n} a_{ik}b_{kj}であったため、[XTX]ij=s=1mxsixsj[X^{\mathsf{T}}X]_{ij} = \sum\limits_{s=1}^{m} x_{si} x_{sj}である。次が成り立つ。

aTXTXb=k,=1nak(s=1mxskxs)b=k,=1ns=1makxskxsb \mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}}\mathbf{X} \mathbf{b} = \sum\limits_{k, \ell = 1}^{n} a_{k} \left(\sum\limits_{s=1}^{m} x_{sk} x_{s\ell}\right) b_{\ell} = \sum\limits_{k,\ell = 1}^{n} \sum\limits_{s=1}^{m} a_{k} x_{sk} x_{s\ell} b_{\ell}

したがって、微分は次のようになる。

(k,=1ns=1makxskxsb)xij==1najxib+k=1nakxikbj \dfrac{\partial \left( \sum\limits_{k,\ell = 1}^{n} \sum\limits_{s=1}^{m} a_{k} x_{sk} x_{s\ell} b_{\ell} \right)}{\partial x_{ij}} = \sum\limits_{\ell=1}^{n}a_{j}x_{i\ell}b_{\ell} + \sum\limits_{k = 1}^{n} a_{k}x_{ik}b_{j}

勾配行列は以下の通りとなる。

X(aTXTXb)=[=1na1x1b+k=1nakx1kb1=1na2x1b+k=1nakx1kb2=1nanx1b+k=1nakx1kbn=1na1x2b+k=1nakx2kb1=1na2x2b+k=1nakx2kb2=1nanx2b+k=1nakx2kbn=1na1xnb+k=1nakxnkb1=1na2xnb+k=1nakxnkb2=1nanxnb+k=1nakxnkbn]=[=1na1x1b=1nanx1b=1na1xnb=1nanxnb]+[k=1nakx1kb1k=1nakx1kbnk=1nakxnkb1k=1nakxnkbn]=[=1nx1ba1=1nx1ban=1nxnba1=1nxnban]+[k=1nx1kakb1k=1nx1kakbnk=1nxnkakb1k=1nxnkakbn] \begin{align*} & \nabla_{\mathbf{X}} (\mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}} \mathbf{X} \mathbf{b}) \\ &= \begin{bmatrix} \sum\limits_{\ell=1}^{n}a_{1}x_{1\ell}b_{\ell} + \sum\limits_{k = 1}^{n} a_{k}x_{1k}b_{1} & \sum\limits_{\ell=1}^{n}a_{2}x_{1\ell}b_{\ell} + \sum\limits_{k = 1}^{n} a_{k}x_{1k}b_{2} & \cdots & \sum\limits_{\ell=1}^{n}a_{n}x_{1\ell}b_{\ell} + \sum\limits_{k = 1}^{n} a_{k}x_{1k}b_{n} \\ \sum\limits_{\ell=1}^{n}a_{1}x_{2\ell}b_{\ell} + \sum\limits_{k = 1}^{n} a_{k}x_{2k}b_{1} & \sum\limits_{\ell=1}^{n}a_{2}x_{2\ell}b_{\ell} + \sum\limits_{k = 1}^{n} a_{k}x_{2k}b_{2} & \cdots & \sum\limits_{\ell=1}^{n}a_{n}x_{2\ell}b_{\ell} + \sum\limits_{k = 1}^{n} a_{k}x_{2k}b_{n} \\ \vdots & \vdots & \ddots & \vdots \\ \sum\limits_{\ell=1}^{n}a_{1}x_{n\ell}b_{\ell} + \sum\limits_{k = 1}^{n} a_{k}x_{nk}b_{1} & \sum\limits_{\ell=1}^{n}a_{2}x_{n\ell}b_{\ell} + \sum\limits_{k = 1}^{n} a_{k}x_{nk}b_{2} & \cdots & \sum\limits_{\ell=1}^{n}a_{n}x_{n\ell}b_{\ell} + \sum\limits_{k = 1}^{n} a_{k}x_{nk}b_{n} \end{bmatrix} \\ &= \begin{bmatrix} \sum\limits_{\ell=1}^{n}a_{1}x_{1\ell}b_{\ell} & \cdots & \sum\limits_{\ell=1}^{n}a_{n}x_{1\ell}b_{\ell} \\ \vdots & \ddots & \vdots \\ \sum\limits_{\ell=1}^{n}a_{1}x_{n\ell}b_{\ell} & \cdots & \sum\limits_{\ell=1}^{n}a_{n}x_{n\ell}b_{\ell} \end{bmatrix} + \begin{bmatrix} \sum\limits_{k = 1}^{n} a_{k}x_{1k}b_{1} & \cdots &\sum\limits_{k = 1}^{n} a_{k}x_{1k}b_{n} \\ \vdots & \ddots & \vdots \\ \sum\limits_{k = 1}^{n} a_{k}x_{nk}b_{1} & \cdots &\sum\limits_{k = 1}^{n} a_{k}x_{nk}b_{n} \end{bmatrix} \\ &= \begin{bmatrix} \sum\limits_{\ell=1}^{n}x_{1\ell}b_{\ell}a_{1} & \cdots & \sum\limits_{\ell=1}^{n}x_{1\ell}b_{\ell}a_{n} \\ \vdots & \ddots & \vdots \\ \sum\limits_{\ell=1}^{n}x_{n\ell}b_{\ell}a_{1} & \cdots & \sum\limits_{\ell=1}^{n}x_{n\ell}b_{\ell}a_{n} \end{bmatrix} + \begin{bmatrix} \sum\limits_{k = 1}^{n} x_{1k}a_{k}b_{1} & \cdots &\sum\limits_{k = 1}^{n} x_{1k}a_{k}b_{n} \\ \vdots & \ddots & \vdots \\ \sum\limits_{k = 1}^{n} x_{nk}a_{k}b_{1} & \cdots &\sum\limits_{k = 1}^{n} x_{nk}a_{k}b_{n} \end{bmatrix} \\ \end{align*}

したがって、次を得る。

X(aTXTXb)=XbaT+XabT=X(abT+baT) \nabla_{\mathbf{X}} (\mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}} \mathbf{X} \mathbf{b}) = \mathbf{X}\mathbf{b}\mathbf{a}^{\mathsf{T}} + \mathbf{X}\mathbf{a}\mathbf{b}^{\mathsf{T}} = \mathbf{X} ( \mathbf{a}\mathbf{b}^{\mathsf{T}} + \mathbf{b}\mathbf{a}^{\mathsf{T}} )

トレーストリック

スカラーは1×11\times 1行列と考えることができ、そうすればトレース(跡)は事実上、恒等関数の役割をする。つまり、双線形形式の値はスカラーなので次が成り立つ。

aTXTXb=Tr(aTXTXb) \mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}} \mathbf{X} \mathbf{b} = \Tr(\mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}} \mathbf{X} \mathbf{b})

また、トレースは巡回性 Tr(ABC)=Tr(CBA)\Tr(ABC) = \Tr(CBA) を持つので、

(aTXTXb)X=Tr(aTXTXb)X=Tr(baTXTX)X \dfrac{\partial (\mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}} \mathbf{X} \mathbf{b})}{\partial \mathbf{X}} = \dfrac{\partial \Tr(\mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}} \mathbf{X} \mathbf{b})}{\partial \mathbf{X}} = \dfrac{\partial \Tr(\mathbf{b}\mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}} \mathbf{X})}{\partial \mathbf{X}}

トレース微分公式によって次を得る。

(baTXTX)X=X(baT+abT) \dfrac{\partial (\mathbf{b}\mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}} \mathbf{X})}{\partial \mathbf{X}} = \mathbf{X}(\mathbf{b}\mathbf{a}^{\mathsf{T}} + \mathbf{a}\mathbf{b}^{\mathsf{T}})

このときb=a\mathbf{b} = \mathbf{a}なら、

(aTXTXa)X=2XaaT \dfrac{\partial (\mathbf{a}^{\mathsf{T}} \mathbf{X}^{\mathsf{T}} \mathbf{X} \mathbf{a})}{\partial \mathbf{X}} = 2\mathbf{X}\mathbf{a}\mathbf{a}^{\mathsf{T}}

このようにトレースの性質を利用して行列微分の計算を簡単な行列操作にする方法を、🔒(25/07/17)トレーストリックという。