我正在尝试使用LSQUnivariateSpline拟合曲线。
from scipy.interpolate import LSQUnivariateSpline, UnivariateSpline
import matplotlib.pyplot as plt
x = np.linspace(-3, 3, 50)
y = np.exp(-x**2) + 0.1 * np.random.randn(50)
t = [-2, -1, 0, 1, 2]
spl = LSQUnivariateSpline(x, y, t, k=4 )
xs = np.linspace(-3, 3, 1000)
plt.plot(x, y, 'ro', ms=5)
plt.plot(xs, spl(xs), 'g-', lw=3)
plt.show()
此代码产生下面的图形。
我要声明dy / dx = 0或d2x / dy = 0之类的边界条件。是否有方法或其他函数来拟合具有边界条件的曲线?