Differenze divise

Questa voce è da wikificare
Questa voce o sezione sull'argomento matematica non è ancora formattata secondo gli standard.

In matematica, una differenza divisa è una quantità, definita in modo ricorsivo su punti distinti. Vengono utilizzate ad esempio nell'interpolazione polinomiale, nei metodi di interpolazione di Newton alle differenze divise e interpolazione di Hermite.

Definizione

Dati k + 1 {\displaystyle k+1} punti

( x 0 , y 0 ) , , ( x k , y k ) . {\displaystyle (x_{0},y_{0}),\ldots ,(x_{k},y_{k}).}

Definiamo le differenze divise come:

[ y ν ] := y ν , ν { 0 , , k } , {\displaystyle [y_{\nu }]:=y_{\nu },\qquad \nu \in \{0,\ldots ,k\},}
[ y ν , , y ν + j ] := [ y ν + 1 , , y ν + j ] [ y ν , , y ν + j 1 ] x ν + j x ν , ν { 0 , , k j } ,   j { 1 , , k } . {\displaystyle [y_{\nu },\ldots ,y_{\nu +j}]:={\frac {[y_{\nu +1},\ldots ,y_{\nu +j}]-[y_{\nu },\ldots ,y_{\nu +j-1}]}{x_{\nu +j}-x_{\nu }}},\qquad \nu \in \{0,\ldots ,k-j\},\ j\in \{1,\ldots ,k\}.}

Definiamo le differenze divise all'indietro come:

[ y ν ] := y ν , ν { 0 , , k } {\displaystyle [y_{\nu }]:=y_{\nu },\qquad \nu \in \{0,\ldots ,k\}}
[ y ν , , y ν j ] := [ y ν , , y ν j + 1 ] [ y ν 1 , , y ν j ] x ν x ν j , ν { j , , k } ,   j { 1 , , k } . {\displaystyle [y_{\nu },\ldots ,y_{\nu -j}]:={\frac {[y_{\nu },\ldots ,y_{\nu -j+1}]-[y_{\nu -1},\ldots ,y_{\nu -j}]}{x_{\nu }-x_{\nu -j}}},\qquad \nu \in \{j,\ldots ,k\},\ j\in \{1,\ldots ,k\}.}

dove j {\displaystyle j} è l'ordine della differenza divisa.

Notazione, differenze divise sui punti di una funzione

Se i punti { x 0 , x 1 , , x k } {\textstyle \{x_{0},x_{1},\dots ,x_{k}\}} vengono dati come valori di una funzione f {\displaystyle f} :

( x 0 , f ( x 0 ) ) , , ( x k , f ( x k ) ) , {\displaystyle (x_{0},f(x_{0})),\ldots ,(x_{k},f(x_{k})),}

si può trovare la notazione

f [ x ν ] := f ( x ν ) , ν { 0 , , k } , {\displaystyle f[x_{\nu }]:=f(x_{\nu }),\qquad \nu \in \{0,\ldots ,k\},}
f [ x ν , , x ν + j ] := f [ x ν + 1 , , x ν + j ] f [ x ν , , x ν + j 1 ] x ν + j x ν , ν { 0 , , k j } ,   j { 1 , , k } . {\displaystyle f[x_{\nu },\ldots ,x_{\nu +j}]:={\frac {f[x_{\nu +1},\ldots ,x_{\nu +j}]-f[x_{\nu },\ldots ,x_{\nu +j-1}]}{x_{\nu +j}-x_{\nu }}},\qquad \nu \in \{0,\ldots ,k-j\},\ j\in \{1,\ldots ,k\}.}

Altre scritture equivalenti sono:

[ x 0 , , x n ] f ; {\displaystyle [x_{0},\ldots ,x_{n}]f;}
f n [ x 0 , , x n ] ; {\displaystyle f_{n}[x_{0},\ldots ,x_{n}];}
[ x 0 , , x n ; f ] ; {\displaystyle [x_{0},\ldots ,x_{n};f];}
D [ x 0 , , x n ] f . {\displaystyle D[x_{0},\ldots ,x_{n}]f.}

Rapporto con le derivate di f(x)

Quando due argomenti risultano coincidenti possiamo ugualmente dare un significato alla corrispondente differenza divisa di ordine 1 {\displaystyle 1} , purché f ( x ) {\displaystyle f'(x)} esista in quel punto[1]:

f [ x 0 , x 0 ] = lim x x 0 f [ x 0 , x ] = lim x x 0 f ( x ) f ( x 0 ) x x 0 = f ( x 0 ) . {\displaystyle f[x_{0},x_{0}]=\lim _{x\to x_{0}}f[x_{0},x]=\lim _{x\to x_{0}}{\frac {f(x)-f(x_{0})}{x-x_{0}}}=f'(x_{0}).}

Più in generale, definiamo

f [ x 0 , x 0 , , x 0 k + 1 ] = f ( k ) ( x 0 ) k ! , {\displaystyle f[\underbrace {x_{0},x_{0},\dots ,x_{0}} _{k+1}]={\frac {f^{(k)}(x_{0})}{k!}},}

la cui esistenza è dimostrabile[2].

Esempi

Differenze divise per ν = 0 {\displaystyle \nu =0} e i primi valori di j {\displaystyle j} :

[ y 0 ] = y 0 ; [ y 0 , y 1 ] = y 1 y 0 x 1 x 0 ; [ y 0 , y 1 , y 2 ] = [ y 1 , y 2 ] [ y 0 , y 1 ] x 2 x 0 = y 2 y 1 x 2 x 1 y 1 y 0 x 1 x 0 x 2 x 0 = y 2 y 1 ( x 2 x 1 ) ( x 2 x 0 ) y 1 y 0 ( x 1 x 0 ) ( x 2 x 0 ) ; [ y 0 , y 1 , y 2 , y 3 ] = [ y 1 , y 2 , y 3 ] [ y 0 , y 1 , y 2 ] x 3 x 0 . {\displaystyle {\begin{aligned}{\mathopen {[}}y_{0}]&=y_{0};\\{\mathopen {[}}y_{0},y_{1}]&={\frac {y_{1}-y_{0}}{x_{1}-x_{0}}};\\{\mathopen {[}}y_{0},y_{1},y_{2}]&={\frac {[y_{1},y_{2}]-[y_{0},y_{1}]}{x_{2}-x_{0}}}={\frac {{\frac {y_{2}-y_{1}}{x_{2}-x_{1}}}-{\frac {y_{1}-y_{0}}{x_{1}-x_{0}}}}{x_{2}-x_{0}}}={\frac {y_{2}-y_{1}}{(x_{2}-x_{1})(x_{2}-x_{0})}}-{\frac {y_{1}-y_{0}}{(x_{1}-x_{0})(x_{2}-x_{0})}};\\{\mathopen {[}}y_{0},y_{1},y_{2},y_{3}]&={\frac {[y_{1},y_{2},y_{3}]-[y_{0},y_{1},y_{2}]}{x_{3}-x_{0}}}.\end{aligned}}}

Per evidenziare il processo ricorsivo, le differenze divise possono essere messe in forma tabellare

x 0 y 0 = [ y 0 ] [ y 0 , y 1 ] x 1 y 1 = [ y 1 ] [ y 0 , y 1 , y 2 ] [ y 1 , y 2 ] [ y 0 , y 1 , y 2 , y 3 ] x 2 y 2 = [ y 2 ] [ y 1 , y 2 , y 3 ] [ y 2 , y 3 ] x 3 y 3 = [ y 3 ] {\displaystyle {\begin{matrix}x_{0}&y_{0}=[y_{0}]&&&\\&&[y_{0},y_{1}]&&\\x_{1}&y_{1}=[y_{1}]&&[y_{0},y_{1},y_{2}]&\\&&[y_{1},y_{2}]&&[y_{0},y_{1},y_{2},y_{3}]\\x_{2}&y_{2}=[y_{2}]&&[y_{1},y_{2},y_{3}]&\\&&[y_{2},y_{3}]&&\\x_{3}&y_{3}=[y_{3}]&&&\\\end{matrix}}}

Rapporto incrementale

Lo stesso argomento in dettaglio: Rapporto incrementale.

Data una funzione f {\displaystyle f} , presi due punti ( x 0 , f ( x 0 ) ) , ( x 1 , f ( x 1 ) ) {\displaystyle (x_{0},f(x_{0})),(x_{1},f(x_{1}))} , la differenza divisa di ordine 1 {\displaystyle 1} :

A 1 = f [ x 0 , x 1 ] = y 1 y 0 x 1 x 0 = f ( x 1 ) f ( x 0 ) x 1 x 0 = Δ f Δ x {\displaystyle A_{1}=f[x_{0},x_{1}]={\frac {y_{1}-y_{0}}{x_{1}-x_{0}}}={\frac {f(x_{1})-f(x_{0})}{x_{1}-x_{0}}}={\frac {\Delta f}{\Delta x}}}

è il rapporto incrementale costruito su due punti per la quantità h = x 1 x 0 {\displaystyle h=x_{1}-x_{0}} .

Invarianza per permutazione

Lo stesso argomento in dettaglio: Funzione simmetrica.

Per induzione matematica, non è difficile dimostrare che

f [ x 0 , x 1 , , x n ] = k = 0 n f ( x k ) j = 0 ,   j k n ( x k x j ) . {\displaystyle f[x_{0},x_{1},\dots ,x_{n}]=\sum _{k=0}^{n}{\frac {f(x_{k})}{\prod _{j=0,\ j\neq k}^{n}{(x_{k}-x_{j})}}}.}

Questa espressione ci permette di affermare che f [ x 0 , x 1 , , x n ] {\displaystyle f[x_{0},x_{1},\dots ,x_{n}]} è una funzione invariante a permutazione dei suoi argomenti, cioè

f [ x 0 , x 1 , , x n ] = f [ x i 0 , x i 1 , , x i n ] , {\displaystyle f[x_{0},x_{1},\dots ,x_{n}]=f[x_{i_{0}},x_{i_{1}},\dots ,x_{i_{n}}],}

dove ( i 0 , i 1 , , i n ) {\textstyle (i_{0},i_{1},\dots ,i_{n})} denota una qualsiasi permutazione di ( 0 , 1 , , n ) {\textstyle (0,1,\dots ,n)} [1].

Note

  1. ^ a b Monegato, Giovanni., Metodi e algoritmi per il calcolo numerico, Clut, [2008], ISBN 9788879922654, OCLC 956017867. URL consultato il 29 aprile 2019.
  2. ^ Isaacson, Eugene., Analysis of numerical methods, Dover Publications, 1994, p. 252, ISBN 0486680290, OCLC 30032279. URL consultato il 29 aprile 2019.

Bibliografia

  • Monegato, Giovanni., Metodi e algoritmi per il calcolo numerico, Clut, [2008], ISBN 9788879922654, OCLC 956017867. URL consultato il 29 aprile 2019.

Voci correlate

  • Interpolazione polinomiale
  • Interpolazione di Hermite
  Portale Matematica: accedi alle voci di Wikipedia che trattano di matematica