如何实现兼具插值与逼近的全局平滑三次B样条曲线拟合？

要实现同时插值部分点、逼近其余点的全局平滑三次B样条曲线，核心思路是给两类点分配差异化权重，通过构建加权线性方程组求解控制顶点，强制插值点被曲线严格经过，同时让曲线尽可能贴合逼近点，保持B样条的全局平滑性。

核心步骤 #

1. 分类数据点与参数化 #

首先将输入点明确分为两类：

• 插值点集：( {P_0^I, P_1^I, ..., P_k^I} )，要求曲线必须严格经过这些点
• 逼近点集：( {P_0^A, P_1^A, ..., P_l^A} )，仅要求曲线尽可能接近这些点

为每个点分配对应的参数值 ( u )，推荐使用弦长参数化（更符合点分布的几何特性）：

• 计算所有点的累积弦长 ( L_i = \sum_{t=0}^{i-1} ||P_t - P_{t+1}|| )
• 归一化后得到参数 ( u_i = L_i / L_{total} )，其中 ( L_{total} ) 是总弦长

2. 构造节点向量 ( U ) #

对于三次B样条（次数 ( d=3 )），节点向量长度为 ( \text{控制点数} + d + 1 )。假设期望控制点数为 ( m+1 )，则节点向量 ( U = {u_0, u_1, ..., u_{m+4}} )：

• 若采用 clamped B样条（首尾端点固定），首尾节点重复 ( d+1=4 ) 次，即 ( u_0=u_1=u_2=u_3=0 )，( u_{m+1}=u_{m+2}=u_{m+3}=u_{m+4}=1 )
• 中间节点可根据插值点的参数值均匀分布，或直接使用插值点的参数值作为节点（提升插值约束的稳定性）

3. 构建加权线性方程组 #

对每个点生成对应的方程，并为两类点分配差异化权重：

• 插值点：分配极大权重 ( W_I )（如 ( 10^6 )，确保方程被严格满足）
• 逼近点：分配正常权重 ( W_A=1 )（仅作为拟合约束）

每个点对应的方程形式为：
[ W \cdot \sum_{t=0}^m N_t(u) \cdot Q_t = W \cdot P ]
其中 ( N_t(u) ) 是第 ( t ) 个三次B样条基函数在参数 ( u ) 处的值，( Q_t ) 是待求解的控制顶点。

将所有方程整理为矩阵形式 ( A \cdot Q = B )：

• 矩阵 ( A ) 的每一行：( [W \cdot N_0(u), W \cdot N_1(u), ..., W \cdot N_m(u)] )
• 向量 ( B ) 的每一行：( W \cdot P )

4. 求解控制顶点 #

根据方程数量与控制点数的关系选择求解方式：

• 若方程数 > 控制点数：用加权最小二乘法求解，既满足插值约束，又最小化逼近点的拟合误差
• 若方程数 = 控制点数：用高斯消元直接求解，此时曲线严格经过所有插值点，同时满足逼近点的约束

关键细节说明 #

• 权重选择：插值点的权重无需无穷大，只需比逼近点权重大3-6个数量级，既能保证插值约束，又避免数值不稳定
• 基函数计算：用德布尔-考克斯（De Boor-Cox）递归公式计算B样条基函数值，确保计算效率与准确性
• 平滑性控制：三次B样条本身具有 ( C^2 ) 全局平滑性，无需额外调整，节点向量的构造会直接影响局部平滑度

内容的提问来源于stack exchange，提问作者Matt-Heun Hong