我正在尝试使用scipy中的几种方法从一组x(分位数)和y(值)样本中插入quantile function(逆CDF)。 由于它是分位数函数,因此值有时会重复。例如,CDF最终在1处变平,因此对于几个递增的y值,x = 1重复:
x = [0, 0.19026078648166053, 0.5364188373245662, 0.9627927389184123, 0.9997059472175255, 0.9997059472175255, 0.9999999999999999, 0.9999999999999999, 0.9999999999999999]
y = [0, 468, 1171, 4918, 10072, 20066, 29982, 45207, 59964]
似乎为函数构建了一些插值方法,并且不喜欢重复x。有些人甚至认为重复x是衍生物(e.g. Krogh)。
我知道如何解决这个问题?
答案 0 :(得分:2)
我刚遇到同样的问题(我有一条重复x值的曲线,想要插入它以在x-y平面的数据点之间获得相同的弧长),并找到了适用于我的以下解决方案: 将x-y函数视为具有x(t)和y(t)的参数函数,其中t是一些递增参数(简单地是x或y数组的索引,或弧长)。 然后,您可以分别进行x(t)和y(t)的插值。 这样您就无法直接访问y(x),但您可以扫描t值以获得您正在寻找的内插x和y值的组合。 在某种2D插值中也可以做同样的事情。