18.01.2015 Views

Appunti di Calcolo Numerico - Esercizi e Dispense - Università degli ...

Appunti di Calcolo Numerico - Esercizi e Dispense - Università degli ...

Appunti di Calcolo Numerico - Esercizi e Dispense - Università degli ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

12.9. Applicazioni <strong>di</strong> MATLAB® nel <strong>Calcolo</strong> <strong>Numerico</strong><br />

centering<br />

and scaling as described in HELP POLYFIT.<br />

Questo ci <strong>di</strong>ce che i risultati che avremo non saranno buoni e, effettivamente, se facciamo il grafico dei dati<br />

del problema e del polinomio <strong>di</strong> interpolazione ottenuto con la polyfit, si nota subito che il polinomio è<br />

completamente errato (ve<strong>di</strong> Figura 12.8) Anche nell’applicare la function interpmonom si ha il messaggio<br />

Figura 12.6: Uso della function polyfit per l’esempio malcon<strong>di</strong>zionato.<br />

<strong>di</strong> avvertimento<br />

>> pmon=interpmonom(x,y)<br />

Warning: Matrix is close to singular or badly scaled.<br />

Results may be inaccurate. RCOND = 5.537667e-31.<br />

Se aggiungiamo al grafico precendente la curva corrispondente al polinomio ottenuto con la interpmonom<br />

osserviamo come le due curve siano tra loro <strong>di</strong>stinte e non interpolano i dati.<br />

Figura 12.7: Uso della function interpmonom per l’esempio malcon<strong>di</strong>zionato.<br />

Proviamo invece ad applicare l’algoritmo <strong>di</strong> Lagrange o delle <strong>di</strong>fferenze <strong>di</strong>vise <strong>di</strong> Newton e aggiungiamo<br />

le nuove curve sul precedente grafico:<br />

>> ylagr=interplagrange(xx,x,y);<br />

>> table=<strong>di</strong>v<strong>di</strong>f(x,y);<br />

>> ynewt=interp<strong>di</strong>v<strong>di</strong>f(xx,x,table);<br />

>> plot(xx,ylagr,xx,ynewt)<br />

Si nota subito che le due curve sono tra loro coincidenti e interpolatorie!<br />

201

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!