對數刻度上的插值

Question

我想插入一些數據，並在日志范圍內繪制結果（ pyplot.loglog ）。 問題是得到的插值看起來很奇怪，並且在對數刻度上繪制時顯示不連續性。 插入日志縮放數據的最佳方法是什么？

pyplot.loglog(x, y, '+')
pyplot.hold(True)
s = scipy.interpolate.InterpolatedUnivariateSpline(x, y)
xs = numpy.logspace(numpy.log10(numpy.min(x)), numpy.log10(numpy.max(x)))
pyplot.loglog(xs, s(xs)) # This looks very strange because of the log scale!

實際上，我通過插入數據的日志來成功實現它，但我想知道是否有更簡單的方法來實現相同的結果？

pyplot.loglog(x, y, '+')
pyplot.hold(True)
s = scipy.interpolate.InterpolatedUnivariateSpline(numpy.log10(x), numpy.log10(y))
xs = numpy.logspace(numpy.log10(numpy.min(x)), numpy.log10(numpy.max(x)))
pyplot.loglog(xs, numpy.power(10, s(numpy.log10(xs)))

Answer 1

看起來首先采用數據的對數，然后擬合是執行此操作的常規方法。 請參閱擬合冪律分布 。