 |
|
 |
|
|
|
situazione con maple
In Maple sono disponibili due librerie per l'interpolazione polinomiale e per il fit di dati, che si caricano rispettivamente con i comandi with(CurveFitting) e with(stats).
Si noti che l'interpolazione polinomiale può essere fatta anche con punti simbolici,
cosa assolutamente impossibile con matlab/octave.
Per esempio:
PolynomialInterpolation([0,1,2,3],[1,0,0,a],z);
Per i fit abbiamo i comandi fit e leastsquare, relativamente facili da usare.
Il problema con Maple sta forse nell'importare i dati: Maple è un po' limitato nei
formati che può importare.
Per esempio si provi a leggere questi dati numerici e a
fittarli meglio possibile con una rette e/o con una parabola.
T:=importdata(`nomefile`,2);
Una volta letti i dati in T, i comandi utilizzabili sono per esempio:
eq_fit1:=fit[leastsquare[[x,y],y=a*x+b,{a,b}]]([T]);
plot(rhs(eq_fit1),x=0..12);
eq_fit2:=fit[leastsquare[[x,y],y=a*x^2+b*x+c,{a,b,c}]]([T]);
plot(rhs(eq_fit2),x=0..12);
|
|
|