Рассчитайте гладкость кубического сплайна SciPy

Question

Я пытаюсь вычислить гладкость кубических сплайнов, сгенерированных с использованием scipy.interpolate.CubicSpline. Я хочу рассчитать гладкость, взяв интеграл от квадрата второй производной:

SciPy CubicSpline принимает integrate и derivativeметоды, но их недостаточно для расчета выше. Что я могу легко сделать, это .derivative(2).integrate(0,20)**2, но это не совсем то, что мне нужно (я беру квадрат после , вычисляется интеграл).

Есть ли простой способвыполнить расчет выше? Я что-то упустил глупо?

Answer 1

Вы почти у цели: вам нужно приготовить отзыв для второго квадрата производной:

>>> spl = CubicSpline(x, y)
>>> der = spl.derivative(2)
>>> from scipy.integrate import quad
>>> quad(lambda x: der(x)**2, 0, 1)     # <--- here