ElectroYou

Buona sera, nell'ambito dell'interpolazione polinomiale, come si calcola l' errore di discretizzazione di Newton?

Grazie

Nel caso generale, leggi qui:
viewtopic.php?f=7&t=52796

Per l'interpolazione con i polinomi di Newton, ne esistono di due tipi: interpolazione alle differenze divise ed alle differenze finite.

Quale ti interessa? Entrambi?

sono interessato alle differenze divise,grazie

Dunque, le differenze divise sono definite come:





e così via all'aumentare dei nodi.

Ora, applicando ripetutamente la definizione (ogni volta dal passo n+1, ottieni ciò che ti serve per il passo n) otteniamo (prova a farlo per i primi 2 nodi):




L'ultimo addendo non è calcolabile numericamente in quanto non si conosce il valore del generico nodo .

Trascurando nei calcoli quel termine otteniamo che l'errore di discretizzazione o troncamento è dato proprio da:



Confrontando questa espressione con quella del caso generale nel link in [2], otteniamo l'uguaglianza:



Supponendo infine che la derivata n-esima sia poco variabile nell'intervallo di interpolazione possiamo approssimare e quindi stimare l'errore di troncamento:



Dall'ultimo passaggio sappiamo quindi che più la derivata è variabile, più la stima sarà peggiore.
Spero di essere stato chiaro, altrimenti dimmi pure

Ti ringrazio infinitamente, sei stato chiarissimo.

Volevo chiederti solo un'altra cosa, perché al denominatore appare il fattoriale di n (dove n= numero di nodi) ?

thx

Ti ringrazio infinitamente, sei stato chiarissimo

Prego, mi fa piacere

perché al denominatore appare il fattoriale di n?

Perché non ti sei fatto la stessa domanda del polinomio nodale e della derivata n-esima?
Dimostriamolo va

Allora..
Dobbiamo stimare l'errore di troncamento nel caso più generale possibile quando si approssima una funzione con un polinomio di grado .
Sappiamo per certo che (a meno di errore di propagazione sui valori di ingresso, che trascuriamo) l'errore sarà nullo nei nodi, allora possiamo scrivere che in generale l'errore di discretizzazione sarà della forma:



Ora dobbiamo determinare .

Per farlo consideriamo una funzione ausiliaria definita così:



(nota che è proprio la definizione di errore di troncamento).
Per come è definita questa funzione ammette zeri, ed essendo polinomiale possiamo dire che la derivata n-esima si annullerà in un punto che chiamiamo .
Allora abbiamo che:



da cui finalmente:



Conoscendo siamo giunti allora alla conclusione che: