Voir la source de Spline

{{Voir paronymes|Spleen}}

[[Fichier:Quadratic spline six segments.svg|vignette|Exemple de spline quadratique.]]

En [[mathématiques appliquées]] et en [[analyse numérique]], une '''spline''' est une [[fonction (mathématiques)|fonction]] définie par morceaux par des [[polynôme]]s.

''Spline'' est un terme anglais qui, lorsqu'il est utilisé en français, est généralement prononcé {{API|[splin]}}, à la française. Il désigne une réglette de bois souple appelée [[wikt:cerce|cerce]] en français. Toutefois, dans l'usage des mathématiques appliquées, le terme anglais ''spline''  est généralisé et le mot français ''cerce'' ignoré.

Dans les problèmes d'[[interpolation numérique|interpolation]], la méthode des splines est très souvent préférée à l'[[interpolation polynomiale]]. Les splines sont également utilisées dans les problèmes de [[Lissage (mathématiques)|lissage]] de données expérimentales ou de statistiques.

Les splines sont utilisées pour représenter numériquement des contours complexes. Leur mise en œuvre est simple. Elles sont fréquemment employées dans les logiciels de dessin ou de conception graphique ; leur usage y a été généralisé par [[Pierre Bézier]] avec les [[B-spline]]s.

== Historique ==

L'origine historique des splines se trouve dans la conception industrielle.

[[Fichier:Spline (PSF).png|vignette|Spline ou ''cerce''.]]
:Autrefois, avant l'existence des outils numériques, la conception industrielle était fondée sur la géométrie et le dessin. Les concepteurs d'une pièce de machine, par exemple, définissaient graphiquement les quelques points, ou ''nœuds'', par lesquels devait nécessairement passer le contour de la pièce ([[dessin technique]]). Le travail du dessinateur industriel consistait à relier ces nœuds par des courbes aussi lisses que possible afin d'éviter les discontinuités, sources de faiblesses mécaniques et de ruptures. Ces courbes devaient être « agréables à l'œil » : il y a un rapport direct entre la notion subjective d'« agrément » et la notion mathématique de [[classe de régularité]] d'une fonction, c'est-à-dire de ses propriétés de continuité et de dérivation.

Pour le tracé graphique des contours, on utilisait une réglette de bois élastique, appelée ''spline'' en anglais et ''cerce'' en français. À l'aide de tenons, on faisait passer la réglette par les nœuds et on traçait entre ceux-ci la courbe prise par la réglette. Les propriétés mécanique d'une réglette de bois sont telles que celle-ci prend, sous déformation, une forme qui a des dérivées continues à un ordre élevé (la dérivée première est continue, la [[dérivée seconde]] aussi, etc.)

=== Sur le terme ''spline'' ===

==== Prononciation ====

''Spline'' est un mot anglais qui, dans cette langue, se prononce {{API|[splaɪn]}}, ce qui le distingue de « {{page h'|Spleen|spleen}} » {{API|[spli:n]}}. En français, on prononce le mot tel qu'il s'écrit : {{API|[splin]}}.

==== Étymologie ====

[[Pierre Bézier]] écrit<ref>{{Article
 | auteur1                = Pierre Bézier
 | titre                  = Courbes et surfaces pour la CFAO
 | périodique             = Techniques de l'ingénieur Applications des mathématiques
 | volume                 = TIB102DUO
 | numéro                 = a1440
 | date                   = 1992
 | lire en ligne          = https://www.techniques-ingenieur.fr/base-documentaire/sciences-fondamentales-th8/applications-des-mathematiques-42102210/courbes-et-surfaces-pour-la-cfao-a1440/
 | consulté le            = 25 décembre 2017
}}.</ref> :
{{Citation bloc|Le mot anglais ''spline'' a malencontreusement prévalu, alors qu'il se traduit en français par les substantifs ''latte'' ou ''baguette''.}}

''The Oxford English Dictionnary 1989'' donne les définitions suivantes pour le mot ''spline'' :
{{Citation étrangère bloc|langue=en|A long narrow, and relatively thin piece or strip of wood, metal, etc.; a slat. A flexible strip of wood or hard rubber used by draftsmen in laying out broad sweeping curves, especially in railroad work.}}

''{{Citation bloc|Une pièce longue, étroite et relativement mince de bois, de métal, etc., une latte. Une bande flexible de bois ou de caoutchouc dur utilisée par les dessinateurs industriels dans le tracé de larges courbes, en particulier dans les travaux ferroviaires.}}''

Les menuisiers, charpentiers et escaliéteurs nommaient cet outil '''cerce''', latte souple de bois bien de fil servant à tracer des courbes « quelconques » (courbes harmonieuses passant par un certain nombre de points et qui ne peuvent être tracées à l'aide du compas).

=== Approximation mathématique des splines ===
La description de la forme prise par la réglette de bois, fondée sur les lois de l'élasticité, est mathématiquement compliquée. On obtient toutefois une approximation acceptable pour les usages pratiques en représentant la forme de la réglette par une ''spline'', c'est-à-dire une fonction définie par morceaux par un polynôme sur chaque intervalle entre nœuds.

Ces polynômes peuvent être calculés en posant des conditions de continuité : ils passent par les nœuds, la [[Tangente_(géométrie)|pente du contour]] est continue, c'est-à-dire que les polynômes à droite et à gauche du nœud ont même dérivée première ; de même, en chaque nœud, la [[courbure]] est continue, c'est-à-dire que les polynômes à droite et à gauche ont même dérivée seconde. Une solution est donnée par des morceaux de polynômes du troisième degré. C'est le degré minimum des polynômes pour que la fonction, sa pente et sa courbure soient continues à chaque nœud. Cette courbe par morceaux s'appelle ''spline''.

Les splines cubiques sont les plus utilisées dans la pratique. Elles sont plus faciles à utiliser que l'interpolation polynomiale et n'ont pas en général de comportements aberrants. Dans quelques cas particuliers, on peut être conduit à utiliser des polynômes du quatrième degré ; toutefois cet usage résulte de la pratique des analystes mais nullement de considérations théoriques.

D'autres familles de fonctions splines (obtenues par minimisation d'une fonctionnelle d'énergie) ont également été développées dans le cadre de l'approximation de données.

== Définition ==
Soit une série de points de donnée, ou ''nœuds'', de coordonnées {{math|(''X{{ind|i}}'', ''Y{{ind|i}}'')}}. On appelle ''spline'' d'interpolation la fonction par morceaux constituée d'un polynôme sur chaque intervalle entre nœuds.

Le degré de la spline est défini par le polynôme de plus haut degré utilisé : si on joint simplement les nœuds par des droites, c'est-à-dire si la spline est constituée de polynômes du premier degré, la spline est de degré 1, etc. Si tous les polynômes ont même degré, on parle de ''spline uniforme''. Si tous les polynômes sont du troisième degré, on parle de ''spline cubique''.

=== Exemples ===
Voici l'exemple de splines d'interpolation de degré 1, d'une part, et de degré 3, d'autre part, fondées sur les quatre mêmes nœuds (les points rouges des graphiques ci-dessous.

À gauche, la spline de degré 1 relie les quatre nœuds par des polynômes de degré 1, c'est-à-dire par des droites. La spline est continue, mais la pente n'est pas continue (elle passe brutalement d'une valeur à une autre). Les polynômes, des fonctions linéaires, permettent, entre les nœuds, des interpolations linéaires.

À droite, une spline cubique dans laquelle les nœuds sont reliés par des polynômes du {{3e}} degré auxquels on a imposé des conditions pour rendre la spline continue au premier degré (les pentes varient de façon continue) et au second degré (la courbure de la spline, plus précisément sa dérivée seconde, varie de façon continue).
{{Gallery
 | title = Exemples de splines
 | lines = 5
 | width = 200
 | height= 250
 | align = center
 | File:Exemple de spline de degré 1.png
 |Exemple de spline de degré 1 pour quatre nœuds.
 | File:Exemple de spline de degré 3.png
 |Exemple de spline de degré 3 pour quatre nœuds.
}}

=== Continuité de la spline ===
Comme les polynômes sont continus, indéfiniment dérivables et que leurs dérivées successives sont continues, la continuité de la spline dépend des conditions de continuité qu'on impose en chaque nœud. La continuité définit les caractéristiques de la jonction entre chaque intervalle. Cela correspond au degré de correspondance entre deux polynômes successifs aux points de jonction.
* {{math|''C''{{ind|0}}}} est la continuité minimum : les polynômes successifs passent bien par les points de jonction ;
* {{math|''C''{{ind|1}}}} est la continuité des tangentes : les polynômes successifs ont des dérivées premières (c'est-à-dire des pentes) égales aux points jonction ;
* {{math|''C''{{ind|2}}}} est la continuité de la courbure :  les polynômes successifs ont de plus des dérivées secondes (des courbures) égales aux points jonction ;
:L'application des conditions {{math|''C''{{ind|0}}}}, {{math|''C''{{ind|1}}}} et {{math|''C''{{ind|2}}}} conduit aux splines cubiques.

== Spline d'interpolation ==
Lorsqu'on a déterminé les éléments de la spline, c'est-à-dire les paramètres des polynômes qui joignent chacun des points donnés (les ''nœuds''), on peut calculer l'ordonnée d'un point situé entre deux nœuds : c'est la valeur du polynôme pour l'abscisse correspondante. La spline permet donc d'interpoler les valeurs comprises entre les nœuds.

=== Inconvénients des splines ===
Quoiqu'elles soient plus stables que d'autres méthodes d'interpolation, les splines n'en présentent pas moins quelques inconvénients:

* l'évaluation des paramètres de la spline peuvent diverger si la pente entre deux nœuds varie fortement d'un intervalle au suivant ;
* la spline dépend du choix du système de coordonnées, elle ne possède donc pas de propriété d’invariance géométrique (cas d'une rotation par exemple) ; toutefois, elle reste insensible à une [[homothétie]] sur l'un des axes ;
* on ne peut utiliser l'interpolation par spline si la courbe à représenter n'est pas une fonction, c'est-à-dire s'il y a plus d'une valeur pour une abscisse donnée (courbes qui « reviennent en arrière » comme les ellipses, les cercles, etc.) Dans ce cas, on utilise des généralisations des splines ([[B-spline]]s, [[NURBS]]).

=== Calcul des éléments d'une spline cubique ===

Considérons ''n'' nœuds de coordonnées {{math|(''X''{{ind|1}}, ''Y''{{ind|1}}), ..., (''X{{ind|n}}'', ''Y{{ind|n}}'')}}. Il y a ''n''-1 intervalles. Nous voulons donc déterminer les paramètres de ''n''-1 polynômes du troisième degré, {{math|''S''{{ind|1}}, ..., ''S''{{ind|''n''–1}}}}.

Chaque polynôme {{mvar|S}} a une équation de la forme :
:<math>S = a_0 x^3 + a_1 x^2 + a_2 x + a_3</math>
Il faut donc déterminer quatre coefficients par polynôme, soit au total 4''n''-4 inconnues, pour lesquelles on doit poser autant de conditions.

Les conditions {{math|''C''{{ind|0}}}}, {{math|''C''{{ind|1}}}} et {{math|''C''{{ind|2}}}} s'expriment de la façon suivante :
* {{math|''C''{{ind|0}}}} : chacun des ''n''-1 polynômes passe par deux nœuds, soit 2''n''-2 conditions ;
* {{math|''C''{{ind|1}}}} : les polynômes ont les mêmes dérivées premières aux ''n''-2 nœuds de jonction, soit ''n''-2 conditions ;
* {{math|''C''{{ind|2}}}} : et ces polynômes ont la même dérivée seconde aux ''n''-2 nœuds de jonction, soit ''n''-2 conditions.

Au total on a ainsi 4''n''-6 conditions. Il manque donc deux conditions pour déterminer les 4''n''-4 inconnues. On adjoint deux conditions supplémentaires en imposant la valeur des dérivées secondes à chaque nœud extrême. Si on impose que les dérivées secondes soient nulles à chaque extrémité, on obtient une ''spline cubique naturelle''. On peut disposer toutefois d'informations supplémentaires qui permettent de fixer de façon plus précise la dérivée seconde aux extrémités.

On obtient ainsi un système de 4''n''-4 équations linéaires à 4''n''-4 inconnues.

=== Algorithme de calcul ===
On trouvera ci-dessous un algorithme de calcul d'une spline cubique d'interpolation, naturelle ou non. Cet algorithme est présenté en langage naturel et mathématique, sans référence à quelque [[langage de programmation]] que ce soit. Il peut être mis en œuvre avec un simple [[tableur]].

{{Boîte déroulante/début|titre=Algorithme de calcul d'une spline d'interpolation}}

[[File:Éléments de calcul d'une spline d'interpolation.png|vignette|Éléments de calcul d'une spline d'interpolation naturelle.]]

:Soit {{mvar|X{{ind|i}}}} et {{mvar|Y{{ind|i}}}} les abscisses et ordonnées de {{mvar|n}} nœuds, les abscisses étant classées par ordre croissant.

:1. Calculer la matrice colonne {{math|'''h'''}} des largeurs des {{math|''n''–1}} intervalles entre nœuds :

:::<math>\mathbf{h}_i = X_{i+1} - X_i</math>, pour {{mvar|i}} de 1 à {{math|''n''–1}}.

:2. Calculer la matrice colonne {{math|'''F'''}} des {{math|''n''–2}} dérivées secondes aux nœuds :

::2.1. Pour les nœuds 2 à {{math|''n''–1}} :

:::<math>\mathbf{F}_i = \frac {Y_{i+1} - Y_i}{\mathbf{h}_i} - \frac {Y_i - Y_{i-1}}{\mathbf{h}_{i-1}}</math>, pour {{mvar|i}} de 2 à {{math|''n''–1}}.

::2.2. Il faut de plus valoriser les dérivées secondes aux extrémités, {{math|'''F'''{{ind|1}}}} et {{math|'''F'''{{ind|''n''}}}}.
:::Si l'on pose {{math|'''F'''{{ind|1}} {{=}} 0}} et {{math|'''F'''{{ind|''n''}} {{=}} 0}}, on obtient une spline naturelle ;
:::si l'on possède des informations sur les extrémités, on affecte à {{math|'''F'''{{ind|1}}}} et/ou à {{math|'''F'''{{ind|''n''}}}} la valeur ''ad hoc''.

:3. Construire la matrice carrée {{math|'''R'''{{ind|''i'',''j''}}}}, de dimensions {{math|''n''×''n''}}, telle que, {{mvar|i}} désigne l'indice des lignes et {{mvar|j}} celui des colonnes :
::3.1 Mettre tous les éléments de la matrice {{math|'''R'''}} à 0.
::3.2 Sur la première ligne, seul le premier élément n'est pas nul : {{math|'''R'''{{ind|1,1}} {{=}} 1}}.
::3.3 Sur la dernière ligne, seul le dernier élément n'est pas nul : {{math|'''R'''{{ind|''n'',''n''}} {{=}} 1}}.

::3.4 La sous-matrice de {{math|'''R'''}} de dimension {{math|(''n''–2) × ''n''}} en excluant la première et la dernière ligne est tridiagonale. Tous ses éléments sont nuls, sauf :

::* La diagonale centrale de la sous-matrice vaut :
:::<math>\mathbf{R}_{i,i} = \frac{h_{i-1} + h_i}{3}</math> pour {{math|''i'' {{=}} 2, ..., ''n''–1}}.

::*La diagonale supérieure vaut :
:::<math>\mathbf{R}_{i,{i+1}} =\frac{h_i}{6}</math> pour {{math|''i'' {{=}} 2, ..., ''n''–1}}.

::*La diagonale inférieure vaut :
:::<math>\mathbf{R}_{i,{i-1}} =\frac{h_{i-1}}{6}</math> pour {{math|''i'' {{=}} 2, ..., ''n''–1}}.

:4. Calculer la matrice colonne {{math|'''M'''}} de dimension {{mvar|n}}, telle que :

:::<math>\mathbf{M} =  \mathbf{R}^{-1} \times \mathbf{F}</math>.

:5. Calculer les matrices colonnes {{math|'''C'''}} et {{math|'''C{{'}}'''}}, de dimensions {{math|''n''–1}}, contenant respectivement les constantes d'intégration d'ordre 1 et 2 des {{math|''n''–1}} polynômes, de la façon suivante :

::Pour {{mvar|i}} de 1 à {{math|''n''–1}} :

:::<math>\mathbf{C}_i = \frac{Y_{i+1} - Y_i}{\mathbf{h}_i} - \frac{\mathbf{h}_i}{6}(\mathbf{M}_{i+1} - \mathbf{M}_i)</math>
:::<math>\mathbf{C'}_i = Y_i  - \mathrm{M}_i \frac{\mathbf{h}_i^2}{6}</math>

:6. Calcul d'une valeur d’interpolation pour une abscisse {{mvar|x}} située dans l'intervalle {{mvar|k}}, entre les nœuds d'abscisses {{mvar|X}}{{ind|k}} et {{math|''X''{{ind|''k''+1}}}} (avec {{mvar|k}} de 1 à {{math|''n''–1}}) :

:::<math>\begin{array}{lcl}y & = & \mathbf{M}_k \frac{(X_{k+1} - x)^3}{6\mathbf{h}_k} \\
                             & + & \mathbf{M}_{k+1} \frac{(x - X_k)^3}{6\mathbf{h}_k} \\
                             & + & \mathbf{C}_k (x - X_k) \\
                             & + & \mathbf{C'}_k
         \end{array}</math>

'''Exemple'''
[[File:Interpolation de la fonction y = sin x par une spline cubique naturelle.png|thumb|Interpolation de la fonction y = sin x par une spline cubique naturelle à l'aide de cinq nœuds. La courbe {{math|''y'' {{=}} sin ''x''}} est en rouge et la spline d'interpolation en bleu.]]

Cet exemple présente l'approximation de la courbe {{math|''y'' {{=}} sin ''x''}} par une spline cubique naturelle pour laquelle on donne les cinq valeurs suivantes (5 nœuds) de la fonction :

:::<math>n=5</math>

:::<math>\begin{array}{rrr} i & X_i   & Y_i      \\
                            1 & 0{,}0 & 0{,}0000 \\
                            2 & 0{,}5 & 0{,}4794 \\
                            3 & 1{,}0 & 0{,}8415 \\
                            4 & 1{,}5 & 0{,}9975 \\
                            5 & 2{,}0 & 0{,}9093 \\
\end{array}</math>

:::<math>\mathbf{h} = \begin{bmatrix} 0{,}5 \\
                                      0{,}5 \\
                                      0{,}5 \\
                                      0{,}5 \\
\end{bmatrix}
\text{ ; }</math><math>
         \mathbf{F} = \left[ \begin{array}{r} 0{,}0000 \\
                                             -0{,}2348 \\
                                             -0{,}4120 \\
                                             -0{,}4884 \\
                                              0{,}0000 \\
\end{array} \right]
\text{ ; }</math><math>
         \mathbf{R} = \left[ \begin{array}{rrrrr} 1{,}0000 & 0        & 0        & 0        & 0        \\
                                                  0{,}0833 & 0{,}3333 & 0{,}0833 & 0        & 0        \\
                                                  0        & 0{,}0833 & 0{,}3333 & 0{,}0833 & 0        \\
                                                  0        & 0        & 0{,}0833 & 0{,}3333 & 0{,}0833 \\
                                                  0        & 0        & 0        & 0        & 1{,}0000 \\ \end{array} \right] \text{ ; }</math>

:::<math>\mathbf{M} = \left[ \begin{array}{r}  0{,}0000 \\
                                              -0{,}5061 \\
                                              -0{,}7928 \\
                                              -1{,}2671 \\
                                               0{,}0000 \\
\end{array} \right] \text{ ; }</math><math>
         \mathbf{C} = \left[ \begin{array}{r}  1{,}0010 \\
                                               0{,}7480 \\
                                               0{,}3516 \\
                                              -0{,}2820\\
\end{array} \right] \text{ ; }</math><math>
         \mathbf{C'}= \left[ \begin{array}{r}  0{,}0000 \\
                                               0{,}5005 \\
                                               0{,}8745 \\
                                               1{,}0503 \\
\end{array} \right]</math>

''Commentaire sur l'exemple'' : La fonction <math>sin(i)</math> est bien représentée par son approximation par une spline naturelle à gauche du graphique parce que la dérivée seconde de <math>sin(i)</math> est bien nulle pour <math>i = 0</math>, comme le suppose le calcul de la spline naturelle. En revanche, à droite du graphique, l'approximation par une spline naturelle est moins
exacte, en effet, la dérivée seconde de <math>sin(i)</math> n'est pas nulle pour <math>X_5 = 2{,}0</math> (elle vaut <math>-0,035</math>), valeur différente de l'hypothèse de nullité posée pour la détermination de la spline naturelle.

{{Boîte déroulante/fin}}

=== Forme explicite d'une spline cubique ===
Pour une fonction <math>f()</math> dont on connaît les valeurs aux nœuds <math>X_i</math>, l'expression générale pour la spline cubique de classe ''C''<sup>2</sup> en un point <math>x</math> dans <math>[X_{i-1}, X_i]</math> est la suivante :

:<math>S_i(x)=\frac{z_i(x-X_{i-1})^3}{6h_i}+\frac{z_{i-1}(X_i-x)^3}{6h_i}+\left[\frac{f(X_i)}{h_i}-\frac{z_ih_i}{6}\right](x-X_{i-1})+\left[\frac{f(X_{i-1})}{h_i}-\frac{z_{i-1}h_i}{6}\right](X_i-x)</math>

où
* <math>f(X_i^{})</math> est la valeur de <math>f()</math> au nœud i,
* <math>h_i^{}=X_i\;\text{-}\;X_{i-1}</math>,
* <math>z_i=f^{\prime\prime}(X_i)</math> est la dérivée seconde de <math>f()</math> au nœud i.

On vérifie facilement que <math>S_i(X_{i-1})=f(X_{i-1})</math> et <math>S_i(X_i)=f(X_i)</math>.

=== Mise en œuvre dans des logiciels ===

[[Fichier:InkscapeBezierHandles.png|vignette|B-splines (courbes de Bézier) utilisées dans Inkscape.]]

Les splines sont souvent intégrées dans des logiciels de dessin, que ce soient des logiciels généralistes, comme les fonctions de dessin des suites de logiciels de [[bureautique]] ([[Microsoft Office]], [[LibreOffice]]…), des logiciels de dessin et de traitement d'image ([[Inkscape]], [[GIMP]]…) ou bien des logiciels de [[dessin assisté par ordinateur|dessin]] ou de [[Conception assistée par ordinateur|conception]] assistés par ordinateur.
Ces logiciels utilisent le plus souvent des [[B-spline]]s.

== Spline de lissage ==
=== Utilisation de splines pour le lissage ===
Les données de mesure sont généralement affectées d'aléas, que ce soit dans les sciences physiques ou dans les sciences humaines (économie, démographie, sociologie, etc.) Pour approcher empiriquement une fonction inconnue dont les mesures sont perturbées par des aléas, on pratique souvent le [[Lissage (mathématiques)|lissage]]. L'idée {{page h'|Heuristique|heuristique}} sous-jacente est que les fonctions qui représentent les phénomènes physique ou en sciences humaines sont généralement très continues : le lissage a donc pour objet de substituer des valeurs continues à des valeurs de mesure dispersées à cause des aléas. Ce problème se pose en particulier pour des mesures d'une grandeur évoluant dans le temps ([[série temporelle|série chronologique]]).

[[File:Spline Smoothing.png|thumb|Principe du lissage par spline : des ressorts lient les nœuds à la réglette qui décrit la courbe lissée ; la forme de celle-ci minimise l'énergie totale de déformation.]]
Il existe de nombreuses méthodes empiriques de lissage, telles que les [[moyenne mobile|moyennes mobiles]], le [[lissage exponentiel]] simple ou double, etc. Les splines de lissage font partie de ces méthodes empiriques.

Considérons des données de mesure dont les abscisses sont strictement croissantes (ce sont souvent des données chronologiques).

L'idée conduisant aux splines de lissage repose sur l'utilisation de la même réglette flexible que celle qui est à l'origine des splines d'interpolation. Toutefois, au lieu d'imposer à la réglette de passer par les nœuds (les valeurs de mesure), on suppose qu'on relie ceux-ci à la cerce par des ressorts. La réglette décrit alors une fonction très continue qui minimise l'énergie totale de déformation des ressorts et de la réglette.

=== Forme d'équilibre d'une fonction de lissage ===
Soit {{mvar|n}} valeurs connues (nœuds), dont les abscisses sont strictement croissantes, que l'on souhaite lisser. On relie ces {{mvar|n}} nœuds à une réglette flexible de raideur {{mvar|r}} par des ressorts ayant tous la même raideur {{mvar|k}}. La forme d'équilibre est celle qui minimise l'énergie totale d'allongement des ressorts et de courbure de la réglette.

==== Énergie d'allongement des ressorts ====
Chaque ressort {{mvar|i}} s'allonge d'une longueur {{mvar|e{{ind|i}}}}. L'énergie d'allongement d'un ressort élastique est proportionnelle au carré de l'allongement ; l'énergie d'allongement du ressort {{numéro|{{mvar|i}}}} vaut :
:<math>w_i = ke_i^2</math> où {{mvar|k}} est la raideur du ressort.

L'énergie totale d'allongement des ressorts est la somme des énergies d'allongement des {{mvar|n}} ressorts :
:<math>W_R =  \sum w_i= \sum k e_i^2</math> (en supposant que la raideur de tous les ressorts est identique et égale à {{mvar|k}}).

==== Énergie de déformation de la réglette ====
L'énergie élastique de courbure d'une cerce élastique est proportionnelle en tout point au carré de la courbure en ce point. On approxime la courbure par la dérivée seconde de la fonction {{math|''s''(''x'')}} représentant la forme de la réglette. Cette approximation est admissible en pratique si la pente de la fonction n'est pas trop grande (>>1). Alors l'énergie de courbure entre {{mvar|x}} et {{math|''x'' + ''dx''}} vaut :
:<math>dW_C = r s''(x)^2 dx</math> où {{mvar|r}} est la raideur de courbure élastique de la cerce.

L'énergie de déformation pour toute la réglette vaut :
:<math>W_C = \int r s''(x)^2\,\mathrm{d}x</math> si l'on suppose qu'elle présente une raideur constante sur toute sa longueur.

==== Condition sur la forme prise par la réglette ====
La forme prise par la réglette est celle qui minimise l'énergie totale de déformation, c'est-à-dire l'énergie d'allongement des ressorts et l'énergie de courbure de la réglette. C'est la fonction qui minimise l'expression :

:<math>W(s, \lambda) = \sum e_i^2 + \lambda \int s''(x)^2dx</math>
:où {{mvar|λ}} est le rapport entre la raideur {{mvar|r}} de la réglette et la raideur {{mvar|k}} des ressorts.

=== Spline cubique de lissage ===

Trouver la fonction {{mvar|s}} qui minimise {{math|''W''(''s'',''λ'')}} pour {{mvar|λ}} donné se simplifie si l'on suppose que cette fonction est une spline cubique : courbe par morceaux constituée de polynômes du troisième degré. Si l'on doit lisser {{mvar|n}} valeurs (ou nœuds), il y a {{math|''n''–1}} intervalles et donc {{math|''n''–1}} polynômes du troisième degré définis chacun par quatre paramètres, soit au total {{math|4''n''–4}} paramètres à déterminer.

On impose à ces [[Fonction polynomiale|fonctions polynomiales]] d'être continues à chaque jonction, d'y avoir une pente (dérivée première) et une dérivée seconde continue, soit :
* {{math|''C''{{ind|0}}}} : chacun des {{math|''n''–1}} polynôme est jointif au suivant, soit {{math|''n''–2}} conditions ;
* {{math|''C''{{ind|1}}}} : les polynômes ont même dérivée première aux {{math|''n''–2}} jonctions, soit {{math|''n''–2}} conditions ;
* {{math|''C''{{ind|2}}}} : et ces polynômes ont même dérivée seconde aux {{math|''n''–2}} jonctions, soit {{math|''n''–2}} conditions.

Soit {{math|3''n''–6}} conditions.

La minimisation de {{math|''W''(''s'',''λ'')}}, où {{mvar|λ}} est fixé, correspond à {{mvar|n}} conditions. Soit donc au total {{math|4''n''–6}} conditions. Il manque deux conditions pour déterminer les {{math|4''n''–4}} inconnues. On impose que les dérivées secondes aux extrémités soient nulles. L'ensemble de ces conditions conduit à {{math|4''n''–4}} équations linéaires à {{math|4''n''–4}} inconnues. (Si l'on dispose d'informations supplémentaires pour les courbures aux points extrêmes, on peut utiliser ces valeurs connues au lieu de la valeur 0.)

=== Algorithme de calcul ===
On trouvera ci-dessous un algorithme de calcul d'une spline cubique de lissage. Cet algorithme est présenté en langage naturel et mathématique, sans référence à quelque langage de programmation que ce soit. Il peut être mis en œuvre avec un simple tableur.

{{Boîte déroulante/début|titre=Algorithme de calcul d'une spline cubique de lissage}}

'''Algorithme de calcul d'une spline cubique de lissage'''<ref name=Moulines>{{Lien web
 | format           = PDF
 | auteur1          = E. Moulines
 | url              = http://www.math.univ-angers.fr/~loustau/rnp.pdf
 | titre            = Étude de cas de régression non paramétrique
}}.</ref>.

[[File:Éléments de calcul d'une spline d'interpolation.png|vignette|Éléments de calcul d'une spline de lissage.]]

: Soit {{mvar|X{{ind|i}}}} et {{mvar|Y{{ind|i}}}} les abscisses et ordonnées de {{mvar|n}} nœuds, les abscisses étant classées par ordre croissant {{math||§=(''n'' > 2)}}.

:1. Calculer la matrice colonne '''h''' {{math||§=des largeurs des}}  {{math||§=''n''–1}} intervalles entre nœuds :
::<math>\mathbf{h}_i = X_{i+1} - X_i</math>, pour {{mvar|i}} de 1 à {{math||§=''n''–1}}.

:2. Construire la matrice tridiagonale {{math||§='''Q'''}} ({{math||§=''n'' × (''n''–2)}}) telle que toutes les valeurs de la matrice sont nulles, sauf :
:*La diagonale centrale :
::<math>\mathbf{Q}_{i,{i-1}}=\frac{1}{h_{i-1}}+\frac{1}{h_i}</math>, pour {{math||§=''i'' {{=}} 2}} à {{math||§=''n''–1}} ;
:*La diagonale supérieure :
::<math>\mathbf{Q}_{i,i}=-\frac{1}{h_i}</math>, pour {{math||§=''i'' {{=}} 1}} à {{math||§=''n''–2}} ;
:*La diagonale inférieure :
::<math>\mathbf{Q}_{i,{i-2}}=-\frac{1}{h_{i-1}}</math>, pour {{math||§=''i'' {{=}} 3}} à {{mvar|n}} ;

:3. Construire la matrice carrée tridiagonale {{math||§='''R'''}} ({{math||§=(''n''–2)×(''n''–2)}}) telle que toutes les valeurs de la matrice sont nulles, sauf :
:*La diagonale principale :
::<math>\mathbf{R}_{i,i}=\frac{h_i+h_{i+1}}{3}</math>, pour {{math||§=''i'' {{=}} 1}} à {{math||§=''n''–2}} ;
:*La diagonale supérieure :
::<math>\mathbf{R}{i,{i+1}}=-\frac{h_i}{6}</math> pour {{math||§=''i'' {{=}} 1}} à {{math||§=''n''–3}} ;
:*La diagonale inférieure :
::<math>\mathbf{R}{i,{i-1}}=-\frac{h_{i-1}}{6}</math> pour {{math||§=''i'' {{=}} 2}} à {{math||§=''n''–2}} ;

:4. Calculer la matrice carrée {{math||§='''K'''}} ({{math||§=''n'' × ''n''}}) :
::<math>\mathbf{K} = \mathbf{Q}\mathbf{R}^{-1}\mathbf{Q}^t</math>;

:5. Soit <math>\mathbf{I}</math> la matrice unitaire {{math||§=''n'' × ''n''}} et <math>\lambda</math> le paramètre de lissage.

:On obtient les ordonnées {{math||§='''s'''}} de la spline de lissage en chacun des {{mvar|n}} nœud d'abscisse {{math||§='''X'''}} et d'ordonnée {{math||§='''Y'''}} pour la valeur <math>\lambda</math> du paramètre de lissage par :
::<math>\mathbf{s} = (\mathbf{I}+\lambda\mathbf{K})^{-1}\mathbf{Y}</math>.

:6. On peut interpoler la valeur de la spline de lissage pour toute abscisse comprise entre les nœuds extrêmes {{math||§='''X'''{{ind|1}}}} et {{math||§='''X'''{{ind|''n''}}}} en appliquant l'algorithme d'interpolation des splines aux points d'abscisse {{math||§='''X'''}} et d'ordonnée {{math||§='''s'''}}.

'''Exemple''' :
[[File:Spline de lissage cubique de quatre points avec lambda {{=}}1.5.png|thumb|center|Lissage des quatre nœuds figurés en rouge avec une spline cubique de lissage de paramètre λ {{=}} 1,5.]]

Cet exemple présente le lissage des quatre nœuds de coordonnées {{math||§=(''X{{ind|i}}'', ''Y{{ind|i}}'')}}, figurés en rouge dans le graphique ci-dessus, par une spline cubique naturelle avec le paramètre de lissage {{math||§=''λ'' {{=}} 1,5}}.

:::<math>\begin{array}{rrr} i & X_i   & Y_i      \\
                            1 & 0{,}0 & 0{,}0 \\
                            2 & 1{,}0 & 6{,}0 \\
                            3 & 3{,}0 & 5{,}0 \\
                            4 & 5{,}0 & 0{,}0 \\
\end{array}</math> ;

:::<math>\mathbf{h} = \begin{bmatrix} 1{,}0 \\
                                      2{,}0 \\
                                      2{,}0 \\
\end{bmatrix}
\text{ ; }</math><math>
\mathbf{Q} = \begin{bmatrix}\begin{array}{rr} -1{,}0000 &  0{,}0000 \\
                                               1{,}5000 & -0{,}5000 \\
                                              -0{,}5000 &  1{,}0000 \\
                                               0{,}0000 & -0{,}5000 \\
\end{array}
\end{bmatrix}
\text{ ; }</math><math>
\mathbf{R} = \begin{bmatrix}\begin{array}{rr}
              1{,}0000 & -0{,}1667 \\
             -0{,}1667 &  1{,}3333 \\
\end{array}
\end{bmatrix}</math> ;

:::<math>\mathbf{K} = \mathbf{Q}\mathbf{R}^{-1}\mathbf{Q}^t = 
\begin{bmatrix}\begin{array}{rrrr}
                      1{,}0213 & -1{,}4681 &  0{,}3830 &  0{,}0638 \\
                     -1{,}4681 &  2{,}2979 & -0{,}9255 &  0{,}0957 \\
                      0{,}3830 & -0{,}9255 &  0{,}8936 & -0{,}3511 \\
                      0{,}0638 &  0{,}0957 & -0{,}3511 &  0{,}1915 \\
\end{array}
\end{bmatrix}</math> ;

:::{{math||§='''I'''}} étant la matrice unité {{math||§=''n'' × ''n''}} et avec le paramètre {{mvar|λ}}, on calcule :
:::<math>(\mathbf{I}+ \lambda\mathbf{K}) =
\begin{bmatrix}\begin{array}{rrrr}
                      2.53191 & -2.20213 &  0.57447 &  0.09574 \\
                     -2.20213 &  4.44681 & -1.38830 &  0.14362 \\
                      0.57447 & -1.38830 &  2.34043 & -0.52660 \\
                      0.09574 &  0.14362 & -0.52660 &  1.28723 \\
\end{array}
\end{bmatrix}
\text{, et : }</math><math>(\mathbf{I}+ \lambda\mathbf{K})^{-1} = 
\begin{bmatrix}\begin{array}{rrrr}
                      0.7051 &  0.3583 &  0.0206 & -0.0840 \\
                      0.3583 &  0.4599 &  0.1843 & -0.0026 \\
                      0.0206 &  0.1843 &  0.5800 &  0.2152 \\
                     -0.0840 & -0.0026 &  0.2152 &  0.8714 \\
\end{array}
\end{bmatrix}</math>.

:::d'où l'on calcule <math>\mathbf{s} = (\mathbf{I} + \lambda \mathbf{K})^{-1}\mathbf{Y}</math>

{{Boîte déroulante/fin}}

=== Paramètre de lissage ===
On suppose généralement que les ressorts qui tirent la cerce ont tous la même raideur {{mvar|k}}. Leur donner des raideurs différentes revient à pondérer les nœuds en considérant que certains sont plus ou moins importants, ou plus ou moins "attractifs" que d'autres. C'est l'expérience des analystes qui permet d'établir ces pondérations.

Si tous les ressorts ont même coefficient de raideur {{mvar|k}}, Le paramètre de lissage, noté souvent {{mvar|λ}}, exprime le rapport entre la raideur élastique de la réglette et celle des ressorts. {{mvar|λ}} peut varier entre 0 et {{math|+∞}} :
* si {{math|''λ'' {{=}} 0}}, les ressorts sont "infiniment" durs par rapport à la cerce. Dans ce cas, la spline passe par les nœuds et l'on obtient une ''spline naturelle'' ;
* si {{mvar|λ}} tend vers {{math|+∞}}, cela revient à supposer que la cerce devient "infiniment" raide par rapport aux ressorts. Dans ce cas, la spline de lissage tend vers la ''droite de régression'' des nœuds.

{{Gallery
 | title = Exemples de splines de lissage cubiques
 | lines = 3
 | width = 200
 | height= 200
 | align = center
 | File:Spline de lissage cubique de quatre points avec lambda {{=}}0.3.png
 |Lissage de quatre nœuds avec une spline cubique de lissage de paramètre λ {{=}} 0,3.
 | File:Spline de lissage cubique de quatre points avec lambda {{=}}1.5.png
 |Lissage de quatre nœuds avec une spline cubique de lissage de paramètre λ {{=}} 1,5.
}}

Pour mieux apprécier l'intensité du lissage, on utilise souvent le paramètre {{mvar|p}}, tel que :

::<math>p=\frac{1}{\lambda+1}</math>, c'est-à-dire :
::<math>\lambda=\frac{1-p}{p}</math>, où {{mvar|p}} varie entre 0 (exclu) et 1.

::{{math|''p'' {{=}} 1}} correspond à {{math|''λ'' {{=}} 0}} : on obtient la spline naturelle passant par les nœuds ;
::{{math|''p''→ 0{{ind|+}}}} correspond à {{math|''λ'' → +∞}} : la spline obtenue tend vers la droite de régression.

=== Optimisation du paramètre de lissage ===
Lorsqu'on étudie de nombreuses données de mesure bruitées, l'objectif du lissage est de fournir une estimation non bruitée de la mesure. Si l'on peut estimer la variance du bruit, la variance des écarts entre données brutes et les données lissées comparée à la variance du bruit peut guider le choix du paramètre de lissage.

On applique ce procédé de façon plus systématique de la façon suivante :

Le paramètre de lissage étant donné, supprimons le nœud {{numéro|{{mvar|i}}}} de la série des mesures et déterminons la spline de lissage de paramètre {{mvar|λ}} avec ce nœud exclu du calcul. L'écart {{mvar|ε{{ind|i}}}} pour le nœud {{numéro|{{mvar|i}}}} entre la valeur brute et la valeur de la spline de lissage mesure l'erreur de prévision pour ce nœud. Réitérons cette opération, avec le même paramètre {{mvar|λ}}, pour les {{mvar|n}} nœuds :
:: l'indicateur <math>\sum \frac{\epsilon_i^2}{n}</math>
est une approximation de la variance de la prévision par spline de lissage sous condition du paramètre {{mvar|λ}}. On appelle cet indicateur ''score de validation croisée'' {{math|CV(''λ'')}}.

En réitérant cette procédure pour plusieurs valeurs de {{mvar|λ}}, on peut rechercher la valeur du paramètre de lissage qui minimise ce score de [[validation croisée]]. Il n'est pas garanti que {{math|CV}} n'aie pas des minimums locaux. Les calculs nécessaires peuvent donc être volumineux et la recherche laborieuse. Heureusement, on dispose d'un algorithme qui permet de déterminer {{math|CV(''λ'')}} sans effectuer le calcul des splines cubiques de lissage. 

== Calcul pratique des splines ==
Que ce soit pour l'interpolation ou le lissage, le calcul met en jeu, si les nœuds sont nombreux, de très grosses matrices qui peuvent poser problème lors de la mise en œuvre du calcul. Les matrices utilisées sont très ''creuses'', c'est-à-dire que tous leurs éléments, sauf ceux de trois diagonales principales, sont nuls. Il existe des logiciels mettant en œuvre des méthodes efficaces de stockage de matrices creuses. Néanmoins, les calculs peuvent être longs.

{{Lien|lang=de|trad=Christian Reinsch|fr=Christian Reinsch|texte=Reinsch}} donne un algorithme qui utilise la [[factorisation de Cholesky]] pour réduire la taille des calculs.

== Notes et références ==
{{Références}}

== Voir aussi ==
{{Autres projets
|wiktionary=spline
}}

=== Articles connexes ===
* [[Courbe de Bézier]]
* [[B-spline]]
* [[NURBS]]
* [[Spline cubique d'Hermite]]
*[[Compression de courbe]]

=== Bibliographie (splines, B-splines, Dm-splines...) ===

* Outre les travaux initiaux de '''De Casteljau''' et '''Bézier''', les ouvrages de référence sont ceux de '''Pierre Jean Laurent''' (''Approximation et optimisation, Hermann, Paris, 1972''), '''[[Carl R. de Boor|Carl de Boor]]''' (''A Practical Guide to Splines, Springer Verlag, Berlin-Heidelberg, 1978''). 
*D.G. Schweikert, An interpolation curve using a spline in tension, J. Math. and Phys., vol 45, pp. 312-317, 1966.
*L.L. Schumaker, Fitting surfaces to scattered data, Approximation Théorie II, G.G. Lorentz, C.K. Chui et L.L. Schumaker (éditeurs), Academic Press, 203-269, (1976).
*M. Atteia, Fonctions "spline" définies sur un ensemble convexe, Numer. Math., 12, 192-210, 1968.
*R. Arcangéli, Some Applications of Discrete D^{m} Splines, Mathematical Methods in Computer Aided Geometric Design, T. Lyche et L.L. Schumaker (Editeurs), Academic Press, 35-44, 1989.
*R. Arcangéli, M. C. López-de-Silanes, and J. J. Torrens, Multidimensional minimizing splines. Theory and applications, Grenoble Sciences, pp.1-4020, 2004.
*J. Duchon, Splines Minimizing Rotation- Invariant Semi- Norms in Sobolev Spaces, Lecture Notes in Math., 571, 85-100, Springer, 1977.
*C. Gout, Z. Lambert and D. Apprato, ''Data approximation : mathematical modelling and numerical simulations,''  {{ISBN|978-2-7598-2367-3}}, [[EDP Sciences]], 2019.
*{{Ouvrage|
prénom1=G.|
nom1=Demengel|
prénom2=J. P.|
nom2=Pouget|
titre=Modèles de Bézier, des B-splines et des NURBS|
éditeur=[[Éditions Ellipses|Ellipses]]|
année=1998|
pages totales=286|
isbn=978-2-7298-9806-9
}}
* {{Ouvrage
 | auteur1       = Michelle Schatzman
 | titre         = Analyse numérique :
 | sous-titre    = Approche mathématique
 | lieu          = Paris
 | éditeur       = [[Éditions Dunod|Dunod]]
 | collection    = Sciences Sup
 | année         = 2004
 | mois          = Juillet
 | pages totales = 480
 | isbn          = 978-2-10-048732-5
}}.

=== Liens externes ===
* {{en}} Richard Bartels, John C. Beatty, Brian A. Barsky, [https://cs.uwaterloo.ca/research/tr/1983/CS-83-09.pdf An Introduction to the ''Use of Splines in Computer Graphics'', TR CS-83-09, ''Introduction à l'usage des splines en dessin par informatique''], [[Université de Waterloo]], Canada ;
* {{Article
 | langue                 = en
 | auteur1                = Christian H. Reinsch
 | titre                  = Smoothing by Spline Functions
 | périodique             = [[Numerische Mathematik]]
 | volume                 = 10
 | date                   = 1967
 | pages                  = 177-183
 | lire en ligne          = https://pdfs.semanticscholar.org/3d59/0e6f0fdf32dbd1aece65b9b9fc58953d5890.pdf
 | consulté le            = 14 janvier 2018
}} ;
*{{Lien web
 | format           = PDF
 | auteur1          = J-P Croisille
 | url              = http://www.math.univ-metz.fr/~croisil/M1-0809/2
 | titre            = INTERPOLATION SPLINE
 | série            = master de mathématiques M1, année 2008-2009
 | année            = 2008
 | site             = iecl.univ-lorraine.fr
 | éditeur          = Institut Élie Cartan de Lorraine
 | consulté le      = 14 janvier 2018
}} ;
* {{Lien web
 | format           = PDF
 | auteur1          = Frédéric PLANCHET
 | url              = http://www.ressources-actuarielles.net/C1256F13006585B2/0/1430AD6748CE3AFFC1256F130067B88E/$FILE/Seance6.pdf?OpenElement
 | titre            = Méthodes de lissage et d’ajustement
 | série            = MODÈLES DE DURÉE - Support de cours 2017-2018
 | mois             = Août
 | année            = 2017
 | site             = ressources-actuarielles.net
 | éditeur          = Institut de science financière et d'assurance
 | page             = 4
 | consulté le      = 14 janvier 2018
}}.

{{Palette|Analyse Numérique}}

{{Portail|analyse}}

[[Catégorie:Spline]]
[[Catégorie:Courbe]]
[[Catégorie:Interpolation polynomiale]]