Spline cubica di Hermite

In analisi numerica la spline cubica di Hermite (chiamata anche cspline), in onore del matematico Charles Hermite, è una funzione spline di 3º grado dove ogni polinomio della spline è nella forma di Hermite (da non confondere con i polinomi di Hermite). La forma di Hermite consiste di due punti di controllo e di due tangenti di controllo per ogni polinomio.

Su una griglia composta dai punti ${\displaystyle x_{k}}$ per ${\displaystyle k=1,...,n}$, l'interpolazione è effettuata su ogni sottointervallo ${\displaystyle (x_{k},x_{k+1})}$ alla volta (dato che i valori della tangente sono predeterminati). Il sottointervallo ${\displaystyle (x_{k},x_{k+1})}$ è normalizzato all'intervallo ${\displaystyle (0,1)}$ tramite la funzione ${\displaystyle t=(x-x_{k})/(x_{k+1}-x_{k})}$.

Con l'intervallo ${\displaystyle (0,1)}$, dato il punto iniziale p₀ con ${\displaystyle t=0}$ ed il punto finale p₁ con ${\displaystyle t=1}$ con la tangente iniziale m₀ con ${\displaystyle t=0}$ e tangente finale m₁ con ${\displaystyle t=1}$, il polinomio è definito da

${\displaystyle {\boldsymbol {p}}(t)=(2t^{3}-3t^{2}+1){\boldsymbol {p}}_{0}+(t^{3}-2t^{2}+t){\boldsymbol {m}}_{0}+(-2t^{3}+3t^{2}){\boldsymbol {p}}_{1}+(t^{3}-t^{2}){\boldsymbol {m}}_{1}}$

dove t ∈ [0, 1].

Le 4 funzioni base di Hermite sono definite come:

${\displaystyle h_{00}(t)=2t^{3}-3t^{2}+1=(1+2t)(1-t)^{2}}$

${\displaystyle h_{10}(t)=t^{3}-2t^{2}+t=t(1-t)^{2}}$

${\displaystyle h_{01}(t)=-2t^{3}+3t^{2}=t^{2}(3-2t)}$

${\displaystyle h_{11}(t)=t^{3}-t^{2}=t^{2}(t-1)}$

Il polinomio diventa: ${\displaystyle {\boldsymbol {p}}(t)=h_{00}(t){\boldsymbol {p}}_{0}+h_{10}(t){\boldsymbol {m}}_{0}+h_{01}(t){\boldsymbol {p}}_{1}+h_{11}(t){\boldsymbol {m}}_{1}.}$

L'interpolazione di ${\displaystyle x}$ nell'intervallo ${\displaystyle (x_{k},x_{k+1})}$ viene ora effettuata con la formula

${\displaystyle {\boldsymbol {p}}(x)=h_{00}(t){\boldsymbol {p}}_{0}+h_{10}(t)h{\boldsymbol {m}}_{0}+h_{01}(t){\boldsymbol {p}}_{1}+h_{11}(t)h{\boldsymbol {m}}_{1}.}$

con ${\displaystyle h=x_{k+1}-x_{k}}$ e ${\displaystyle t=(x-x_{k})/h}$. Si noti che i valori della tangente sono stati scalati di ${\displaystyle h}$ in raffronto alla equazione presentata sull'intervallo unitario.

La formula garantisce un unico percorso tra i due punti iniziale e finale.

Un set di dati, ${\displaystyle (x_{k},{\boldsymbol {p}}_{k})}$ con ${\displaystyle k=1,...,n}$, viene interpolato applicando la procedura sopra esposta su ogni sottointervallo, dove le tangenti sono scelte in modo opportuno. Le tangenti per gli intervalli che condividono gli stessi punti finali sono uguali.

La scelta delle tangente non è univoca e ci sono diversi metodi che si possono applicare.

La scelta più semplice è data dalla differenza su 3 punti e non richiede una lunghezza di intervallo costante,

${\displaystyle {\boldsymbol {m}}_{k}={\frac {{\boldsymbol {p}}_{k+1}-{\boldsymbol {p}}_{k}}{2(x_{k+1}-x_{k})}}+{\frac {{\boldsymbol {p}}_{k}-{\boldsymbol {p}}_{k-1}}{2(x_{k}-x_{k-1})}}}$

per i punti interni ${\displaystyle k=2,...,n-1}$.

Una cardinal spline si ottiene^[1] se

${\displaystyle {\boldsymbol {m}}_{k}=(1-c){\frac {{\boldsymbol {p}}_{k+1}-{\boldsymbol {p}}_{k-1}}{2}}}$

usata per calcolare le tangenti. Il parametro ${\displaystyle c}$ detto tensione è nell'intervallo ${\displaystyle (0,1)}$. Rappresenta la "lunghezza" della tangente. ${\displaystyle c=1}$ comporta tangenti lunghe pari a zero, e ${\displaystyle c=0}$ porta al caso delle spline di tipo Catmull-Rom.

Per le tangenti si applica

${\displaystyle {\boldsymbol {m}}_{k}={\frac {{\boldsymbol {p}}_{k+1}-{\boldsymbol {p}}_{k-1}}{2}}}$

Una Catmull–Rom spline si ottiene come caso speciale delle cardinal spline.

La curva prende il nome da Edwin Catmull e Raphael (Raphie) Rom. In grafica computerizzata, le Catmull–Rom splines sono usate per ottenere un effetto morbido di interpolazione tra i keyframes di un video.

Una spline di Kochanek-Bartels è una generalizzazione su un modo di calcolare le tangenti partendo dai punti ${\displaystyle {\boldsymbol {p}}_{k-1}}$, ${\displaystyle {\boldsymbol {p}}_{k}}$ e ${\displaystyle {\boldsymbol {p}}_{k+1}}$, dove sono possibili 3 parametri detti tensione, bias e continuità.

• Catmull, Edwin and Rom, Raphael, A class of local interpolating splines, in R.E. Barnhill and R.F. Riesenfed (eds.) Computer Aided Geometric Design, Academic Press, New York, 1974, 317-326.

• Introduction to Catmull-Rom Splines, MVPs.org
• Catmull-Rom spline - Indopedia, the Indological Encyclopedia, su indopedia.org.
• Interpolating Cardinal and Catmull-Rom splines, su ibiblio.org.
• Interpolation methods: linear, cosine, cubic and hermite (with C sources), su local.wasp.uwa.edu.au. URL consultato il 23 agosto 2009 (archiviato dall'url originale il 15 settembre 2008).
• Rapahel Rom's Homepage, su webee.technion.ac.il.

Portale Matematica: accedi alle voci di Wikipedia che trattano di Matematica