B-样条基函数:定义 (B-spline Basis Functions: Definition)

B-样条基函数:定义

B-spline Basis Functions: Definition


 

上一页动机 回目录 下一页重要性质  

 

 

    贝塞尔基函数用作权重。B-样条基函数也一样;但更复杂。但是它有两条贝塞尔基函数所没有的特性,即(1)定义域被节点细分(subdivided); (2) 基函数不是在整个区间非零。实际上,每个B样条基函数在附近一个子区间非零,因此,B-样条基函数相当局部

    设U m + 1个非递减数的集合,u0 <= u2 <= u3 <= … <= umui称为节点knots, 集合U 称为节点向量knot vector, 半开区间[ui, ui+1是第i个节点区间(knot span。注意某些ui可能相等,某些节点区间会不存在。如果一个节点 ui 出现 k (即,ui = ui+1 = … = ui+k-1), 其中 k > 1, ui 是一个重复度(multiplicityk 多重节点,写为 ui(k)。 否则,如果ui只出现一次,它是一个简单节点。如果节点等间距(即, ui+1ui 是一个常数,对 0 <= i <= m – 1),节点向量或节点序列称为均匀的;否则它是非均匀的

    节点可认为是分隔点,将区间[u0, um]细分为节点区间。所有B-样条基函数被假设定义域在[u0, um]上。在本文中,我们经常使用u0 = 0um = 1,所以定义域是闭区间[0,1]

    为了定义B-样条基函数,我们还需要一个参数,基函数的次数(degreep,第ipB-样条基函数,写为Ni,p(u),递归定义如下:

 

 

 

       上述公式通常称为Cox-de Boor递归公式。 这个定义看起来很复杂;但是不难理解。如果次数(degree)为零(即, p = 0),这些基函数都是阶梯函数,这也是第一个表达式所表明的。即,如果u是在第i个节点区间[ui, ui+1)上基函数Ni,0(u)1。 例如,如果我们有四个节点u0 = 0, u1 = 1, u2 = 2u3 = 3, 节点区间 0, 1 2[0,1), [1,2), [2,3)0次基函数是N0,0(u) = 1 [0,1) ,在其它区间是0N1,0(u) = 1  [1,2)上,在其它区间是0N2,0(u) = 1[2,3)上,其它区间是0。如下图所示:

 

 

       为了理解p大于0时计算Ni,p(u)的方法,我们使用三角计算格式。所有节点区间列在左边(第一)列,所有零次基函数在第二列。见下图。

 

        为了计算Ni,1(u),需要Ni,0(u)Ni+1,0(u)。因此,我们可以计算N0,1(u), N1,1(u), N2,1(u), N3,1(u) 等等。所有这些Ni,1(u)写在第三列。一旦所有Ni,1(u)计算完毕,我们可以计算Ni,2(u)并将其放在第四列。继续这个过程直到所有需要的Ni,p(u)的计算完毕。

        上面我们获得了针对节点向量 U = { 0, 1, 2, 3 }N0,0(u), N1,0(u)N2,0(u) 。现在计算N0,1(u)N1,1(u)。要计算N0,1(u),因为i = 0p = 1,从定义出发有

 

因为u0 = 0, u1 = 1u2 = 2,上式变为

 

因为N0,0(u)[0,1)上非零且N1,0(u)[1,2)上非零,如果u[0,1) (resp., [1,2) ), 只有N0,0(u) (resp., N1,0(u) )N0,1(u)有贡献。因此,如果u [0,1), N0,1(u) uN0,0(u) = u而如果u [1,2), N0,1(u) (2 – u)N1,0(u) = (2 – u)。相似的计算得到N1,1(u) = u – 1如果u [1,2), N1,1(u) = 3 – u 如果 u [2,3)上。下图中,黑色和红色线分别是N0,1(u)  N1,1(u)。注意N0,1(u) (resp., N1,1(u))[0,1) [1,2) (resp., [1,2)  [2,3))是非零的。

 

一旦获得N0,1(u)N1,1(u),可以计算N0,2(u)。由定义得到下式:

 

代入节点值得到

 

注意 N0,1(u) [0,1)[1,2)上非零而N1,1(u) [1,2) [2,3)上非零。因此,我们有三种情况要考虑:

 

1. u [0,1):
    这种情况,只有N0,1(u)N0,2(u)的值有贡献。因此,N0,1(u)u, 得到

 

2. u [1,2):
    这种情况, N0,1(u)N1,1(u)都对 N0,2(u)有贡献。因此N0,1(u) = 2 – u N1,1(u) = u – 1 [1,2),得到 

 

   【译注:上式中间的式子的第二项应为:0.5(3-u)(u-1)

3. u [2,3):
    这种情况,只有N1,1(u) N0,2(u)有贡献。因此N1,1(u) = 3 – u [2,3)上,得到,

 

 

    如果我们画出上述三种情况的曲线段,我们会看到两个相邻曲线段连接起来形成了在节点上的曲线。更确切地,第一种和第二种情况的曲线段在u = 1处连接起来,而第二种和第三种情况的曲线段在u = 2处连接起来。注意合成曲线是光滑的,但是如果节点向量包含多重节点通常就不是这样的。

 

 

 

两个重要的观察

       因为 Ni,1(u) 是从 Ni,0(u) Ni+1,0(u)计算的而 因为Ni,0(u)Ni+1,0(u) 在区间[ui, ui+1)[ui+1, ui+2)分别是非零的,Ni,1(u) 在这两个区间都是非零的。换句话说,Ni,1(u)[ui, ui+2)上是非零的。相似地,因为 Ni,2(u) 依赖于Ni,1(u) Ni+1,1(u)且因为这两个基函数在[ui, ui+2)[ui+1, ui+3)分别是非零的,Ni,2(u)[ui, ui+3)上非零。总之,为确定基函数Ni,p(u), 的非零定义域,可以追溯到三角计算格式直到回到第一列。例如,假设我们想找到 N1,3(u)的非零定义域。基于上述讨论,我们可从西北和西南方向追溯直到第一列为止,如下图中蓝色虚线所示。因此 N1,3(u) [u1, u2), [u2, u3), [u3, u4) [u4, u5)上是非零的。或,相等地,它在[u1, u5)上非零。.

 

总之,我们有下列观察:

    基函数 Ni,p(u) [ui, ui+p+1)上非零。或,相等地,Ni,p(u) p+1个节点区间[ui, ui+1), [ui+1, ui+2), …, [ui+p, ui+p+1)上非零。

      接着,我们看相反的方向。给定一个节点区间[ui, ui+1),我们想知道哪个基函数会在计算中使用这个区间。我们可以以这个节点区间开始并画一个西北界限箭头和一个西南界限的箭头。所有封闭在楔形里的基函数使用 Ni,0(u)(为什么?)因此在该区间是非零的。因此,所有在[ui, ui+1)上非零的p 次基函数是这个楔形和包含所有Ni,p(u) 的列的交集。实际上,这一列和两个箭头形成一个等边三角形,而这一列是垂直边。 从 Ni,0(u) 数到 Ni,p(u) p+1列。因此,等边三角形的垂直边至多有p+1 项,即 Ni,p(u), Ni-1,p(u), Ni-2,p(u), …, Ni-p+2,p(u), Ni-p+1,p(u) Ni-p,p(u)

 

      让我们看上图。为了找到所有3次在 [u4, u5) 上非零的基函数,画出两个箭头和所有在垂直边的函数是我们想要的。这个例子,是N1,3(u), N2,3(u), N3,3(u), N4,3(u).用黄色三角表示。蓝色 (resp., 红色) 三角显示的是在[u3, u4) (resp., [u2, u3) )上非零的3次基函数。注意在[u2, u3)上只有33次基多项式。.

      总之,我们观察到下列特性:

    在任何一个节点区间 [ui, ui+1), 最多有 p+1p 次基函数非零,即:Ni-p,p(u), Ni-p+1,p(u), Ni-p+2,p(u), …, Ni-1,p(u) Ni,p(u),

 

 

系数的意义是什么?

    最后,让我们研究下Ni,p(u)定义中系数的意义。当计算 Ni,p(u) 时,它使用Ni,p-1(u)Ni+1,p-1(u)。前者在 [ui, ui+p)上非零。如果u 是在这个半开区间,那么uui u 和这个区间左端之间的距离,区间长度是ui+pui, ,而(uui) / (ui+pui) 是上述距离的比且在01之间。见下图。第二项,Ni,p-1(u),在[ui+1, ui+p+1)上非零。如果u 在该区间,那么ui+p+1u u 到该区间右端的距离,ui+p+1ui+1 是区间长度,而(ui+p+1u) / (ui+p+1ui+1) 是这两个距离的比且值在01之间。因此, Ni,p(u) Ni,p-1(u) Ni+1,p-1(u)的线性组合,有两个系数,都在 u上是线性的,在01之间。

 

 

 

 

 


上一页动机 回目录 下一页重要性质


 

译注:

  1.  本文翻译是“B-样条曲线(B-spline Curves)教程”中的一部分,其余翻译部分见“B-样条曲线(B-spline Curves)教程目录”。
  2.  “B-样条曲线(B-spline Curves)教程”是翻译自C.-K. Shene博士的CS3621 Introduction to Computing with Geometry Notes的第6部分“B-spline Curves”。
  3. 本文原文地址:B-spline Basis Functions: Definition   。
  4. 本文首发“博士数学家园

 

© 版权声明
THE END
喜欢就支持一下吧
点赞0 分享
评论 抢沙发
头像
欢迎您留下宝贵的见解!
提交
头像

昵称

取消
昵称表情代码图片