%
%************三次均匀B样条基函数************
function y=Nfun(i,t)
%参数i,
y=0;
for j=03-i)
y=y+(-1)^j*nchoosek(4,j)*(t+3-i-j)^3;
end
y=y/6;
%************K次非均匀B样条基函数************
function result = Bbase(i,k,u,t)
%第i段k次B样条基,Deboor递推递归算法
%t为变量,u(i)<=t<u(i+1),k=0时result=1;
if k==0
if u(i)<=t && t<u(i+1) %注意1=u(i)<=t<u(i+1)=1时的情况,这里要用t<=u(i+1);
result=1;
return;
3-i)