Régression linéaire Matlab

[seb__d]

Bonjour

Après une experimentation avec relevé de point, j'aimerai effectuer une régression linéaire sur ces point à l'aide de matlab. Tout d'abord j'aimerai une regression normal d'ordre 1 (une doite) et deuxiement j'aimerai que cette droite par 0. En cherchant dans Matlab je n'arrive pas a trouver la fonction qui me permet de faire cette operation. Est ce que vous pourriez m'aider s'il vous plait.

Merci d'avance.
-----

[invite24357dd0]

Bonjour,
Tu peux utiliser la fonction polyfit (help polyfit). Voici un exemple de code (fichier .m ; a partir de donnees vraiment stupides) :

clear;

% Donnees originales
a = [[0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20];
     [2 5 3 5 7 9 11 13 15 17 19 20 22 24 25 26 28 30 33 36 39]];
plot(a(1,:), a(2,:), 'd');

% Regression lineaire de degre 0 : f(x) = m*x
% Coupe donc l'axe en (0, 0)
p = polyfit(a(:,1), a(:,2), 0);
b(1,:) = a(1,:);
b(2,:) = a(1,:) * p;

% Regression lineaire d'ordre 1 : f(x) = m*x + n
p = polyfit(a(:,1), a(:,2), 1);
c(1,:) = a(1,:);
c(2,:) = (p(2) + a(1,:) * p(1));
plot(a(1,:), a(2,:), 'd', b(1,:), b(2,:), c(1,:), c(2,:));
title('Donnees et regressions lineaires');
legend('donnees', 'ordre 0', 'ordre 1');

Cela donne ce graphique-ci :

[invite7f2944e0]

Le graphe de cette régression est faux, quasiment tous les points sont en dessous de la courbe!
Et pour cause, ce n'est pas :
p = polyfit(a(:,1), a(:,2), 1);
mais
p = polyfit(a(1,, a(2,, 1);

On obtient une courbe qui fit beaucoup mieux ainsi...

[lou_ibmix_xi]

Il était temps que tu répondes... Après 4 ans...

[A voir en vidéo sur Futura]