Python Least-Squares Natural Splines

Question

가중치가 최소 인 사각형을 최소화하는 자연스러운 스플라인에 맞는 숫자 패키지를 찾으려고합니다.

부 자연스러운 스플라인을 원한다면 scipy에 패키지가 있습니다.

import numpy as np 
import matplotlib.pyplot as plt 
from scipy import interpolate, randn 

x = np.arange(0,5,1.0/6) 
xs = np.arange(0,5,1.0/500) 

y = np.sin(x+1) + .2*np.random.rand(len(x)) -.1 

knots = np.array([1,2,3,4]) 
tck = interpolate.splrep(x,y,s=0,k=3,t=knots,task=-1) 
ys = interpolate.splev(xs,tck,der=0) 

plt.figure() 
plt.plot(xs,ys,x,y,'x') 

Answer 1

this page에서이 tar file의 내부 spline.py 파일은 기본적으로 자연 스플라인 맞춤을한다. this page에도 대부분 원하는대로 소유권을 주장하는 코드가 있습니다. pyD3D package도 pyDataUtils 모듈에 natural spline function이 있습니다. 이 마지막 하나가 나에게 가장 유망 해 보인다. 그러나, 그것은 자신의 노트를 설정할 수있는 옵션을 가지고 있지 않습니다. 어쩌면 소스를 보면이를 수정하는 방법을 찾을 수 있습니다.

또한 Scipy 메일 목록에 s = 0.0을 사용하면 메시지 작성자에 따라 위 절차를 사용하여 스플라인을 자연 스럽 게 만든다는 메시지가 나타납니다 (this message). 나는 this splmake 함수에 자연스러운 스플라인을 적용 할 수있는 옵션이 있지만 소스를 보면 아직 구현되지 않은 것으로 나타났습니다.