带形状参数的五次Bezier曲线的光滑拼接.pdf

第39卷第4期2023年8月山西大同大学学报（自然科学版）Journal of Shanxi Datong University(Natural Science Edition)Vol.39 No.4Aug.2023带形状参数的五次Bezier曲线的光滑拼接孙明灿，师晶（闽南理工学院 信息管理学院，福建 石狮 362700）摘要：针对复杂曲线的造型研究了带形状参数的五次Bezier曲线的光滑拼接问题。利用基函数性质和曲线性质，得到了曲线间的G1、G2、C1、C2光滑拼接定理。通过计算实例验证了曲线间拼接条件的有效性。关键词：五次Bezier曲线；光滑拼接；曲线造型中图分类号：TP391文献标识码：Adoi：10.3969/j.issn.1674-0874.2023.04.008曲线光滑拼接是计算机辅助几何设计中的重要研究内容1，它不仅能满足工业曲线复杂化、多样化的需求，而且广泛应用于电子信息、工业制造及流体力学等领域。Bezier曲线是曲线造型的常用方法，它具有图形直观、方便调控等优点2-3。由于改变Bezier曲线的控制顶点，会影响曲线的形状，所以其不能构造较复杂的工程曲线和曲面4-5。另外，随着曲线求导次数的增加，计算量逐步增大，并且容易出现一些难以察觉的不可导的点，这些缺点限制了Bezier曲线在曲线造型中的应用6-8。研究了带形状参数的五次Bezier曲线的光滑拼接问题，并给出曲线基函数和曲线的性质，得到了曲线间不同条件下的光滑拼接定理。通过计算实例说明曲线的光滑拼接定理不但能构造较复杂的工程曲线，而且还可为高次曲线的拼接提供参考。1 曲线定义及性质定义1 设 -2,1,t 0,1，称关于t的多项式 b0,5(t)=(1-t)(1-t)4b1,5(t)=(4+-2t)(1-t)3tb2,5(t)=(6+)(1-t)2t2b3,5(t)=(4-+2t)(1-t)t3b4,5(t)=(1-+t)t4（1）为带形状参数的五次Bezier曲线基函数。式（1）中的基函数满足以下性质：（1）非负性：bi,5(t)0,t 0,1,(i=0,1,2,3,4)。（2）规范性：i=04bi,5(t)=1,t 0,1。（3）对称性：bi,5(t;)=b4-i,5(1-t;)(i=0,1,2,3,4)。（4）端点性质：b0,5(0)=b4,5(1)=1,bi,5(0)=0,(i=1,2,3,4),bi,5(0)=0,(i=1,2,3,4),bi,5(0)=0,(i=2,3,4),bi,5(1)=0,(i=0,1,2,3),bi,5(1)=0,(i=0,1,2)。（5）退化性：当形状参数=0时，五次Bezier曲线基函数退化为传统的五次Bernstein基函数。（6）线性无关性：基函数 bi,5(t)(i=0,1,2,3,4)是线性无关的。（7）最大值：基函数bi,5(t)(i=0,1,2,3,4)在区间0,1上有一个最大值。当=1时，式（1）中的基函数的图形如图1。图1 基函数图形收稿日期：2022-12-19基金项目：国家科技重大专项资助项目2016ZX05045-004-05；福建省中青年教师教育科研项目JAT200761作者简介：孙明灿（1983-），河南平顶山人，讲师，硕士，研究方向：计算机辅助几何设计。E-mail:文章编号：1674-0874（2023）04-0036-042023年定 义 2 给 定R3空 间 中 5 个 控 制 顶 点Pi(i=0,1,2,3,4)，定义带形状参数的五次 Bezier 曲线为p(t)=i=04bi,5(t)Pi，其中，t 0,1。带形状参数的五次 Bezier 曲线p(t)具有如下性质：（1）端点性质：曲线p(t)自首端点P0开始，至末端点P4结束，并且在首末端点的切矢模长分别等于控制多边形首末边边长的（+4）倍，即p(0)=P0,p(1)=P4,p(0)=(+4)(P1-P0)，p(1)=(+4)(P4-P3)。证明 由式(1)及定义2可得：p(t)=i=04bi,5(t)Pi=P0(1-t)(1-t)4+P1(4+-2t)(1-t)3t+P2(6+)(1-t)2t2+P3(4-+2t)(1-t)t3+P4(1-+t)t4。将t=0,t=1分别代入上式可得p(0)=P0,p(1)=P4，故曲线p(t)分别插值于首末两端点P0,P4。又因为p(t)=P0(5t-4)(1-t)3+P1-2t(8+4-5t)+4(1-t)2+P2 2t(6+)(1-t)(1-2t)+P3-2t(8-6+5t)-3+12 t2+P4(5t-4+4)t3。（2）将t=0,t=1分 别 代 入 上 式 后 整 理 可 得p(0)=(+4)(P1-P0)p(1)=(+4)(P4-P3)，故曲线p(t)分别与直线P0P1,P3P4相切，且在首末端点的切矢模长分别等于控制多边形首末边边长的（+4）倍。（2）几何不变性：曲线p(t)的形状仅与控制顶点Pi(i=0,1,2,3,4)有关，而与坐标系的方向和位置无关。证 明 对 控 制 顶 点 为Pi(i=0,1,2,3,4)的 曲 线p1(t)进行线性变换M和平移变换N后，得到控制顶点为Qi(i=0,1,2,3,4)的曲线p2(t)，即：p2(t)=Mp1(t)+N=Mi=04bi,5(t)Pi+Ni=04bi,5(t)Pi=i=04Mbi,5(t)Pi+i=04Nbi,5(t)Pi=i=04(M+N)bi,5(t)Pi=i=04bi,5(t)Qi,故曲线具有几何不变性。（3）凸包性：曲线p(t)在其控制顶点Pi(i=0,1,2,3,4)构成的凸包内。证明 因为基函数bi,5(t)0，且i=04bi,5(t)=1，(0 t 1;i=0,1,2,3,4)，所以当t 0,1时，p(t)是特征多边形各顶点的加权平均，且权因子为基函数，即曲线p(t)是控制顶点Pi(i=0,1,2,3,4)的凸线性组合，故曲线p(t)落在控制顶点Pi(i=0,1,2,3,4)构成的凸包内。（4）对称性：将曲线p(t)的控制多边形的次序取反后，控制顶点Pi(i=0,1,2,3,4)定义的曲线与控制顶点P4-j(j=0,1,2,3,4)定义的曲线相同，仅方向相反，即p(t;P0,P1,P2,P3,P4)=p(1-t;P4,P3,P2,P1,P0)-2,1,t 0,1。证明 由定义2得p(t;P0,P1,P2,P3,P4)=i=04bi,5(t)Pi=j=04b4-j,5(1-t)P4-j=p(1-t;P4,P3,P2,P1,P0)。故曲线满足对称性。（5）变差缩减性：曲线p(t)比其控制多边形的波动小，即曲线p(t)的光滑性不低于其控制多边形的光滑性。（6）保凸性：当曲线p(t)的控制多边形是凸多边形时，曲线p(t)是凸曲线。2 曲线拼接设：p1(t;1)=i=04bi,5(t)Pi,0 t 1,-2 1 1p2(t;2)=i=04bi,5(t)Qi,0 t 1,-2 2 1分别为控制顶点Pi,Qi(i=0,1,2,3,4)定义的带形状参数的五次Bezier曲线。定理 两条带形状参数的五次 Bezier 曲线p1(t;1)和p2(t;2)在连接点P4=Q0处G1光滑拼接的充要条件是：1+42+4=Q1-Q0Q0-P3，（3）其中：0。特别地，当=1时，两曲线在连接点P4=Q0处达到C1光滑拼接。证明 若两曲线p1(t;1)和p2(t;2)在连接点P4=Q0处G1光滑拼接，则它们在连接点处具有相同的一阶切矢方向，即孙明灿等：带形状参数的五次Bezier曲线的光滑拼接37山西大同大学学报(自然科学版)2023年p1(1;1)=p2(0;2),0。（4）由曲线的端点性质得p1(1;1)=(1+4)(P4-P3)p2(0;2)=(2+4)(Q1-Q0)将上式代入式（4）化简得(1+4)(P4-P3)=(2+4)(Q1-Q0)。又因为P4=Q0，整理后可得式(3)是两曲线p1(t;1)和p2(t;2)在连接点P4=Q0处G1光滑拼接的充要条件。当=1时，即1+42+4=Q1-Q0Q0-P3成立时，两曲线在连接点P4=Q0处具有相同的一阶切矢，此时这两条曲线达到C1光滑拼接。为方便讨论两条曲线在连接点P4=Q0处G2和C2光滑拼接的情况，设P3P4与-2(12+51)P3+4(21+3)P4组成的平行四边形的面积为S1，且P3P4与4(22+3)Q0-2(12+52)Q1+2(6+2)Q2组成的平行四边形的面积为S2。定理 2 两条带形状参数的五次 Bezier 曲线p1(t;1)和p2(t;2)在连接点P4=Q0处G2光滑拼接的充要条件是：两曲线满足G1光滑拼接，并且有S2=S1（5）其中：=(1+4)2(2+4)2。证明 两曲线p1(t;1)和p2(t;2)在连接点P4=Q0处G2光滑拼接的充要条件是：两曲线除满足G1光滑拼接外，还需在连接点处具有相同的曲率矢:1 p1(1;1)=2 p2(0;2)，即|p1(1;1)p1(1;1)|p1(1;1)3=|p2(0;2)p2(0;2)|p2(0;2)3（6）由曲线的端点性质及式（2）得p1(1;1)=(1+4)(P4-P3)p2(0;2)=(2+4)(Q1-Q0)（7）p1(1;1)=-2(12+51)P3+4(21+3)P4p2(0;2)=4(22+3)Q0-2(12+52)Q1+2(6+2)Q2（8）将式（7）和式（8）代入式（6）整理后即可得出式（5），得证。定理3 当形状参数1=2=1时，两条带形状参数的五次 Bezier曲线p1(t;1)和p2(t;2)在连接点P4=Q0处C2光滑拼接的充要条件是：两曲线满足C1光滑拼接及-3P3+2P4=2Q0-3Q1+Q2。证明 两曲线p1(t;1)和p2(t;2)在连接点P4=Q0处C2光滑拼接的充要条件是：两曲线除满足C1光滑拼接外，还需在连接点处具有相同的二阶导数，即p1(1;1)=p2(0;2)，（9）将C1光滑拼接及式（8）代入式（9）整理后即可得出结论，得证。3 计算实例在曲线造型中，适当选取控制顶点，构造带形状参数的五次Bezier曲线。通过调整参数取值可以灵活改变Bezier曲线的形状。根据曲线间的光滑拼接定理，可设计较复杂的工程曲线和曲面。3.1 实例一在构造高脚酒杯曲面时，母线需要两条五次Bezier曲线拼接而成。首先，根据控制顶Pi,Qi(i=0,1,2,3,4)分别计算出两条五次Bezier曲线p1(t;1)和p2(t;2)。然后，在控制多边形P0P1P2P3P4和Q0Q1Q2Q3Q4中，分别调整两个形状参数1,2以获得较满意的曲线。当1=0.6，2=0.2时，根据定理 2将两条曲线在连接点P4(Q0)处进行光滑拼接，拼接后生成的高脚酒杯曲面的母线，如图2。最后，将母线旋转即可得到高脚酒杯曲面，如图3。图2 高脚酒杯曲面的母线图3 高脚酒杯曲面382023年3.2 实例二在构造烛台曲面时，母线需要4条五次Bezier曲线拼接而成。首先，根据控制顶点分别计算出4条五次Bezier曲线p1(t;1)，p2(t;2)，p3(t;3)和p4(t;4)。然后，在 4 个 控 制 多 边 形 中，分 别 调 整 4 个 形 状 参 数1,2,3,4以获得较满意的曲线。当1=0.89,2=0.05,3=0.13,4=0.51时，根据定理1和定理3将四条曲线分别在连接点处进行光滑拼接，拼接后生成的烛台曲面的母线如图4。最后，将母线旋转可得到烛台曲面，如图5。图4 烛台曲面的母线图5 烛台曲面4 结语研究带形状参数的五次 Bezier 曲线的光滑拼接问题，得到了曲线基函数性质、曲线性质及曲线间的 G1、G2、C1、C2连续性条件。通过计算实例验证了曲线间光滑拼接的有效性。在今后的工作中将对高次 Bezier曲面的光滑拼接进行探索。参考文献1 施法中.计算机辅助几何设计与非均匀有理 B样条 M.北京:高等教育出版社，2001:306-454.2 查东东，刘华勇，王曾珍.带形状参数的三次三角域 Bzier曲面 J.计算机工程与科学，2021，43（11）:1994-2002.3 王海波，杨当福，佘卫勤，等.带2个形状参数的多项式可展曲面造型 J.浙江大学学报（理学版），2021，48（2）:131-142.4 吴丽娟，李博，ABEYSINGHE A S S，等.B样条曲面拼接算法的设计与实现 J.沈阳师范大学学报(自然科学版)，2019，37（6）:549-553.5 王崇.基于近似光滑的样条曲面拼接方法研究 D.长春：吉林大学，2022.6 师晶.一种基于几何约束的插值曲线的参数连续性 J.沈阳大学学报(自然科学版)，2019，31（1）:78-83.7 孙明灿，师晶.一种代数曲线的 C3连续性条件 J.新乡学院学报，2018，35（12）:9-12.8 江卯，喻德生.三次 T-Bzier曲线间的混合延拓 J.浙江大学学报(理学版)，2015，42（6）:696-703.Smooth Stitching of Quintic Bezier Curve with a Shape ParameterSUN Ming-can,SHI Jing(School of Information Management,Minnan University of Science and Technology,Shishi Fujian,362700)Abstract:The problem of smooth concatenation of quintic Bezier curves with shape parameters was studied for the modelingof complex curves.By utilizing the properties of basis functions and curves,the G1,G2,C1,C2smooth concatenation theorems be-tween curves were obtained.The effectiveness of the splicing conditions between curves was verified through calculation exam-ples.Key words:quintic Bezier curve;smooth stitching;curve modeling责任编辑 高彩云孙明灿等：带形状参数的五次Bezier曲线的光滑拼接39