微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

如何在 Math.NET 中对多元正态分布进行采样?

如何解决如何在 Math.NET 中对多元正态分布进行采样?

在这里遗漏了什么明显的东西吗? Math.NET 有 a wealth 个概率分布抽样类,但没有多元正态分布。它有 normalMatrixnormal 类——是否有一种简单的方法可以采用其中任何一个来对由均值向量和协方差矩阵定义的多元正态分布进行采样?

解决方法

根据@robert-dodier 的建议,MatrixNormal 分布成为 p = 1 处的多元正态分布。这比有原生 Multinormal 分发类更冗长,但不多:

using MathNet.Numerics.LinearAlgebra;

Vector Sample<T>(System.Random random,Vector<T> mean,Matrix<T> cov)
{
    return MatrixNormal.Sample(
        random,mean.ToColumnMatrix(),cov,Matrix<T>.Build.DenseIdentity(1)).Column(0);

}

但是,只允许正定协方差矩阵,因为分布执行 Cholesky 分解。

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。