Interpolazione per mezzo di polinomi

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]
[Graphics:../Images/index_gr_191.gif]
[Graphics:../Images/index_gr_192.gif]

[Graphics:../Images/index_gr_193.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]
[Graphics:../Images/index_gr_195.gif]
[Graphics:../Images/index_gr_196.gif]
[Graphics:../Images/index_gr_197.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_198.gif]) = [Graphics:../Images/index_gr_199.gif])

[Graphics:../Images/index_gr_200.gif]
[Graphics:../Images/index_gr_201.gif]
[Graphics:../Images/index_gr_202.gif]

Disegnamolo ... "a occhio"  pare proprio il pol interpolante

[Graphics:../Images/index_gr_203.gif]

[Graphics:../Images/index_gr_204.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]
[Graphics:../Images/index_gr_206.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]
[Graphics:../Images/index_gr_208.gif]
[Graphics:../Images/index_gr_209.gif]
[Graphics:../Images/index_gr_210.gif]
[Graphics:../Images/index_gr_211.gif]
[Graphics:../Images/index_gr_212.gif]
[Graphics:../Images/index_gr_213.gif]

poi con tre punti

[Graphics:../Images/index_gr_214.gif]
[Graphics:../Images/index_gr_215.gif]
[Graphics:../Images/index_gr_216.gif]
[Graphics:../Images/index_gr_217.gif]
[Graphics:../Images/index_gr_218.gif]
[Graphics:../Images/index_gr_219.gif]
[Graphics:../Images/index_gr_220.gif]

con quattro

[Graphics:../Images/index_gr_221.gif]
[Graphics:../Images/index_gr_222.gif]
[Graphics:../Images/index_gr_223.gif]
[Graphics:../Images/index_gr_224.gif]
[Graphics:../Images/index_gr_225.gif]
[Graphics:../Images/index_gr_226.gif]
[Graphics:../Images/index_gr_227.gif]

Ora è chiaro ...

Partiamo nuovamente da due punti

[Graphics:../Images/index_gr_228.gif]
[Graphics:../Images/index_gr_229.gif]

[Graphics:../Images/index_gr_230.gif]

Dalla geometria analitica sappiamo che vale la seguente equazione

                             [Graphics:../Images/index_gr_231.gif].

Con  Lagrange

[Graphics:../Images/index_gr_232.gif]
[Graphics:../Images/index_gr_233.gif]
[Graphics:../Images/index_gr_234.gif]
[Graphics:../Images/index_gr_235.gif]
[Graphics:../Images/index_gr_236.gif]

[Graphics:../Images/index_gr_237.gif]

Ripetiamo con i seguenti tre punti

[Graphics:../Images/index_gr_238.gif]
[Graphics:../Images/index_gr_239.gif]

[Graphics:../Images/index_gr_240.gif]

La curva deve essere una parabola

[Graphics:../Images/index_gr_241.gif]
[Graphics:../Images/index_gr_242.gif]
[Graphics:../Images/index_gr_243.gif]
[Graphics:../Images/index_gr_244.gif]
[Graphics:../Images/index_gr_245.gif]

[Graphics:../Images/index_gr_246.gif]

Proviamo ancora con i seguenti quattro punti

[Graphics:../Images/index_gr_247.gif]
[Graphics:../Images/index_gr_248.gif]

[Graphics:../Images/index_gr_249.gif]

[Graphics:../Images/index_gr_250.gif]
[Graphics:../Images/index_gr_251.gif]
[Graphics:../Images/index_gr_252.gif]
[Graphics:../Images/index_gr_253.gif]
[Graphics:../Images/index_gr_254.gif]

[Graphics:../Images/index_gr_255.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]

[Graphics:../Images/index_gr_257.gif]

Tabuliamo la funzione per ottenere una lista di cinque punti

[Graphics:../Images/index_gr_258.gif]

[Graphics:../Images/index_gr_259.gif]

[Graphics:../Images/index_gr_260.gif]
[Graphics:../Images/index_gr_261.gif]
[Graphics:../Images/index_gr_262.gif]
[Graphics:../Images/index_gr_263.gif]
[Graphics:../Images/index_gr_264.gif]
[Graphics:../Images/index_gr_265.gif]
[Graphics:../Images/index_gr_266.gif]
[Graphics:../Images/index_gr_267.gif]
[Graphics:../Images/index_gr_268.gif]
[Graphics:../Images/index_gr_269.gif]
[Graphics:../Images/index_gr_270.gif]
[Graphics:../Images/index_gr_271.gif]
[Graphics:../Images/index_gr_272.gif]

[Graphics:../Images/index_gr_273.gif]

È soddisfacente ? Proviamo ora con nove punti

[Graphics:../Images/index_gr_274.gif]

[Graphics:../Images/index_gr_275.gif]

[Graphics:../Images/index_gr_276.gif]
[Graphics:../Images/index_gr_277.gif]
[Graphics:../Images/index_gr_278.gif]
[Graphics:../Images/index_gr_279.gif]
[Graphics:../Images/index_gr_280.gif]
[Graphics:../Images/index_gr_281.gif]
[Graphics:../Images/index_gr_282.gif]
[Graphics:../Images/index_gr_283.gif]
[Graphics:../Images/index_gr_284.gif]

[Graphics:../Images/index_gr_285.gif]

E con 17?

[Graphics:../Images/index_gr_286.gif]

[Graphics:../Images/index_gr_287.gif]

[Graphics:../Images/index_gr_288.gif]
[Graphics:../Images/index_gr_289.gif]
[Graphics:../Images/index_gr_290.gif]
[Graphics:../Images/index_gr_291.gif]
[Graphics:../Images/index_gr_292.gif]
[Graphics:../Images/index_gr_293.gif]
[Graphics:../Images/index_gr_294.gif]
[Graphics:../Images/index_gr_295.gif]
[Graphics:../Images/index_gr_296.gif]

[Graphics:../Images/index_gr_297.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_298.gif]

[Graphics:../Images/index_gr_299.gif]
[Graphics:../Images/index_gr_300.gif]

[Graphics:../Images/index_gr_301.gif]

Cominciamo con solo tre punti

[Graphics:../Images/index_gr_302.gif]
[Graphics:../Images/index_gr_303.gif]
[Graphics:../Images/index_gr_304.gif]
[Graphics:../Images/index_gr_305.gif]
[Graphics:../Images/index_gr_306.gif]

[Graphics:../Images/index_gr_307.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]
[Graphics:../Images/index_gr_309.gif]
[Graphics:../Images/index_gr_310.gif]
[Graphics:../Images/index_gr_311.gif]
[Graphics:../Images/index_gr_312.gif]

[Graphics:../Images/index_gr_313.gif]

[Graphics:../Images/index_gr_314.gif]
[Graphics:../Images/index_gr_315.gif]
[Graphics:../Images/index_gr_316.gif]
[Graphics:../Images/index_gr_317.gif]

[Graphics:../Images/index_gr_318.gif]

[Graphics:../Images/index_gr_319.gif]
[Graphics:../Images/index_gr_320.gif]
[Graphics:../Images/index_gr_321.gif]
[Graphics:../Images/index_gr_322.gif]

[Graphics:../Images/index_gr_323.gif]

Come si può vedere, la rappresentazione alla Lagrange è accettabile solo in un intervallo ridotto.


Converted by Mathematica      May 26, 2003