Bézier曲线降阶矩阵方法.pdfVIP

维普资讯 第24卷 第6期 工 程 数 学 学 报 Vo1．24No．6 2007年12月 CHINESEJOURNALOFENGINEERING MATHEMATICS Dec．2007 文章编~-：1005—3o85(2oo7)o6—1007—08 B6zier曲线降阶的矩阵方法术 任水利，张凯院， 叶正麟 (西北工业大学理学院，西安 710072) 摘 要：本文根据升阶的逆过程并结合矩阵代数知识，给出了B6zier曲线降阶的矩阵向量表达式，并且 通过分析得到降阶矩阵实际上是升阶矩阵的广义逆矩阵，同时给出了降阶曲线的误差分析。在某 些情形下，我们得到降阶曲线更好的误差估计。 关键词：B6zier曲线；降阶：广义逆矩阵；降阶矩阵 分类号：AMS(2000165D17 中图分类号：TP391 文献标识码：A 1 引言 在 CAD 和 CAM 系统 中，为实现不 同阶数的曲线 曲面 的数据转换 ，提高计算效率 和稳定性 ，经常 需要对 曲线 曲面进行 降阶处理 。B6zier曲线是几何造型 中常用 的工具 之一 ，关于 B6zier曲线 曲面 的降阶方法 ， 目前相对而 言 已有较为成熟 的理论和方法 。 早在 1972年，Forrest[1】曾对 B6zier曲线的降阶做过研究；进入80年代后，有更多的学者 对 B6zier曲线曲面的降阶算法作出了贡献，如 Farin等[2-s】从各个方面对该问题进行了探讨， 都取得了一定的研究成果。国内的学者对这一领域的研究也作出了相应的贡献，主要成果有； 胡事 民[99]等从 B6zier曲线的退化条件出发，利用扰动控制顶点和约束最小二乘法来实现降阶， 并将这种方法推广到其他曲线 曲面降阶 问题；陈国栋和王国瑾 1【O,ll】利用 B6zier曲线本身的 几何性质并结合广义逆矩阵、最佳一致逼近和最小二乘理论，给 出了3种不同的降阶方法， 同时讨论了保端点插值 问题；康宝生 [12】等应用 B6zier曲线的几何性质和 B6zier曲线的升阶 公式，基于遗传算法，给出了有理 B6zier曲线降阶的新算法。本文 以矩阵向量的形式给出 了B6zier曲线降阶的表达式，通过分析得到降阶矩阵实际上是升阶矩阵的广义逆矩阵，同时 给出了降阶曲线的误差分析。 2 预备知识 为了用矩阵向量形式研究降阶过程，我们引入广义逆矩阵的概念 定义1 设矩阵A ∈Cm ，若矩阵 ∈C 满足如下四个 Penrose方程 AXA=A，XAX= ， (AX)H=Ax， (XA)H= ， 则称 为 的Moore-Penrose逆，记为 +。 定义2 对任意A∈C ，若 ∈C 满足Penrose方程中的 ()，(J)，… ，(f)等方程， 则称 为 的 {i，J，… ，f}一逆，记为 (，’，)，其全体记为A{i，J，… ，f)。 收稿 日期：2005—12—22．作者简介：任水利(1975年7月生)，男，博士研究生．研究方向：CAGD，数值代数 基金项目：国家自然科学基金． 维普资讯 1008 工 程 数 学 学 报 第24卷 特别地，引入如下定义 定义3 对任意 A∈C ，若 X ∈C 满足 Penrose方程中的第一个方程，即A = ， 则称 为 的广义逆，记为 一。 我们主要用到定义3所定义的广义逆。 一 般地，B6zier曲线的降阶问题可描述如下： 设 {只)0为给定的一组控制顶点，则可以确定一条n次B6zier曲线 pn(t)=∑ ()，t∈0[1] i=0