六次PH曲线C^1 Hermite插值.pdf

方林聪等：六次 PH 曲线 C Hermite插值 下曲线参数的符号选择与可行解的性质之间的关系，提出了一种不需要对形状特征进行附加计算的方 法，选择最优的可行解．Yong和 Zheng[11]给出了五次 PH 曲线 Hermite插值的几何方法，并讨论了 最优解的选取标准．Choi和Kwonl【2]比较了符合插值数据条件得到的三次和五次 PH 曲线的相似度 后 进一步提出了用曲线的拓扑结构作为度量标准来选择最优解的方法 [13]．Jiittler[14】讨论了七次PH 曲线的Hermite插值方法，该方法需要对一元四次方程进行求解．Kong等人 [15】采用重新参数化方法 研究了C 连续的PH 曲线的 Hermite插值 问题．Pelosi等人 [16]提出了通过指定控制多边形和节点 序列来构造 C0连续的平面五次PH样条曲线的方法．Wang和 Fang[17]讨论了平面四次PH 曲线的 控制多边形的几何性质，并给出了平面四次PH 曲线的G Hermite插值方法，该方法可以有效地避免 在曲线段上出现奇异点． 本文采用平面 B6zier曲线的复数形式，讨论六次PH 曲线的C 连续的Hermite插值问题，并指 出在构造六次PH曲线的过程中，可以有效地避免奇异点在曲线段上出现．在给定C 连续的Hermite 条件下，可以得到四条五次 PH 曲线 [8_11]，而本文将给出六次 PH 曲线的Hermite插值方法，并指出 作为初始条件之一 可以提供一个 自由参数对最终的曲线进行调整，因此在 曲线构造上将更具灵活性． 2 六次 PH 曲线 一 条平面参数曲线P(t)= (()，(t))的复数形式为P(t)=x(t)+iy(t)．文中以粗斜体字母表示 复变量，而普通斜体字母表示实变量．一条 B~zier曲线 P(t)是PH 曲线当且仅当存在实多项式 n(t) 和复多项式Q(t)满足 (参见文献 [17]) P )=n(t)-Q ()． (2．1) 记 B~zier曲线 P(t)的控制定点为 Pj，J=0… ．，6，APJ=PJ+1一PJ为控制顶点的一阶向前差分， Lj=lIzxP~ll为控制多边形的边长， 为对应控制多边形的边 APj的辐角，B ()为n次Bernstein 多项式， B2(t)=()(1一t)n-Jtj， 对于六次B~zier曲线P()，其一阶导数P )是五次复多项式，且可 以表示为 5 P )=6∑APj (t)． (2．2) j=0 六次PH曲线可以根据 (2．1)中多项式R(t)和Q(t)的次数不同分为两类，第一类曲线中，n(t)和Q(t) 分别为如下的三次和一次多项式： n(t)=aoBa(t)+alB3(t)+a2B3(t)+Ba(t)，Q(t)：zoB~(t)+ZlB (￡)； 而对于第二类 曲线，R(t)和 Q(￡)则分别为如下的一次和二次多项式： R(t)=aoB (t)+B (t)，Q(t)=zoBo2(t)+zlB~(t)+z2B~(t)， 其中aj和 j为待定系数．在这两类六次PH 曲线中，根据代数基本定理可知，实多项式n(t)均至少 存在一个实根，因此，六次PH 曲线上至少存在一个奇异点，为非正则 曲线． 在采用 n次B~zier曲线对给定的Hermite数据进行插值时，由B~zier控制多边形的性质可知， 连续的Hermite条件等价于给定了曲线两端的四个控制顶点，即Po，P1，P 一1和 P ．下文对第一类 和第二类六次PH 曲线的C 连续的Hermite插值分别进行阐述． 800 中国科学：数学 第 44卷 第 7期 3 六次 PH 曲线的 Hermite插值 3．1 第一类曲线 将 (2．1)和 (2．2)展开，利用比较系数法可以得到 6AP0= 0o ， 30AP1=2aoz