论文部分内容阅读
曲面的光滑拼接和有理参数化是计算机辅助几何设计中的两个基本问题.构造过渡曲面来光滑地连接两个或者多个实体模型这一过程称为拼接.由曲面的隐式代数表示转换成有理参数表示这一过程称为有理参数化.本文主要研究含参代数曲面族的光滑拼接和有理参数化.所谓含参代数曲面族是指由含参数的多项式的零点集定义的代数曲面族.令R表示实数域,X：={x,y,z}是三个未定元构成的集合,（?）：={∈1,,εm}是有限个参数构成的集合.记R[（?）][X]：=（R[（?）]）[X]是多项式环R[（?）]上的多项式环.（?）E∈Rm,定义由∈导出的标准特定化同态σ∈：R[（?）][X]→R[X].R上所有闭区间构成的集合记为IR,m个IR的笛卡尔乘积记为IRm.定义1（?）f∈R[（?）][X],如果每个∈i均在各自的定义区间内连续变动,即∈i∈Ii∈IR,i=1,,m,或者将其写成向量形式,即（?）∈（?）：=I1×I2×…×Im∈IRm,那么我们可以得到一族多项式{σ∈（?）∈R[X]：E∈3}.令V（?）：={V（σ∈（f））：E∈3},这里V（σ∈（f））表示R[X]中的多项式σ∈（f）的零点集,我们称V（?）是由f定义的含参代数曲面族.含参代数曲面族模型可以用来描述计算机辅助几何设计中的很多问题.例如：一方面,在众多工程应用中,人们获取的曲面数据通常是含误差的.然而,对于很多奇异性问题而言,微少的扰动亦能使问题的解发生本质性的变化.这时,我们可以使用含参代数曲面族表示含误差的曲面,那么我们只要讨论这个新的含参问题的关于其参数稳定/连续的解即可.另一方面,在几何造型中,人们经常会遇到类似于贝壳结构的物体,即物体是由带厚度的曲面片构成的.含参代数曲面族亦可以表示带厚度的曲面,其厚度可以通过调整各参数的定义区间来控制.与代数区间曲面模型相比,含参代数曲面族模型的范围更广,并且不会出现由于区间运算所导致的“膨胀”；与等距代数曲面模型相比,含参代数曲面族模型的计算更加简单,但是如何自由地控制其厚度是一个比较困难的问题.本文的主要工作可以分为下面三个部分.（1）含参代数曲面族的光滑拼接.2010年,林群和]R,okne讨论了代数区间曲面的拼接问题.然而,他们未能给出代数区间曲面的几何连续性的精确定义,并且由于区间运算会导致运算结果发生“膨胀”,他们构造的拼接曲面比输入曲面“厚”.本文的模型含参代数曲面族—可以有效地改进上述问题.首先,我们将Warren在1989年提出的代数曲面的rescaling连续的定义推广到含参代数曲面族情形.定义2假设含参代数曲面族V（（?））和V（（?））相交于一不可约的含参代数曲线族C,这里f,g∈R[（?）][X].称V（（?））和V（（?））沿着C Gk rescaling连续,如果1.（?）E∈（?）,存在v∈￡,使得V（σ∈（?））和V（σv（g））沿着其公共交线Gk rescal-ing连续；2.（?）v∈￡,存在E∈（?）,使得V（σv（g））和V（σ∈（f））沿着其公共交线Gk rescal-ing连续.上述定义保证了我们构造的拼接曲面族和输入曲面族两者的厚度恰好“匹配”.接下来,我们将拼接理论中的理想论方法推广到含参代数曲面族情形.定理1假设含参代数曲面族V（g（?））和代数曲面V（h）相交于一含参代数曲线族C,这里g∈R[（?）][X],h∈R[X].（?）p1,p2∈R[（?）][X],我们令f=p1g,p2hk+1∈<g,hk+1>（?）[（?）][X],那么含参代数曲面族V（f（?））和V（g（?））沿着C Gk rescaling连续.人们在处理曲面拼接问题时,通常尽可能选择低次数的拼接曲面.然而,上述定理并未给出计算低次数拼接曲面族的方法.所以,下面我们使用全局Grobner系统方法,提出计算多个含参多项式理想的交的算法IMPI,这里算法SS是计算全局Grobner系统的Suzuki-Sato算法.算法1IMPI（J1,,Jr）输入：R[（?）][X]中的r个理想Ji=,这里Fi是R[（?）][X]中的有限个多项式构成的集合,i=1,,r.输出：R[（?）][X]中理想（?）Ji的一个全局Grobner系统9.IMPRI1:令P←{1-（?）ti,t1F1,,trFr},再将P视为R[（?）][t1,,tr,X]的子集.IMPI2:令（?）←SS（P,（?））,这里（?）是满足{l1,lr}》X的分块序.IMPI3:令（?）←（?）.对于每个（Si,Ti,Gi’）∈g’,令Gi←{g:g∈Gi’∩R[（?）][X]}, g←g∪（Si,Ti,Gi）.IMPI4:返回9.最后,我们将算法IMPI应用到含参代数曲面族的拼接问题中.假设g1,gr是R[（?）][X]中互异的多项式,h1,.,hr是R[X]中的多项式,且满足V（gi,（?）i）和V（hi）相交于一含参代数曲线族Ci,i=1,r.根据定理1,如果f∈R[（?）][X]满足f∈<g1,h1k+1>∩…∩<gr,hrk+1>（?）R[（?）]+[X],那么含参代数曲面族V（f（?））和V（gi（?））沿着Ci Gk rescaling连续,i=1,.,r.我们计算IMPI（<g1,h1k+1>,,<gr,hrk+1>）,得到(g1,h1k+1>∩…∩<gr,hrk+1>的一个全局Grobner系统g={（S1,T1,G1）,,（Sl,Tl,Gl）}这时,参数空间Rm被分割成l个部分V（Si）\V（Ti）,i=1,,l.（?）∈∈V（Si）\V（Ti）和.f∈<g1,h1k+1>∩…∩<gr,hrk+1>,σ∈（f）均可表示成σ∈（Gi）的多项式组合.所以,这时我们必须对（?）再做一次判断,即是否存在某个V（Si0）\V（Tiv）,满足（?）1∪…∪（?）r（?）V（Si0）\V（Ti0）.如果上述条件成立,那么由任一f∈(Gi。>定义的含参代数曲面族V（f（?）..（?）r）均为所求；否则,算法失败.（2）单参二次曲面族的有理参数化.称一个含参二次曲面族是单参二次曲面族,如果其定义多项式中仅含一个参数并且该参数仅出现在多项式的某一项中.假设V（fI）是由关于x,y,z,ω的二次齐次多项式f（x,y,z,w;∈）定义的单参二次曲面族,这里X：=（x,y,z,ω）T是曲面族上点的齐次坐标,∈是参数.我们将全体单参二次曲面族分成两类：∈出现在f的ω2项中的和∈出现在f的zω项中的.显然,对于∈出现在f的其他项中的单参二次曲面族,我们可以通过简单的变量替换将其转化成上面形式之一基于经典的球极投影法,我们提出了参数化V（fI）的标准型方法算法CM.算法2CM（f（x,y,z,w;∈））输入：V（fI）的定义多项式f（x,y,z,w;∈）.输出：V（fI）的一个有理参数表示.CM1：将f（x,y,z,w;∈）写成形如XTA∈X的向量形式.CM2：（?）∈I,计算相应的射影变换X=P∈X’:=Pe（x’,y’,z’,w’）T,这里P∈是含∈的非奇异矩阵,使得X’T（P∈TA∈P∈）X是下面三种标准型之一c1x’2+y’2+z’2-w’2; c2x’2+y’2+z’2-w’2; c3x’2+y’2-z’2CM3：使用球极投影法,分别计算标准型C1,C2,C3的有理参数表示PC1, PC2,PC3（不唯一）CM4：（?）∈I,将其对应的标准型Ci的有理参数表示,PCi代入到相应的射影变换X=P∈X’中即可.然而我们发现,使用算法CM计算出的有理参数表示有时会在某一点处出现“不连续”现象.这种“不连续”现象之所以出现,是由于我们在化简系数矩阵A。时,将∈与定义多项式中其他项的系数“一视同仁”.实际上,由于∈是参数,我们应当将∈与其他项的系数区别开来,对其单独处理.基于上述想法,我们改进了标准型方法,改进的标准型方法可以修复这种“不连续”现象.下面我们仅给出参数化∈出现在f的ω2项中的V（fI）的算法SCM.算法3SCM（f（x,y,z,w;∈））输入：V（fI）的定义多项式f（x,y,z,w;∈）.输出：V（fI）的一个关于∈连续的有理参数表示.SCM1：将f（x,y,z,w;∈）写成形如XTA∈X的向量形式.SCM2：计算射影变换X=PX’:=P（x’,y’,z’,w’）T,这里P是不含∈的非奇异矩阵,使得X’T（PTA∈P）X’是下面五种标准型之一SC1x’2+y’2+z’2+l（∈）w’2; SC2x’2+y’2-z’2+l（∈）w’2; SC3x’2+y’2+2a34z’w’+l（∈）w’2; SC4x’2-y’2+2a34z’w’+l（∈）w’2; SC5x’2+2a24y’w’+2a34z’w’+l（∈）w’2.这里l（∈）是关于∈的线性多项式.不失一般性,我们假设a34是不等于零的实数.SCM3：使用球极投影法,计算V（fI）对应的标准型SCi的有理参数表示SPCi（不唯一）.SCM4：将V（fI）对应的SPCi代入到射影变换X=PX’中即可.（3）非奇异三次拼接曲面的几何信息和有理参数化.非奇异三次拼接曲面亦可视为特殊的含参三次曲面族.2003年,王相海讨论了如何参数化非奇异三次拼接曲面,然而其计算的参数表示不是有理的.2006年,伍铁如和程宏路讨论了由f=（b1（y-d2）+b2（x-d1））（x2+y2+z2-r2）-（x-d1）（y-d2）定义的特殊非奇异三次拼接曲面的有理参数化,这里b1,b2,d1,d2,r是参数.上述工作在本文中得到了进一步地推广和深化.我们将非奇异三次拼接曲面分成两类：特殊形式曲面和一般形式曲面.对于特殊形式曲面,我们分析了其几何信息,并且证明了定理2特殊形式曲面只能是Fi（i=3,4,5）曲面.为了进一步判断特殊形式曲面的曲面类型,我们使用杨路提出的多项式判别系统方法得到了下面定理.定理3特殊形式曲面的D（λ）的判别式序列[D1,D2,D3,D4]为[1,-8b4b2+3632,16b42b2b0-18b24b12+14b4b3b3b1-6b4b32b0-4b4b23-3b33b1+b32b22,256b43b03-192b24b3b1b02+144b25b2b12b0128b24b22b02-27b42b14+144b4b32b2b02-6b4b32b12b0-80b4b3b22b1b0+18b4b3b2b13+16b4b24b0-4b4b23b12-27b34b02-18b33b2b1b0-4b33b13-4b32b23b0+b32b22b12],这里bi是D（λ）的λi项的系数,j=0,1,,4.1.如果下面条件之一成立：D2<0∧D3<0∧D4>0;D2≥0∧D3≤0∧D4>0;D2<0∧D3≥0;D2=0∧D3>0,那么特殊形式曲面是F3曲面.2.如果下面条件之一成立：D2≤0∧D3<0∧D4≤0;D2=0∧D3=0∧D4<0;D2>0∧D3<0∧D4=0;D2>0∧D4<0,那么特殊形式曲面是F4曲面.3.如果下面条件之一成立：D2>0∧D3>0∧D4≥0;D2≥0∧D3=0∧D4=0,那么特殊形式曲面是F5曲面.这里“八”表示逻辑交,即A八B成立当且仅当A和B同时成立.根据Berry和Patterson在2001年提出的非奇异三次曲面的有理参数化算法,我们得到了特殊形式曲面的统—Hilbert-Burch矩阵、相关的3×3实矩阵和有理参数表示.对于一般形式曲面,我们亦有类似的结论成立.