种插值与逼近运动物体活动标架的新方法.docVIP

一种插值与逼近运动物体活动标架的新方法 姜忠鼎1 马利庄1,2 1(浙江大学CADCG国家重点实验室 杭州 310027) 2(Centre for Advanced Media Technology Nanyang Technological University 639798 Singapore) Email: {zdjiang, mlz}@cad.zju.edu.cn 摘 要 插值运动物体给定活动标架（原点位置和三个坐标轴朝向）是计算机图形学、机器人等领域的一个重要问题.本文提出一种采用B样条曲线插值与逼近运动物体活动标架的新方法.采用4个欧拉参数对正交活动标架的旋转变换矩阵进行参数化,得到了一个简单的优化方程.采用迭代法求解该优化方程来逼近运动轨迹任意处伪活动标架的旋转变换矩阵.本文也证明了插值和逼近所引起的误差是可控的.由于活动标架的计算只涉及2阶或3阶线性方程组,该方法具有很高的运行效率. 关键词 运动设计, 插值与逼近, 活动标架, B样条 中图法分类号 TP391 在计算机图形学、机器人等领域,人们通常设定运动物体的若干关键位置即一组物体所在的局部坐标向量,对这些关键位置进行适当的插值或逼近,以确定运动坐标系的一般位置,进而可确定物体的一般位置.物体的朝向一般可用欧拉角来表示,因此朝向的插值问题可简单地转化为3个欧拉角的插值问题.但欧拉角又有它的局限性.因旋转矩阵是不可交换的,欧拉角的旋转一定要按某个特定的次序进行；等量的欧拉角变化不一定引起等量的旋转变化,导致了旋转的不均匀性；欧拉角还有可能导致自由度的丧失,即所谓的gimble lock现象.Shoemake为解决因采用欧拉角表示引起的麻烦,最早将四元数引入计算机动画领域,并提出用单位四元数空间上的Bezier样条来插值四元数[1].Barr等人提出了一个采用四元数对带有角速度约束的景物的朝向进行光滑插值方法,他们的方法允许对轨迹端点处的角速度进行约束.首先把关键帧处的旋转角转化为四元数,然后把非欧几里德空间四元数路径的切向加速度取极小,并用有限差分和最优化方法求解得到的能量方程[2].Kim通过构造一组新的基,提出了把R3空间曲线变换到单位四元数的一般性方法[3]. 在CAD/CAM、机器人等领域常要求插值方法必须保证所得运动坐标系既是单位正交,又为NURBS曲线表示形式.最近这一问题引起学者们的广泛兴趣.Roschel于1998年对运动物体的有理运动表示方法作了综述[4].许多商用机器人控制系统采用分段线性或分段园弧插值方式.它简单实用,但也存在明显缺点,如物体的运动在连接处很不自然[5].为此Wang和Yang引入了一种近似弧长参数化的5次样条曲线[6]；Farouki和shah则引入了一种所谓的PH曲线[7]；但这两种方法中曲线活动标架的表达式并非真正的NURBS表示形式.为使活动标架具有精确的NURBS表示形式,Ge和Ravani引入了刚体运动群SE（3）的单位四元数表示方法[8],本质是将R3中正交变换表示成四元数,然后将给定运动位置的插值问题转化为四元数的插值问题.此后又有作者发展了相应四元数插值方法[9,10].采用四元数,活动标架最终可表示成NURBS形式,但表示形式仍复杂.例如,若利用次NURBS曲线对局部坐标原点路径采用四元数方法插值,则最终的活动标架为次NURBS形式.为降低表示次数,提高计算效率,有必要对此作进一步研究,以满足CAD虚拟环境中快速计算的需要.此外,仍缺乏一个好的误差估计公式,而这又是实际应用中必需的. 本文提出了一种用B样条曲线统一插值和逼近活动标架坐标原点和三个坐标轴朝向的新方法.该方法并不精确计算物体在运动时所对应的正交活动标架,而是先得到一个伪活动标架,然后求得最佳逼近该伪活动标架的正交活动标架.同时我们还证明了逼近误差是可控的.采用了逐步求精的迭代技术,可根据需要得到更精确的正交活动标架.迭代求解计算只涉及2阶或3阶线性方程组,该方法具有很高的运行效率. 本文内容安排如下.第1节对活动标架插值问题进行了描述并给出我们用B样条曲线统一插值活动标架原点和三个坐标轴朝向的方法；第2节证明了活动标架的逼近误差是可控的；第3节给出逼近伪活动标架的方法和该算法的一个运行实例；最后是本文的结论. 插值运动物体活动标架 令表示三维欧几里德空间.采用Cartesian坐标架和来分别表示固定标架和活动标架.固定标架中的点表示为,活动标架中的点表示为. 设 为活动标架通过变换（1.1）得到的 个相关位置. (1.1) 其中为平移向量,为3(3 正交矩阵.插值问题可以描述为寻找一参数运动 满足