简介一套基于Matlab的NURBS非均匀有理B样条工具箱面向需要处理复杂曲线曲面建模、数据拟合与几何造型的 CAD、计算机图形学及工程计算领域用户。压缩包仅 1.47MB目前已有 103 人学习使用。工具箱围绕 NURBS 核心知识体系组织包括控制点网格、knot 向量与权重因子的定义与操作内置 nurbscreate、nurbsedit、nurbsplot、nurbsfit 等常用函数可支持从创建、编辑、可视化到数据点拟合的完整流程。通过调整控制点与 knot 向量还能实现多个曲线曲面的无缝拼接部分目录提供交互式 GUI 工具便于直观调节参数并观察结果。此外该工具箱支持与 CAD 软件交换 IGES/STEP 格式数据并结合其他 Matlab 工具实现逆向工程、有限元前处理、动画制作等高级应用。掌握这套资源可大幅提升 NURBS 建模与分析的效率。1. 从 CAD 回到 Matlab为什么还要研究 Nurbs Toolbox for Matlab当你在做逆向工程、船体外形、叶片曲面拟合或机器人轨迹规划时CAD 工具内部用 NURBS 定义精确曲线曲面但要把几何算法放入 Matlab 做参数优化、仿真或与其他工具箱串联Matlab 自带样条相关函数只处理多项式形式表达不了圆弧这类有理几何。于是出现了一批以“Nurbs Toolbox for Matlab.zip”为名的第三方工具箱它们把控制点、权重、节点向量封装成统一结构解压后 addpath 即用。这类工具箱不是官方产品却因为贴近工业几何格式而在航空航天、船舶和逆向工程团队里流传。适合需要快速验证 NURBS 算法、不想从底层实现 B 样条递归公式的工程师。下面按最可靠的常用方案讲安装、数据结构、基本操作和排错。2. 把 Nurbs Toolbox 装进 Matlabzip 解压、路径设置和最小验证2.1 不要只 unzipMATLAB 搜索路径决定你能否调用拿到 zip 后很多人第一反应是解压到临时目录然后直接构造coefs调用nrbmak结果 Matlab 报 Undefined function。原因不是代码错而是工具箱路径没进path。与从 Add-On Explorer 安装不同zip 分发的外部工具箱不会自动注册重启后更不会保留。常见做法是解压到固定位置再手动执行addpath(genpath(...))。下面这段负责解压、加路径并保存% 解压到用户默认路径下的 toolboxes 目录 unzip(NurbsToolbox.zip, fullfile(userpath, toolboxes)); % 确认解压出的根目录名常见是 NurbsToolbox 或 nurbs_toolbox tb_path fullfile(userpath, toolboxes, NurbsToolbox); % genpath 会包含所有子目录避免漏掉 demo 或 utilities addpath(genpath(tb_path)); % 保存到 pathdef.m否则重启 Matlab 后失效 status savepath; if status 0 disp(路径已保存); else disp(当前用户无权限写入 pathdef.m请在 startup.m 中补 addpath); end逻辑说明unzip的第二参数指定目标位置userpath在不同操作系统下的默认值不同但都是用户可写目录。genpath返回递归子目录适合工具箱内部存在nurbs类或utils目录的情况。savepath失败时Matlab 不会给出明显错误只看返回状态码此时把addpath写进startup.m是更稳妥的替代方案。不要直接解压到 Program Files 下的 Matlab 安装目录。那个目录权限受限savepath很容易失败而且 Matlab 升级重装时会被清掉。路径中也不要带中文或首尾空格某些版本的which命令会因此找不到函数。2.2 有 install.m 时先读一下它再决定是否执行很多 NURBS 工具箱 zip 会附带install.m。这个脚本通常做的事是向当前目录写入startup.m或在pathdef.m中追加路径。但它的运行环境未必与你的系统一致尤其当 zip 来自不同维护者时install.m里的路径分隔符、子目录名可能写死。我一般会先打开install.m看前几行再决定是否运行。常见执行方式如下cd(fullfile(userpath, toolboxes, NurbsToolbox)); % 存在 install.m 时运行它 if exist(install.m, file) run(install.m); end % 验证主函数是否可见 which nrbmakwhich nrbmak如果返回完整路径说明安装完成如果返回空继续执行path命令检查目标目录是否在列表中。一个容易混淆的情况是Matlab 自带的spline、ppmak属于 Curve Fitting Toolbox函数路径在 matlabroot 下所以which nrbmak找不到不代表系统坏了只是第三方路径缺失。这与 Robotics Toolbox 等第三方 Robotic 工具箱很像这类独立工具箱的安装成败往往不是代码问题而是工具所有权和路径管理。install.m执行时若提示无法写入pathdef.m就退回手动addpath。在 R2023b 之后的新版 Matlab 上MathWorks 收紧了安装目录写入权限所以这一情况更常见。下面这张表可以帮你快速决定用哪种安装方式安装方式适用情况主要风险Add-On Explorer官方来源自动管理路径第三方 zip 不一定能搜到手动addpath本地 zip控制力最强重启后失效需 savepathinstall.m脚本维护者提供的快速安装可能写入错误路径未知行为多2.3 最小验证不用 demo用一条直线确认安装很多教程会让你运行 demo 脚本但 demo 脚本依赖绘图首次运行可能因为显卡或旧函数名卡住。更快的最小验证是直接构造一条三次 NURBS 直线并画出它% 用 nrbmak 构造一条从 (0,0) 到 (1,1) 的三次曲线 coefs [0 1/3 2/3 1; % x 坐标 0 1/3 2/3 1; % y 坐标 0 0 0 0; % z 坐标 1 1 1 1]; % 权重全 1 knots [0 0 0 0 1 1 1 1]; % 4 个控制点order4节点数8 crv nrbmak(coefs, knots); figure; nrbplot(crv, 30);这里coefs是 4×n 齐次坐标矩阵前三行是 XYZ第四行是权重。knots长度必须等于控制点数加 order。order4对应三次曲线即 degree3。整条曲线几何上是一条对角线但控制点会把曲线拉向内部折线显示时能看到控制多边形与曲线的差异。如果这条线能画出来说明nrbmak和nrbplot都已可用安装闭环完成。3. Nurbs 曲线核心对象coefs、order 与 knots 的三方关系3.1 权重和齐次坐标NURBS 与 B 样条的唯一差异NURBS 全称 Non-Uniform Rational B-Splines和普通 B 样条的唯一区别是每个控制点带权重。曲线上某个点的计算方法是先对权重乘控制点的齐次坐标求和再除以权重之和。这个除法让曲线变成有理函数能表达 B 样条无法精确表达的圆弧、椭圆和抛物线。工具箱内部把这类对象统一存储为结构体而不是一个个独立变量。常见的字段包括form、number、order、knots和coefs。其中coefs使用齐次坐标第四行是权重而不是可选的独立数组。权重为 1 时NURBS 退化为 B 样条权重大于 1 时曲线被拉向对应控制点小于 1 则被推开。这个“拉近推开”的效果就是 CAD 软件里权重滑块的物理意义。字段名维度或类型含义formchar对象类型曲线为 B-NURB曲面为 B-NURBSnumber1×2 或标量控制点数量曲面为 [nu nv]order1×2 或标量degree1knotscell 或 vector节点向量曲面为两个 cellcoefs4×n 或 4×nu×nv齐次控制点坐标3.2 用 nrbmak 手工构造一条二次 NURBS以半径 1 的四分之一圆弧为例这是 NURBS 工具箱里最常见的起点示例。圆弧需要三个控制点中间控制点权重取根号二分之一节点向量在两端各重复三次% 四分之一圆弧P0(1,0), P1(1,1), P2(0,1) x [1 1 0]; y [0 1 1]; z [0 0 0]; w [1 sqrt(2)/2 1]; % 中间权重让曲线呈圆弧 coefs [w.*x; w.*y; w.*z; w]; knots [0 0 0 1 1 1]; % 二次曲线order3节点长度控制点数order arc nrbmak(coefs, knots); % 求参数 0 和 1 处的点验证端点 [pnt, ~] nrbdeval(arc, [0 1]); disp(pnt(1:2,:));代码要点coefs的第四行必须和前三行同步乘权重不能单独跳过。knots一定要单调不减两端重复次数等于 order。这里order3表示二次曲线但内部节点少所以退化为一段有理贝塞尔曲线。nrbmak输入顺序是coefs, knots两参数不可颠倒。如果少写权重那一行size(coefs,1)变成 3nrbmak会直接报错而不是自动补权重为 1。很多新手在这里卡住总是下意识把 NURBS 当成普通 CAD 点列。3.3 曲面是两条方向的笛卡尔积曲面对象是曲线概念的二维扩展。控制点不再是 n× 数组而是 nu×nv 网格coefs变成 4×nu×nv 的三维数组。knots是两个 cell分别存放 U 方向和 V 方向的节点向量。一个最简单的 2×2 平面曲面可以这样构造% 四个角点每个都是 4×1 齐次坐标 c1 [0; 0; 0; 1]; c2 [1; 0; 0; 1]; c3 [0; 1; 0; 1]; c4 [1; 1; 0; 1]; % cat 的第三维对应 V 方向 coefs cat(3, c1, c2, c3, c4); % 结果维度 4×2×2 knots {[0 0 1 1], [0 0 1 1]}; % 每个方向都是线性阶数 srf nrbmak(coefs, knots); nrbplot(srf, [10 10]);这个例子两个方向都是一次所以曲面是平面上的双线性面片。真正有工程价值的曲面通常从一个截面曲线开始用nrbextrude或nrbrevolve拉伸旋转得到不需要逐点排列控制网格。理解coefs的三维维度是手工构造曲面的基础后面做控制点投影、曲面拼接时都要依赖这个布局。4. Nurbs 实操节点插入、升阶、求导和插值拟合的边界4.1 精化不改变几何nrbkntins 和 nrbdegelev曲线外观由控制点、权重和节点向量共同决定局部加密不能只增加控制点还要同步修改节点向量。工具箱提供nrbkntins插入节点插入后控制点数量增加但原曲线几何完全不变。这是后续做拟合、局部调整的第一步。% 在参数 0.25 和 0.75 处插入节点 arc2 nrbkntins(arc, [0.25 0.75]); % 查看新节点向量 knots_new arc2.knots;nrbkntins第二个参数是待插入节点位置必须落在定义域区间内。插入后控制点数会相应增加但不能期望它改变曲线的局部形状。要改变形状需要挪移控制点或调整权重。升阶是另一个不改变几何的精化操作arc3 nrbdegelev(arc, 1); % 一次提升一个阶数升阶后权重会被重新计算不能沿用原有权重数组。nrbkntins和nrbdegelev是 NURBS 里最常用的两个精化命令它们不逼近、不迭代结果在数学上是精确等价的。4.2 求值和求导用 nrbdeval 拿到曲线上的点和切线方向轨迹规划和曲面法向量计算都需要求导。nrbdeval同时返回参数点处的坐标和一阶导使用方式如下tt [0 0.5 1]; [pnts, ders] nrbdeval(arc, tt); % pnts 是 3×N 坐标矩阵 % ders 是 3×N 切向量矩阵tt是参数向量单位区间对应节点定义域。ders的列是对应参数处的 dP/du单位与曲线坐标一致。如果tt超出节点向量定义域结果可能为空或 NaN。另外一阶导不等于切线方向本身工程使用前最好归一化。曲面对象也可以调用nrbdeval但参数要传 cell 数组例如{[0.5], [0.3]}返回的导数会有两个分量。这里就不再展平为单个矩阵了。4.3 插值与拟合阶数不是越高越好参数化方式决定成败如果要用 NURBS 通过一组型值点常见做法是先做参数化再求解线性方程组得到控制点。许多 zip 版本提供nrbinterp接口但参数化方式需要自己控制。假设有 5 个二维点pts [0 0; 1 1; 2 1.5; 3 0.5; 4 0]; % 均匀参数化第二个参数传空 fit1 nrbinterp(pts, []); % 弦长参数化结果更符合点列分布 seg sqrt(sum(diff(pts, 1, 1).^2, 2)); tt cumsum([0; seg]); tt tt / tt(end); fit2 nrbinterp(pts, tt);均匀参数化适合点距相近的输入一旦点距差异大会出现明显过冲。弦长参数化把参数间隔与几何距离关联起来能有效抑制振荡。两种方式对比下来大部分工程点列都应优先采用弦长参数化。插值阶数没有必要设得很高。一般order4即三次已经能保证二阶连续性视觉上足够光顺。超过五次容易产生 Runge 现象控制点坐标出现量级跳动。若曲线仍无法达到误差要求正确做法是先用nrbkntins增加内部节点再做最小二乘拟合而不是继续升阶。下面这张表总结常用函数的边界函数作用是否改变几何关键注意点nrbmak新建 Nurbs 对象建立coefs 和 knots 顺序固定nrbkntins节点精化否节点必须在定义域内nrbdegelev升阶否阶数不宜超过 5 次nrbdeval求值和求导否参数越界返回 NaNnrbinterp型值点插值是参数化方式影响极大验证一个操作是否改变几何最直接的方法是在同样参数点上比较采样坐标。比如插入节点后用nrbdeval在 50 个点处比较误差应在浮点精度内。5. Nurbs 可视化验证nrbplot 采样、报错定位和小技巧5.1 nrbplot 的第二个参数采样数不是控制点数量nrbplot(crv, 80)里的 80 只决定等参数采样点数不改变 Nurbs 对象本身。80 大于控制点数不会增加控制点小于控制点数也不会降低精度。但显示效果对采样数敏感节点多时采样数不足会看到明显折线。figure; nrbplot(fit2, 80); % 80 个采样点绘制曲线 nrbctrlplot(fit2); % 叠加控制多边形 grid on; exportgraphics(gcf, fit_curve.eps, ContentType, vector);exportgraphics在 R2020a 之后的版本可用输出矢量 eps 适合论文插图。这个命令搭配nrbctrlplot能一次得到曲线和控制多边形比起反复调整nrbplot参数更直观。5.2 三个高频报错knot 排序、尺寸不匹配和 NaNnrbmak报错时先查三件事节点向量必须单调不减coefs第一维必须是 4节点向量长度必须等于控制点数加 order。这三条是 Nurbs 数据结构的基本约束可以在建模前写一个断言函数拦截function ok check_nurbs(nrb) k nrb.knots; assert(all(diff(k) -1e-14), knot 必须单调不减); assert(size(nrb.coefs, 1) 4, coefs 行数必须为 4); assert(length(k) nrb.number nrb.order, ... 节点数必须等于控制点数order); ok true; end若曲线显示正常但计算输出 NaN多半是权重出现了 0 或负数。权重 0 会让分母求和失效负权重虽然能显示但会产生反直觉参数化不建议在生产代码里使用。还有一个隐蔽错误Matlab 内置函数名与工具箱函数冲突。which -all nrbmak可以查看是否有多个同名文件如果有Matlab 会优先选择搜索路径里靠前的那一个。常见的是用户自己写过同名nrbmak.m结果一直调用旧版本。5.3 一个小技巧用圆弧权重验证工具箱是否符合预期最后分享一个自检技巧。四分之一圆弧的权重和节点向量都是确定值如果工具箱正确实现了有理求值那么参数中点 0.5 处应该精确落在 (√2/2, √2/2) 上。用几行代码就能验证arc nrbmak([1 1 0; 0 1 1; 0 0 0; 1 sqrt(2)/2 1], ... [0 0 0 1 1 1]); [pnt, ~] nrbdeval(arc, [0 0.5 1]); disp(pnt(1:2, :));输出第二列如果接近 0.7071说明工具箱的实现正确偏差超过 1e-6 就要怀疑路径里掺入了旧版本。这个小技巧比运行完整 demo 更接近底层因为它只测了最小单元二次有理贝塞尔圆弧。以后拿到任何新的 Nurbs Toolbox for Matlab 压缩包都可以用这段代码作为第一道验收门槛。本文还有配套的精品资源点击获取