您当前的位置: 首页 >  VT LI

三维曲线

VT LI 发布时间:2021-01-28 00:10:03 ,浏览量:4

 三维曲线的表示方式有几种方式:包括参数多项式曲线,艾特肯多项式算法,拉格朗日基多项式,埃尔米特形式曲线,贝塞尔曲线等。最后我们会讨论样条曲线的方式。

 

一:参数多项式曲线:

说道曲线,我们很自然的想到一种应用最广的曲线函数,就是圆的参数描述方式:

而多项式参数的表示形式是类似如下形式:

其中n被称为多项式的次数。这也是一个单项式形式或叫幂形式。

 

二:艾特肯算法:

他是利用一种分而治之的方式实现的额曲线,也就是一条仅插入前n-1个点的曲线,忽略最后的点;另一条插入最后n-1个点的曲线,二不用担心第一个点。然后,将这两条曲线混合在一起即可。

艾特肯算法可以将二次曲线创建为两条线段的混合。


相应的两个控制点之间的额线性差值公式是:

,

将连个线段再做一次差值运算得到二次曲线为:

也就是根据你分的段先计算每段的结果,再结合每两端的结果,不断下去。相当于归并的思想。

艾特肯的算法是有限的,如下图左侧展示的是三条线段混合形成两个二次段。右图展示了两条二次曲线的混合产生了我们一直在寻找的最终结果:一个插值所有四个控制点的三次样条。

但是这个多项式插值病没有真正为我们提供集合设置中曲线设计所需的控制类型。

三:拉格朗日基多项式:

关于拉格朗日多项式的公式为:

用于节点矢量t1,L,tn的拉格朗日基多项式看起来像:

如果我们带入均匀矢量(0,1/3,2/3,1)的话。相应对的值为:

具体计算为:

通过上式可以看出,拉格朗日基多项式是根据当前计算的i来确定其他项对自己的差值。比如就是把第一项0作为基础,跟其他项(除了自己)相减。分号上方是t跟除了0之外的相减,下方是0减去除了自己外的值。最终得到下式。

所以p(t)为:

下面展示了多项式分解的曲线:

总和得到一条我们最终的曲线

 

四:多项式插值的汇总

艾特肯的算法是一种基于重复线性插值的几何方法,通过它可以在不知道曲线的多项式的情况下计算给定t的曲线上的点。拉格朗日插值通过创建仅依赖于节点矢量的基函数来工作。但是多项式插值不理想在于不能控制越线(指超过了两个中间控制点的y值)的情况。

 

五:埃尔米特曲线:

它是通过端点处的切线控制曲线的形状,它通过列出起始位置和结束位置以及导数来制定曲线。三次曲线只有四个系数,这允许仅指定一阶导数,即端点处的速度。因此,描述埃尔米特形式的三次曲线可归结为以下4条信息:

t=0时的起始位置。

t=0时的一阶导数(初始速度)

t=1时的结束位置。

t=1时的一阶导数(最终速度)

 

y以下是一些三次埃尔米特曲线的例子:

假设三次多项式函数为P(t) = c0 + c1*x + c2*x² + c3*x³。

然后我们知道p(t)的导数是:。而导数就是他的速度

可得埃尔米特条件的方程组是:

p(0)=p0=>c0=p0                       (这里是相当于把0带入P(t))

v(0)=v0=>c1=v0                       (这里是相当于把0带入)

v(1)=v1=>c1+2c2+3c3=v1       (这里是相当于把1带入)

p(1)=p1=>c0+c1+c2+c3=p1    (这里是相当于把1带入P(t))

前后两个公式指定端点,中间两个指定速度。

然后求解上面的方程组,产生一种通过埃尔米特位置和导数来计算单项式系数的方法:

上面解出了四个未知项。

 

使用矩阵的方式表示为:

然后我们合并后两个矩阵得:

也就得到三次埃尔米特基函数:

这个是三次埃尔米特曲线的设计,当然也可以得到更高次数的埃尔米曲线设计。在使用三次样条的情况下,可以指定终点的位置(“0阶”导数)和速度(一阶导数)。当同时还指定了加速度(二阶导数)时,会出现5次埃尔米特曲线。

 

六:贝塞尔曲线:

贝塞尔曲线采用的是近似方法而不是插值:虽然他们确实通过了第一个和最后一个控制点,但他们只能通过内部点附近。所以贝塞尔控制点被称为“控制点”而不是“节点”。

艾特肯的算法在贝塞尔曲线算法中的对应物是de Casteljau算法,他是一种通过重复线性插值构造贝塞尔曲线的递归几何级数。

6.1关于de Casteljau算法:

 他可以通过原来的4个控制点开始,可以推导出3个新点,和,然后这三个点的每一对之间的另一轮插值将给出两个点和,最终的插值将产生点,这个点就是我们想要的曲线的点。

得到de Casteljau的递归关系:

下面这张图展示了t在0.25,0.5,0.75时的点。

 

下面是de Casteljau算法中的三次曲线的层级关系

Vector3 deCasteljau(
    int n,                //曲线的阶,点的数量
    Vector3 points[],     //点的数组,将被覆盖
                          //因为该算法是在原地执行的
    float t               //要评估的参数值
){
    //原地执行转换
    while (n > 1)
    {
        --n;
        //执行下一轮的插值
        //曲线的次数减1
        for (int i = 0; i < n; ++i)
        {
            points[i]=points[i] * (1.0f - t) + points[i+1] * t;
        }
    }
    //最后剩下的一个就是结果
    return points[0];
    
}

线性情况直接来自递归关系,无需任何实际操作,具体如下:

二次多项式为:

使用矩阵表示法的二次贝塞尔曲线

单项式形式的三次贝塞尔曲线:

通过单项式的三次贝塞尔曲线可以得到:

单项式形式的三次贝塞尔曲线矩阵:

6.2:伯恩斯坦基多项式

de Casteljau是基于系数t的幂实现的多项式,我们还可以通过收集控制点上的项而不是t的幂来写出贝塞尔形式的多项式。

一次多项式:

二次多项式:

三次多项式:

这个模式会更清晰,每个项都是有一个常量系数,一个(1-t)的幂和一个t的幂。t的幂以递增的顺序编号。常量系数如下图帕斯卡三角形:

也就是任意次数的贝塞尔曲线可以定义为:

伯恩斯多项式有以下属性:

1.总和为1:

也就是对弈t的所有值总和为1,不然无法定义中心坐标

2.凸包属性:

曲线有界限以保持在控制点的凸包内。可以将其与拉格朗日基多项式进行比较,后者不会保持在区间内,也导致了拉格朗日基多项式插值不遵守凸包属性。

3.端点插值:

根据需要,可以让第一个和最后一个多项式达到1,他们两个端点就会接触。因为且。

4.全局支持:

所有多项式的开区间(0,1)上都是非零的,也就是说,除了端点之外的整个曲线都是非零的。控制点的混合权重非零的区域称为控制点的支持。只要控制点有支持,它就会对曲线产生一些影响。

当移动一个特定的控制点并且只有该控制点附近的曲线部分受到影响时,称为局部支持。

5.局部最大值:

虽然每个控制点会对整个曲线产生影响,但每个控制点在曲线的一个特定点上施加的影响最大。

 

6.3 贝塞尔导数及其余埃尔米特形式的关系:

我们先来看看贝塞尔曲线的导数。一阶导是速度,二阶导加速度。速度跟曲线的切线方向有关,加速度与曲率有关。

然后我们已知:

将他带入上式得:

然后我们考虑端点t=0和t=1的速度:

可以看到,b1-b0给出的是从第一个控制点到第二个控制点的矢量;而b3-b2是第三个控制点到最后一个控制点的矢量。t=0时曲线起点处的切线瞄准第一个控制点,t=1时曲线末端的切线瞄准第三个控制点。

现在我们知道贝塞尔控制点和曲线速度之间的关系了,很容易从贝塞尔转换为埃尔米特形式。

我们可以看到一阶导跟两个端点相关,而二阶导则是与三个端点相关。我们可以推导下:

三次贝塞尔曲线的加速度:

带入:

得:

在端点处,加速度是:

可以看到就是由三个控制点决定的。

一般来说,如果移动控制点bi,则会影响i阶导数和在曲线的开始出的更高阶的导数,但是不会影响编号较低的导数。

 

 

 

 

关注
打赏
查看更多评论

VT LI

暂无认证

  • 4浏览

    0关注

    46博文

    0收益

  • 0浏览

    0点赞

    0打赏

    0留言

私信
关注
热门博文
立即登录/注册

微信扫码登录