Fonction de tests

src/main.c

@@ -38,13 +38,26 @@ int main() {
 
   point p = findBestDroite(v);
   printf("Sol Analytique : a: %f b: %f\n",p.x,p.y);
+
+
   point drt;
   drt.x = 6;
   drt.y = 6;
-  p = descenteGradient(0.0001,drt,v);
+  p = descenteGradient(0.000001,drt,v);
   printf("Sol Numérique : a: %f b: %f\n",p.x,p.y);
 
+  //Tests
+  point drtTest;
+  drtTest.x = 0;
+  drtTest.y = 2;
+  vector vecTest = nuagePoints(drtTest, 30, 0, 20);
+  p = descenteGradient(0.000001,drtTest,v);
+  //Très proche de la sol numérique donc peut-être augmenter le bruit
+  printf("Sol Numérique : a: %f b: %f\n",p.x,p.y);
+
+
   vector_free(&v);
+  vector_free(&vecTest);
 
   return EXIT_SUCCESS;
 }

src/regressionLin.c

@@ -37,7 +37,7 @@ point descenteGradient(double precision, point drt, vector pts){
   double oldptx = 0;
   double oldpty = 0;
   double testSortie=1;
-  double lambda = 0.0001;
+  double lambda = 0.000001;
   while(testSortie > precision){
     gradientx = 0;
     gradienty = 0;
@@ -56,15 +56,28 @@ point descenteGradient(double precision, point drt, vector pts){
     oldpty = pt.y;
 
     pt.x = oldptx - (lambda * gradientx);
-//    printf("ptx : %f\n",pt.x);
     pt.y = oldpty - (lambda * gradienty);
-  //  printf("pty : %f\n",pt.y);
 
 
     //Sortie
     testSortie = sqrt(pow(pt.x - oldptx,2) + pow(pt.y - oldpty,2));
-//    printf("%f\n",testSortie);
   }
 
   return pt;
 }
+
+vector nuagePoints(point drt, int len, int minP, int maxP){
+  vector vec;
+  vector_init(&vec);
+  srand( time( NULL ) );
+  double bruit;
+  point add;
+  // Il faut prend drt proche de la solution Analytique
+  for(int i = 0; i < len; i++){
+      bruit = ((double) rand()*(1-(-1))/(double)RAND_MAX-(-1));
+      add.x = ((double) rand()*(maxP-minP)/(double)RAND_MAX-minP);
+      add.y = drt.x * add.x + drt.y + bruit;
+      vector_push(&vec,add);
+  }
+  return vec;
+}

src/regressionLin.h

@@ -6,6 +6,7 @@
 #include <stdlib.h>
 #include <string.h>
 #include <math.h>
+#include <time.h>
 
 #include "vector.h"
 
@@ -18,4 +19,7 @@ point findBestDroite(vector pts);
 // Numérique
 point descenteGradient(double precision, point drt, vector pts);
 
+// Tests
+vector nuagePoints(point drt, int len, int minP, int maxP);
+
 #endif