Changeset 782 in ETALON

Timestamp:

Jul 11, 2018, 12:25:10 PM (6 years ago)

Author:

moutardier

Message:

correction of issu in print_data.py and add comment

Location:

BPM

Files:

• motor_control.py (modified) (2 diffs)
• print_datas.py (modified) (6 diffs)

BPM/motor_control.py

@@ -10 +10 @@
 # moteur 2 -> verticale
 
+
+#motor_reset(2)
+#motor_disable(2)
+
+motor_enable(1)
+move_motor_absolute(1,5000,-5000)
+move_motor_absolute(1,5000,0)
+motor_disable(1)
+
 """
-motor_enable(2)
-move_motor_absolute(2,5000,-50000)
-while motor_get(1, 'status') != 'READY':
-    time.sleep(1)
-time.sleep(5)
+for i in range(5):
+    motor_enable(2)
+    move_motor_absolute(2,5000,-5000)
+    while motor_get(2, 'status') != 'READY':
+        time.sleep(1)
+    time.sleep(10)
+
+
+    move_motor_absolute(2,5000,0)
+    while motor_get(2, 'status') != 'READY':
+        time.sleep(1)
+    motor_disable(2)
+    time.sleep(10)
 """
-""" 
-move_motor_absolute(2,5000,0)
-while motor_get(1, 'status') != 'READY':
-    time.sleep(1)
-"""
+
 #controller_hello()
 #move_motor_relative(1,1000,-1000)
@@ -37 +50 @@
 #horizontal_acquisition(4000,0,-4000,"test")    
 
-acquisition(-100000,100000,20000,-50000,50000,10000,"test")
+#acquisition(-100000,100000,20000,-50000,50000,10000,"test")
 
 #acquisition(-10000,10000,1000,-50000,0,5000,"test")

BPM/print_datas.py

@@ -3 +3 @@
 from scipy import stats
 
-def float_list(l):
+def float_list(l): # convert a list of char in a list of float 
     L = []
     for i in l:
@@ -9 +9 @@
     return(L)
 
-def product_list(L,x, b = None):
-    if b is None:
-        Lx = []
-        for i in L:
-            Lx.append(i*x)
-        return(Lx)
-    else:
-        Lx = []
-        for i in L:
-            Lx.append(i*x +b)
-        return(Lx)
+def product_list(L,x, b = 0): # linear convertion of a list return L*x + b (b is 0 as default) 
+    Lx = []
+    for i in L:
+        Lx.append(i*x +b)
+    return(Lx)
 
 
-def read_data(fichier):
+def read_data(fichier): # program to read data from a file of data with position of motor and them Va, Vb, Vc and Vd
     fichier = open(fichier, "r")
-    A = fichier.read()
-    L = [[],[],[],[],[]]
-    a = ""
-    l = 0
-    for i in A:
+    Data = fichier.read()
+    Data_list = [[],[],[],[],[]]
+    value = ""
+    index = 0
+    for i in Data:
         if (i == " ") :
-            if l == 0: 
-                L[0].append(a)
-                l +=1
-            elif l == 1:
-                L[1].append(a)
-                #print(a)
-                l += 1
-            elif l == 2:
-                L[2].append(a)
-                l += 1
-            elif l == 3:
-                L[3].append(a)
-                l += 1
+            if index == 0: 
+                Data_list[0].append(value)
+                index +=1
+            elif index == 1:
+                Data_list[1].append(value)
+                #print(value)
+                index += 1
+            elif index == 2:
+                Data_list[2].append(value)
+                index += 1
+            elif index == 3:
+                Data_list[3].append(value)
+                index += 1
             else:
-                L[4].append(a)
-                l = 0
-            a = ""
+                Data_list[4].append(value)
+                index = 0
+            value = ""
         elif (i == "\n") :
-            a = ""
+            value = ""
         else:
-            a += i
+            value += i
     fichier.close()
     for i in range(5):
-        L[i] = float_list(L[i])
-    L[0] = product_list(L[0],10**-6)
+        Data_list[i] = float_list(Data_list[i])
+    Data_list[0] = product_list(Data_list[0],10**-6)
     for i in range(1,5):
-        L[i] = product_list(L[i],1000)
-    return(L)
+        Data_list[i] = product_list(Data_list[i],1000)
+    return(Data_list)
 
 
 
-def print_graph(L):
+def print_graph(L): # program to print graph of data set in a list with [[position][Va][Vb][Vc][Vd]], also do a linear regretion and print the equation
 #    a = np.array([L[0][i],L[1][i]] for i in range(len(L[0])))
     fig = plt.figure(1)
@@ -86 +80 @@
     plt.show()
 
-def residu(L):
+def residu(L): # program to do graph of residu of data give in a list with [[position][Va][Vb][Vc][Vd]]
     fig = plt.figure(1)
     plt.clf()
@@ -122 +116 @@
 """
 
-def min_max(L,mini = None ,maxi = None):
+def min_max(L,mini = None ,maxi = None): #program to cut a list L = [[position][Va][Vb][Vc][Vd]] with minimum value of position = mini and maximum value of position = maxi 
     if mini is None :
-        mini = -1
+        mini = L[0][0]
     if maxi is None :
-        maxi = len(L[0])
-    i_min = 0
-    i_max = 0
+        maxi = L[0][len(L[0])-1]
+    if mini > maxi:
+        mini,maxi = maxi,mini
+    index_list = []
     for i in range(len(L[0])):
-        if  L[0][i] <= mini :
-            i_min = i
-        if  L[0][i] <= maxi :
-            i_max = i
-    if i_max < len(L[0]) :
-        i_max += 1
-    return([L[0][i_min:i_max],L[1][i_min:i_max],L[2][i_min:i_max],L[3][i_min:i_max],L[4][i_min:i_max]])
+        if  L[0][i] >= mini and L[0][i] <= maxi :
+            index_list.append(i)
+    print(index_list)
+    return([[L[0][i] for i in index_list],[L[1][i] for i in index_list],[L[2][i] for i in index_list],[L[3][i] for i in index_list],[L[4][i] for i in index_list]])
 
 
-def convertion_step_mm(L,begin_mm,end_mm):
+def convertion_step_mm(L,begin_mm,end_mm): # length convertion 
     begin_step = L[0]
     end_step = L[len(L)-1]
@@ -155 +147 @@
     
 def print_data(fichier, min = None, max = None):
-    L = min_max(read_data(fichier),min,max)
+    L = min_max(read_data(fichier), mini = min, maxi = max)
     print_graph(L)
     #print(L)
@@ -163 +155 @@
     residu(L)
 
-#print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_160153_0.txt",0.1,0.3)
-#print_residu("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_160153_0.txt",0.1,0.3)
+print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_160153_0.txt",0.1,0.3)
+print_residu("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_160153_0.txt",0.1,0.3)
 
 #print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_163036_0.txt")