该错误代码为% 给定控制点数组P P = [0 0; 1 2; 3 3; 4 1; 6 2]; % 构造节点向量 k = 4; % B样条的次数 n = size(P,1) - 1; % 控制点的个数 t = [zeros(1,k), linspace(0,1,n-k+2), ones(1,k)]; % 计算B样条曲线的系数矩阵和节点向量 c = spline(t, [zeros(k,2); P; zeros(k,2)]); % 求出每一段的起始点、终止点和方程 for i = 1 : n % 起始点和终止点 t1 = t(i+k); t2 = t(i+k+1); P1 = c(i,:); P2 = c(i+1,:); disp(['第', num2str(i), '段：起始点(', num2str(P1(1)), ',', num2str(P1(2)), ')，终止点(', num2str(P2(1)), ',', num2str(P2(2)), ')']); % 方程 syms x y; f = c(i,1) * (1 - (x - t1)/(t2 - t1))^3 + 3 * c(i+1,1) * ((x - t1)/(t2 - t1) - 1) * (1 - (x - t1)/(t2 - t1))^2 + 3 * c(i+2,1) * ((x - t1)/(t2 - t1))^2 * (1 - (x - t1)/(t2 - t1)) + c(i+3,1) * ((x - t1)/(t2 - t1))^3 - y; disp(['第', num2str(i), '段方程：y = ', char(simplify(solve(f, y)))]); end % 绘制B样条曲线 ts = linspace(0,1,100); xs = ppval(c(:,1), ts); ys = ppval(c(:,2), ts); plot(xs, ys, P(:,1), P(:,2), 'o');

这个错误是因为在调用 spline 函数时，输入的 x 和 y 数组长度不一致导致的。在你的代码中，x 数组是节点向量 t，y 数组是控制点数组 P，它们的长度不相等。为了解决这个问题，你可以将 P 数组修改为一个列向量，然后将其转置，这样得到的 y 数组长度就与 x 数组相等了。具体地，可以将 spline 函数的调用改为：

c = spline(t, [zeros(k,2); P(:).'; zeros(k,2)]);

这里使用了 P(:).' 将 P 数组转化为一个列向量并转置，然后再将其与两个零向量拼接起来，得到一个与 t 长度相等的数组。