三次插值ClampedB样条曲线Matlab代码实现
本文代码根据 插值样条曲线原理实现。
在代码实现过程中与原博客中矩阵计算有误差,主要体现在边缘型值点对控制点的系数求解上。教程中系数求解为0,1/6,4/6,1/6,0…,代码与本人手算求解为0,1/4,7/12,1/6,0…
计算繁杂,如有有误望指正。
结果如下:
代码如下:
clc;clear;close;
%输入型值点 l+1个
points=[[4,4];[5,6];[7,7];[8,8.5];[9,9]];
plot(points(:,1),points(:,2),'or','MarkerSize',10);
hold on
%三次样条曲线 k=m-n-1=3
%准均匀B样条曲线时间节点会比型值点多6个 m+1=l+1+6
tArr=getTimeArray(points);
%根据准均匀节点与插值节点计算控制节点 n+1=m-3=l+1+2个
controlPoint=getControlPoint(tArr,points);
%根据控制点与时间数组进行绘制插值样条曲线
for times=0:0.002:1
b=getB(times,tArr);
p=(b.')*controlPoint;
plot(p(1,1),p(1,2),'.b');
hold on
end
%获取准均匀时间节点数组
function tArr=getTimeArray(points)
num=size(points,1);
tArr=zeros(1,num+6);
for i=1:1:num+6
if i<=4
tArr(1,i)=0;
end
if i>=num+3
tArr(1,i)=1;
end
if i>4 && i<num+3
tArr(1,i)=1/(num-1)*(i-4);
end
end
end
%%根据时间节点与型值点获取控制点 l+1个
%根据非纽结边界条件 第一个点与第二个点的三次求导相同 最后一个点与倒数第二个点三次求导相同