三维曲线拟合代码

以下是一个基于Python的三维曲线拟合代码示例，使用了scipy库中的interpolate模块和numpy库：

import numpy as np
from scipy import interpolate

# 构造数据
x = np.linspace(0, 1, 20)
y = np.linspace(0, 1, 20)
z = np.sin(x * y)

# 三维插值
interp_func = interpolate.interp2d(x, y, z, kind='cubic')

# 构造新的网格
x_new = np.linspace(0, 1, 100)
y_new = np.linspace(0, 1, 100)

# 生成新的插值结果
z_new = interp_func(x_new, y_new)

# 绘制插值结果
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

fig = plt.figure()
ax = Axes3D(fig)
X, Y = np.meshgrid(x_new, y_new)
ax.plot_surface(X, Y, z_new, cmap='rainbow')
plt.show()

这个代码示例中，我们首先构造了一个由20个点组成的三维曲线，然后使用`interp2d`函数将它进行了三维插值，得到了一个新的网格和插值结果。最后使用`plot_surface`函数将插值结果绘制出来，以便于可视化观察。