贝塞尔曲线插值实为误解：标准贝塞尔曲线仅端点插值（必过首尾控制点），不强制经过中间控制点；其本质是通过控制点定义光滑形状的参数曲线，而非传统意义的插值算法。
什么是贝塞尔曲线插值？别被名字骗了
贝塞尔曲线本身不是插值算法——它不强制经过控制点，只用控制点定义形状。如果你想要“经过给定点”的曲线，实际需要的是**插值型贝塞尔样条（如三次B样条插值）或分段贝塞尔拼接**，但标准贝塞尔曲线（BezierCurve）本身不插值端点以外的点。
常见误解：把 std::vector<point></point> 丢进一个叫“贝塞尔插值”的函数，期待它穿过所有点。这不会发生。真正的贝塞尔曲线仅保证起点和终点在首尾控制点上，中间点只是“引导方向”。
怎么用C++实现标准贝塞尔曲线求值（最常用场景）
绝大多数工程需求是：给定一组控制点，计算曲线上某参数 t ∈ [0,1] 对应的坐标。这是德卡斯特里奥（de Casteljau）算法或伯恩斯坦多项式直接求值。
推荐用迭代版 de Casteljao，数值更稳定、易理解、无阶乘/组合数溢出风险：
立即学习“C++免费学习笔记（深入）”；
Point bezierEvaluate(const std::vector<Point>& points, double t) {
auto pts = points;
while (pts.size() > 1) {
for (size_t i = 0; i < pts.size() - 1; ++i) {
pts[i].x = (1 - t) * pts[i].x + t * pts[i + 1].x;
pts[i].y = (1 - t) * pts[i].y + t * pts[i + 1].y;
}
pts.pop_back();
}
return pts[0];
}-
t = 0→ 返回第一个控制点；t = 1→ 返回最后一个控制点 - 控制点数量决定阶数：4个点 = 三次贝塞尔（最常用）
- 避免对
t做等距采样后直接连直线——曲率大时视觉断层明显；建议用自适应步长或固定高密度（如t += 0.005） - 若需高性能批量求值（如GPU上传前预计算），可展开为显式三次多项式：
B(t) = a·t³ + b·t² + c·t + d，但需手算系数，且失去通用性
如果真要“插值已知点”，该选什么替代方案？
想让曲线严格穿过输入点集 {P₀, P₁, ..., Pₙ}，贝塞尔曲线不是第一选择。更合理的技术路径是：
- 用三次B样条（
BSpline）并设节点向量为“非均匀周期性”+端点插值条件（clamped）→ 曲线经过首尾点，内部近似 - 构造分段三次贝塞尔：每段 4 个控制点，强制满足 C¹ 连续（位置+一阶导连续），通过解线性方程组反推中间控制点 —— 这才是所谓“贝塞尔插值”的实际含义
- 更简单粗暴：用
std::spline（C++26 中的<splines>）或第三方库如ceres::Spline/libigl::fit_bspline - 警告：
std::vector<Point>直接喂给任意“贝塞尔插值”函数却没说明其是否做控制点反解，大概率只是画了一条无关的贝塞尔曲线
容易踩的坑：坐标类型、浮点精度与内存布局
写完发现曲线抖动、端点偏移、性能差？先检查这几个硬伤：
- 用
float存控制点坐标？在 t 接近 0 或 1 时，(1-t)和t的浮点误差会放大，尤其当点坐标量级大（如屏幕坐标 > 10000）。改用double或定点缩放 - 控制点存成
std::vector<std::pair<double,double>>？不如自定义struct Point { double x,y; }，避免结构体拆包开销和对齐问题 - 每次调用
bezierEvaluate都拷贝整个points？对动画帧率敏感场景，应传 const 引用 + 栈上临时数组（std::array<Point,16>）代替std::vector - 没验证
t范围！t < 0或t > 1会导致外推，结果完全不可控。加 assert 或 clamp：t = std::clamp(t, 0.0, 1.0)
真正难的从来不是公式，而是搞清你要的到底是“画一条光滑过渡的控制线”，还是“让曲线精确穿过测量数据点”——这两个目标对应的数学工具和代码结构完全不同。
