4-8.m

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

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

x0=[0 2 2 4 6 8 10];          %%赋初始值
y0=[0 2 -1 0 -2 -3 0];
z0=[0 -2 -1 0 1 2 1];
k=3;

n0=length(x0);               %%反算插值曲线控制顶点
u0=LJD(n0,k,x0,y0,z0);
x1=CHZH(n0,k,x0,u0);
y1=CHZH(n0,k,y0,u0);
z1=CHZH(n0,k,z0,u0);

n1=length(x1);               %%使用deBoor求值算法求出样条曲线
u1=HJ(n1,k,x1,y1,z1);
x=deBoor(n1,k,x1,u1);
y=deBoor(n1,k,y1,u1);
z=deBoor(n1,k,z1,u1);

plot3(x,y,z,'LineWidth',2)      %%画图
hold on
plot3(x0,y0,z0,'m:o')            %%型值点比较
plot3(x1,y1,z1,'k-.s')           %%所求控制点比较
title('三次B样条插值曲线','FontSize',12,'FontName','宋体')
xlabel('x','FontSize',12)
ylabel('y','FontSize',12)
zlabel('z','FontSize',12)
axis vis3d
hold off