Uno dei modi più semplici per interpolare una lista di punti è basato sui polinomi algebrici.
Mettiamo alla prova il metodo trovato da lagrange con cinque punti la cui ordinata è ottenuta
in modo casuale con il comando Random
![[Graphics:../Images/index_gr_190.gif]](../Images/index_gr_190.gif)
![[Graphics:../Images/index_gr_191.gif]](../Images/index_gr_191.gif)
Strisciamo, sara utile in seguito
{{0,1.4323274350876178`},{1,3.5342120026828856`},
{2,6.8216251767071086`},{3,7.672220834210907`},
{4,1.9366138720527097`}}
![[Graphics:../Images/index_gr_194.gif]](../Images/index_gr_194.gif)
![[Graphics:../Images/index_gr_196.gif]](../Images/index_gr_196.gif)
La funzione (polinomio) interpolante di Lagrange si scrive
(vedi anche interp-tremetodi.nb in cui si vedo anche come le
singole funzioni di Lagrange soddisfino la regola di Kronecker ) =
)
![[Graphics:../Images/index_gr_200.gif]](../Images/index_gr_200.gif)
![[Graphics:../Images/index_gr_201.gif]](../Images/index_gr_201.gif)
Disegnamolo ... "a occhio" pare proprio il pol interpolante
![[Graphics:../Images/index_gr_203.gif]](../Images/index_gr_203.gif)
Se eseguiamo anche il test di introdurre la lista delle ascisse x in p1[t],
ricaviamo una lista contenente le ordinate fornite in ingresso !
![[Graphics:../Images/index_gr_205.gif]](../Images/index_gr_205.gif)
Fino all'ultima cifra ...
1.4323274350876178`,3.534212002682885`,
6.8216251767071086`,7.672220834210907`,
1.9366138720527097`}
Esaminiamo alcuni casi molto semplici per cercare di comprendere l'idea di Lagrange;
iniziamo con due punti
![[Graphics:../Images/index_gr_207.gif]](../Images/index_gr_207.gif)
![[Graphics:../Images/index_gr_210.gif]](../Images/index_gr_210.gif)
![[Graphics:../Images/index_gr_212.gif]](../Images/index_gr_212.gif)
poi con tre punti
![[Graphics:../Images/index_gr_214.gif]](../Images/index_gr_214.gif)
![[Graphics:../Images/index_gr_217.gif]](../Images/index_gr_217.gif)
![[Graphics:../Images/index_gr_219.gif]](../Images/index_gr_219.gif)
con quattro
![[Graphics:../Images/index_gr_221.gif]](../Images/index_gr_221.gif)
![[Graphics:../Images/index_gr_224.gif]](../Images/index_gr_224.gif)
![[Graphics:../Images/index_gr_226.gif]](../Images/index_gr_226.gif)
Ora è chiaro ...
Partiamo nuovamente da due punti
![[Graphics:../Images/index_gr_228.gif]](../Images/index_gr_228.gif)
Dalla geometria analitica sappiamo che vale la seguente equazione
.
Con Lagrange
![[Graphics:../Images/index_gr_232.gif]](../Images/index_gr_232.gif)
![[Graphics:../Images/index_gr_234.gif]](../Images/index_gr_234.gif)
![[Graphics:../Images/index_gr_236.gif]](../Images/index_gr_236.gif)
Ripetiamo con i seguenti tre punti
![[Graphics:../Images/index_gr_238.gif]](../Images/index_gr_238.gif)
La curva deve essere una parabola
![[Graphics:../Images/index_gr_241.gif]](../Images/index_gr_241.gif)
![[Graphics:../Images/index_gr_243.gif]](../Images/index_gr_243.gif)
![[Graphics:../Images/index_gr_245.gif]](../Images/index_gr_245.gif)
Proviamo ancora con i seguenti quattro punti
![[Graphics:../Images/index_gr_247.gif]](../Images/index_gr_247.gif)
![[Graphics:../Images/index_gr_250.gif]](../Images/index_gr_250.gif)
![[Graphics:../Images/index_gr_252.gif]](../Images/index_gr_252.gif)
![[Graphics:../Images/index_gr_254.gif]](../Images/index_gr_254.gif)
Tutto sembra procedere senza intoppi ma, per testare ora il grado di applicabilità
dell'idea di Lagrange, consideriamo una funzione apparentemente innocua,
la funzione di Runge
![[Graphics:../Images/index_gr_256.gif]](../Images/index_gr_256.gif)
Tabuliamo la funzione per ottenere una lista di cinque punti
![[Graphics:../Images/index_gr_258.gif]](../Images/index_gr_258.gif)
![[Graphics:../Images/index_gr_260.gif]](../Images/index_gr_260.gif)
![[Graphics:../Images/index_gr_262.gif]](../Images/index_gr_262.gif)
![[Graphics:../Images/index_gr_264.gif]](../Images/index_gr_264.gif)
![[Graphics:../Images/index_gr_266.gif]](../Images/index_gr_266.gif)
![[Graphics:../Images/index_gr_268.gif]](../Images/index_gr_268.gif)
![[Graphics:../Images/index_gr_270.gif]](../Images/index_gr_270.gif)
![[Graphics:../Images/index_gr_272.gif]](../Images/index_gr_272.gif)
È soddisfacente ? Proviamo ora con nove punti
![[Graphics:../Images/index_gr_274.gif]](../Images/index_gr_274.gif)
![[Graphics:../Images/index_gr_276.gif]](../Images/index_gr_276.gif)
![[Graphics:../Images/index_gr_278.gif]](../Images/index_gr_278.gif)
![[Graphics:../Images/index_gr_280.gif]](../Images/index_gr_280.gif)
![[Graphics:../Images/index_gr_282.gif]](../Images/index_gr_282.gif)
![[Graphics:../Images/index_gr_284.gif]](../Images/index_gr_284.gif)
E con 17?
![[Graphics:../Images/index_gr_286.gif]](../Images/index_gr_286.gif)
![[Graphics:../Images/index_gr_288.gif]](../Images/index_gr_288.gif)
![[Graphics:../Images/index_gr_290.gif]](../Images/index_gr_290.gif)
![[Graphics:../Images/index_gr_292.gif]](../Images/index_gr_292.gif)
![[Graphics:../Images/index_gr_294.gif]](../Images/index_gr_294.gif)
![[Graphics:../Images/index_gr_296.gif]](../Images/index_gr_296.gif)
Come si può osservare, all'aumentare dei punti della curva, nonostante l'aumento del grado del polinomio, non si ottiene una rappresentazione fedele. La funzione trovata passa sì attraverso tutti i punti, ma presenta delle oscillazioni che si
discostano in modo significativo
Proviamo ancora con un altra funzione:
![[Graphics:../Images/index_gr_299.gif]](../Images/index_gr_299.gif)
![[Graphics:../Images/index_gr_300.gif]](../Images/index_gr_300.gif)
Cominciamo con solo tre punti
![[Graphics:../Images/index_gr_302.gif]](../Images/index_gr_302.gif)
![[Graphics:../Images/index_gr_305.gif]](../Images/index_gr_305.gif)
Vista la natura della curva era prevedibile che con soli tre punti il metodo di Lagrange fallisse. Ma vediamo come si comporta con cinque punti
![[Graphics:../Images/index_gr_308.gif]](../Images/index_gr_308.gif)
![[Graphics:../Images/index_gr_311.gif]](../Images/index_gr_311.gif)
![[Graphics:../Images/index_gr_314.gif]](../Images/index_gr_314.gif)
![[Graphics:../Images/index_gr_317.gif]](../Images/index_gr_317.gif)
![[Graphics:../Images/index_gr_319.gif]](../Images/index_gr_319.gif)
![[Graphics:../Images/index_gr_322.gif]](../Images/index_gr_322.gif)
Come si può vedere, la rappresentazione alla Lagrange è accettabile solo in un intervallo ridotto.