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 平面のデータポイント間で同じ弧の長さを取得したい)、次の解決策にたどり着きました。 x(t) と y(t) を持つパラメトリック関数の、ここで t は増加するパラメーター (単純に x または y 配列のインデックス、または弧の長さ) です。次に、x(t) と y(t) の補間を別々に行うことができます。この方法では、y(x) に直接アクセスすることはできませんが、t 値をスキャンして、探している補間された x 値と y 値の組み合わせを取得できます。おそらく、ある種の 2D 補間でも同じことができます。

于 2016-02-23T10:00:50.680 に答える