关于矩阵函数$(I + x\Sigma)^{1/2}$对$x$的求导问题咨询
嗨,我来帮你梳理这个矩阵求导的问题~咱们一步步拆解推导,把思路理清楚:
第一步:明确矩阵平方根的定义
因为$\Sigma$是正定矩阵,只要$x$的取值能保证$I + x\Sigma$正定(比如$x \geq 0$,或者更宽泛的范围只要矩阵所有特征值为正),那么$(I + x\Sigma)^{1/2}$就是唯一的正定矩阵$A(x)$,满足核心等式:
$$A(x)^2 = I + x\Sigma$$
第二步:对核心等式两边求导
你提到的矩阵乘积求导法则$\frac{dUV}{dx} = \frac{dU}{dx}V + U\frac{dV}{dx}$完全正确,咱们用它对等式两边关于$x$求导:
- 左边展开:$\frac{d}{dx}[A(x)A(x)] = A'(x)A(x) + A(x)A'(x)$
- 右边求导:$\frac{d}{dx}(I + x\Sigma) = \Sigma$
由此得到关键方程:
$$A'(x)A(x) + A(x)A'(x) = \Sigma$$
第三步:利用可交换性简化方程
这里有个核心性质:$A(x)$和$\Sigma$是可交换的。原因很简单:$\Sigma$和$I+x\Sigma$显然可交换,而矩阵平方根$A(x)$可以表示为$\Sigma$的幂级数(正定矩阵的平方根可通过泰勒展开定义):
$$A(x) = (I + x\Sigma)^{1/2} = \sum_{k=0}^\infty \binom{1/2}{k} (x\Sigma)^k$$
每一项都是$\Sigma$的幂,自然和$\Sigma$可交换。进一步,$A'(x)$是这个级数逐项求导的结果,同样是$\Sigma$的幂的组合,所以$A'(x)$和$A(x)$也可交换。
既然$A'(x)$和$A(x)$可交换,左边的式子就能合并:
$$2A(x)A'(x) = \Sigma$$
第四步:推导最终导数表达式
两边同时左乘$A(x)^{-1}$,得到:
$$A'(x) = \frac{1}{2}A(x)^{-1}\Sigma$$
又因为$A(x)^2 = I + x\Sigma$,所以$A(x)^{-1} = A(x)/(I + x\Sigma)$,代入后也可以写成更直观的形式:
$$\frac{d}{dx}(I + x\Sigma)^{1/2} = \frac{1}{2}(I + x\Sigma)^{-1/2}\Sigma$$
由于可交换性,也等价于$\frac{1}{2}\Sigma(I + x\Sigma)^{-1/2}$
对你之前尝试的补充
你用到的矩阵乘积求导法则是推导的关键起点,之所以一开始没直接得到结果,是因为缺少了矩阵可交换性这个重要细节,这也是矩阵求导和标量求导的核心区别之一~
备注:内容来源于stack exchange,提问作者jason 1

