Différences

no way to compare when less than two revisions

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

@@ 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>