Outils pour utilisateurs

Outils du site


teaching:progappchim:regression_lineaire_2013

Différences

Ci-dessous, les différences entre deux révisions de la page.

Lien vers cette vue comparative

teaching:progappchim:regression_lineaire_2013 [2013/11/28 13:38] (Version actuelle)
villersd créée
Ligne 1: Ligne 1:
 +====== Régression linéaire ======
 +Entrée de couples, calcul et affichage de la [[http://​fr.wikipedia.org/​wiki/​M%C3%A9thode_des_moindres_carr%C3%A9s#​R.C3.A9gression_lin.C3.A9aire|droite de moindres carrés]]
  
 +<sxh python; title : fit_linear.py>​
 +#​!/​usr/​bin/​env python
 +# -*- coding: utf-8 -*-
 +# version un peu aménagée du travail de BD et EH, ba2 chimie 2012-2013
 +
 +import matplotlib.pyplot as plt
 +import pylab
 +import numpy
 +
 +def droite_moindre_carre(x):​
 +    return a * x + b
 +
 +tabx = []
 +taby = []
 +# Avoir le jeu de data
 +
 +while True:
 +    x = raw_input("​Enter X ")
 +    y = raw_input("​Enter Y ")
 +    tabx.append(float(x))
 +    taby.append(float(y))
 +    another = raw_input("​Continue ? y/n \n")
 +    if another != "​y":​
 +        break
 +# le remplir tant que l'​utilisateur n'a pas mis autre chose que y comme réponse
 +
 +temp_min = 0.1
 +temp_max = 20
 +#axes horizontaux
 +
 +a,b = pylab.polyfit(tabx,​ taby, 1)
 +# calcul du moindre carre
 +
 +x_moindre_carre = numpy.linspace(min(tabx),​ max(tabx), 100)
 +#​linspace(a,​b,​n) => fournit une liste de valeurs de a à b, en n étapes
 +y_moindre_carre = map(droite_moindre_carre,​ x_moindre_carre)
 +# calcul des valeurs de la courbe de moidnres carrés
 +# map => Utilise une fonction sur chaque element d'une liste
 +
 +plt.axis([temp_min,​ temp_max,​0,​100])
 +plt.plot(tabx,​ taby, '​or'​)
 +# o => point, r => red
 +plt.plot(x_moindre_carre,​ y_moindre_carre,​ '​b'​)
 +# b = blue, pas de o => line
 +plt.show()
 +# trace le graphique
 +</​sxh>​
teaching/progappchim/regression_lineaire_2013.txt · Dernière modification: 2013/11/28 13:38 par villersd