基于路径跟踪的同伦迭代法的几何约束求解方法地研究.pdfVIP

第五届中国计算机图形学大会，2004年9月23日～26日，西安 基于路径跟踪的同伦迭代法的几何约束求解方法的研究木 曹春红卢奕南李文辉 (吉林大学计算机科学与技术学院，长春 130012) 摘 要：几何约束求解问题等价于求解一系列的非线性方程组的问题。非线性方程组可以通过经典的 Newton·Raph∞n方法来求解。路径跟踪是牛顿算法的迭代应用。一般的同伦法求解非线性方程组，要比 一个合适的初值来保证其局部收敛，而一般对初值的选取带有很大的盲目性，这就使方法的有效性受到制约。现在 借助于同伦函数，很好的解决了这个问题。基于跟踪方案的同伦迭代法适用于求解一切多项式方程组。由于在路径 跟踪的每一步都去掉了预估这个环节t而且其划分的等分数较少，所以计算效率比一般连续法高，计算复杂度也比 一般连续法小。 关键词：几何约束求解；欠约束；同伦迭代方法；路径跟踪 1引言 L J 几何约束求解问题是当前基于约束设计研究 定初值础后按照迭代公式逐次计算新值墨1， 中的热点问题。一个约束描述了一个应该被满足的 II。 ¨ 关系，一旦用户已经定义了一系列的关系，那么在 修改参数之后，系统会自动选择合适的状态来满足 bⅣ。。一般来讲，这种方法的求解速度快，但是 约束，这种思想方法叫做基于约束的模型。系统解 需要一个好的初值而且对初值的选取比较敏感。当 决约束的部分叫约束求解器。Ⅲ国内外很多学者运 初值发生微小的变动的时候，用这种方法求解可能 用数值计算理论、人工智能理论、图论、自由度分 会发散或者收敛到一个用户不想要的解。文献【4] 析理论等对约束求解进行了深入的研究，几何约束 提出的区间分析法虽然能在较大范围内收敛，可判 问题求解的方法发展到现在，大致可以分为三类： 断在给定区间内有无解，并对有解区间求出其中的 基于代数的几何约束求解方法、基于规则的几何约 全部解，但它仍存在选择合适的初值区间的问题。 束求解方法和基于图论的几何约束求解方法。 同伦方法在一定程度上可以弥补这个缺点口“’7】。同 基于代数方法的几何约束求解方法，简单的 伦方法是一种有效的数值迭代方法，它有较强的整 说，是指首先建立一个坐标系将几何约束问题转化 体收敛性，能可靠的求出方程组的全部或一批孤立 为一个代数方程组，然后在用数值计算的方法和符 解。但是，几何约束问题得到的方程组往往是欠约 号代数的方法求解这个代数方程组。相应地，我们 束的，为此在应用同伦方法求解时，要采用齐次同 将基于代数方法的几何约束求解方法分为数值的几 伦法或系数同伦法，需要对问题进行一些前处理， 何约束求解方法和符号的几何约束求解方法。对非 这不仅使求解过程复杂化，而且要个别问题个别处 线性方程组的求解，可采用符号和数值求解两种方 理，降低了求解过程的通用性。对较大规模和高度 式。符号求解可以求出问题的封闭形式的解析解， 亏欠的问题，解题的计算效率也限制了方法的实际 当然是最理想的，但其难度往往也是很大的。随着 应用。 问题数学模型的增大，消元过程变得愈加复杂【2】。 基于同伦连续法的思路，我们提出了一种将同 为进行数值求解，最普遍的是采用牛顿迭代法， 伦函数与传统迭代法结合起来的新数值迭代法，我 们称之为同伦迭代法。类似于常规的同伦方法，同 比如GossardD]提出的方法。设约束方程组为^(xD