三次插值ClampedB样条曲线Matlab代码实现_matlab三次样条函数原始实现-CSDN博客

本文链接：https://blog.csdn.net/tao123123q/article/details/125609090

本文介绍了如何在Matlab中实现三次Clamped B样条曲线的插值，特别关注了边缘型值点对控制点系数的计算。在实际编程过程中，发现教程中的系数与计算结果存在差异，提出了修正后的系数公式，并提供了相应的代码实现。

摘要由CSDN通过智能技术生成

三次插值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个
%根据非纽结边界条件 第一个点与第二个点的三次求导相同 最后一个点与倒数第二个点三次求导相同