如何计算具有给定隐函数的贝塞尔曲线的控制点

Question

我正试图编写一个程序来绘制bezier曲线段，以实现我课本上的知识。我想尝试制作一个绘图程序，用隐式函数绘制一些东西。我希望我的程序按以下方式运行。

• 用户输入一个隐含函数，比如y=2x^3
• 计算控制点（不知道怎么做）
• 绘制贝塞尔曲线（完成）
• 下一个函数

理论上，我可以通过直接替换来绘制曲线，但我可能会对曲线进行一些修改，我想实现我所学到的，所以，我想明智地做这件事（直接替换似乎很愚蠢）。Thz^^

编辑1：假设边界是由用户提供的

Answer 1

第一步是为曲线生成参数化表达式。给定的例子可以很容易地转换：

c(t) = (t, 2 * t^3)^T

现在用单调基础表示这条曲线：

c(t) = / 0  1  0  0  * (1, t, t^2, t^3)^T
        0  0  0  2 /
     = C * M(t)

在该表达式中，第一矩阵C是系数矩阵。我们所需要做的就是将这个矩阵转换为Bernstein基。将单调基转换为Bernstein基的矩阵是：

       / 1 - 3t + 3t^2 -  t^3    / 1 -3  3 -1     /  1  
B(t) = |     3t - 6t^2 + 3t^3 | = | 0  3 -6  3 | *  |  t  |
       |          3t^2 - 3t^3 |   | 0  0  3 -3 |    | t^2 |
                         t^3 /    0  0  0  1 /     t^3 /

这个方程可以反过来得到：

       / 1   1   1   1 
M(t) = | 0  1/3 2/3  1 | * B(t)
       | 0   0  1/3  1 |
        0   0   0   1 /

将其代入曲线方程，我们得到：

c(t) = C * M(t)
           / 1   1   1   1 
     = C * | 0  1/3 2/3  1 | * B(t)
           | 0   0  1/3  1 |
            0   0   0   1 /

可以计算第一个矩阵乘积：

c(t) = / 0  1/3  2/3  1  * B(t)
        0   0    0   2 / 

这给了你贝塞尔曲线的控制点：

p0 = (0, 0)^T
p1 = (1/3, 0)^T
p2 = (2/3, 0)^T
p3 = (1, 2)^T

这个过程可以应用于任何多项式曲线。

形式方程的通解

y = a + b * x + c * x^2 + d * x^3

是：

p0 = (0, a)^T
p1 = (1/3, a + b/3)^T
p2 = (2/3, a + 2b/3 + c/3)^T
p3 = (1, a + b + c + d)^T