关于控制多边gerbs曲线拐点的研究.docxVIP

关于控制多边gerbs曲线拐点的研究 nurbs作为cal形式的基本工具，自出版以来，它引起了人们的广泛关注和研究，尤其是对转折点的分析和计算。1992年,Manocha与Canny 提出了利用曲率公式计算真参数化(Proper Parametrization)参数曲线及有理参数曲线拐点的数值方法;1995年,Goodman等人研究了NURBS曲线变差缩减性质,得到曲线拐点个数与其控制多边形拐点个数之间的一般关系;1997年,Li与Robert利用代数曲线的切线与曲线交点的重数关系,研究了判定齐次坐标空间曲线拐点的方法,得到了计算有理参数曲线拐点的代数方程法;1999年,Sakai又借助于Mathematica,得到了利用端点斜率判定三次有理曲线段上拐点的方法;2001年,Liu Chaoyang又对NURBS保凸性问题进行了研究,得到了NURBS曲线的凸性定理。然而不难发现,为满足CAD工程实践的具体需求,尚需要对NURBS拐点的分析判断作更加深入的研究。 CAD工程师主要通过交互调节曲线的控制顶点或者权因子来控制曲线的形状,经常需要判断所设计曲线的凸凹性。根据Liu Chaoyang的凸性定理,凸的控制多边形将产生凸的NURBS曲线。但是实践表明,非凸的控制多边形也可以产生凸的NURBS曲线。将控制多边形(权因子)调节到什么程度曲线会出现拐点是交互设计中非常关心的问题。单凭设计者的经验来控制显然是不科学的。因此,需要研究拐点出现的条件以设计适合计算机处理的自动判断算法,这对于CAD工程十分重要。 作者在文献通过引入转向矢量,给出了平面曲线拐点的一个数学定义,并演绎推导出曲线在拐点处的若干性质。本文在该文献基础上,首先从变形的观点出发,研究了一般平面参数曲线变形中的几何特征变化规律,得到了曲线形状变化过程中产生拐点的充分必要条件以及确定拐点的一个泛函方程;同时,结合NURBS的特殊要求,研究了非凸多边形NURBS曲线出现拐点的条件和判据。本文的结果可应用于调节与控制曲线形状的算法设计,也可指导CAD造型工程师的交互设计。 1 曲线变形及其几何特征的变化 1.1 平面正则曲线的深度 定义1设有向平面π的单位法矢量为B,a、b是π上两个矢量。a在π上绕其起点按照π的定向法则正向转动到与b的方向相同时扫过的有向角为a,b,则a,b的矢量积(叉积)为: a×b=B|a||b|sina,b 容易验证,定义1是经典矢量积在平面上的归一化表示,它将经典平面矢量叉积的方向与数量关系统一在一个完整的数学式中。 定义2平面曲线的拐点、转向矢量。见文献。 下面的引理1～3亦可由文献的相关理论直接推出。 引理1设平面连续曲线r(t)具有分段C1连续的单位切矢量,K(t)是与r(t)对应的转向矢量;那么r(t)在某区间的拐点个数等于K(t) 在该区间的方向改变次数。 引理2如果两条同源平面正则曲线r1(t)、r2(t)的切矢、主法矢分别为T1(t)、T2(t)、N1(t)、N2(t)。若在t*点满足0≤(T1(t*),T2(t*))π,N2(t*)在N1(t*)上的投影与N1(t*)的方向相反,则r1(t)、r2(t) 在t*点的回转方向是相反的。 引理3假设平面正则曲线r(t)在起点与终点的K矢量有相同的方向,那么r(t)的拐点个数为偶数;若r(t)在起点与终点的K矢量有不同的方向,则r(t)的拐点个数为奇数。 注意:引理1中r(t)具有分段C1连续的单位切矢量是一个非常重要的条件,故在后文中如无特殊说明,都假定所讨论的曲线满足该条件。 1.2 曲线变形与法矢、转向矢的变化 参数曲线r1=r1(t)在变形位移(偏移量)r2=r2(t)的作用下,变形为r(t),其数学表示为: r(t)=r1(t)+r2(t)t∈[a,b] (1) 变形往往导致曲线的凸凹性、连续性等发生变化。依据上节分析,不难得到以下结论: 定理1设平面凸曲线r的两个端点位置固定不变,r上一点P变形前后的切矢与主法矢分别为T1、T2与N1、N2;如果T1、T2满足0≤(T1,T2)π,而N2在N1上的投影与N1的方向相反,则变形后r失去原有凸性(出现拐点或者没有拐点但改变原有回转方向)。 证明r的两个端点固定意味着曲线不能有整体转动,这就避免了变形过程中r绕一个端点刚性旋转而导致N1,N2大于π的情形。将变形前与变形后r的两种状态看成是两条同源曲线,根据引理2知,在P点附近变形前后的转向矢量方向是相反的。此时若r上另外一点Q变形前后的转向矢量方向不变,则根据引理1知r变形后出现拐点;若r上所有点变形后的转向矢量方向都不同于变形前,则曲线整体转向发生改变。故定理成立。[证毕] 上述基于法矢方向变化的方法,只能定性地分析曲线变形前后凸凹性的变化,无法与位移量之间建立数学上解析