三次B样条曲线反求控制点，在matlab里编程

以下是一个示例程序，可以通过三次B样条曲线的控制点来反求曲线上的点：

% 定义三次B样条曲线的控制点
P = [0 0; 1 1; 2 -1; 3 0; 4 2; 5 1];

% 计算三次B样条曲线的参数向量
m = size(P, 1) - 1;
t = linspace(0, 1, m + 3);
t = t(2:end-1);

% 定义要求解的曲线上的点
t0 = 0.5;

% 计算基函数
N = zeros(m+3, 4);
for i = 1:m+3
    N(i, :) = bSplineBasis(i, 3, t, t0);
end

% 计算曲线上的点
x = N * P(:,1);
y = N * P(:,2);

% 绘制曲线和控制点
plot(x, y, 'b', P(:,1), P(:,2), 'ro');
axis equal;

function N = bSplineBasis(i, k, t, x)
% 计算三次B样条基函数
if k == 1
    if t(i) <= x && x < t(i+1)
        N = 1;
    else
        N = 0;
    end
else
    N = ((x-t(i))/(t(i+k-1)-t(i))) * bSplineBasis(i, k-1, t, x) + ...
        ((t(i+k)-x)/(t(i+k)-t(i+1))) * bSplineBasis(i+1, k-1, t, x);
end
end

此程序输出一个三次B样条曲线，并在控制点上显示红色圆圈。要求解曲线上的点，请更改`t0`的值。