3

我正在尝试使用 scipy 中的几种方法从一组 x（分位数）和 y（值）样本中插入分位数函数（逆 CDF）。由于它是一个分位数函数，因此这些值有时会重复。例如，CDF 最终在 1 处变平，因此 x=1 重复几个递增的 y 值：

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 是导数（例如 Krogh）。

知道如何解决这个问题吗？

4

1 回答 1

1

我刚刚遇到了同样的问题（我有一条重复 x 值的曲线，并希望对其进行插值以在 xy 平面中的数据点之间获得相同的弧长），并得出以下对我有用的解决方案：将您的 xy 函数视为具有 x(t) 和 y(t) 的参数函数，其中 t 是一些增加的参数（只是 x 或 y 数组的索引，或弧长）。然后您可以分别对 x(t) 和 y(t) 进行插值。这样您就不能直接访问 y(x)，但您可以扫描 t 值以获得您正在寻找的插值 x 和 y 值的组合。同样的可能也可以在某种 2D 插值中完成。

于 2016-02-23T10:00:50.680 回答