简介这份资源面向计算机图形学、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_{i0}^{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_{i1} \ P_{i2} \ P_{i3} \end{bmatrix} ]其中 ( t \in [0,1] )。这个矩阵形式的好处是每段曲线独立计算四个控制点确定一段段与段之间自然C2连续。VC里只需要把控制点存成std::vectorCPoint2D然后按四个一组滑动窗口计算即可。注意这里的均匀B样条默认节点向量是均匀分布的如果你需要准均匀或非均匀矩阵会变但工程上八成场景均匀就够。2.2 控制点、节点向量与曲线点容器怎么定义在VC里我一般定义两个结构体一个存二维点一个存曲线计算结果。控制点用double而不是int因为后续可能要做缩放和变换。节点向量如果均匀其实不需要显式存储直接用段索引映射即可。但为了扩展性还是留一个std::vectordouble knots均匀情况下可以按i/(n4)生成。// 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::vectorPoint2D ctrlPts; // 控制点 std::vectorPoint2D curvePts; // 采样后的曲线点 int samplesPerSeg; // 每段采样数默认20 BSpline3() : samplesPerSeg(20) {} void SetControlPoints(const std::vectorPoint2D pts) { ctrlPts pts; } void Generate(); // 核心计算 };这里samplesPerSeg控制每段曲线的离散精度。20 段在普通屏幕上已经看不出折线如果要做高精度导出可以调到 50 或 100。控制点数量至少 4 个否则无法构成一段三次B样条。如果少于 4 个常见做法是重复端点或直接退化为直线但我不建议在算法层做这种兼容最好在调用前就保证输入合法。2.3 均匀三次B样条的段数与参数映射假设有 ( n ) 个控制点均匀三次B样条的有效段数是 ( n-3 )。每段对应四个连续控制点第 ( i ) 段用 ( P_i, P_{i1}, P_{i2}, P_{i3} )。参数 ( 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_{i1} (-3t^3 3t^2 3t 1) x_{i2} t^3 x_{i3} \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包含两端所以每段会生成samplesPerSeg1个点相邻段之间会有一个重复点上一段的终点和下一段的起点相同。如果不想重复可以在s循环里从 1 开始但那样第一段的起点会丢所以更常见的做法是保留重复点后续画线时直接按顺序连重复点不影响视觉。参数samplesPerSeg越大曲线越平滑但点也越多导出DXF或SVG时要注意文件大小。3.2 控制点输入与边界处理端点重复与准均匀方案前面提到均匀三次B样条不经过端点。如果你希望曲线从第一个控制点附近开始、到最后一个控制点附近结束最简单的办法是在控制点数组前后各加一个重复点。比如原始控制点是P0, P1, P2, P3, P4变成P0, P0, P1, P2, P3, P4, P4这样曲线会靠近P0和P4。但注意这样会改变段数从n-3变成n1-3 n-2。下面是一个封装函数std::vectorPoint2D BuildClampedCtrlPts(const std::vectorPoint2D raw) { if (raw.size() 2) return raw; std::vectorPoint2D 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[i1].x, (int)m_spline.ctrlPts[i1].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[i3]时越界或者曲线点为空。原因控制点数量小于4但代码没有做前置检查。解决在Generate()开头加if (n 4) return;同时在UI层限制用户至少输入4个点。如果业务上确实可能少于4个那就退化为直线或二次B样条但那是另一套逻辑。我一般会在SetControlPoints里就做一次检查少于4个直接清空并返回错误码。4.4 浮点累积误差导致曲线抖动现象曲线在局部出现微小抖动放大后能看到点不在平滑轨迹上。原因t用float累加或者基函数计算时用了低精度类型。解决所有中间变量用doublet用(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::vectorPoint2D 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::vectordouble BuildChordKnots(const std::vectorPoint2D pts) { std::vectordouble 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::vectorPoint2D的连续内存布局对缓存友好但每段计算仍然要访问四个点。优化点在于把基函数系数预先算成表避免每次重复计算t2、t3。另外curvePts的push_back在点数多时会频繁扩容可以先用reserve预留(n-3)*(samplesPerSeg1)的空间。如果还要更快可以把Generate()改成多线程每段独立计算最后按段索引合并。但说实话十万级控制点的三次B样条在普通PC上也就几十毫秒除非你在做实时渲染否则没必要过度优化。从那以后我每次写B样条相关代码都会先把控制点数量检查、端点重复策略和采样数这三个参数在纸上过一遍再动手敲键盘。很多翻车不是算法错而是输入假设没对齐。希望帮到你。本文还有配套的精品资源点击获取