基于刘徽割圆术的等距线逼近算法-浙江大学数学系.pdfVIP

基于刘徽割圆术的等距线逼近算法* 寿华好 刘利刚 王国瑾 (浙江大学数学系 杭州 310027) 摘 要 本文给出了一种基于刘徽割圆术的平面 NURBS 曲线的等距线的逼近算法．利用正多边 形代替圆所扫掠出的区域边界来近似等距曲线，所得到的逼近曲线是与基曲线同次的 NURBS 曲 线，并且可以达到任意的精度． 关键词 NURBS 曲线，等距曲线，曲线逼近． 分类号：TP391 1 引 言 等距曲线曲面(Offset curve/surface)作为 CAD/CAM 系统中的一个重要几何运算在数控领域 有直接的应用，此外在机器人路径规划，实体造型，计算机图形学，形位公差学，公路铁路设计， 带厚度物体如箱包的设计，以及艺术花纹设计等领域中有广泛的应用．然而等距运算是一个非常 困难的几何运算，与基曲线曲面相比，等距曲线曲面有非常高的代数次数，比如一条3 次 NURBS [1] 曲线的等距曲线的代数次数高达 10 次 ． 给定平面正则参数曲线C(t) = (x (t), y (t)) ，距离为d 的等距曲线Cd (t) 为 Cd (t) = C(t) ± d ⋅ N(t), ′ ′ (y (t),−x (t)) 其中N(t) = 为C(t) 的单位法向量．由于法向量N(t) 的分母出现根式，一条有 ′ 2 ′ 2 x (t) + y (t) 理 NURBS 曲线的等距曲线一般不再是一条有理曲线，从而无法被通用的CAD/CAM 系统进行有 效的处理．鉴于以上困难，目前用较简单的曲线曲面去逼近等距曲线曲面仍是十分重要且切实可 行的办法． 关于等距曲线的逼近算法已有不少工作，它们大致上可以分为三类: (1) 基于移动控制多边形 或控制顶点的方法．如 Tiller 和 Hanson[2]是用平移控制多边形的方法；而 Cobb[3]是用移动控制顶 点的方法；Coquillart[4] 的方法则是 Cobb 的方法的改进．(2) 基于插值或拟合的方法．如 Klass[5] 用 3 次Hermite 曲线逼近等距线； Hoschek[6]提出了在最小二乘意义下的等距线逼近方法； Pham[7] 用 3 次均匀 B 样条插值的方法产生逼近等距线；Sederberg 和 Buehler[8]用区间算术和 Hermite 插 值的方法逼近 Bézier 曲线的等距线；最近 Piegl 和 Tiller[9]又提出了一种基于插值的 NURBS 等距 线逼近方法．(3) 基于逼近基圆的方法．Lee 等[10]提出了一种先用 2 次 Bézier 曲线段逼近基圆再 与基曲线卷积产生逼近等距线的方法． 平面曲线的等距线问题也可以看成是基曲线与基圆的 Minkowski 和的边界问题．Lee 等人[10] 的方法很有特色地解决了整体误差问题，然而该方法按 Piegl 和 Tiller[9] 的评论有两个主要缺点： (1)等距逼近曲线有大量的控制顶点；(2)等距逼近曲线的次数很高，即对r 次多项式曲线而言，逼 近曲线是3r − 2 次的有理曲线；对r 次有理曲线而言，逼近曲线是5r − 4 次的有理曲线，逼近曲 线的构造比较复杂．鉴于目前的计算机具有处理大量控制顶点的能力，所以最关键的问题是如何 把逼近曲线的次数降下来． 本文提出了一种基于刘徽割圆术的等距曲线逼近算法．首先用正多边形来逼近等距圆，将正 多边形沿着基曲线的扫掠区域边界作为等距曲线的逼近曲线，所得到的逼近曲线是与基曲线同次 * 本文受国家自然科学基金(No.,国家重点基础研究“973 ”项目基金(No. G1998030600), 浙江省自然 科学基金(No.698025)和浙江省教委科研基金资助. 数的 NURBS 曲线，逼近曲线的构造十分简单，并且可以达到任意的精度． 2 NURBS 等距曲线的逼近算法 我国古代数学家刘徽[11] 曾利用割圆术(即用内接正6n 边形来逼近圆) 的方法来计算圆周率π ， 祖冲之将其推算到小数点后 6 位，