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 .
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.
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.}
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] .
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}}}
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}}
.
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] .
^ a b Monegato, Giovanni., Metodi e algoritmi per il calcolo numerico , Clut, [2008], ISBN 9788879922654 , OCLC 956017867 . URL consultato il 29 aprile 2019 .
^ Isaacson, Eugene., Analysis of numerical methods , Dover Publications, 1994, p. 252, ISBN 0486680290 , OCLC 30032279 . URL consultato il 29 aprile 2019 .