Changeset 784 in ETALON

Timestamp:

Jul 18, 2018, 4:52:30 PM (6 years ago)

Author:

moutardier

Message:

new error of wago

Location:

BPM

Files:

• motor_error.txt (modified) (1 diff)
• print_datas.py (modified) (4 diffs)

BPM/motor_error.txt

@@ -72 +72 @@
 
 
+Error No.: 1812
+Designation : TMS_OVERCURRENT_HW
+Description : Overcurrent output stage(750-673)
+Possible cause/Solution : Overcurrent output stage (750-673): Short circuit or ground leak on the motor circuit
+What we have done : change motor (because the frist one had died)
+
 
 

BPM/print_datas.py

@@ -18 +18 @@
     return(Lx)
 
-
-def read_data(fichier): # program to read data from a file of data with position of motor and them Va, Vb, Vc and Vd
+def read_data_xmotor_ymotor_Va_Vb_Vc_Vd_sum_x_y(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")
+    Data = fichier.read()
+    Data_list = [[[] for i in range(9)] for j in range(3)]
+    value = ""
+    index = 0
+    for i in Data:
+        if (i == " ") :
+            if index == -1: 
+                bpm_index = int(value)
+                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:
+                Data_list[4].append(value)
+                index = 0
+            value = ""
+        elif (i == "\n") :
+            value = ""
+            index = -1
+        else:
+            value += i
+    fichier.close()
+    for i in range(5):
+        Data_list[i] = float_list(Data_list[i])
+    Data_list[0] = product_list(Data_list[0],10**-6)
+    for i in range(1,5):
+        Data_list[i] = product_list(Data_list[i],1000)
+    return(Data_list)
+
+
+
+def read_data_position_Va_Vb_Vc_Vd(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")
     Data = fichier.read()
@@ -58 +98 @@
 
 
-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])))
-    global global_index_figure
+
+
+def print_graph(L): # program to print graph of X_libera vs X_motor with a list L = [[[x_motor_position][y_motor_position][Va][Vb][Vc][Vd][sum][x_libera_position][y_libera_position]][same for BPM1][same for BPM2]], also do a linear regretion and print the equation
+    global global_index_figure # to creat different figure
     fig1 = plt.figure(global_index_figure)
     global_index_figure += 1
     plt.clf()
-    lr1 = stats.linregress(L[0],L[1]) # return tuple (pente,ordonnee a l'origine, coef de correlation, p-value, erreur standard de l'estimation)
+    for bpm_data in range(L):
+        lr1 = stats.linregress(bpm_data[0],bpm_data[7]) # return tuple (pente,ordonnee a l'origine, coef de correlation, p-value, erreur standard de l'estimation)
+        plt.plot(bpm_data[0],bpm_data[7], "r", label = "BMP_ref")
+        plt.plot(bpm_data[0],[lr1.slope*i + lr1.intercept for i in bpm_data[0]], "r--", label = "BMP_ref regression line, \nerror = "+str(round(lr1.stderr,3))+",\ncorrelation coefficient = "+str(round(lr1.rvalue,3)))
+
+
+
+
+
+
+
+    """
+    lr1 = stats.linregress(L[1],L[1]) # return tuple (pente,ordonnee a l'origine, coef de correlation, p-value, erreur standard de l'estimation)
     plt.plot(L[0],L[1], "r", label = "BMP_ref")
     plt.plot(L[0],[lr1.slope*i + lr1.intercept for i in L[0]], "r--", label = "BMP_ref regression line, \nerror = "+str(round(lr1.stderr,3))+",\ncorrelation coefficient = "+str(round(lr1.rvalue,3)))
@@ -85 +138 @@
     plt.gca().yaxis.set_tick_params(labelsize = 8)
 #    plt.show()
-
-def residu(L): # program to do graph of residu of data give in a list with [[position][Va][Vb][Vc][Vd]]
+"""
+def residu_position_Va_Vb_Vc_Vd(L): # program to do graph of residu of data give in a list with [[position][Va][Vb][Vc][Vd]]
     global global_index_figure
     fig1 = plt.figure(global_index_figure)
@@ -153 +206 @@
 
     
-def print_data(fichier, min = None, max = None):
-    L = min_max(read_data(fichier), mini = min, maxi = max)
+def print_data_position_Va_Vb_Vc_Vd(fichier, min = None, max = None):
+    L = min_max(read_data_position_Va_Vb_Vc_Vd(fichier), mini = min, maxi = max)
     print_graph(L)
     #print(L)
 
-def print_residu(fichier, min = None, max = None):
-    L = min_max(read_data(fichier), mini = min, maxi = max)
-    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_163036_0.txt")
-
-#print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_vertical_acquisition_20180612_0.txt")
-#print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_vertical_acquisition_20180612_0.txt")
+def print_residu_position_Va_Vb_Vc_Vd(fichier, min = None, max = None):
+    L = min_max(read_data_position_Va_Vb_Vc_Vd(fichier), mini = min, maxi = max)
+    residu_position_Va_Vb_Vc_Vd(L)
+
+print_data_position_Va_Vb_Vc_Vd("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_160153_0.txt",0.1,0.3)
+print_residu_position_Va_Vb_Vc_Vd("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_160153_0.txt",0.1,0.3)
+
+#print_data_position_Va_Vb_Vc_Vd("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_163036_0.txt")
+
+#print_data_position_Va_Vb_Vc_Vd("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_vertical_acquisition_20180612_0.txt")
+#print_data_position_Va_Vb_Vc_Vd("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_vertical_acquisition_20180612_0.txt")
 plt.show()