B样条插值的matlab程序

function B_Spline_3u（P_xP_y）
%3次均匀B样条差值函数
%P_x：控制点的x坐标；P_y:控制点的y坐标；k:次幂
P_size=size（P_x）;
P_num=P_size（2）;
xx=zeros（1（P_num-3）*11）;
yy=xx;
for i=1:（P_num-3）
    k=1;
    for u=0:0.1:1         
        xx（k+（i-1）*11）=P_x（i）*（1/6）*（1-u）^3 + P_x（i+1）*（1/6）*