Differenze divise: differenze tra le versioni

Versione delle 22:38, 11 apr 2019

Versione attuale delle 14:32, 17 mag 2022

(20 versioni intermedie di 6 utenti non mostrate)

Riga 1:

{{W|matematica|aprile 2019}}

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 ==

Siano <math display="inline">\{x_0, x_1,\dots, x_n\}</math>, <math display="inline">n+1</math> punti assegnati, che inizialmente supponiamo distinti.

Dati <math>k+1</math>punti

:<math>(x_0, y_0),\ldots,(x_{k}, y_{k}).</math>

Definiamo la '''differenza divisa di ordine <math display="inline">0</math> di di <math display="inline">f(x)</math>''':

Definiamo le '''differenze divise''' come:

<math>

f[x_0] := y_0 = f(x_0)

</math>

:<math>[y_\nu] := y_\nu, \qquad \nu \in \{ 0,\ldots,k\},</math>

dove una scrittura equivalente per <math>f[x_0]</math>è <math display="inline">A_0</math>.Definiamo la '''differenza divisa di ordine <math display="inline">1</math>''':

:<math>[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\}.</math>

Definiamo le '''differenze divise all'indietro''' come:

<math>

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} = f[x_1,x_0]

</math>

:<math>[y_\nu] := y_{\nu},\qquad \nu \in \{ 0,\ldots,k\}</math>

che è il [[rapporto incrementale]] <math display="inline">\frac{\Delta y}{\Delta x}</math> costruito su due punti.

:<math>[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\}.</math>

dove <math>j</math> è l''''ordine''' della differenza divisa.

Definiamo la '''differenza divisa di ordine <math display="inline">2</math>''':

== Notazione, differenze divise sui punti di una funzione ==

<math>

Se i punti <math display="inline">\{x_0, x_1,\dots, x_k\}</math> vengono dati come valori di una funzione <math>f</math>:

A_2 = f[x_0,x_1, x_2] = \frac{f[x_1, x_2]-f[x_0,x_1]}{x_2-x_0}

</math>

:<math>(x_0, f(x_0)),\ldots,(x_{k}, f(x_{k})),</math>

E in generale la '''differenza divisa di ordine <math display="inline">n</math>''':

si può trovare la notazione

<math>

A_n = f[x_0,x_1,\dots, x_n] =

:<math>f[x_\nu] := f(x_{\nu}), \qquad \nu \in \{ 0,\ldots,k \},</math>

\frac

:<math>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\}.</math>

{f[x_1,x_2,\dots, x_n]-f[x_0,x_1,\dots, x_{n-1}]}

{x_n-x_0}

Altre scritture equivalenti sono:

</math>

:<math>[x_0,\ldots,x_n]f;</math>

:<math>f_n[x_0,\ldots,x_n];</math>

:<math>[x_0,\ldots,x_n;f];</math>

:<math>D[x_0,\ldots,x_n]f.</math>

== Rapporto con le derivate di ''f''(''x'') ==

Quando due argomenti risultano coincidenti possiamo ugualmente dare un significato alla corrispondente differenza divisa di ordine <math>1</math>, purché <math>f'(x)</math> esista in quel punto<ref name=":0">{{Cita libro|cognome=Monegato, Giovanni.|titolo=Metodi e algoritmi per il calcolo numerico|url=https://www.worldcat.org/oclc/956017867|accesso=2019-04-29|data=[2008]|editore=Clut|OCLC=956017867|ISBN=9788879922654}}</ref>:

:<math>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).</math>

Più in generale, definiamo

:<math>f[\underbrace{x_0,x_0,\dots,x_0}_{k+1}] = \frac{f^{(k)}(x_0)}{k!},</math>

la cui esistenza è dimostrabile<ref>{{Cita libro|cognome=Isaacson, Eugene.|titolo=Analysis of numerical methods|url=https://www.worldcat.org/oclc/30032279|accesso=2019-04-29|data=1994|editore=Dover Publications|p=252|OCLC=30032279|ISBN=0486680290}}</ref>.

== Esempi ==

Differenze divise per <math>\nu=0</math> e i primi valori di <math>j</math>:

:<math>\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}</math><!-- the \mathopen command is there because latex otherwise thinks that [...] denotes an optional argument -->

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

:<math>

\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}

</math>

=== Rapporto incrementale ===

{{vedi anche|Rapporto incrementale}}

Data una funzione <math>f</math>, presi due punti <math>(x_0, f(x_0)),(x_1, f(x_1))</math>, la differenza divisa di ordine '''<math>1</math>''':

:<math>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}</math>

è il [[rapporto incrementale]] costruito su due punti per la quantità <math>h = x_1 - x_0</math>.

== Invarianza per permutazione ==

{{vedi anche|Funzione simmetrica}}

Per induzione matematica, non è difficile dimostrare che

:<math>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)}}}.</math>

Questa espressione ci permette di affermare che <math>f[x_0,x_1,\dots, x_n]</math> è una funzione '''[[Funzione simmetrica|invariante a permutazione]]''' dei suoi argomenti, cioè

:<math>f[x_0,x_1,\dots, x_n] = f[x_{i_0},x_{i_1},\dots, x_{i_n}],</math>

dove <math display="inline">(i_0,i_1,\dots, i_n)</math> denota una qualsiasi [[permutazione]] di <math display="inline">(0, 1, \dots, n)</math><ref name=":0" />.

Per induzione matematica non è difficile dimostrare che

== Note ==

<math>

<references />

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)}}

}

== Bibliografia ==

</math>

*{{Cita libro|cognome=Monegato, Giovanni.|titolo=Metodi e algoritmi per il calcolo numerico|url=https://www.worldcat.org/oclc/956017867|accesso=2019-04-29|data=[2008]|editore=Clut|OCLC=956017867|ISBN=9788879922654}}

== Voci correlate ==

Questa espressione ci permette di affermare che <math display="inline">f[x_0,x_1,\dots, x_n]</math> è una funzione '''invariante a [[permutazione]]''' dei suoi argomenti, cioè

{{Div col}}

* [[Interpolazione polinomiale]]

* [[Interpolazione di Hermite]]{{Div col end}}

{{Portale|matematica}}

<math>

f[x_0,x_1,\dots, x_n] = f[x_{i_0},x_{i_1},\dots, x_{i_n}]

</math>

[[Categoria:Analisi numerica]]

dove <math display="inline">(i_0,i_1,\dots, i_n)</math> denota una qualsiasi permutazione di <math display="inline">(0, 1, \dots, n)</math>.