4-1.m

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

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

x=[0 1 2 3 4 5 6];                   %%附初值
y=[0 1 3 5 3 2 0];
z=[0 0 0 0 0 0 0];                   %%便于调用函数
k=2;

n=length(x);    

u=HJ(n,k,x,y,z);                     %%调用自定义函数确定节点矢量
x1=deBoor(n,k,x,u);                  %%调用自定义deBoor算法函数
y1=deBoor(n,k,y,u);

plot(x1,y1,'LineWidth',2)            %%画图
hold on
plot(x,y,'m:o')                      %%画原始点比较图
xlabel('x','FontSize',16)
ylabel('y','FontSize',16)
title('二次B样条曲线','FontSize',12,'FontName','宋体')
hold off