根据给定坐标点进行B样条全局插值的方法

B样条曲线全局插值的目标是构造一条B样条曲线，使其精确地通过一组给定的型值点（数据点）。与局部插值（如分段插值）不同，全局插值考虑所有点来确定整条曲线的控制顶点。

1. B样条曲线基础

一条k次（degree = k）的B样条曲线定义为： C(u)=∑i=0n​Ni,k​(u)Pi​ 其中：

C(u) 是曲线上的点。

Pi​ 是控制顶点（i=0,1,...,n）。

Ni,k​(u) 是k次B样条基函数。

u 是参数，属于节点向量U={u0​,u1​,...,um​}，其中m=n+k+1。

2. 全局插值步骤

给定d+1 个型值点Q0​,Q1​,...,Qd​，我们希望构造一条k次B样条曲线通过这些点。

步骤1：确定次数k 选择B样条曲线的次数k（通常k=2, 3, 4）。次数越高，曲线越光滑，但计算越复杂。

步骤2：确定控制顶点数n+1 对于通过d+1 个点的插值，控制顶点的数量n+1 必须满足 n+1≥d+1。为了精确插值，通常设置n=d，即控制顶点数等于数据点数。如果 n>d，问题会变成超定或欠定，需要最小二乘等方法。这里我们讨论n=d 的情况。

步骤3：构造节点向量U 节点向量U={u0​,u1​,...,um​} 的长度 m=n+k+1=d+k+1。

• 边界节点：为了使曲线通过端点，通常将前k+1个和后k+1个节点设置为重复的（clamped B-spline）：
  • u0​=u1​=...=uk​=0
  • um−k​=um−k+1​=...=um​=1
• 内部节点：uk+1​,uk+2​,...,um−k−1​ 的值。这通常通过参数化来完成。

步骤4：参数化（Parameterization） 为每个型值点Qj​ 分配一个参数值tj​ (j=0,1,...,d)。常用方法有：

• 均匀参数化：tj​=j/d
• 弦长参数化：t0​=0, tj​=tj−1​+∥Qj​−Qj−1​∥/L，其中L=∑i=1d​∥Qi​−Qi−1​∥ 是总弦长。
• 中心参数化：tj​=(tj−1​+tj+1​)/2，但通常使用弦长法。

步骤5：构建节点向量（基于参数化） 使用参数值tj​ 构造内部节点。一种常用方法是平均法（Averaging Method）： 对于i=1,2,...,n−k：uk+i​=k1​∑j=ii+k−1​tj​ 这个公式将参数值平均分配到内部节点上。

步骤6：建立线性方程组 要求曲线在参数 tj​ 处通过点 Qj​： for  C(tj​)=∑i=0n​Ni,k​(tj​)Pi​=Qj​这可以写成矩阵形式： NP=Q 其中：

• N 是(d+1)×(n+1) 的系数矩阵，元素 Nji​=Ni,k​(tj​)。
• P (n+1)×2 的控制顶点矩阵（假设是2D点）。
• Q 是 (d+1)×2 的型值点矩阵。

由于 n=d，N 是方阵。对于clamped B样条，N 通常是可逆的。

步骤7：求解控制顶点 解线性方程组：P=N−1Q

步骤8：绘制B样条曲线 使用求得的控制顶点 Pi​ 和节点向量U，通过de Boor算法或其他方法绘制B样条曲线。