Bonjour à tous,
Je dois écrire un code matlab qui permet de faire un fitting entre une courbe expérimentale et autres analytiques .( courbes non linéaire)
j'aimerais trouver la courbe la plus proche de celle de l’expérience (avec R la plus proche de 1), est ce que vous pouvez me renseigner sur la méthode de moindre carrée SVP.
Je travaille sur matlab version R2011b .
Merci d'avance
Salut
Est-ce tu as déjà une idée de la forme générale qu'a ta courbe ?
C'est une parabole, une exponentielle, une gaussienne ...?
salut,
c'est une parabole.
voila code :
clear all
close all
clc
x=1:100;
%z=16.18*x.^2+89;
w=[10 21 28 37 43 55 68 78 82 94 90 100 ];%valeurs exprimentales
a= [2000 8000 1.2e4 2.2e4 2.6e4 5.2e4 6.8e4 8.7e4 8.8e4 14e4 11.9e4 15e4];%valeurs exprimentales
figure,hold on;
for i=10:1:20;
y=i*x.^2+5;
plot(x,y,w,a,'*')
end
T'as plusieurs façons :
t'as la fonction polyfit sur Matlab qui fait tout le travail pour toi. T'as juste à faire polyfit(w,a,2) et ça te donne les coefficient m,n,p de ta parabole mw^2+nw+p
Sinon tu passes par l'algèbre linéaire, matlab est fait pour ça (et c'est d'ailleurs comme ça que polyfit fonctionne) :
Tu crées une matrice W avec en première ligne les w au carré, en 2e ligne les w et en 3e lignes des 1. Si M=[m,n,p] est la matrice que tu cherches, alors:
MW = a (tu peux voir que ça marche au niveau des dimensions des matrices: (1x3)*(3x12)=(1*12) )
donc M = a*pinv(W)
où pinv est la fonction matlab qui implémente la pseudo-inverse de W car W n'est pas carrée.
bonjour ,
merci pour votre réponse j'ai essayée avec polyfit F=polyfit(w,a,2) et ca donne
F =
1.0e+003 *
0.0154 -0.0733 1.8857
est ce que vous pouvez m'expliquer ces chiffres là !!!!
Je l'ai déjà fait ci-dessus. Et en tapant "help polyfit" dans matlab, tu as tout le détail de la fonction
