c++ - 使用 Boost 的 pdf 功能

Question

愚蠢的问题，但我花了太多时间在网上寻找答案，但没有成功。我有一个 boost::random::gamma_distribution 对象和一个浮点值，我想为其计算 pdf。

我究竟应该包含哪些 Boost 模块以及如何调用计算 gamma 的 pdf 的函数？

谢谢

Answer 1

我已经查看了 random/gamma_distribution.hpp 并且没有返回 pdf 的方法，因此 gamma_distribution 的实例对您没有帮助。然而，boost::math::gamma_distribution提供实现说明和公式(底部的表格)以使用库函数 gamma_p_derivative 定义 pdf。

现在你可以自己拼一个pdf函数了:

#include <boost/math/special_functions/gamma.hpp>

// Makes sense for k, theta, x greater than 0.
double gamma_pdf(double k, double theta, double x) {
    return boost::math::gamma_p_derivative(k, x / theta) / theta;
}

基本上就是这样。由于 gamma.hpp 包含所需的定义，因此您无需在编译期间链接任何其他库。