基于几何区间裁剪求交算法-计算数学专业论文.docxVIP

插 插 图 插 图 2.1 Be′zier曲线完全包含在凸包内 5 2.2 de Causteljau算法 5 2.3 第二分割算法,其中有三个patch被标记 8 将三角Be′zier曲面分割成5个子patch 8 将三角Be′zier曲面分割成三个子patch 9 2.6 第三分割算法,其中只有中间的子patch被标记 9 2.7 将三角Be′zier曲面分割成七个子patch 10 2.8 将三角Be′zier曲面分割成四个子patch 10 2.9 Be′zier Clipping:控制多边形的凸包和t轴的交 11 2.10 二次裁剪算法的一次迭代过程 . . . . . . . . . . . . . . . . . . . . . 12 3.1 几何区间裁剪曲线求根收敛阶证明 . . . . . . . . . . . . . . . . . . 19 3.2 Tτ,g , Y 和Z之间的角 20 3.3 例1: 计算时间(10?5s)和精度e 23 3.4 例2: 计算时间(10?5s)和精度e 25 3.5 例3: 计算时间(10?5s)和精度e 25 4.1 三角域上参数区域 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30 4.2 三角曲面向直线Lr 投影 32 4.3 3次三角曲面沿着Lr 投影后控制顶点的分布图 33 4.4 所有的最大值和最小值分别确定n次的Be′zier曲线 34 4.5 光线/三角曲面分割示意图 . . . . . . . . . . . . . . . . . . . . . . . 34 XI 中国科学技术大学学位论文原创性和授权使用声明 本人声明所呈交的学位论文，是本人在导师指导下进行研究工 作所取得的成果。除已特别加以标注和致谢的地方外，论文中不包 含任何他人已经发表或撰写过的研究成果。与我一同工作的同志对 本研究所做的贡献均已在论文中作了明确的说明。 本人授权中国科学技术大学拥有学位论文的部分使用权，即： 学校有权按有关规定向国家有关部门或机构送交论文的复印件和电 子版，允许论文被查阅和借阅，可以将学位论文编入有关数据库进 行检索，可以采用影印、缩印或扫描等复制手段保存、汇编学位论 文。 保密的学位论文在解密后也遵守此规定。 作者签名： 年 月 日 第 1 第 1 章 绪论 第 1 第 1 章 绪论 PAGE PAGE 11 PAGE PAGE 10 第 1 章 绪论 曲线曲面的求交一直以来都是计算机辅助设计和制造领域中最重要的问题 之一。 寻找快速而且稳定的求交算法成为了很多研究者的关注课题。 本章主要 回顾了曲线、 曲面求交算法的发展。 首先回顾了方程求根和曲线求交算法，然 后列出了光线/曲面求交算法及其应用技术，最后给出了本文的内容结构。 1.1 方程求根和曲线求交 下面回顾一些曲线求交算法，诸如拟牛顿法，隐式化法，细分法和Be′zier Clipping，二次裁剪，几何区间裁剪算法。 拟Newton法是为了克服Newton法需要求Hesse矩阵的缺点发展起来的求解 根的有效方法。 拟Newton法就是利用目标函数值及其一阶导数信息，给出目标 函数曲率的近似值。 设F (x) = 0，假定F ′(x)在点xk+1处连续，从而对任意的ε 0，存在δ 0， 使得当∥x ? xk+1∥ δ时，有 ∥F ′(x) ? F ′(xk+1) ? Gk+1(x ? xk+1)∥ ≤ ∥s ? xk+1∥ k+1令x = xk，Hk+1是G? k+1 的近似，从而Hk+1(F ′(xk+1) ? F ′(xk )) = xk+1 ? xk 是合 理的。 则拟Newton迭代算法为： [ Hk+1 = Hk + 1 +  kyT Hk yk ] k sTk yk · sT  kksk sT sT yk k k  ksk yT Hk k ksT yk k  (1.1) 初始矩阵H0一般取为单位矩阵I，然后通过公式1.1修正Hk。综上所述，拟Newton 法较Newton法有以下优点：仅需要一阶导数，Hk 保持正定，使得方法具有下降 性质，而且每次迭代仅需O(n2)次乘法。 隐式化法是用于求曲线交点的代数方法。 首先将一条曲线的参数方程隐式 化，然后将另一条曲线的参数方程代入到隐式方程中得到一个多项式。 从而曲 线的求交问题转化为多项式的求根问题。 实验结果表明，对于低次多项式，隐 式法是快速的求交方法。 但是对于大于5次的曲线在实际中隐式法无效且数值不 稳定。 分割算法充分利用Be′zier曲线的凸包性质。 如果两