简介：这份资源面向计算机图形学、CAD开发及数值计算方向的VC++程序员与学习者，提供三次B样条曲线算法的完整工程实现，帮助理解曲线平滑建模与局部控制的核心思路。压缩包共38个文件，约2.29MB，以h头文件与cpp源文件为主体，配合obj、pdb、exe等编译产物，另有ico、bmp、rc等界面资源及dsp、dsw工程文件，构成可直接用VC6打开编译的MFC单文档项目。内容围绕B样条基础理论、基函数构造、Knot向量与参数u、控制点调整、De Boor算法等关键知识点展开，并借助GDI完成曲线绘制与可视化输出，便于读者对照代码理解算法流程。目前已有528人学习下载，适合希望从理论推导过渡到工程落地、掌握曲线插值与交互式界面实现的开发者参考。

1. 三次B样条曲线算法在VC++里到底怎么落地：从控制点到可运行曲线

如果你手头有一组离散控制点，想生成一条光滑、局部可调、不产生龙格现象的曲线，三次B样条基本是绕不开的选择。它不像贝塞尔曲线那样改一个点整条线都跟着动，也不像插值样条那样每个点都必须穿过。三次B样条的核心价值在于：用四个控制点决定一段曲线，段与段之间自动保持C2连续，工程上做轨迹规划、字体轮廓、相机路径、数据平滑都用它。但很多资料只给公式不给代码，或者给一堆数学推导却不说VC++里怎么组织数据结构、怎么分段、边界怎么处理。这篇笔记就按我实际拆解一个VC++实现的思路来写，从控制点输入到曲线点输出，把参数、代码和踩过的坑都摊开。适合有C++基础、需要在MFC或控制台里快速出曲线的人。

2. 三次B样条数学形式与VC++数据结构映射

2.1 从基函数到分段矩阵：为什么工程上更爱用矩阵形式

三次B样条的标准定义是 ( C(u)=\sum_{i=0}^{n} N_{i,3}(u) P_i )，其中 ( N_{i,3} ) 是三次B样条基函数，通常用Cox-de Boor递推得到。但实际写代码时，递推公式每算一个点都要重新计算基函数，效率低且容易在节点区间边界上出问题。更常见的做法是把均匀三次B样条写成矩阵形式，每一段曲线只和四个控制点有关：

[ C_i(t)=\frac{1}{6} \begin{bmatrix} 1 & t & t^2 & t^3 \end{bmatrix} \begin{bmatrix} 1 & 4 & 1 & 0 \ -3 & 0 & 3 & 0 \ 3 & -6 & 3 & 0 \ -1 & 3 & -3 & 1 \end{bmatrix} \begin{bmatrix} P_{i} \ P_{i+1} \ P_{i+2} \ P_{i+3} \end{bmatrix} ]

其中 ( t \in [0,1] )。这个矩阵形式的好处是：每段曲线独立计算，四个控制点确定一段，段与段之间自然C2连续。VC++里只需要把控制点存成std::vector<CPoint2D>，然后按四个一组滑动窗口计算即可。注意这里的均匀B样条默认节点向量是均匀分布的，如果你需要准均匀或非均匀，矩阵会变，但工程上八成场景均匀就够。

2.2 控制点、节点向量与曲线点容器怎么定义

在VC++里我一般定义两个结构体：一个存二维点，一个存曲线计算结果。控制点用double而不是int，因为后续可能要做缩放和变换。节点向量如果均匀，其实不需要显式存储，直接用段索引映射即可。但为了扩展性，还是留一个std::vector<double> knots，均匀情况下可以按i/(n+4)生成。

// Point2D.h struct Point2D { double x; double y; Point2D(double _x = 0, double _y = 0) : x(_x), y(_y) {} }; // BSpline3.h class BSpline3 { public: std::vector<Point2D> ctrlPts; // 控制点 std::vector<Point2D> curvePts; // 采样后的曲线点 int samplesPerSeg; // 每段采样数，默认20 BSpline3() : samplesPerSeg(20) {} void SetControlPoints(const std::vector<Point2D>& pts) { ctrlPts = pts; } void Generate(); // 核心计算 };

这里samplesPerSeg控制每段曲线的离散精度。20 段在普通屏幕上已经看不出折线，如果要做高精度导出可以调到 50 或 100。控制点数量至少 4 个，否则无法构成一段三次B样条。如果少于 4 个，常见做法是重复端点或直接退化为直线，但我不建议在算法层做这种兼容，最好在调用前就保证输入合法。

2.3 均匀三次B样条的段数与参数映射

假设有 ( n ) 个控制点，均匀三次B样条的有效段数是 ( n-3 )。每段对应四个连续控制点：第 ( i ) 段用 ( P_i, P_{i+1}, P_{i+2}, P_{i+3} )。参数 ( t ) 从 0 到 1 遍历，代入矩阵公式得到曲线点。这里有个容易翻车的地方：很多人以为曲线会经过第一个和最后一个控制点，实际上均匀三次B样条不经过端点，曲线起点在 ( P_0P_1P_2P_3 ) 决定的段上，离 ( P_0 ) 还有一段距离。如果你需要曲线经过端点，要么用准均匀B样条（重复端点），要么在两端各补一个虚拟控制点。我一般是在输入控制点前后各复制一个端点，这样曲线会靠近端点但不一定完全重合，视觉上更自然。

3. VC++实现三次B样条：从矩阵乘法到采样输出

3.1 核心生成函数与矩阵系数展开

直接写矩阵乘法在C++里不算复杂，但为了减少循环开销，我习惯把矩阵展开成显式系数。对于均匀三次B样条，每段曲线点的坐标可以写成：

[ x(t) = \frac{1}{6} \left[ (1-t)^3 x_i + (3t^3 - 6t^2 + 4) x_{i+1} + (-3t^3 + 3t^2 + 3t + 1) x_{i+2} + t^3 x_{i+3} \right] ]

这个展开式比矩阵乘法少一层循环，而且系数可以预先算好。下面是对应的VC++实现：

void BSpline3::Generate() { curvePts.clear(); int n = (int)ctrlPts.size(); if (n < 4) return; // 控制点不足，直接返回 int segCount = n - 3; // 有效段数 for (int i = 0; i < segCount; ++i) { const Point2D& P0 = ctrlPts[i]; const Point2D& P1 = ctrlPts[i + 1]; const Point2D& P2 = ctrlPts[i + 2]; const Point2D& P3 = ctrlPts[i + 3]; for (int s = 0; s <= samplesPerSeg; ++s) { double t = (double)s / samplesPerSeg; double t2 = t * t; double t3 = t2 * t; // 三次B样条基函数展开 double b0 = (1 - 3*t + 3*t2 - t3) / 6.0; double b1 = (4 - 6*t2 + 3*t3) / 6.0; double b2 = (1 + 3*t + 3*t2 - 3*t3) / 6.0; double b3 = t3 / 6.0; double x = b0*P0.x + b1*P1.x + b2*P2.x + b3*P3.x; double y = b0*P0.y + b1*P1.y + b2*P2.y + b3*P3.y; curvePts.push_back(Point2D(x, y)); } } }

这段代码里b0到b3就是三次B样条的四个基函数。注意s从 0 到samplesPerSeg包含两端，所以每段会生成samplesPerSeg+1个点，相邻段之间会有一个重复点（上一段的终点和下一段的起点相同）。如果不想重复，可以在s循环里从 1 开始，但那样第一段的起点会丢，所以更常见的做法是保留重复点，后续画线时直接按顺序连，重复点不影响视觉。参数samplesPerSeg越大曲线越平滑，但点也越多，导出DXF或SVG时要注意文件大小。

3.2 控制点输入与边界处理：端点重复与准均匀方案

前面提到均匀三次B样条不经过端点。如果你希望曲线从第一个控制点附近开始、到最后一个控制点附近结束，最简单的办法是在控制点数组前后各加一个重复点。比如原始控制点是P0, P1, P2, P3, P4，变成P0, P0, P1, P2, P3, P4, P4，这样曲线会靠近P0和P4。但注意这样会改变段数，从n-3变成n+1-3 = n-2。下面是一个封装函数：

std::vector<Point2D> BuildClampedCtrlPts(const std::vector<Point2D>& raw) { if (raw.size() < 2) return raw; std::vector<Point2D> out; out.push_back(raw.front()); // 重复起点 out.insert(out.end(), raw.begin(), raw.end()); out.push_back(raw.back()); // 重复终点 return out; }

这样处理之后，曲线起点和终点会分别落在第一个和最后一个控制点上吗？严格说不会完全重合，但距离非常近，视觉上可以接受。如果你需要严格通过端点，那就得用准均匀B样条，节点向量两端重复度设为4。准均匀的矩阵和均匀不同，代码里需要单独处理节点区间。我一般不在基础版本里做准均匀，因为大部分可视化场景均匀已经够用，而且准均匀的边界段计算容易出索引越界。

3.3 在MFC中绘制曲线：CDC与双缓冲

VC++里如果是在MFC对话框或视图里画曲线，直接CDC::Polyline或者MoveTo/LineTo就行。但控制点多、采样密的时候，频繁刷新会闪。我习惯用双缓冲：先在内存DC里画好，再一次性贴到屏幕。下面是一个简化的绘制函数：

void CMyView::DrawBSpline(CDC* pDC) { if (m_spline.curvePts.size() < 2) return; CRect rect; GetClientRect(&rect); CDC memDC; memDC.CreateCompatibleDC(pDC); CBitmap bmp; bmp.CreateCompatibleBitmap(pDC, rect.Width(), rect.Height()); memDC.SelectObject(&bmp); memDC.FillSolidRect(&rect, RGB(255, 255, 255)); // 画控制点 CPen penCtrl(PS_SOLID, 1, RGB(180, 180, 180)); memDC.SelectObject(&penCtrl); for (size_t i = 0; i + 1 < m_spline.ctrlPts.size(); ++i) { memDC.MoveTo((int)m_spline.ctrlPts[i].x, (int)m_spline.ctrlPts[i].y); memDC.LineTo((int)m_spline.ctrlPts[i+1].x, (int)m_spline.ctrlPts[i+1].y); } // 画曲线 CPen penCurve(PS_SOLID, 2, RGB(0, 120, 215)); memDC.SelectObject(&penCurve); memDC.MoveTo((int)m_spline.curvePts[0].x, (int)m_spline.curvePts[0].y); for (size_t i = 1; i < m_spline.curvePts.size(); ++i) { memDC.LineTo((int)m_spline.curvePts[i].x, (int)m_spline.curvePts[i].y); } pDC->BitBlt(0, 0, rect.Width(), rect.Height(), &memDC, 0, 0, SRCCOPY); }

这里控制点用灰色细线连，曲线用蓝色粗线。注意坐标转换：如果控制点坐标是数学坐标（y轴向上），而屏幕坐标y轴向下，需要在输入前做一次翻转，否则曲线会上下颠倒。我一般是在SetControlPoints之前统一把y坐标取反，或者用rect.Height() - y转换。双缓冲的代价是每次重绘都要创建位图，如果曲线点特别多（比如几万个），可以考虑只创建一次位图重复使用。

4. 避坑与排查：三次B样条在VC++里最容易翻车的五个点

4.1 曲线不光滑或出现折线：采样数不够还是控制点顺序错了

现象：生成的曲线看起来一段一段的折线，尤其是曲率大的地方。原因通常有两个：一是samplesPerSeg太小，比如设成5，每段只有6个点，连起来就是折线；二是控制点顺序不是按路径顺序排列的，导致段与段之间方向突变。解决：先把samplesPerSeg调到20以上，如果还不光滑，检查控制点数组是不是有交叉或跳跃。我遇到过有人把控制点按x坐标排序后传入，结果曲线自交，这就是顺序错了。控制点必须按你期望的路径顺序排列，不能排序。

4.2 端点不贴合：均匀B样条不经过首末控制点

现象：曲线起点离第一个控制点很远，终点也离最后一个控制点很远。原因：均匀三次B样条本身不插值端点，这是数学性质，不是bug。解决：如果必须贴合，用前面说的端点重复法，或者改用准均匀B样条。注意端点重复后，曲线只是靠近端点，不是精确通过。要精确通过，得用插值样条或者把端点作为约束反算控制点。工程上大部分场景“靠近”就够了，不必追求精确通过。

4.3 控制点少于4个导致崩溃或空结果

现象：程序在Generate()里访问ctrlPts[i+3]时越界，或者曲线点为空。原因：控制点数量小于4，但代码没有做前置检查。解决：在Generate()开头加if (n < 4) return;，同时在UI层限制用户至少输入4个点。如果业务上确实可能少于4个，那就退化为直线或二次B样条，但那是另一套逻辑。我一般会在SetControlPoints里就做一次检查，少于4个直接清空并返回错误码。

4.4 浮点累积误差导致曲线抖动

现象：曲线在局部出现微小抖动，放大后能看到点不在平滑轨迹上。原因：t用float累加，或者基函数计算时用了低精度类型。解决：所有中间变量用double，t用(double)s / samplesPerSeg直接算，不要用t += step累加。累加会引入舍入误差，段数多了之后误差累积，曲线就会抖。另外，基函数展开式里的系数最好一次性算完，不要分步近似。

4.5 在MFC中刷新闪烁或曲线残影

现象：拖动控制点时曲线闪烁，或者旧曲线没擦掉。原因：没有用双缓冲，或者Invalidate之后没有UpdateWindow。解决：用内存DC双缓冲，如3.3节所示。如果还有残影，检查是不是在OnDraw里重复创建了画笔和位图而没有释放。GDI对象泄漏也会导致绘制异常。我习惯把双缓冲的位图和内存DC作为成员变量，只在窗口大小变化时重建，而不是每次重绘都创建。

5. 进阶：把三次B样条封装成可复用类与参数化技巧

5.1 支持任意维度与导出SVG路径

基础版本只处理二维点，但三次B样条的公式对任意维度都成立。如果你要做三维轨迹，只需要把Point2D换成Point3D，基函数部分完全不用改。更进一步，可以把Generate()写成模板函数，维度作为模板参数。导出方面，曲线点可以直接转成SVG的path指令，用M开头，后面跟一串L，虽然不如贝塞尔指令紧凑，但兼容性好。下面是一个简单的SVG导出片段：

void ExportSVG(const std::vector<Point2D>& pts, const std::string& file) { std::ofstream ofs(file); ofs << "<svg xmlns=\"http://www.w3.org/2000/svg\" width=\"800\" height=\"600\">\n"; ofs << "<path d=\"M " << pts[0].x << " " << pts[0].y; for (size_t i = 1; i < pts.size(); ++i) { ofs << " L " << pts[i].x << " " << pts[i].y; } ofs << "\" stroke=\"blue\" fill=\"none\" stroke-width=\"2\"/>\n"; ofs << "</svg>\n"; }

注意SVG的y轴向下，如果你的控制点是数学坐标，导出前要翻转y。另外点数太多时SVG文件会很大，可以在导出前做一次道格拉斯-普克抽稀，但那是另一个话题了。

5.2 参数化：均匀、弦长与向心

均匀参数化就是每段t等距，实现最简单，但控制点分布不均匀时曲线会出现“过冲”或“扁平”。弦长参数化按控制点之间的距离分配参数，向心参数化按距离的平方根分配，后者在控制点急转弯时更稳定。VC++里实现弦长参数化只需要在生成节点向量时累加距离：

std::vector<double> BuildChordKnots(const std::vector<Point2D>& pts) { std::vector<double> knots; knots.push_back(0.0); double total = 0.0; for (size_t i = 1; i < pts.size(); ++i) { double dx = pts[i].x - pts[i-1].x; double dy = pts[i].y - pts[i-1].y; total += sqrt(dx*dx + dy*dy); knots.push_back(total); } for (auto& k : knots) k /= total; // 归一化到[0,1] return knots; }

然后基函数计算时用节点向量查表，而不是直接用t。这样曲线在控制点密集的地方会更贴合。但注意弦长参数化后，每段的t范围不再是[0,1]，需要做区间映射。我一般只在控制点分布极不均匀时才用，普通场景均匀参数化足够。

5.3 反算控制点：让曲线通过指定点

有时候需求是“曲线必须经过这几个点”，但三次B样条是逼近不是插值。这时候需要反算控制点：给定数据点 ( Q_i ) 和节点向量，解一个三对角线性方程组得到控制点 ( P_i )。VC++里可以用追赶法解三对角方程，代码不长但容易在边界条件上出错。我一般用现成的线性代数库，比如Eigen，但如果你不想引入第三方库，手写追赶法也就三十行。反算之后再用Generate()生成曲线，就能保证曲线经过指定点。这个技巧在轨迹规划里很实用，但注意反算对节点向量敏感，节点选不好方程会病态。

5.4 性能：十万级控制点的分段计算与内存布局

如果控制点上万，std::vector<Point2D>的连续内存布局对缓存友好，但每段计算仍然要访问四个点。优化点在于：把基函数系数预先算成表，避免每次重复计算t2、t3。另外，curvePts的push_back在点数多时会频繁扩容，可以先用reserve预留(n-3)*(samplesPerSeg+1)的空间。如果还要更快，可以把Generate()改成多线程，每段独立计算，最后按段索引合并。但说实话，十万级控制点的三次B样条在普通PC上也就几十毫秒，除非你在做实时渲染，否则没必要过度优化。

从那以后我每次写B样条相关代码，都会先把控制点数量检查、端点重复策略和采样数这三个参数在纸上过一遍，再动手敲键盘。很多翻车不是算法错，而是输入假设没对齐。希望帮到你。

本文还有配套的精品资源，点击获取