nurbs曲线拟合钢轨廓形代码

以下是一个使用NURBS曲线拟合钢轨廓形的示例代码，代码使用Python语言编写，使用了NumPy和matplotlib库：

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

# 读取钢轨边缘点的数据
data = np.loadtxt('rail_data.txt')

# 将钢轨边缘点数据转换为NURBS曲线所需的格式
x, y = data[:, 0], data[:, 1]
tck, u = splprep([x, y], s=0)

# 生成NURBS曲线的参数化点
u_new = np.linspace(u.min(), u.max(), 1000)
x_new, y_new = splev(u_new, tck)

# 绘制原始数据点和拟合曲线
plt.plot(x, y, 'o', label='原始数据点')
plt.plot(x_new, y_new, label='NURBS曲线')
plt.legend()
plt.show()

在代码中，我们首先从文件中读取钢轨边缘点的数据，然后使用`splprep`函数将数据转换为NURBS曲线所需的格式。`splprep`函数返回一个三元组`tck`，其中包含了NURBS曲线的控制点、节点向量和次数。我们使用`tck`和生成的参数化点`u_new`，使用`splev`函数生成NURBS曲线的坐标点`x_new`和`y_new`。最后，我们使用`matplotlib`库绘制原始数据点和NURBS曲线。