简介MATLAB NURBS工具箱是一套面向MATLAB环境的NURBS建模扩展库适用于需要处理非均匀有理B样条曲线的工程师与科研人员常用于CAD、CAM、CAE及计算机图形学中的复杂几何建模与数据拟合。压缩包共60个文件其中59个为.m函数源码覆盖曲线曲面创建、控制点编辑、参数化、插值拟合、修剪旋转等完整操作流程另有1个txt文件提供配套说明整体仅44KB轻量易部署。目前已有452人学习下载。借助工具箱中的nrbcreate、nrbplot、nrbinterpolate、nrbrevolve等核心函数用户可以快速构建并可视化NURBS曲线与曲面完成从基础绘制到复杂变换的各类几何运算同时支持IGES/STEP等格式导出便于与外部软件衔接。这份资源适合希望系统掌握NURBS理论并需要在MATLAB中落地应用的初学者和进阶开发者。1. 一个 zip 包背后的 MATLAB NURBS 工具箱到底解决什么问题当你从网络或同事手里拿到一个名为 MATLAB NURBS工具箱.zip 的压缩包时很容易产生两个错觉一是解压出来就能用二是它安装后会和 MATLAB 里自带的优化工具箱那样出现在 Add-On 管理器里。实际这两条都不成立。这个 zip 包通常是把第三方编写的 NURBS 建模函数集合打包分发没有安装向导也不会自动注册到 MATLAB 的路径系统。解开后你可能看到 nrbline、nrbcirc、nrbkntins 这些看似神秘的 m 文件如果直接把路径指到某个子目录很多函数会因为找不到对应文件而报错。本文会从 NURBS 的基本数据结构讲起然后带你完成解压、路径配置、曲线曲面的构建与求值最后给出一套验证工具箱是否真正生效的方法。适合需要用 MATLAB 做曲线曲面建模、样条拟合或机器人路径规划的工程师也适合正在做几何设计相关课程作业的学生。2. 先立住理论NURBS 的数据结构就是 MATLAB 代码的地图2.1 B 样条基函数与节点向量最容易读错的参数NURBS 是 Non-Uniform Rational B-Splines 的缩写核心是在 B 样条基础上引入权因子和齐次坐标。要理解工具箱里的每一个函数在做什么必须先看它的三个基本输入控制点矩阵、节点向量和阶次。控制点用坐标列出比如二维曲线的控制点就是 N×2 矩阵三维则是 N×3节点向量knot vector是一组非减实数例如[0 0 0 0.5 1 1 1]它决定基函数在哪里拼接也决定参数域的范围。阶次order在工具箱里习惯用p表示实际次数是p-1比如 order 为 3 意味着二次 B 样条。很多人在第一次写代码时会把节点向量和控制点顺序搞混或者把 order 当成控制点个数这是最典型的错误。因为单元格数组形式的节点向量允许在不同维度的方向上定义不同的次数和拼接方式比如曲面 UV 方向的节点向量可以不同。还要注意同一条曲线同样的控制点如果节点向量分布不一样形状会差别很大。B 样条基函数由 Cox-de Boor 递推公式计算工具箱内部已经封装好你不需要自己写但必须理解参数 t和几何坐标之间的关系。2.2 工具箱用一个 struct 装下整个曲线或曲面NURBS 工具箱的数据结构很有设计感一条曲线或一个曲面就是一个 MATLAB struct。字段通常包含form、dim、order、knots、coefs。form是描述类型的字符串比如line、circ、surfdim是几何空间的维数2 表示二维3 表示三维order是阶数向量knots是 cell 数组每个单元格对应一个参数方向的节点向量coefs是控制点齐次坐标矩阵曲线是 4×(n1)曲面是 4×(n1)×(m1)其中第四行存放的是权重系数。这个第四行的存在正是区分普通 B 样条和 NURBS 的关键。理解coefs的行布局对后续改权重很重要。前三行是 x、y、z 坐标但注意它们是权重乘以原始坐标之后的齐次坐标值不是直接可见的控制点位置。如果你直接改coefs(1,i)你改变的是齐次坐标的 x 分量而最终显示出来的实际控制点坐标会被第四行的权重重新除回来。所以正常修改权重的方法是同时调整前三行或者用工具箱提供的变换函数。下面的代码创建一个直线段 NURBS 曲线然后打印出它的结构字段让你直观看到每个部分长什么样p nrbline([0 0 0], [1 1 0]); disp(p);nrbline用两个端点构造线段返回一个 struct。运行后可以在命令行看到form是linedim是 3order是 2knots是{ [0 0 1 1] }coefs是 4 行 2 列。这个简单的例子能让你验证路径是否配置正确同时熟悉数据结构的形态。如果你在命令行能正常显示这个 struct 而不是报错说明工具箱的核心文件已经可用了。3. 解锁 zip 包用 addpath 正确把工具箱装进 MATLAB3.1 解压与安置该把文件放到哪个目录拿到 MATLAB NURBS工具箱.zip 后右键解压不要直接双击打开看两眼就关掉更不要直接在压缩包内运行 m 文件。很多第三方工具箱都会带多层子目录比如nurbs/下面还有nurbs/或者分成inst、examples、tests等。把它解压到 MATLAB 工作路径之外的一个固定位置比较好比如D:\MatlabTools\nurbs_toolbox或 Linux 下的~/MatlabTools/nurbs_toolbox这样即使切换当前工作目录也不会影响调用。如果你解压到临时目录每次重启 MATLAB 都得重新配置浪费时间还容易忘。解压后用文件管理器确认根目录下存在名为nrbline.m、nrbeval.m、nrbkntins.m等文件。有些压缩包会把这些核心文件放在更深的子目录里比如notoolbox/或src/此时你最好不要挨个手动添加每一个文件而是使用genpath把根目录下所有子文件夹一次性加入搜索路径这才是完整安装的常见做法。3.2 安全添加路径genpath、savepath 与 which 验证打开 MATLAB 后把当前目录切换到解压后的根目录执行以下代码完成路径注册% 假设你解压到了 D:\MatlabTools\nurbs_toolbox toolbox_path D:\MatlabTools\nurbs_toolbox; % 递归获取该目录下所有子文件夹路径 allpaths genpath(toolbox_path); % 逐个添加到MATLAB搜索路径 addpath(allpaths); % 保存到路径定义文件下次启动自动生效 status savepath; disp(status); % 验证核心函数是否可见 which nrbline这段代码的关键点在于genpath它会生成一个以分号分隔的路径字符串把所有子目录都包含进来。addpath可以直接接收这个多路径字符串不需要循环。最后savepath把路径写入pathdef.m如果返回1表示写入失败通常是因为 MATLAB 安装目录没有写权限这时可以改用savepath(C:\Users\你的用户名\Documents\MATLAB\pathdef.m)保存到用户目录。which nrbline是验证安装是否成功的最直接手段正常会输出D:\MatlabTools\nurbs_toolbox\nrbline.m这样的完整路径。如果在which nrbline后仍然显示nrbline not found不要急着怀疑文件损坏先检查解压目录内是否真的存在该文件再确认你设置的toolbox_path是否指向了包含nrbline.m的那一层。常见错误是把路径指到了nurbs_toolbox/的上一级或指到了内含nrbline.m的子文件夹的下一级。3.3 和官方工具箱不同这个工具箱不是从 Add-On 管理器装的很多用户把 MATLAB 自带工具包的安装方式套到这个第三方库上结果在 Home 页面的 Add-On 资源管理器里怎么也搜不到它。这不奇怪NURBS 工具箱不是 MathWorks 官方发布的工具包不能通过 MATLAB 的集成化安装流程来处理。它本质就是一堆 m 文件、p 文件和演示脚本的集合所以安装这个动作在 MATLAB 里对应的就是添加搜索路径。正因如此卸载也很简单直接从路径中移除相关目录即可和matlab下载安装教程里讲的完整安装完全不冲突。同样地这个工具箱和 matlab优化工具箱 也没有任何概念上的重叠。优化工具箱解决的是数学规划问题NURBS 工具箱解决的是几何表达和样条计算问题二者目录结构、函数命名和数据格式都截然不同。如果在使用 NURBS 相关函数时报错说找不到确定类型的输入参数不要去检查官方优化工具箱那说明你的路径配置还没有做到位。4. 动手造曲线曲面从线条、圆到蒙皮样条4.1 基础图元nrbline、nrbcirc 与 nrb4surf配置好之后先用几个基础图元把工具箱跑通。下面这段代码首先生成一条直线段、一段圆弧和一个由四点定义的平面曲面然后绘制出来% 直线段起点(0,0,0)终点(1,1,0) line_nrb nrbline([0 0 0], [1 1 0]); % 圆弧圆心(0,0,0)半径1从0度到90度 circ_nrb nrbcirc(1, [0 0 0], 0, pi/2); % 四点曲面四个角点形成一个平面矩形 surf_nrb nrb4surf([0 0 0], [1 0 0], [0 1 0], [1 1 0]); % 将三条几何对象绘制在一张图上 figure; nrbplot(line_nrb, 20); hold on; nrbplot(circ_nrb, 30); nrbplot(surf_nrb, [20 10]); axis equal; grid on;nrbcirc的第一个参数是半径第二个参数是圆心坐标第三和第四个参数是起始和终止角度单位为弧度。注意如果你写nrbcirc(1, [0 0 0], 0, 90)那就会画出一个 90 弧度的圆弧因为工具箱不会自动把度数转成弧度。nrb4surf要求四个角点按对角顺序传入它会自动生成一片双线性曲面非常适合作为后续拉伸或蒙皮操作的初始片体。nrbplot的第二个参数是采样点数这里给了直线 20 个点、圆弧 30 个点采样越多绘制越平滑但也会增加计算量。4.2 改动几何移动控制点、修改权因子、插入节点NURBS 的一大优势是局部支撑性移动一个控制点只会影响附近区域的形状。下面的例子以一条二阶 B 样条曲线为对象展示三种常见修改操作% 原始控制点经过三个点的二次曲线 ctrl_pts [0 0 0; 1 2 0; 2 0 0]; curve nrbmak(ctrl_pts, [0 0 0 1 1 1]); % 在参数中点插入一个新节点 new_curve nrbkntins(curve, 0.5); % 移动原始曲线的第二个控制点到新位置 curve.coefs(1,2) 1.5; curve.coefs(2,2) 0.5; % 将权因子改为0.2让曲线远离该控制点 curve.coefs(4,2) 0.2; % 绘制对比 nrbplot(curve, 50); hold on; plot(ctrl_pts(:,1), ctrl_pts(:,2), ro);这里nrbmak是直接利用控制点矩阵和节点向量构造曲线返回的curve里coefs前三行默认已经包含了权重为 1 时的坐标值。当你把第二列第四行的权因子改成 0.2 时曲线会被推离那个控制点因为除数变小了。修改coefs(1,2)和coefs(2,2)则是在改变齐次坐标最终的可视化坐标会自动折算。nrbkntins返回新的曲线节点插入不会改变几何形状但会让该参数位置拥有局部修改能力。这个特性在后续做精细调整时非常有用。4.3 求值、求导和可视化从数据到一阶导数曲线建模完成后往往需要把曲线转成离散点供机器人路径规划或加工代码使用工具箱里的nrbeval专门做这件事。下面代码对圆弧在参数点 0 到 1 之间采样并同时计算一阶导数值% 在参数域[0,1]上取100个等距点 t linspace(0, 1, 100); % 对圆弧求值得到欧氏坐标点集 pts nrbeval(circ_nrb, t); % 求一阶导数 dn nrbderiv(circ_nrb); [~, dpts] nrbeval(dn, t); % 绘制圆弧和导数箭头示意 figure; plot(pts(1,:), pts(2,:), b-, LineWidth, 2); hold on; step 10; % 每隔10个点画一条导数线 for i 1:step:100 quiver(pts(1,i), pts(2,i), dpts(1,i), dpts(2,i), 0.1, r); end axis equal; grid on;nrbeval的第一个返回值是曲线上的欧氏坐标第二个返回值是对应的权重值如果你不需要权重可以写成pts nrbeval(circ_nrb, t);但要注意你得到的是一个 3×N 的矩阵。nrbderiv返回一个新的 NURBS 结构体这个结构体的coefs存储的是导数的控制点后续再用nrbeval就可以得到各参数点的切向量。上面代码里用quiver把每个离散点的导数方向画成箭头你会发现切向量方向始终与圆弧在该点的切线方向一致这是验证求导正确性的一个直观方法。下表列出几个最常用函数的作用和参数格式方便你在写脚本时快速查阅函数名典型调用作用返回值nrblinenrbline(p0, p1)创建直线段NURBS 结构体nrbcircnrbcirc(r, center, ang0, ang1)创建圆弧NURBS 结构体nrbmaknrbmak(coefs, knots)由控制点和节点向量构造曲线或曲面NURBS 结构体nrbkntinsnrbkntins(nrb, knots)插入节点或细化控制网格NURBS 结构体nrbdegelevnrbdegelev(nrb, num)提升一个方向的阶次NURBS 结构体nrbeval[pts, w] nrbeval(nrb, t)对参数值求几何坐标坐标矩阵和权重nrbderivdn nrbderiv(nrb)生成导数几何体NURBS 结构体nrbtformnrbtform(nrb, mat)应用变换矩阵移动/旋转/缩放NURBS 结构体nrbdegelev和nrbtform是和nrbkntins并列的高频操作。提升阶次能让曲线在保持形状不变的情况下拥有更多可用自由度而nrbtform接收一个 4×4 的齐次变换矩阵可以一次性把曲线整体平移到新位置比逐点改coefs要稳妥。5. 验证和排错怎样确定工具箱真的跑对了5.1 用单位圆验证权因子是否生效NURBS 与普通样条的最大区别在权因子而权因子最容易在数值上被弄错。下面这段代码用半圆弧验证几何精度如果权因子没生效或者被错误地归一化圆半径不会稳定在 1% 创建从0度到180度的半圆弧圆心原点半径1 circ nrbcirc(1, [0 0 0], 0, pi); % 在参数点0到1上取5个点 t [0 0.25 0.5 0.75 1]; pts nrbeval(circ, t); % 计算这些点到原点的距离 r sqrt(pts(1,:).^2 pts(2,:).^2); % 显示半径 disp(r);正常输出应该是五个都接近 1例如1.0000 1.0000 1.0000 1.0000 1.0000。如果发现某个点半径偏离明显说明coefs中权重行可能被意外覆盖。注意nrbeval返回的坐标已经是把权重除回去的欧氏坐标所以半径计算就是直接取 x² y² 开根号不需要再额外处理。这个方法同时可以验证你在 4.2 节中手动修改coefs(4,:)是否造成了不可预料的形状突变因为单位圆的权重值不是常数半圆弧中间控制点的权重是cos(π/4)量级而不是 1。5.2 三个容易踩的坑路径没生效、函数冲突、节点向量范围第一坑是路径没生效。如果你用addpath添加之后当前绘图能用但重启 MATLAB 后又报Undefined function or variable nrbline说明savepath没有真正保存成功。解决方法是检查savepath返回值如果为 1就在用户目录下创建一个pathdef.m或者直接使用matlabrc开机脚本。最简单的定位办法是每次打开 MATLAB 后先执行which nrbline如果找不到立刻重新addpath(genpath(你的工具箱路径))不要继续往下跑其他的模型代码。第二坑是函数冲突。你可能同时安装了多个版本的 NURBS 工具箱或者 MATLAB 自带的其他工具箱里恰好有一个同名函数。用which nrbline -all会列出所有匹配路径如果不止一个需要把你要用的版本的上级目录调到第一位。不需要卸载另一个工具箱只需要调整搜索顺序。注意新版 MATLAB 里 m 文件与小写的差异不会导致冲突但旧版本对大写敏感如果你下载的 zip 内文件名为nrbLINE.m而调用时写小写Windows 系统下可能同样可运行在 Linux/macOS 下则必报错。第三坑是节点向量范围。大多数 NURBS 工具箱默认参数域是[0,1]你在调用nrbeval时如果传入参数t1.2函数不会自动裁剪而是尝试在节点向量之外求值结果会产生NaN或发散值。解决办法是先在调用前用min(t)和max(t)检查范围。如果确实需要循环曲线或周期样条不要自己强行算使用nrbrevolve这类专用函数生成旋转体它们内部已经处理了节点向量的闭合条件。在正式把 NURBS 几何输出到加工代码或仿真环境之前多花 30 秒执行which nrbline -all确认只找到工具箱自带的路径再用 5.1 的单位圆脚本做一次精度校验这两步能帮你挡掉大多数路径和数值层面的低级错误。本文还有配套的精品资源点击获取 MATLAB NURBS工具箱安装与实战：从数据结构到曲线曲面建模
简介MATLAB NURBS工具箱是一套面向MATLAB环境的NURBS建模扩展库适用于需要处理非均匀有理B样条曲线的工程师与科研人员常用于CAD、CAM、CAE及计算机图形学中的复杂几何建模与数据拟合。压缩包共60个文件其中59个为.m函数源码覆盖曲线曲面创建、控制点编辑、参数化、插值拟合、修剪旋转等完整操作流程另有1个txt文件提供配套说明整体仅44KB轻量易部署。目前已有452人学习下载。借助工具箱中的nrbcreate、nrbplot、nrbinterpolate、nrbrevolve等核心函数用户可以快速构建并可视化NURBS曲线与曲面完成从基础绘制到复杂变换的各类几何运算同时支持IGES/STEP等格式导出便于与外部软件衔接。这份资源适合希望系统掌握NURBS理论并需要在MATLAB中落地应用的初学者和进阶开发者。1. 一个 zip 包背后的 MATLAB NURBS 工具箱到底解决什么问题当你从网络或同事手里拿到一个名为 MATLAB NURBS工具箱.zip 的压缩包时很容易产生两个错觉一是解压出来就能用二是它安装后会和 MATLAB 里自带的优化工具箱那样出现在 Add-On 管理器里。实际这两条都不成立。这个 zip 包通常是把第三方编写的 NURBS 建模函数集合打包分发没有安装向导也不会自动注册到 MATLAB 的路径系统。解开后你可能看到 nrbline、nrbcirc、nrbkntins 这些看似神秘的 m 文件如果直接把路径指到某个子目录很多函数会因为找不到对应文件而报错。本文会从 NURBS 的基本数据结构讲起然后带你完成解压、路径配置、曲线曲面的构建与求值最后给出一套验证工具箱是否真正生效的方法。适合需要用 MATLAB 做曲线曲面建模、样条拟合或机器人路径规划的工程师也适合正在做几何设计相关课程作业的学生。2. 先立住理论NURBS 的数据结构就是 MATLAB 代码的地图2.1 B 样条基函数与节点向量最容易读错的参数NURBS 是 Non-Uniform Rational B-Splines 的缩写核心是在 B 样条基础上引入权因子和齐次坐标。要理解工具箱里的每一个函数在做什么必须先看它的三个基本输入控制点矩阵、节点向量和阶次。控制点用坐标列出比如二维曲线的控制点就是 N×2 矩阵三维则是 N×3节点向量knot vector是一组非减实数例如[0 0 0 0.5 1 1 1]它决定基函数在哪里拼接也决定参数域的范围。阶次order在工具箱里习惯用p表示实际次数是p-1比如 order 为 3 意味着二次 B 样条。很多人在第一次写代码时会把节点向量和控制点顺序搞混或者把 order 当成控制点个数这是最典型的错误。因为单元格数组形式的节点向量允许在不同维度的方向上定义不同的次数和拼接方式比如曲面 UV 方向的节点向量可以不同。还要注意同一条曲线同样的控制点如果节点向量分布不一样形状会差别很大。B 样条基函数由 Cox-de Boor 递推公式计算工具箱内部已经封装好你不需要自己写但必须理解参数 t和几何坐标之间的关系。2.2 工具箱用一个 struct 装下整个曲线或曲面NURBS 工具箱的数据结构很有设计感一条曲线或一个曲面就是一个 MATLAB struct。字段通常包含form、dim、order、knots、coefs。form是描述类型的字符串比如line、circ、surfdim是几何空间的维数2 表示二维3 表示三维order是阶数向量knots是 cell 数组每个单元格对应一个参数方向的节点向量coefs是控制点齐次坐标矩阵曲线是 4×(n1)曲面是 4×(n1)×(m1)其中第四行存放的是权重系数。这个第四行的存在正是区分普通 B 样条和 NURBS 的关键。理解coefs的行布局对后续改权重很重要。前三行是 x、y、z 坐标但注意它们是权重乘以原始坐标之后的齐次坐标值不是直接可见的控制点位置。如果你直接改coefs(1,i)你改变的是齐次坐标的 x 分量而最终显示出来的实际控制点坐标会被第四行的权重重新除回来。所以正常修改权重的方法是同时调整前三行或者用工具箱提供的变换函数。下面的代码创建一个直线段 NURBS 曲线然后打印出它的结构字段让你直观看到每个部分长什么样p nrbline([0 0 0], [1 1 0]); disp(p);nrbline用两个端点构造线段返回一个 struct。运行后可以在命令行看到form是linedim是 3order是 2knots是{ [0 0 1 1] }coefs是 4 行 2 列。这个简单的例子能让你验证路径是否配置正确同时熟悉数据结构的形态。如果你在命令行能正常显示这个 struct 而不是报错说明工具箱的核心文件已经可用了。3. 解锁 zip 包用 addpath 正确把工具箱装进 MATLAB3.1 解压与安置该把文件放到哪个目录拿到 MATLAB NURBS工具箱.zip 后右键解压不要直接双击打开看两眼就关掉更不要直接在压缩包内运行 m 文件。很多第三方工具箱都会带多层子目录比如nurbs/下面还有nurbs/或者分成inst、examples、tests等。把它解压到 MATLAB 工作路径之外的一个固定位置比较好比如D:\MatlabTools\nurbs_toolbox或 Linux 下的~/MatlabTools/nurbs_toolbox这样即使切换当前工作目录也不会影响调用。如果你解压到临时目录每次重启 MATLAB 都得重新配置浪费时间还容易忘。解压后用文件管理器确认根目录下存在名为nrbline.m、nrbeval.m、nrbkntins.m等文件。有些压缩包会把这些核心文件放在更深的子目录里比如notoolbox/或src/此时你最好不要挨个手动添加每一个文件而是使用genpath把根目录下所有子文件夹一次性加入搜索路径这才是完整安装的常见做法。3.2 安全添加路径genpath、savepath 与 which 验证打开 MATLAB 后把当前目录切换到解压后的根目录执行以下代码完成路径注册% 假设你解压到了 D:\MatlabTools\nurbs_toolbox toolbox_path D:\MatlabTools\nurbs_toolbox; % 递归获取该目录下所有子文件夹路径 allpaths genpath(toolbox_path); % 逐个添加到MATLAB搜索路径 addpath(allpaths); % 保存到路径定义文件下次启动自动生效 status savepath; disp(status); % 验证核心函数是否可见 which nrbline这段代码的关键点在于genpath它会生成一个以分号分隔的路径字符串把所有子目录都包含进来。addpath可以直接接收这个多路径字符串不需要循环。最后savepath把路径写入pathdef.m如果返回1表示写入失败通常是因为 MATLAB 安装目录没有写权限这时可以改用savepath(C:\Users\你的用户名\Documents\MATLAB\pathdef.m)保存到用户目录。which nrbline是验证安装是否成功的最直接手段正常会输出D:\MatlabTools\nurbs_toolbox\nrbline.m这样的完整路径。如果在which nrbline后仍然显示nrbline not found不要急着怀疑文件损坏先检查解压目录内是否真的存在该文件再确认你设置的toolbox_path是否指向了包含nrbline.m的那一层。常见错误是把路径指到了nurbs_toolbox/的上一级或指到了内含nrbline.m的子文件夹的下一级。3.3 和官方工具箱不同这个工具箱不是从 Add-On 管理器装的很多用户把 MATLAB 自带工具包的安装方式套到这个第三方库上结果在 Home 页面的 Add-On 资源管理器里怎么也搜不到它。这不奇怪NURBS 工具箱不是 MathWorks 官方发布的工具包不能通过 MATLAB 的集成化安装流程来处理。它本质就是一堆 m 文件、p 文件和演示脚本的集合所以安装这个动作在 MATLAB 里对应的就是添加搜索路径。正因如此卸载也很简单直接从路径中移除相关目录即可和matlab下载安装教程里讲的完整安装完全不冲突。同样地这个工具箱和 matlab优化工具箱 也没有任何概念上的重叠。优化工具箱解决的是数学规划问题NURBS 工具箱解决的是几何表达和样条计算问题二者目录结构、函数命名和数据格式都截然不同。如果在使用 NURBS 相关函数时报错说找不到确定类型的输入参数不要去检查官方优化工具箱那说明你的路径配置还没有做到位。4. 动手造曲线曲面从线条、圆到蒙皮样条4.1 基础图元nrbline、nrbcirc 与 nrb4surf配置好之后先用几个基础图元把工具箱跑通。下面这段代码首先生成一条直线段、一段圆弧和一个由四点定义的平面曲面然后绘制出来% 直线段起点(0,0,0)终点(1,1,0) line_nrb nrbline([0 0 0], [1 1 0]); % 圆弧圆心(0,0,0)半径1从0度到90度 circ_nrb nrbcirc(1, [0 0 0], 0, pi/2); % 四点曲面四个角点形成一个平面矩形 surf_nrb nrb4surf([0 0 0], [1 0 0], [0 1 0], [1 1 0]); % 将三条几何对象绘制在一张图上 figure; nrbplot(line_nrb, 20); hold on; nrbplot(circ_nrb, 30); nrbplot(surf_nrb, [20 10]); axis equal; grid on;nrbcirc的第一个参数是半径第二个参数是圆心坐标第三和第四个参数是起始和终止角度单位为弧度。注意如果你写nrbcirc(1, [0 0 0], 0, 90)那就会画出一个 90 弧度的圆弧因为工具箱不会自动把度数转成弧度。nrb4surf要求四个角点按对角顺序传入它会自动生成一片双线性曲面非常适合作为后续拉伸或蒙皮操作的初始片体。nrbplot的第二个参数是采样点数这里给了直线 20 个点、圆弧 30 个点采样越多绘制越平滑但也会增加计算量。4.2 改动几何移动控制点、修改权因子、插入节点NURBS 的一大优势是局部支撑性移动一个控制点只会影响附近区域的形状。下面的例子以一条二阶 B 样条曲线为对象展示三种常见修改操作% 原始控制点经过三个点的二次曲线 ctrl_pts [0 0 0; 1 2 0; 2 0 0]; curve nrbmak(ctrl_pts, [0 0 0 1 1 1]); % 在参数中点插入一个新节点 new_curve nrbkntins(curve, 0.5); % 移动原始曲线的第二个控制点到新位置 curve.coefs(1,2) 1.5; curve.coefs(2,2) 0.5; % 将权因子改为0.2让曲线远离该控制点 curve.coefs(4,2) 0.2; % 绘制对比 nrbplot(curve, 50); hold on; plot(ctrl_pts(:,1), ctrl_pts(:,2), ro);这里nrbmak是直接利用控制点矩阵和节点向量构造曲线返回的curve里coefs前三行默认已经包含了权重为 1 时的坐标值。当你把第二列第四行的权因子改成 0.2 时曲线会被推离那个控制点因为除数变小了。修改coefs(1,2)和coefs(2,2)则是在改变齐次坐标最终的可视化坐标会自动折算。nrbkntins返回新的曲线节点插入不会改变几何形状但会让该参数位置拥有局部修改能力。这个特性在后续做精细调整时非常有用。4.3 求值、求导和可视化从数据到一阶导数曲线建模完成后往往需要把曲线转成离散点供机器人路径规划或加工代码使用工具箱里的nrbeval专门做这件事。下面代码对圆弧在参数点 0 到 1 之间采样并同时计算一阶导数值% 在参数域[0,1]上取100个等距点 t linspace(0, 1, 100); % 对圆弧求值得到欧氏坐标点集 pts nrbeval(circ_nrb, t); % 求一阶导数 dn nrbderiv(circ_nrb); [~, dpts] nrbeval(dn, t); % 绘制圆弧和导数箭头示意 figure; plot(pts(1,:), pts(2,:), b-, LineWidth, 2); hold on; step 10; % 每隔10个点画一条导数线 for i 1:step:100 quiver(pts(1,i), pts(2,i), dpts(1,i), dpts(2,i), 0.1, r); end axis equal; grid on;nrbeval的第一个返回值是曲线上的欧氏坐标第二个返回值是对应的权重值如果你不需要权重可以写成pts nrbeval(circ_nrb, t);但要注意你得到的是一个 3×N 的矩阵。nrbderiv返回一个新的 NURBS 结构体这个结构体的coefs存储的是导数的控制点后续再用nrbeval就可以得到各参数点的切向量。上面代码里用quiver把每个离散点的导数方向画成箭头你会发现切向量方向始终与圆弧在该点的切线方向一致这是验证求导正确性的一个直观方法。下表列出几个最常用函数的作用和参数格式方便你在写脚本时快速查阅函数名典型调用作用返回值nrblinenrbline(p0, p1)创建直线段NURBS 结构体nrbcircnrbcirc(r, center, ang0, ang1)创建圆弧NURBS 结构体nrbmaknrbmak(coefs, knots)由控制点和节点向量构造曲线或曲面NURBS 结构体nrbkntinsnrbkntins(nrb, knots)插入节点或细化控制网格NURBS 结构体nrbdegelevnrbdegelev(nrb, num)提升一个方向的阶次NURBS 结构体nrbeval[pts, w] nrbeval(nrb, t)对参数值求几何坐标坐标矩阵和权重nrbderivdn nrbderiv(nrb)生成导数几何体NURBS 结构体nrbtformnrbtform(nrb, mat)应用变换矩阵移动/旋转/缩放NURBS 结构体nrbdegelev和nrbtform是和nrbkntins并列的高频操作。提升阶次能让曲线在保持形状不变的情况下拥有更多可用自由度而nrbtform接收一个 4×4 的齐次变换矩阵可以一次性把曲线整体平移到新位置比逐点改coefs要稳妥。5. 验证和排错怎样确定工具箱真的跑对了5.1 用单位圆验证权因子是否生效NURBS 与普通样条的最大区别在权因子而权因子最容易在数值上被弄错。下面这段代码用半圆弧验证几何精度如果权因子没生效或者被错误地归一化圆半径不会稳定在 1% 创建从0度到180度的半圆弧圆心原点半径1 circ nrbcirc(1, [0 0 0], 0, pi); % 在参数点0到1上取5个点 t [0 0.25 0.5 0.75 1]; pts nrbeval(circ, t); % 计算这些点到原点的距离 r sqrt(pts(1,:).^2 pts(2,:).^2); % 显示半径 disp(r);正常输出应该是五个都接近 1例如1.0000 1.0000 1.0000 1.0000 1.0000。如果发现某个点半径偏离明显说明coefs中权重行可能被意外覆盖。注意nrbeval返回的坐标已经是把权重除回去的欧氏坐标所以半径计算就是直接取 x² y² 开根号不需要再额外处理。这个方法同时可以验证你在 4.2 节中手动修改coefs(4,:)是否造成了不可预料的形状突变因为单位圆的权重值不是常数半圆弧中间控制点的权重是cos(π/4)量级而不是 1。5.2 三个容易踩的坑路径没生效、函数冲突、节点向量范围第一坑是路径没生效。如果你用addpath添加之后当前绘图能用但重启 MATLAB 后又报Undefined function or variable nrbline说明savepath没有真正保存成功。解决方法是检查savepath返回值如果为 1就在用户目录下创建一个pathdef.m或者直接使用matlabrc开机脚本。最简单的定位办法是每次打开 MATLAB 后先执行which nrbline如果找不到立刻重新addpath(genpath(你的工具箱路径))不要继续往下跑其他的模型代码。第二坑是函数冲突。你可能同时安装了多个版本的 NURBS 工具箱或者 MATLAB 自带的其他工具箱里恰好有一个同名函数。用which nrbline -all会列出所有匹配路径如果不止一个需要把你要用的版本的上级目录调到第一位。不需要卸载另一个工具箱只需要调整搜索顺序。注意新版 MATLAB 里 m 文件与小写的差异不会导致冲突但旧版本对大写敏感如果你下载的 zip 内文件名为nrbLINE.m而调用时写小写Windows 系统下可能同样可运行在 Linux/macOS 下则必报错。第三坑是节点向量范围。大多数 NURBS 工具箱默认参数域是[0,1]你在调用nrbeval时如果传入参数t1.2函数不会自动裁剪而是尝试在节点向量之外求值结果会产生NaN或发散值。解决办法是先在调用前用min(t)和max(t)检查范围。如果确实需要循环曲线或周期样条不要自己强行算使用nrbrevolve这类专用函数生成旋转体它们内部已经处理了节点向量的闭合条件。在正式把 NURBS 几何输出到加工代码或仿真环境之前多花 30 秒执行which nrbline -all确认只找到工具箱自带的路径再用 5.1 的单位圆脚本做一次精度校验这两步能帮你挡掉大多数路径和数值层面的低级错误。本文还有配套的精品资源点击获取