计算机图形学十一伯恩斯坦多项式

錶哥

已于 2023-08-08 02:21:18 修改

阅读量304

点赞数 1

分类专栏：计算机图形学文章标签：概率论

于 2022-12-21 11:08:32 首次发布

本文链接：https://blog.csdn.net/weixin_44683202/article/details/128392353

版权

计算机图形学专栏收录该内容

23 篇文章 9 订阅

订阅专栏

伯恩斯坦多项式

开礼包的概率

现在在各种游戏里，都有充斥着开礼包的骗钱方法。那天刀为例，有0.1的概率开到高级物品，假设我们有三个这样的礼包，那么我们可以得到如下一些结果：

首先是拿命玩游戏的人，开到三个高级物品，这个概率很简单，是0.1的三次方。

其次是还能活几天的人，开到2个高级物品，那么这个概率是多少呢？我们知道开到的概率是0.1，那么开不到的概率自然是 1-0.1=0.9。那么开到两个的概率是 (0.1)² ，一个没开到的概率是0.9，总的概率是(0.1)²0.9嘛？错！在这里还存在着排列组合的关系，例如我前两个开到后一个没开到，它的概率是(0.1)²0.9，而我第一个没开到后两个开到，它的概率是(0.1)²0.9，当然还有中间没开到前后开到的情况，这些三种情况都属于开到2个高级物品，因此总概率应该是(0.1)²0.9*3。

然后是欧皇，开到一个高级物品，一样存在排列组合的情况，概率应该为30.1(0.9)²。

最后是正常人，一个都开不到，概率是0.9的三次方。

所有的概率相加的结果是1。

伯恩斯坦多项式

通过上面的例子，我们可以假设有个礼包开到东西的概率为t，我们一共开n次，那么开到i次东西的概率应该怎么计算呢？

i=n时很好计算，t的n次方，同样i=0也很好计算，自然是（1-t）的n次方，相对麻烦的是中间的，因为前面说了，要考虑到排列组合的情况，当然不管怎么排列组合，每个组合的概率应该都是t的i次方(1-t)的(n-i)次方。

接下来的问题就是算开到i次时，存在多少种排列组合的情况了(假设值为k)。根据高中知识，我们可以轻松得到公式为
在这里插入图片描述
（注：规定0！=1，因此即使i=0，除数不会为0）

简单的理解下，我们把问题想成，我们有i个开到的礼包，要插入到这n个位置里，那么当我插入第一个礼包的时候，有n个空位给我放，插入第二个的时候，有n-1个空位，知道插到第i个，此时有n-i+1个空位，那么结果就是
在这里插入图片描述
但是这里面还存在重复的情况，例如我第1个礼包放1号位，第2个礼包放2号位，和第1个礼包放2号位，第2个礼包放1号位是一样的，因此我们要去重。i 个礼包一共会有 i! 次情况，因此最终会得到上面的排列组合结果。

当然了，i=0或者i=n同样适用于这个公式，因此开到 i 次东西的概率应该为：

在这里插入图片描述
我们常用或者来代表伯恩斯坦多项式，其中部分我们常用来代替。
因此伯恩斯坦多项式可以表示为：

其中n为非0正整数，i的取值范围为0-n，t的范围为0-1.

在这里插入图片描述
从曲线可以看出，在t=0.3处会有一个极大值，并向两边衰减，并且所有伯恩斯坦曲线都有了类似的性质。如果我们想要从10次开出3次，那么从统计学来说只有概率等于
0.3(i/n)时，开出3次的概率最高。可以用导数对其进行证明。

贝塞尔曲线，就是使用伯恩斯坦多项式来定义。因此它的一些性质同样会影响到贝塞尔曲线的性质，我们接下来看看伯恩斯坦多项式有哪些性质。

非负性

在这里插入图片描述

归一性

伯恩斯坦多项式的求和等于1，在前面开礼包的例子里，我们已经提到了，所有概率加起来等于1。
我们现在在代数的形式上证明一下：
首先在这里插入图片描述正好是一个二项展开式。

对称性

在这里插入图片描述
可以看出，关于n/2对称。
我们也可以从下图曲线看出其对称性：

对称，也就是说，第i个和倒数第i个是对称的。

证明：
在这里插入图片描述
两者相等。

递归性

递归性是什么意思呢，就是说一个n阶的伯恩斯坦多项式在这里插入图片描述，它可以写成两个n-1阶的伯恩斯坦多项式的组合，即：

很多贝塞尔曲线的定理和证明里都是做这个组合公式的应用。
证明：
首先：
上面式子也是排列组合的基本性质，常用于杨辉三角。

有个很聪明的理解方法：在这里插入图片描述即n里面取i个，那么我们就可以拆分成：
从n-1个里面取i个并且从1个里面取0个，那么对应的就是
从n-1个里面取i-1个并且从1个里面取1个，那么对应的就是
这两个步骤合并在一起就是从n里面取i个。
然后我们可以得到：
在这里插入图片描述
然后就变成了
左项提取1-t，右项提取t，即可得到我们最上面的式子。
这里还有一个问题，即i=0怎么办？那i-1不就等于-1了。

两个特殊情况：
1.当i=0时，i-1不就等于-1了。
这种情况下，式子如下：
在这里插入图片描述
2.当i=n时，会出现的情况，和上面一样理解，可以得到

导数及其最大值

在这里插入图片描述
有个导数，那么求最值就很简单了，我们可以直接令导数的值等于0即可。
我们来解一下：

即：
这个解是唯一的，因此伯恩斯坦多项式会在t=i/n处有唯一的局部最大值。

升阶公式

升阶的意思，即把n阶的伯恩斯坦多项式写成n+1阶。即把它乘以(1-t)或者t，如下：
在这里插入图片描述
两者相加，即可得到：

我们不能把的二次多项式写成和x的三次方相关的三次多项式，除非x=0，但是在伯恩斯坦多项式却可以。

积分

在这里插入图片描述因此伯恩斯坦多项式的总面积并不是1，曲线下的总面积为1/(n+1)，与i无关。

三角域的伯恩斯坦多项式

这个在三角域贝塞尔曲面里非常重要！！！！！
什么叫三角域的伯恩斯坦多项式，我们还是拿开礼包的例子来说明，如下：
假如我有一个礼包，可以开到A,B,C三种物体中的一个，并且不存在一个都开不到的情况。其中开到A的概率为u，开到B的概率为v，开到C的概率为w，并且u+v+w =1。
那么我问如果我靠这个礼包n次，开到i个A，j个B，k个C的概率是多少？因为必出，所以i+j+k=n;

我们可以得到这样一个伯恩斯坦多项式：
在这里插入图片描述
其中u，v，w就是我们三角形重心坐标。我们称之为三角域的伯恩斯坦多项式。

在这里插入图片描述有多少项？我们知道中i的取值范围是0到n，因此一共会有n+1项，那么对于三角域的伯恩斯坦多项式有多少项呢？
当i=0时，j的取值范围可以是0到n，而当我j确定了，k也确定了，例如及那么k必须等于n-1,否则不能保证i+j+k = 1。因此i=0时一共有n+1项，当i=1时有n项。因为j的取值范围变为0到n-1。那么我i可以从0到n，对应的项从n+1慢慢递减，即可得到(n+1)+n+(n-1)+…+1，其中一共n+1项，头尾相加乘以项数除以2，
即在这里插入图片描述一共有项。
其中三角域的伯恩斯坦多项式与普通伯恩斯坦多项式没有太大区别。
都有非负性，归一性等。
其中递归性单独说一下，因为递归性要用其来推到de Casteljau算法。

我们用常理来理解，原本要开n个礼包，其中想要i个A，j个B，k个C。但是现在我们发现，妈的钱不够了，只能买的起n-1个礼包了，那么必然只能有下面三种情况：

情况1：i-1个A，j和k不变，那么概率是在这里插入图片描述
情况2：j-1个B，i和k不变，那么概率是

情况3：k-1个C，i和j不变，那么概率是

这个时候，突然有个大佬说他送你一个礼包，那么如果：
1：在情况1后，你用送的礼包在u的概率下开出了A，那种总概率就是在这里插入图片描述
2：在情况2后，你用送的礼包在v的概率下开出了B，那种总概率就是

3：在情况3后，你用送的礼包在w的概率下开出了C，那种总概率就是
上面的结果变成了，n个礼包开到了i个A，j个B，k个C，因此：

这个公式就是我们的递归公式。