B样条插值是通过反解控制点使曲线精确经过给定数据点，需解线性方程组；而普通B样条仅由控制点和节点定义分段多项式曲线，不保证过控制点。
什么是B样条插值，和普通样条有什么区别
B样条插值不是直接让曲线穿过所有控制点（那是插值型B样条的常见误解），而是通过一组控制点 control_points 和节点向量 knots 定义一条分段多项式曲线，它不一定经过控制点；若要实现“插值”，即让曲线精确经过给定数据点 data_points，必须反解控制点——这叫插值型B样条拟合，本质是解线性方程组。
关键判断：标准 std::vector + 手写 De Boor 算法只能做求值（evaluating），不能自动插值。你要的是“给定 P₀…Pₙ，构造一条 B 样条使其满足 C(tᵢ) = Pᵢ”，这就绕不开矩阵求解。
如何构造插值型B样条：节点、次数、参数化三要素缺一不可
插值成败高度依赖三个输入是否合理：
-
degree：通常取 3（三次B样条），太低（如 1）欠光滑，太高（≥5）易振荡且病态 -
parameterization：不能用等距tᵢ = i，推荐 chord length（弦长参数化）：t[0] = 0，t[i] = t[i-1] + distance(data_points[i], data_points[i-1]) -
knots：需满足插值条件数量匹配。对n+1个数据点、次数p，需n+1个控制点 → 节点数应为(n+1) + p + 1；常用开节点（open uniform）：前p+1个为 0，后p+1个为 1，中间均匀分布
常见错误现象：matrix is singular 或插值结果严重偏离——八成是节点向量没按开形式构造，或参数 tᵢ 重复导致基函数矩阵列线性相关。
立即学习“C++免费学习笔记（深入）”；
用Eigen解控制点：核心是构造基函数矩阵 A
对每个数据点 i，计算其参数 tᵢ 处所有非零基函数值 N_{j,p}(tᵢ)，填入矩阵 A[i][j]（j 是控制点索引）。最终解 A × X = P，其中 X 是待求控制点坐标（x/y/z 分开解）。
// 示例：计算单个 N_{j,p}(t) 的伪代码（De Boor 递推）
double basis(int j, int p, double t, const std::vector<double>& knots) {
if (p == 0) return (knots[j] <= t && t < knots[j+1]) ? 1.0 : 0.0;
double left = (t - knots[j]) / (knots[j+p] - knots[j]);
double right = (knots[j+p+1] - t) / (knots[j+p+1] - knots[j+1]);
return left * basis(j, p-1, t, knots) + right * basis(j+1, p-1, t, knots);
}
实操建议：
- 不要手算整个矩阵：循环
i（数据点）、j（控制点），调用basis(j, p, t[i], knots)填A(i,j) - 使用
Eigen::MatrixXd A(n+1, n+1)，Eigen::VectorXd Px(n+1)存 x 坐标，再用A.colPivHouseholderQr().solve(Px) - 若
A条件数 > 1e8，改用FullPivLu或加小正则项（A += 1e-10 * Eigen::MatrixXd::Identity()）
为什么你跑出来的曲线不经过数据点
最常被忽略的点：De Boor 算法求值时，t 必须落在有效区间内，且节点向量必须与插值时一致。
- 插值阶段用了开节点
[0,0,0,0,0.25,0.5,0.75,1,1,1,1]（p=3, n=6），那求值时t必须 ∈ [0,1]，且边界行为要处理：t=0→ 曲线起点 = 第一个控制点？错，只在均匀/特殊节点下成立；实际应严格按 De Boor 在t=0处计算，它可能不等于任一控制点 - 检查
t[i]是否全部 ∈[knots[p], knots[n+1]]（有效定义域），否则基函数返回 0 导致求值失效 - 可视化调试：打印
C(t[0]),C(t[1])… 看是否 ≈data_points[0],data_points[1]；不等就回头查A矩阵第 0 行是否全零或异常小
真正麻烦的从来不是公式，而是节点向量的手动构造和参数化的隐含假设。
