4-2.m

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

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

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

n=length(x)-1;
u=LJD2(n+1,2,x,y,z);       %%调用自定义累加弦长算法函数确定节点矢量 
Dx=JCHZH(n,x,u);           %%调用自定义地基于基运算的插值函数确定控制顶点
Dy=JCHZH(n,y,u);
Dz=zeros(1,n+2);           %%便于调用函数

n1=length(Dx);               %%使用deBoor求值算法求出样条曲线
u1=HJ(n1,2,Dx,Dy,Dz);
x1=deBoor(n1,2,Dx,u1);
y1=deBoor(n1,2,Dy,u1); 

plot(x1,y1,'LineWidth',2)      %%画图
hold on
plot(x,y,'m:o')              %%型值点比较
plot(Dx,Dy,'k-.s')           %%所求控制顶点比较
xlabel('x','FontSize',16)
ylabel('y','FontSize',16)
title('二次B样条插值曲线','FontSize',12,'FontName','宋体')
hold off