Courbes de Bézier

Option agrégation C
Christophe Ritzenthaler
October 28, 2009

Abstract
Résumé : on décrit ici les propriétés des courbes de Bézier et leur application au
lissage.

Il est rappelé que le jury n’exige pas une compréhension exhaustive du texte. Vous
êtes laissé(e) libre d’organiser votre discussion comme vous l’entendez. Des suggestions de
développement, largement indépendantes les unes des autres, sont proposées en fin de texte.
Vous n’êtes pas tenu(e) de les suivre. Il vous est conseillé de mettre en lumière vos connais-
sances à partir du fil conducteur constitué par le texte. Le jury demande que la discussion soit
accompagnée d’exemples traités par ordinateur. Il est souhaitable que vous organisiez votre
présentation comme si le jury n’avait pas connaissance du texte. Le jury aura néanmoins le
texte sous les yeux pendant votre exposé.

1 Premières définitions
1.1 Polynômes de Bernstein
Définition 1.1. Soit n un entier positif. On définit pour i ∈ {0, . . . , n}, les polynômes de
Bernstein par la formule  
n
Bi,n (X) = X i (1 − X)n−i .
i
Dans la suite on considérera uniquement ces polynômes sur [0, 1].

Proposition 1.1. On a les propriétés suivantes:

1. ∀i ∈ {0, . . . , n}, ∀u ∈ [0, 1], Bi,n (u) ≥ 0.
Pn
2. Partition de l’unité : i=0 Bi,n (u) = 1.

3. Bi,n atteint son maximum en u = i/n.

4. Définition récursive : Bi,n (u) = (1 − u)Bi,n−1 (u) + uBi−1,n−1 (u). On a convenu que
Bi,n = 0 si i < 0 ou i > n.

1
1.2 Courbes de Bézier
Définition 1.2. Pour n + 1 points du plan Pi = (x0 , y0 ), . . . , (xn , yn ) ∈ R2 , on définit la
courbe de Bézier C associé à l’ensemble de ces points par
( Pn
x(u) = i=0 Bi,n (u)xi
Pn
y(u) = i=0 Bi,n (u)yi
P
pour u ∈ [0, 1]. En abrégé, C(u) = Bi,n (u)Pi .
Les Pi sont appelés points de contrôle. Un exemple de courbe de Bézier avec 4 points de
controle est donné Fig.1.

Figure 1: Courbe de Bézier

Proposition 1.2 (Propriétés des courbes de Bézier). 1. Contrôle local : le point de contrôle
Pi influence la courbe au voisinage de u = i/n. De plus C(0) = P0 et C(1) = Pn .
2. Si on note Cn (P0 , . . . , Pn ) la courbe associée aux points de contrôle (Pi )0,...,n , on a la
définition récursive

Cn (P0 , . . . , Pn )(u) = (1 − u)Cn−1 (P0 , . . . , Pn−1 )(u) + uCn−1 (P1 , . . . , Pn )(u) (1)

3. La courbe est à l’intérieur de l’enveloppe convexe des points de contrôle.

4. La droite ((x0 , y0 ); (x1 , y1 )) (resp. ((xn−1 , yn−1 ); (xn , yn )) est tangente à la courbe en
(x0 , y0 ) (resp. (xn , yn )).
5. Une courbe de Bézier est infiniment dérivable.
6. Une ligne droite rencontre au plus autant de fois la courbe de Bézier que la ligne polyg-
onale joignant les points de contrôle (ceci signifie que la courbe suit assez fidèlement
cette ligne sans osciller).
Remarque 1. Le fait d’ajouter des points de contrôle permet d’influer sur la forme de la
courbe. Cependant il n’est pas judicieux d’élever trop le degré d’une courbe. On préfère
donc utiliser le cas n = 3 (courbe de Bézier cubique) et les relier entre elles de manière à
garantir une continuité C 1 . La courbe de Bézier cubique se trace en partant du point P0 ,
en se dirigeant vers P1 et en arrivant au point P3 selon la direction (P2 ; P 3). En général,
la courbe ne passe ni par P1 ni par P2 : ces points sont simplement là pour donner une
information de direction. La distance entre P0 et P1 détermine la ‘longueur’ du déplacement
dans la direction de P1 avant de tourner vers P3 .

2
2 Méthodes de calcul
2.1 Algorithme de De Casteljau
La formule (1) permet de construire un algorithme pour déterminer récursivement la po-
sition d’un point C(u) sur la courbe grâce à des interpolations linéaires en partant des
points de contrôle. Bien qu’aussi coûteux que l’algorithme de Horner d’évaluation directe
des polynômes, il est avantageux car moins sujet aux erreurs d’arrondis.

2.2 Algorithme de subdivision

Dans bien des applications, on a besoin de discrétiser la courbe en calculant les positions
de ses points pour des valeurs uk , k = 0, . . . , N . Pour pouvoir appliquer un algorithme de
dichotomie, on suppose N = 2p et uk = k/2p . Dans ce qui suit on va appliquer la méthode
pour une courbe de Bézier cubique mais elle s’étend facilement aux degrés plus élevés.
Tout commence par l’écriture matricielle de l’équation de la courbe
 
P0
3 P1 
 3 2 2
 
C(u) = (1 − t) 3t(1 − t) 3t (1 − t) t  
P2
P3
  
1 0 0 0 P0
3 −3 3 0 0
 2
 P1 
= 1 t t t    .
3 −6 3 0 P2 
1 3 −3 1 P3

Soit  
1 0 0 0
−3 3 0 0
M = .
3 −6 3 0
1 3 −3 1
On cherche à paramétrer la courbe C(t) avec t ∈ [0, 1/2] sous la forme d’une courbe de Bézier.
Il faut donc écrire Cg (t) = C(t/2), t ∈ [0, 1] :
   
1 0 0 0 P0
3 0 1/2 0 0 
 2
 P1 
Cg (t) = C(t/2) = 1 t t t   M  .
0 0 1/4 0  P2 
0 0 0 1/8 P3

Pour déterminer les points de contrôle de la courbe Cg , on doit écrire

 
P0
P1 
Cg (t) = 1 t t2 t3 M S 
 
P2  .


P3

Il faut donc que

   
1 0 0 0 1 0 0 0
0 1/2 0 0  1/2 1/2 0 0 
S = M −1 

M =  .
0 0 1/4 0  1/4 1/2 1/4 0 
0 0 0 1/8 1/8 3/8 3/8 1/8

3
Ainsi si les Pi,j sont les points de contrôle à l’étape j pour une des parties de la courbe, on
obtient, par symétrie, deux nouveaux ensembles de points de contrôle Pi,j+1 (pour la partie
0
[0, t/2]) et Pi,j+1 (pour la partie [t/2, 1])
     0   
P0,j+1 P0,j P0,j+1 1/8 3/8 3/8 1/8 P0,j
P1,j+1  P1,j  0
P1,j+1   0 1/4 1/2
 P1,j  .
1/4  
P2,j+1  = S P2,j  ,
     0 
P2,j+1   0 0 1/2 1/2 P2,j 
P3,j+1 P3,j 0 0 0 0 1 P3,j
P3,j+1
0 0
Rappelons que les points P0,j , P4,j = P0,j , P4,j appartiennent à la courbe. On arrête la
récurrence lorsque la distance entre ces points est plus petite qu’une borne fixée.

3 Application au lissage et généralisation

3.1 Lissage
Les courbes de Bézier composent l’outil de la base du dessin vectoriel qui repose sur la
transcription mathématique des objets. Les courbes de Bézier cubiques, les plus utilisées,
se retrouvent en graphisme et dans de multiples systèmes de synthèse d’images, tels que
PostScript, Metafont et The GIMP, pour dessiner des courbes ‘lisses’ joignant des points ou
des polygones de Bézier.
Ainsi, étant donné un chemin polyngonal Pi , on pourra approximer la tangente au point Pi
−−−−→ −−−−→
par le vecteur Pi Pi+1 renormalisé pour avoir une longueur ρ|Pi Pi+1 |. Le paramêtre ρ contrôle
en quelque sorte l’importance de la courbure au niveau des points de contrôle.
Cependant, on constate que cette méthode conduit à des phénomènes de ‘bosses’ après chaque
point de contrôle. D’où l’idée de calculer la tangente au point Pi comme la moyenne de la
−−−−→ −−−−→
direction précédente Pi−1 Pi et de la direction suivante Pi Pi+1 . De plus on renormalise avec
un coefficient ρ.

3.2 Courbes de Bézier rationnelles

Bien que pratiques et intuitives, les courbes de Bézier souffrent de nombreux problèmes. Le
plus important est qu’elles ne permettent pas de tracer des courbes aussi simples que des
arcs de cercles comme le précise le théorème suivant.
Théorème 3.1. Il n’existe pas de courbe paramétrée polynomiale qui décrive un arc de cercle.
Pour pallier à ce problème, on peut généraliser les courbes de Bézier.

Définition 3.1 (Courbes de Bézier rationnelles). On se donne n+1 points Pi et pour chaque
point un poids wi > 0. On définit alors
Pn
i=0 Bi,n (u)Pi
C(u) = P n , u ∈ [0, 1].
i=0 wi Bi,u

4 Suggestions
Soulignons qu’il s’agit d’un menu à la carte et que vous pouvez choisir d’étudier certains
points, pas tous, pas nécessairement dans l’ordre, et de façon plus ou moins fouillée. Vous

4
pouvez aussi vous poser d’autres questions que celles indiquées plus bas. Il est vivement
souhaité que vos investigations comportent une partie traitée sur ordinateur et, si possible,
des représentations graphiques de vos résultats.

1. On pourra démontrer les résultats de la Prop.1.1.

2. Quelle est l’utilité de Prop.1.1 (4) ?
3. On pourra détailler la remarque 1.
4. Commenter le paragraphe 2.1.
5. Implanter l’algorithme de subdivision.
6. On pourra illustrer la section 3.1 avec quelques exemples de courbes.
7. On pourra démontrer le Th.3.1.
8. Montrer comment tracer un arc de cercle avec une courbe de Bézier rationnelle.

5 Les remarques du jury pour les textes

Introduction et conclusion :
”Le jury n’ayant a priori pas lu le texte, le candidat commencera par présenter celui-ci. Un
plan en début d’exposé est apprécié, annonçant en particulier les propriétés du modèle que
le candidat va dégager. Il est important d’expliquer le problème et le modèle, de l’illustrer,
ainsi que d’y revenir en fin d’exposé. Le modèle mathématique a-t-il les propriétés attendues
? Des propriétés parasites surprenantes ? A-t-on résolu le problème posé ?”

La programmation :
”Il est vivement souhaité que des illustrations informatiques (simulation, résolution numérique
ou formelle, cas particuliers éclairants . . . ) soient présentées, mais il ne s’agit pas d’une
épreuve de programmation. Un programme qui ne fonctionne pas n’est en rien rédhibitoire
et le jury appréciera un regard critique du candidat sur une tentative non aboutie. Une
utilisation raisonnée des fonctions des logiciels disponibles est plus appréciée qu’une repro-
grammation d’algorithmes standards. Bien intégré dans l’exposé, un tel travail peut en
revanche devenir pertinent pour illustrer les insuffisances d’une méthode naı̈ve.”

Les démonstrations :
”S’il est exclu de plaquer une démonstration d’un théorème du programme dans l’exposé,
les démonstrations mathématiques de certaines assertions du texte sont très appréciées.
Lorsqu’une démonstration est ébauchée dans le texte, le candidat peut choisir de la compléter.
Il est alors particulièrement apprécié que le candidat précise les point mathématiques nécessaires
pour une démonstration rigoureuse. Le candidat peut, tout comme le texte, utiliser des ar-
guments heuristiques s’il les signale comme tels. Cependant le candidat ne doit pas oublier
qu’il s’agit d’une épreuve de l’agrégation externe de mathématiques, et qu’un exposé ne com-
portant aucun argument mathématique précis est vivement déconseillé.”

À ne pas faire :
”Le principal travers observé chez les candidats est la répétition linéaire du texte, y compris
des passages non compris en espérant que le jury ne demandera pas de détails. Rappelons
qu’utiliser des notions que l’on ne comprend pas, dans cette épreuve comme dans les autres,

5
est une faute lourdement sanctionnée. Enfin, rappelons qu’aucun développement n’est at-
tendu. Le candidat est libre de proposer des démonstrations de résultats utilisés,mais le jury
peut les refuser, ou demander au candidat d’en donner seulement les grandes lignes.”

À faire :
”Quelques qualités appréciées : prise de distance et d’initiative par rapport au texte, étude
d’un exemple ou d’un cas simple pour comprendre le texte et le faire comprendre au jury,
simplification ou, à l’inverse, généralisation du problème proposé, étude qualitative ou heuris-
tique, critique du modèle.”