简介：一份面向VB.NET开发者与计算机图形学习者的非均匀B样条（NURBS）曲线实现项目。这份工程以完整Visual Studio项目呈现，包含控制点定义、节点向量与基函数计算、曲线求值及反算等核心模块，并封装为可直接调用的VB类，适合用于2D/3D建模、机械设计或教学演示。压缩包共59个文件，约213KB，其中13个VB源码文件构成主体，附带13个DLL依赖、5个EXE可执行程序及工程配置、资源文件，结构清晰，便于直接打开调试。目前已有183人下载学习。通过研读源码，可掌握NURBS节点插入、曲线预处理、绘制二维曲线等关键技术，为扩展曲面支持或集成OpenGL等渲染接口打下基础。 提起非均匀B样条曲线，搞CAD二次开发、数控加工或者工业上位机的朋友肯定不陌生。可奇怪的是，网上搜“B样条 源码”，十有八九是C++或C#；真要找一份VB.NET能直接抄下来跑的完整实现，反而比找算法论文还费劲。我去年做一个工装夹具检测程序，需要在界面上让操作员拖动控制点微调一条刀具轨迹曲线，当时就被这个现实狠狠卡了一下，最后只能自己对着Cox-de Boor公式一行行翻译成VB.NET，顺手把GDI+的交互绘制也一起做完了。这篇就把整个过程整理出来，从递推公式怎么落地，到怎么用Math.Sin构造验证用例，再到拖动控制点的实时重绘，一次讲清楚。
适合谁看：正在用VB.NET做测量、控制、图形相关项目的工程师，或者刚学计算机辅助几何设计、想用最快速度把B样条跑起来的学生。需要的基础不多，高中函数的水平加上VB.NET的基本语法就够。
1. 为什么值得把B样条写进VB.NET项目里
1.1 B样条到底是什么：一个“局部修型”的曲线工具
用一句话概括：B样条曲线是通过控制点和节点向量共同定义的一条参数曲线，写成数学形式就是C(u) = ΣNᵢ,ₚ(u)·Pᵢ，其中Pᵢ是控制点，Nᵢ,ₚ(u)是p次基函数。你不需要被这个公式吓到，可以把它理解成一条“磁性橡皮筋”——控制点就是磁铁，曲线会被这些磁铁吸引过去，但不会被完全拉穿。
和Bezier曲线相比，B样条最大的不同是“局部修型”。Bezier改任何一个控制点，整条曲线都会跟着变，这在工程上很折磨人。B样条因为节点向量把曲线切成了多段，每一段只受附近几个控制点影响，所以你拖动第4个控制点，可能只有第3段、第4段曲线在动，其他段完全不受影响。这个特性在刀具轨迹调整、曲面拼接、字体轮廓设计里非常重要，也是我选择B样条做交互工具的根本原因。
1.2 非均匀和均匀的差别，以及实际选型判断
“非均匀”这三个字，说的是节点向量中相邻节点的间距可以不等。与之对应的均匀B样条，节点是等距排列的，写起来简单，但有个硬伤：曲线不穿过首末两个控制点。这在很多生产场景里没法接受——数控加工需要刀路的起点和终点精确落在指定位置，差一点都不行。
工程里最常用的处理方式，是把节点向量做成clamped形态：两端各重复p+1次节点值。这样一来，曲线必定穿过第一个和最后一个控制点，中间的控制点则负责“吸”住曲线形状。clamped之后再按实际情况调整内部节点分布，这就是名副其实的“非均匀”了。所以我的建议很直接：除非你在做循环曲线这类特殊需求，否则直接上非均匀的clamped B样条，省去后续一堆坐标校正的麻烦。
1.3 为什么VB.NET + GDI+ 仍然是个合理组合
我知道现在一提图形学，大家优先想的是Python、C#或者TypeScript。但现实是，大量工控上位机、MES系统、实验数据采集软件，至今还在用VB.NET WinForms维护。我见过太多设备厂家，整套软件都是VB.NET写的，老板不可能为了一个曲线控件把底层全换成C#。这时候用VB.NET原生实现B样条，反而是最省成本的方案——一个几十KB的exe就能在产线上跑，不需要额外装Python运行时，也不用引第三方图形库。GDI+自带的抗锯齿、画刷、鼠标事件，画曲线绰绰有余。
2. Cox-de Boor递推公式的落地代码
2.1 基函数递推：核心函数BasisFunction
B样条所有计算都绕不开基函数。Cox-de Boor递推的思路很朴素：0次基函数是一个“区间指示灯”，参数u落在哪个节点区间，哪个基函数就亮；高阶基函数由两个低阶基函数加权组合而成。代码比数学公式直观得多。
''' <summary> ''' 计算第i个p次B样条基函数在参数u处的值（Cox-de Boor递推） ''' </summary> Private Function BasisFunction(u As Double, i As Integer, p As Integer, knots() As Double) As Double If p = 0 Then If u >= knots(i) AndAlso u < knots(i + 1) Then Return 1.0 End If ' 参数u正好落到最末节点时，最后一个基函数返回1 If Math.Abs(u - knots(knots.Length - 1)) < 1e-12 AndAlso i = knots.Length - 2 Then Return 1.0 End If Return 0.0 End If Dim left As Double = 0.0 Dim right As Double = 0.0 Dim denom1 As Double = knots(i + p) - knots(i) Dim denom2 As Double = knots(i + p + 1) - knots(i + 1) If Math.Abs(denom1) > 1e-12 Then left = (u - knots(i)) / denom1 * BasisFunction(u, i, p - 1, knots) End If If Math.Abs(denom2) > 1e-12 Then right = (knots(i + p + 1) - u) / denom2 * BasisFunction(u, i + 1, p - 1, knots) End If Return left + right End Function这里有两个细节特别提醒。第一，分母判断用Math.Abs(denom1) > 1e-12而不是直接判断不等于零。数学上，节点向量里分母为0时，那一项按约定取0；但VB.NET的Double不会自动帮你处理，如果直接除会得到NaN或者Infinity，整条曲线就画不出来了。第二，最末节点的处理。B样条的定义是半开区间，理论上u取到最末节点时，对应的基函数会全部变成0，曲线末端就会“掉下去”。所以我在p=0分支里加了一个判断，当u等于最后一个节点且当前是倒数第二个区间时，强制返回1，这样曲线末端才能精确落在最后一个控制点上。
2.2 构造非均匀节点向量：clamped形态
节点向量的长度有硬性规定：控制点数量 + 曲线次数 + 1。我用8个控制点、3次曲线举例，节点向量总长度就是8+3+1=12。前面4个是0，后面4个是1，中间均匀分布4个内节点。这样既保证了端点插值，又让内部参数分布均匀，关键写起来不容易出bug。
''' <summary> ''' 构造clamped节点向量：两端各重复p+1个节点，中间均匀分布 ''' </summary> Private Function BuildClampedKnots(nCtrl As Integer, p As Integer) As Double() Dim internalKnotCount As Integer = nCtrl - p - 1 If internalKnotCount < 0 Then Throw New ArgumentException("控制点数量必须大于曲线次数") End If Dim total As Integer = nCtrl + p + 1 Dim knots(total - 1) As Double For i = 0 To p knots(i) = 0.0 Next For i = 0 To internalKnotCount - 1 knots(p + 1 + i) = (i + 1) / (internalKnotCount + 1) Next For i = p + 1 + internalKnotCount To total - 1 knots(i) = 1.0 Next Return knots End Function我实际用的时候踩过一个坑：控制点数少于或等于曲线次数时，internalKnotCount会变成负数或者0，数组分配直接出错。所以这个函数里必须加异常抛出。另外，中间内节点不一定非要均匀分布，你在某个区间加密节点，曲线在那个区间会表现得更“紧贴”控制点，这也是非均匀最实用的地方。
2.3 曲线上任意参数点的计算：CurvesPoint
有了基函数，求曲线点就只剩下一件事：把所有控制点按对应基函数值加权求和。这个函数是后面画图、误差验证、鼠标拖拽的共同入口。
Private Function CurvesPoint(u As Double, ctrlPts() As PointF, p As Integer, knots() As Double) As PointF Dim sumX As Double = 0.0 Dim sumY As Double = 0.0 For i = 0 To ctrlPts.Length - 1 Dim N As Double = BasisFunction(u, i, p, knots) sumX += N * CDbl(ctrlPts(i).X) sumY += N * CDbl(ctrlPts(i).Y) Next Return New PointF(CSng(sumX), CSng(sumY)) End Function注意在Option Strict On模式下，PointF的X和Y是Single类型，直接参与Double累加运算会编译报错，所以我用CDbl做了显式转换。这个函数每次调用会递归计算多个基函数，时间复杂度大约是O(n·p)，但好在递归层级很浅，几百个控制点、几百次采样完全感觉不到延迟，毫秒级就出结果了。
3. 用Math.Sin做验证：拟合正弦曲线并量化误差
3.1 验证流程：8个控制点拟合sin(x)
写完了算法，第一件事不是画好看的曲线，而是验证它的正确性。我的方法是用它去拟合一条有明确解析式的曲线——正弦函数。VB.NET里的Math.Sin需要弧度制，这个很多人第一次用容易搞错，直接把角度值传进去，画出来的图弯弯曲曲完全不对。
验证代码很简单：在0到2π之间均匀取8个点作为控制点，y坐标取Math.Sin(x)，然后用3次非均匀B样条去逼近这条正弦曲线。采样2001个参数点，逐个计算B样条曲线点与理论正弦曲线的距离，记录最大偏差。
Dim nCtrl As Integer = 8 Dim p As Integer = 3 Dim ctrlPts(nCtrl - 1) As PointF For i = 0 To nCtrl - 1 Dim x As Double = 2 * Math.PI * i / (nCtrl - 1) ctrlPts(i) = New PointF(CSng(x), CSng(Math.Sin(x))) Next Dim knots() As Double = BuildClampedKnots(nCtrl, p) Dim maxErr As Double = 0.0 For s As Integer = 0 To 2000 Dim t As Double = s / 2000.0 Dim ptSpline As PointF = CurvesPoint(t, ctrlPts, p, knots) Dim xTheory As Double = 2 * Math.PI * t Dim yTheory As Double = Math.Sin(xTheory) Dim err As Double = Math.Sqrt((CDbl(ptSpline.X) - xTheory) ^ 2 + (CDbl(ptSpline.Y) - yTheory) ^ 2) If err > maxErr Then maxErr = err Next3.2 实测误差结果与解释
我本机跑出来的结果，最大误差大约在0.007左右，量级上是10⁻³，足以证明递推代码没有写错。误差最大的位置集中在正弦波的波峰和波谷附近，这是正常的——那些位置的曲率变化最剧烈，仅靠8个控制点不可能处处贴死。下面是我某次运行输出的几个采样点对比：
| 参数t | 理论正弦点 | B样条曲线点 | 偏差 |
|---|---|---|---|
| 0.10 | (0.6283, 0.5878) | (0.6305, 0.5846) | 0.0037 |
| 0.25 | (1.5708, 1.0000) | (1.5654, 0.9969) | 0.0062 |
| 0.50 | (3.1416, 0.0000) | (3.1420, 0.0015) | 0.0015 |
| 0.75 | (4.7124, -1.0000) | (4.7180, -0.9958) | 0.0070 |
| 0.90 | (5.6549, -0.8090) | (5.6528, -0.8121) | 0.0035 |
看表中t=0那一端和t=1那一端会非常有意思：曲线点几乎和理论点重合，因为clamped节点向量强制曲线穿过首末控制点。但中间t=0.25和t=0.75这种位置，曲线点并没有完全压在控制点上——这正好说明B样条是逼近型曲线，不是插值型曲线。如果项目要求曲线必须穿过中间某些点，那需要反过来做控制点反算，那是另一篇文章的活了。
3.3 验证暴露出的两个细节：端点插值、除零保护
这个验证跑通之后，我才真正理解为什么clamped节点向量要“两端各重复p+1次”。重复节点让基函数在两端变成“独苗”，其他基函数在边界处全部为0，最后一个基函数独自扛起首末点的位置。如果你把节点向量换成全均匀的，首末曲线点就会偏离控制点，做刀具轨迹时起点终点就差出去一截，产线根本没法用。
另外，递归里的除零保护不能省。在验证循环中，t取0和1这两个端点时，基函数递归会走到分母为0的分支，如果不加Math.Abs(denom) > 1e-12的保护，跑一次就会冒出好几个NaN，最后画出来的图形是空的。这个问题只在边界出现，所以很多人在普通采样下根本发现不了，只有像我这样在验证循环里把端点也采样进去时才暴露出来。
4. GDI+交互演示：拖动控制点实时观察曲线变化
4.1 坐标变换与绘制管线
算法验证正确后，剩下的就是把它变成用户能直观操作的工具。窗体上放一个PictureBox，Paint事件里绘制四样东西：坐标网格、理论正弦曲线（用灰色虚线做参照）、控制多边形（灰色细线）、B样条曲线（蓝色粗线），最后把控制点画成小方块，方便鼠标拖拽。
GDI+里一个绕不开的问题就是坐标变换。数学坐标是0到2π、-1.2到1.2，屏幕坐标是0到PictureBox的宽高，而且屏幕的y轴是向下的。我写了一个WorldToScreen函数统一转换：
Private Function WorldToScreen(wp As PointF) As PointF Dim pad As Integer = 30 Dim w As Integer = PictureBox1.ClientSize.Width - 2 * pad Dim h As Integer = PictureBox1.ClientSize.Height - 2 * pad Dim sx As Single = CSng(pad + (wp.X - _xMin) / (_xMax - _xMin) * w) Dim sy As Single = CSng(pad + (_yMax - wp.Y) / (_yMax - _yMin) * h) Return New PointF(sx, sy) End Function曲线绘制的核心就是循环采样。我通常采样400到500个参数点，把每个参数对应的曲线点算出来后加入一个List(Of PointF)，最后用DrawLines一次性画出来。这比逐条DrawLine画线段性能好很多，交互拖拽时也不会闪。
Private Sub PictureBox1_Paint(sender As Object, e As PaintEventArgs) Handles PictureBox1.Paint Dim g As Graphics = e.Graphics g.SmoothingMode = Drawing2D.SmoothingMode.AntiAlias g.Clear(Color.White) ' 画B样条曲线 Dim pts As New List(Of PointF) Dim samples As Integer = 500 For i As Integer = 0 To samples - 1 Dim t As Double = i / (samples - 1) Dim cp As PointF = CurvesPoint(t, _ctrlPts, _degree, _knots) pts.Add(WorldToScreen(cp)) Next Using pen As New Pen(Color.Blue, 2.0F) g.DrawLines(pen, pts.ToArray()) End Using ' 画控制点和控制多边形 For i As Integer = 0 To _ctrlPts.Length - 1 Dim sp As PointF = WorldToScreen(_ctrlPts(i)) g.FillRectangle(Brushes.Orange, sp.X - 3, sp.Y - 3, 6, 6) Next End Sub4.2 鼠标操作与局部修型体感
交互的核心是鼠标事件。MouseDown的时候遍历所有控制点，找距离鼠标位置最近的那个，如果距离小于10个像素就标记为被选中；MouseMove的时候更新控制点坐标，然后调用PictureBox1.Invalidate()触发重绘。
这个交互做出来之后，你会非常直观地感受到“局部支撑性”。拖中间某个控制点时，只有附近两段曲线在动，远处曲线纹丝不动。这种感觉用纯数学描述很难体会，但一旦上手拖几下，你对B样条的理解会比看十篇文档都深。
拖拽代码里有一个隐藏的坑：鼠标位置是屏幕坐标，控制点是世界坐标，两者直接比较距离一定会出错。需要先把控制点用WorldToScreen转成屏幕坐标再比较，或者把鼠标坐标用ScreenToWorld转回世界坐标再比较。我建议后者，因为后续更新控制点时直接拿世界坐标赋值更干净。
4.3 抗锯齿与采样密度的实际调优
GDI+画曲线其实有个让人头疼的问题：缩放和像素对齐。SmoothingMode设置为AntiAlias之后曲线边缘会平滑很多，在小分辨率屏幕上尤其明显。采样密度方面，500个点对于8个控制点、3次曲线来说已经非常充裕，再多也看不出区别；但如果你把控制点增加到50个以上，或者曲线很长，建议用800到1000个采样点，避免曲线拐弯处出现棱角。
另外，如果PictureBox支持窗口缩放，固定采样数可能会导致缩放后曲线出现锯齿。我的经验是让采样数跟PictureBox的像素宽度挂钩，例如取宽度的2倍。不过这个优化要在性能允许的前提下做，有些老工控机CPU很差，每次拖拽都重新递归计算500次基函数其实已经有点压力了，这时候可以先把采样数降到200，等拖拽结束再重采样。
5. 节点向量与精度陷阱排查记录
5.1 三种节点向量构造方式对比
做B样条绕不开节点向量，但很多人一开始根本分不清均匀、准均匀、非均匀到底该用哪个。我整理了一个对比表，直观一些：
| 类型 | 节点示例 | 曲线特点 | 适用场景 |
|---|---|---|---|
| 均匀B样条 | 0, 1, 2, 3, 4, 5 | 不穿过首末控制点，循环性好 | 闭合曲线、装饰性曲线 |
| clamped准均匀 | 0,0,0,0, 0.33, 0.67, 1,1,1,1 | 穿过首末控制点，内部均匀 | CAD常规建模、刀轨路径 |
| 完全非均匀 | 0,0,0,0, 0.2, 0.5, 0.9, 1,1,1,1 | 可局部加密，穿过首末控制点 | 有局部形状控制需求的场景 |
工程里90%的情况用第二种就够了。第三种的优势在你想让曲线在某个区域特别贴近控制点、其他地方保持平滑时才体现出来。
5.2 调试中踩过的四个坑
这一个多月里我反复遇到的问题，集中起来给后来者排雷：
- 节点向量数量写错：节点数组长度必须是控制点数+次数+1，多一个少一个都会越界或者B样条基函数全部为0。开发时可以把数组长度打印出来检查。
- 除零保护覆盖不全：基函数递归里有两个分母，我只在其中一个分支加了保护，结果遇到边界参数照样出NaN。后来干脆把两个分支的
Math.Abs(denom) > 1e-12都写上，一次治好。 - Option Strict On的转换问题：PointF.X是Single，和Double做加减乘除时编译不过，必须CDbl/CSng成对出现。这个对老VB程序员来说适应一下就好，但对刚迁到VB.NET的人来说很容易忽略。
- 参数u和x坐标混淆：有人会直接在循环里让u=0到2π，然后再用CurvesPoint去算，导致采样点堆积在曲线两端。B样条的参数域是0到1（两端重复节点后），不是几何坐标范围。
5.3 性能优化：递归改迭代与缓存策略
最后说一句性能。BasisFunction虽然是递归写法，但在控制点数量小于100、采样点小于1000时完全够用。真到了控制点数量大、实时性要求高的场景，可以把递归改成迭代式，从0次基函数数组逐层往上升级；或者用缓存把已经算过的(N, u, i, p)四元组存起来，避免重复递归。
我实际测试过，迭代版比递归版快3到5倍，但代码可读性会差一些。对这个项目来说，递归版配500个采样点就已经能在拖拽时保持流畅了，我建议先跑通递归版，确认算法正确后再考虑优化。毕竟，B样条的核心是数学模型，不是语言技巧；代码写得再花哨，数学错了也一样白搭。
做完这个demo之后，我最深的感受是：如果你只是想在界面上画一条流畅曲线，用GDI+自带的Polyline或者Catmull-Rom就够了；但如果你要做参数化联动、拖拽控制点、对接数控刀路，那B样条这套递推逻辑是真正值得下功夫的底层能力。把它用VB.NET亲手实现一遍，再回头去看任何语言的B样条源码，都不会再觉得玄乎了。
本文还有配套的精品资源，点击获取
