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