import numpy as np
import matplotlib.pyplot as plt
import multipolyfit as mpf
xs = np.array([[1,2,3,4,5], [1,-1,1,-1,1]]).T
y = np.array([2,2,4,4,6])
model = mpf.multipolyfit(xs, y, 2, model_out = True)
yhat = np.asarray(map(model, xs[:,0], xs[:,1]))
map(model, [1.5],[1])
X = np.arange(-5, 5, 0.25)
Y = np.arange(-5, 5, 0.25)
X, Y = np.meshgrid(X, Y)
Z = np.asarray(map(model, X.reshape(-1),Y.reshape(-1))).reshape(Z.shape)
fig = plt.figure()
ax = fig.gca(projection='3d')
surf = ax.plot_surface(X, Y, Z, rstride=1, cstride=1, cmap=cm.coolwarm,
linewidth=0, antialiased=False)
plt.show()
Inspiraci hledej zde