matlab - 如何计算多元正态概率密度函数的导数

Question

是否有任何内置函数计算给定点的多元正态概率密度函数的梯度值？

编辑：发现这个如何在matlab中评估函数的导数？但这不是我要找的

Edit2: 好的，这就是我正在使用的http://www.mathworks.co.uk/help/stats/mvnpdf.html案例 3，寻找关于 X 的导数值

Answer 1

我可以建议您查看 Peterson 和 Pedersen 的 The Matrix Cookbook（可在线免费获得 - 只需 google 即可）。您的问题的解析解在 p39，方程 325（2008 年版）。

我们甚至不需要 Matlab 来做这个！

编辑：正如 YBE 所暗示的，也许我应该在我的答案中包含解决方案。因此，设 p(x) 表示多元高斯 pdf，以均值向量 m 和协方差矩阵 S 为特征。然后：

dp(x) / dx = -p(x) * S^(-1) * (x - m)

和

d^2p / dx dx' = p(x) * (S^(-1) (x - m)(x - m)' S^(-1) - S^(-1))

如果你想要一个 Matlab 函数，那么：

function Gradient = MultNormD1(x, Mu, Sigma)
Gradient = -1 * mvnpdf(x, Mu, Sigma) * (Sigma \ (x - Mu));