[MatLab] Malcondizionamento della matrice di Vandermonde

[prince931]

Salve sono nuovo del forum volevo chiedervi gentilemte il vostro aiuto in questo esercizio che non sono in grado di fare.



Grazie attendo una vostra risposta.

[Quinzio]

Questo è quanto ottengo con polyfit (blu) e Vandermonde (verde).

[vict85]

Il grafico non dovrebbe uscire così.. Potresti scrivere i comandi MATLAB che hai usato?

[Quinzio]

Si certo...
ma del resto la matrice di Vandermonde è malcondizionata e anche polyfit lavora male perchè ha dei coefficienti molto alti.
Tutta quella dispersione attorno al valore corretto credo che siano dei roundoff.

x=1200.1:0.1:1200.9;
x([4 6 8])=[];
y=[3 1.5 1.5 1 1 0 ];
c_p=polyfit(x,y, 5);
V=vander(x);
c_v=V\y';
xx=linspace(1200,1201,1000);
subplot(2,2,1)
plot(xx,polyval(c_p,xx))
subplot(2,2,2)
plot(xx,polyval(c_v,xx))
subplot(2,2,3)
plot(x,y,'r-*')
subplot(2,2,4)
plot(xx,polyval(c_p,xx),xx,polyval(c_v,xx),x,y,'r*')