简介一份面向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形态两端各重复p1次节点值。这样一来曲线必定穿过第一个和最后一个控制点中间的控制点则负责“吸”住曲线形状。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 基函数递推核心函数BasisFunctionB样条所有计算都绕不开基函数。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曲线末端就会“掉下去”。所以我在p0分支里加了一个判断当u等于最后一个节点且当前是倒数第二个区间时强制返回1这样曲线末端才能精确落在最后一个控制点上。2.2 构造非均匀节点向量clamped形态节点向量的长度有硬性规定控制点数量 曲线次数 1。我用8个控制点、3次曲线举例节点向量总长度就是83112。前面4个是0后面4个是1中间均匀分布4个内节点。这样既保证了端点插值又让内部参数分布均匀关键写起来不容易出bug。 summary 构造clamped节点向量两端各重复p1个节点中间均匀分布 /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.00370.25(1.5708, 1.0000)(1.5654, 0.9969)0.00620.50(3.1416, 0.0000)(3.1420, 0.0015)0.00150.75(4.7124, -1.0000)(4.7180, -0.9958)0.00700.90(5.6549, -0.8090)(5.6528, -0.8121)0.0035看表中t0那一端和t1那一端会非常有意思曲线点几乎和理论点重合因为clamped节点向量强制曲线穿过首末控制点。但中间t0.25和t0.75这种位置曲线点并没有完全压在控制点上——这正好说明B样条是逼近型曲线不是插值型曲线。如果项目要求曲线必须穿过中间某些点那需要反过来做控制点反算那是另一篇文章的活了。3.3 验证暴露出的两个细节端点插值、除零保护这个验证跑通之后我才真正理解为什么clamped节点向量要“两端各重复p1次”。重复节点让基函数在两端变成“独苗”其他基函数在边界处全部为0最后一个基函数独自扛起首末点的位置。如果你把节点向量换成全均匀的首末曲线点就会偏离控制点做刀具轨迹时起点终点就差出去一截产线根本没法用。另外递归里的除零保护不能省。在验证循环中t取0和1这两个端点时基函数递归会走到分母为0的分支如果不加Math.Abs(denom) 1e-12的保护跑一次就会冒出好几个NaN最后画出来的图形是空的。这个问题只在边界出现所以很多人在普通采样下根本发现不了只有像我这样在验证循环里把端点也采样进去时才暴露出来。4. GDI交互演示拖动控制点实时观察曲线变化4.1 坐标变换与绘制管线算法验证正确后剩下的就是把它变成用户能直观操作的工具。窗体上放一个PictureBoxPaint事件里绘制四样东西坐标网格、理论正弦曲线用灰色虚线做参照、控制多边形灰色细线、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坐标混淆有人会直接在循环里让u0到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样条源码都不会再觉得玄乎了。本文还有配套的精品资源点击获取 VB.NET实现非均匀B样条曲线：Cox-de Boor递推与GDI+交互拖拽实战
简介一份面向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形态两端各重复p1次节点值。这样一来曲线必定穿过第一个和最后一个控制点中间的控制点则负责“吸”住曲线形状。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 基函数递推核心函数BasisFunctionB样条所有计算都绕不开基函数。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曲线末端就会“掉下去”。所以我在p0分支里加了一个判断当u等于最后一个节点且当前是倒数第二个区间时强制返回1这样曲线末端才能精确落在最后一个控制点上。2.2 构造非均匀节点向量clamped形态节点向量的长度有硬性规定控制点数量 曲线次数 1。我用8个控制点、3次曲线举例节点向量总长度就是83112。前面4个是0后面4个是1中间均匀分布4个内节点。这样既保证了端点插值又让内部参数分布均匀关键写起来不容易出bug。 summary 构造clamped节点向量两端各重复p1个节点中间均匀分布 /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.00370.25(1.5708, 1.0000)(1.5654, 0.9969)0.00620.50(3.1416, 0.0000)(3.1420, 0.0015)0.00150.75(4.7124, -1.0000)(4.7180, -0.9958)0.00700.90(5.6549, -0.8090)(5.6528, -0.8121)0.0035看表中t0那一端和t1那一端会非常有意思曲线点几乎和理论点重合因为clamped节点向量强制曲线穿过首末控制点。但中间t0.25和t0.75这种位置曲线点并没有完全压在控制点上——这正好说明B样条是逼近型曲线不是插值型曲线。如果项目要求曲线必须穿过中间某些点那需要反过来做控制点反算那是另一篇文章的活了。3.3 验证暴露出的两个细节端点插值、除零保护这个验证跑通之后我才真正理解为什么clamped节点向量要“两端各重复p1次”。重复节点让基函数在两端变成“独苗”其他基函数在边界处全部为0最后一个基函数独自扛起首末点的位置。如果你把节点向量换成全均匀的首末曲线点就会偏离控制点做刀具轨迹时起点终点就差出去一截产线根本没法用。另外递归里的除零保护不能省。在验证循环中t取0和1这两个端点时基函数递归会走到分母为0的分支如果不加Math.Abs(denom) 1e-12的保护跑一次就会冒出好几个NaN最后画出来的图形是空的。这个问题只在边界出现所以很多人在普通采样下根本发现不了只有像我这样在验证循环里把端点也采样进去时才暴露出来。4. GDI交互演示拖动控制点实时观察曲线变化4.1 坐标变换与绘制管线算法验证正确后剩下的就是把它变成用户能直观操作的工具。窗体上放一个PictureBoxPaint事件里绘制四样东西坐标网格、理论正弦曲线用灰色虚线做参照、控制多边形灰色细线、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坐标混淆有人会直接在循环里让u0到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样条源码都不会再觉得玄乎了。本文还有配套的精品资源点击获取