4-3.m

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

pcode HJ;
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=3;

n=length(x);    

u=HJ(n,k,x,y,z);
x1=deBoor(n,k,x,u);
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