matlab中nurbs曲线拟合钢轨廓形代码

以下是MATLAB代码示例，用于拟合NURBS曲线，以逼近钢轨轮廓形状：

% 钢轨坐标数据
x = [...];
y = [...];

% 定义控制点权重
w = ones(size(x));

% 定义控制点
n = 10; % 控制点数量
P = zeros(n, 3); % 控制点矩阵（每行为一个控制点的坐标）
P(:,1) = linspace(min(x), max(x), n); % 在x方向上均匀分布
P(:,2) = spline(x, y, P(:,1)); % 使用样条插值计算控制点坐标
P(:,3) = ones(n, 1); % 控制点权重

% 定义节点向量
m = n+2; % 节点向量长度
k = linspace(0, 1, m); % 在0到1之间均匀分布

% 定义次数
p = 3; % 曲线次数

% 计算NURBS曲线
N = spcol(k, p, linspace(min(x), max(x), 1000)); % 计算节点矩阵
C = nrbmak(P, w, k); % 构造NURBS曲线
curve = nrbdeval(C, N); % 计算NURBS曲线上的点

% 绘制曲线和原始数据
figure;
plot(curve(1,:), curve(2,:), 'LineWidth', 2);
hold on;
plot(x, y, '.', 'MarkerSize', 10);
axis equal;

请注意，此示例假定您已经有钢轨坐标数据存储在`x`和`y`向量中。为了使用此代码，您需要将`x`和`y`替换为您的实际数据，并相应地更改控制点和节点向量的数量和分布。