4-5.m

de Boor 求值算法求作二次和三次B-样条曲线非插值

pcode HJ;                        %%生成p码文件，加快执行速度
pcode deBoor;

x=[0 2 2 4 6 8 10];               %%赋初值
y=[0 2 -1 0 -2 -3 0];
z=[0 -2 -1 0 1 2 1];
k=2;

n=length(x);

u=HJ(n,k,x,y,z);                  %%使用Hartley-Judd法确定节点矢量
x1=deBoor(n,k,x,u);               %%使用deBoor法求曲线
y1=deBoor(n,k,y,u);
z1=deBoor(n,k,z,u);
    
plot3(x1,y1,z1,'LineWidth',2)     %%画图
hold on
plot3(x,y,z,'m:o')                %%原始点比较
title('二次B样条曲线','FontSize',12,'FontName','宋体')
xlabel('x','FontSize',12)
ylabel('y','FontSize',12)
zlabel('z','FontSize',12)
axis vis3d;
hold off