计算逆 CDF 的最有效和数值稳定的算法是什么的概率函数,假设 PDF和 CDF是已知的,但逆 CDF 不是?
显然,这与寻找非线性函数的根相同, 和对于给定的. 我们知道/假设是:
- 单调不减(即CDF);
- 至少可微分次(它的一阶导数是);
- 连续的(我们假设不是增量函数)。
- 我对 PDF 的形状是高斯分布乘以任意次数多项式的通用混合(在这种情况下,可以解析计算 CDF)的情况特别感兴趣。
有几种标准方法可用于一般非线性函数的求根(例如参见数值配方的第 9 章)。我将使用的方法是 Brent 的方法,最后可能需要几个Newton-Raphson步骤来细化根(如果有的话);但我想知道根据上述假设是否有更好的方法。
另外,请注意,我需要计算大量的逆 CDF() 此类中不同的 CDF;查找表或预计算是不可行的。