张火火和西西弗一起推石头

以执念对抗不完美

贝塞尔曲线原理、推导及Matlab实现

一、简介

贝塞尔曲线提出

在数学的数值分析领域中,贝塞尔曲线(English:Bézier curve)是计算机图形学中相当重要的参数曲线。更高维度的广泛化贝塞尔曲线就称作贝兹曲面,其中贝兹三角是一种特殊的实例。

贝塞尔曲线于1962年,由法国工程师皮埃尔·贝兹Pierre Bézier)所广泛发表,他运用贝塞尔曲线来为汽车的主体进行设计。贝塞尔曲线最初由保尔·德·卡斯特里奥于1959年运用德卡斯特里奥算法开发,以稳定数值的方法求出贝塞尔曲线。

贝塞尔曲线与F1的渊源

说到皮埃尔·贝兹,就不得不提他作为雷诺公司的一名工程师,曾引领了设计和制造业的转型,将它们从单纯使用数学和计算工具带向了计算机辅助设计和三维模型,同时是实体造型、几何模型和物理模型领域的奠定者之一。

而雷诺公司早期于1977年开始,作为车队和发动机供应商在不同时期里参加了F1(一级方程式赛事)。在1977年的最后五站比赛中,雷诺以厂队的身份踏进了F1比赛,并且将涡轮增压器引入其第一款赛车雷诺 RS01中却遭众人嘲笑,然而在法国的分站赢得胜利之后,不到两年时间,所有一流车队都改用涡轮发动机。在1992年雷诺凭英国车手文素和威廉士车队拿下史上第一个年度车队冠军,而接下来五年的年度车队冠军也都是使用雷诺发动机(除1995年是班尼顿,其余皆为威廉士),不过随着其在98年退出F1,雷诺在F1的霸业也宣告结束。雷诺在2000年接手贝纳通车队重返一级方程式赛车。在2002年,车队改名为雷诺一级方程式车队,更在2005年和2006年赛季中同时夺得车队和车手世界冠军。

现今雷诺公司旗下的F1车队为AlpineF1车队(法语:Alpine F1 Team),全名倍世AlpineF1车队(法语:BWT Alpine F1 Team)。本赛季Alpine车队的新赛车A524在赛季伊始便严重缺乏竞争力;赛季揭幕战巴林大奖赛上,排位赛中双车排名垫底在第一节即被淘汰,正赛17、18名带回。在沙特大奖赛,双车17、18位发车,正赛加斯利因变速箱问题退赛,奥康带回13名。位列厂队垫底,心疼他们的两名车手。

跑题了,赶紧回归正题。

二、原理

特性

贝塞尔曲线有如下特性:

  • 通过n个控制点P1,P2,P3,...,Pn控制曲线形状。
  • 曲线只经过起点P1和终点Pn.

构建贝塞尔曲线

  1. 一次曲线(线性):
    一次贝塞尔曲线由两个端点P0P1组成,曲线由线段直接连接这两个点,B(t)描述一条由P0P1的直线,可以视为一个P0P1的连续中间点Q0


一次贝塞尔曲线演示.gif

线性贝塞尔曲线演示动画,t在[0,1]区间
  1. 二次曲线
    二次贝塞尔曲线由定点P0,P1,P2的函数B(t)描述。为构建贝塞尔曲线,引入中间点Q0,Q1
    • P0,P1之间的连续点Q0,描述一条线性贝塞尔曲线;
    • P1,P2之间的连续点Q1,描述一条线性贝塞尔曲线;
    • Q0,Q1之间的连续点B(t),描述一条二次贝塞尔曲线。


二次贝塞尔曲线的结构.png

二次贝塞尔曲线的结构


image

二次贝塞尔曲线演示动画,t在[0,1]区间
  1. 三次曲线
    对于三次曲线,由定点P0,P1,P2,P3的函数B(t)描述。可以加入三个可由线性贝塞尔曲线描述的中间点Q0,Q1,Q2,和两个由二次贝塞尔曲线描述的点R0,R1来描述:
    • P0,P1之间的连续点Q0,描述一条线性贝塞尔曲线;
    • P1,P2之间的连续点Q1,描述一条线性贝塞尔曲线;
    • P2,P3之间的连续点Q2,描述一条线性贝塞尔曲线;
    • Q0,Q1之间的连续点R0,描述一条二次贝塞尔曲线;
    • Q1,Q2之间的连续点R1,描述一条二次贝塞尔曲线;
    • R0,R1之间的连续点B(t),描述一条三次贝塞尔曲线。


三次贝塞尔曲线的结构.png

三次贝塞尔曲线的结构


image

三次贝塞尔曲线演示动画,t在[0,1]区间
  1. 四次曲线
    对于四次曲线,由定点P0,P1,P2,P3,P4的函数B(t)描述。可以加入四个可由线性贝塞尔曲线描述的中间点Q0,Q1,Q2,Q3,三个由二次贝塞尔曲线描述的点R0,R1,R2,和两个由三次贝塞尔曲线描述的点S0,S1来描述:
    • P0,P1之间的连续点Q0,描述一条线性贝塞尔曲线;
    • P1,P2之间的连续点Q1,描述一条线性贝塞尔曲线;
    • P2,P3之间的连续点Q2,描述一条线性贝塞尔曲线;
    • P3,P4之间的连续点Q3,描述一条线性贝塞尔曲线;
    • Q0,Q1之间的连续点R0,描述一条二次贝塞尔曲线;
    • Q1,Q2之间的连续点R1,描述一条二次贝塞尔曲线;
    • Q2,Q3之间的连续点R2,描述一条二次贝塞尔曲线;
    • R0,R1之间的连续点S0,描述一条三次贝塞尔曲线;
    • R1,R2之间的连续点S1,描述一条三次贝塞尔曲线;
    • S0,S1之间的连续点B(t),描述一条四次贝塞尔曲线。


四次贝塞尔曲线的结构.png

四次贝塞尔曲线的结构


image

四次贝塞尔曲线演示动画,t在[0,1]区间
  1. 高阶曲线
    由上文可知,为构建更高阶的曲线,需要引入更多的中间点。同理可推得五次曲线:


image

五次贝塞尔曲线演示动画,t在[0,1]区间

构建贝塞尔曲线公式

  1. 一次曲线(线性):
    P0的坐标为(α,β)P1的坐标为(λ,μ),引入中间点Q0的坐标为(x,y),则有:

xαλx=t1t

可推得:

(1)x=(1t)α+tλ

同理,可得:

yβμy=t1t

可推得:

(2)y=(1t)β+tμ

整理可得:

(3)(xy)=(1t)(αβ)+t(λμ)

可将式(3)简写为:

(4)B(t)=(1t)P0+tP1,t[0,1]

得到一次贝塞尔曲线公式。
2. 二次曲线
P0,P1上的中间点为Q0P1,P2上的中间点为Q1Q0,Q1上的中间点为R0R0即构建曲线的连续点。
带入式(4)可得:

(5)Q0=(1t)P0+tP1Q1=(1t)P1+tP2R0=(1t)Q0+tQ1

将式(5)中的Q0Q1代入R0中,可得二次贝塞尔曲线公式:

(6)B(t)=(1t)2P0+2t(1t)P1+t2P2,t[0,1]

  1. 三次曲线
    同理,通过引入中间点Q0,Q1,Q2,R0,R1,S0,代入使用式(6),可得三次贝塞尔曲线公式:

(7)B(t)=(1t)3P0+3t(1t)2P1+3t2(1t)P2+t3P3,t[0,1]

  1. 高阶曲线
    可以迭代得到四次、五次贝塞尔曲线的公式:

(8)B(t)=(1t)4P0+4t(1t)3P1+6t2(1t)2P2+4t3(1t)P3+t4P4,t[0,1]

(9)B(t)=(1t)5P0+5t(1t)4P1+10t2(1t)3P2+10t3(1t)2P3+5t4(1t)P4+t5P5,t[0,1]

三、推导

一般化

观察上文中15阶贝塞尔曲线,发现其各项中的常数系数满足组合数规律,故对于n阶贝塞尔曲线,给出定点P0,P1,P2,...,Pn,其贝塞尔曲线可用下式表示:

(10)B(t)=(n0)(1t)nt0P0+(n1)(1t)n1t1P1++(nn1)(1t)1tn1Pn1+(nn)(1t)0tnPn,t[0,1]

可化简为:

(11)B(t)=i=0(ni)(1t)nitiPi,t[0,1]

式中组合数:

(ni)=n!i!(ni)!

多项式表达

同时,式(11)也可表示为:

(12)B(t)=i=0bi,n(t)Pi,t[0,1]

其中多项式:

(13)bi,n(t)=(ni)(1t)niti,i=0,1,2,,n

(13)被称作nBernsteain多项式

Bernsteain多项式
假设开展一个实验,其成功的概率为u。若连续进行这个实验n次,成功的次数恰好为i次,那么其成功概率是多少呢?
能够确定成功概率为:bi,n(t)=(ni)(1t)niti,i=0,1,2,,n.
这便是Bernsteain多项式。

四、应用(Matlab实现)

实现目标

  1. 能够进行控制点坐标的输入。
  2. 计算贝塞尔曲线。
  3. 绘制贝塞尔曲线及其控制点。

脚本编写

main.m(主函数,调用输入函数、贝塞尔曲线计算函数及绘制贝塞尔曲线)

%% 初始化
clc
clear

%% 输入
n = input("请输入控制点个数:\n");
P = pos_input(n);

%% 调用贝塞尔曲线计算函数
P = P';  % 坐标矩阵转置,便于曲线运算及绘图
[t, p] = bezier_curve(P, n - 1);  % 阶数为控制点个数 - 1

%% 绘制贝塞尔曲线
plot(t, p, '-b');  % 绘制贝塞尔曲线
hold on;
plot(P(1,:), P(2,:), '-ro');  % 绘制控制点
title(sprintf('贝塞尔曲线 (阶数: %d)', n - 1));
xlabel('X');
ylabel('Y');
grid on;
legend('贝塞尔曲线', '控制点');
hold off;

bezier_curve.m(贝塞尔曲线函数,计算n阶贝塞尔曲线)

function [Px, Py] = bezier_curve(control_points, order)
    % 计算贝塞尔曲线上的点
    t = linspace(0, 1, 1000);  % 在[0, 1]范围内生成点
    curve = zeros(2, numel(t));  % 初始化曲线矩阵,numel函数用于获取数组元素的数目
    for i = 0 : order
        curve = curve + nchoosek(order, i) * ((1 - t) .^ (order - i)) .* ...
        (t .^ i) .* control_points(:, i + 1);  % 该计算公式即上文中式(11),nchoosek函数用于计算组合数
    end
    Px = curve(1, :);
    Py = curve(2, :);
end

pos_input.m(输入函数,输入控制点坐标)

function P = pos_input(I)
   % 输入贝塞尔曲线控制点 
    P = ones(I, 2);  % 初始化坐标矩阵 
        for i = 1 : I
            P(i, :) = input("请输入输入坐标点P" + num2str(i - 1) + "(以一维矩阵形式输入): \n");
        end
end

图像结果

六阶贝塞尔曲线

六阶贝塞尔曲线.png

十四阶贝塞尔曲线

十四阶贝塞尔曲线.png

五、参考

  1. 维基百科_贝塞尔曲线
posted @   zbyisgudi  阅读(1534)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 在鹅厂做java开发是什么体验
· 百万级群聊的设计实践
· WPF到Web的无缝过渡:英雄联盟客户端的OpenSilver迁移实战
· 永远不要相信用户的输入:从 SQL 注入攻防看输入验证的重要性
· 全网最简单!3分钟用满血DeepSeek R1开发一款AI智能客服,零代码轻松接入微信、公众号、小程
点击右上角即可分享
微信分享提示