Changeset 779 in ETALON

Timestamp:

Jun 15, 2018, 2:13:45 PM (6 years ago)

Author:

moutardier

Message:

add fonction to print residu

Location:

BPM

Files:

• initialise_motors.py (modified) (2 diffs)
• motor_control.py (modified) (1 diff)
• print_datas.py (modified) (4 diffs)

BPM/initialise_motors.py

@@ -406 +406 @@
         motor_enable(1)
         motor_enable(2)
+        move_motor_absolute(1,10000,begin1)
+        move_motor_absolute(2,10000,begin2)
+        while (motor_get(1, 'status') != 'READY') or (motor_get(2, 'status') != 'READY'):
+            time.sleep(1)
+        print("bagin position ok ?")
+        time.sleep(15)
 #        L = [[],[],[],[],[],[]]
         for i in range(begin1, end1, pas1):            
@@ -417 +423 @@
             time.sleep(1)
         time.sleep(1)
-
         acqui_column(begin2, end2, pas2, data_name_and_path, end1)
-        
+        print("end position ok ?")
+        time.sleep(15)
         motor_disable(1)
         motor_disable(2)

BPM/motor_control.py

@@ -37 +37 @@
 #horizontal_acquisition(4000,0,-4000,"test")    
 
-acquisition(-10000,10000,10000,-50000,0,25000,"test")
+acquisition(-15000,15000,10000,-50000,50000,25000,"test")
 
 #acquisition(-10000,10000,1000,-50000,0,5000,"test")

BPM/print_datas.py

@@ -16 +16 @@
 
 
-def print_data(fichier):
+def read_data(fichier):
     fichier = open(fichier, "r")
     A = fichier.read()
@@ -52 +52 @@
     for i in range(1,5):
         L[i] = product_list(L[i],1000)
-    print_graph(L)
+    return(L)
+
+
 
 def print_graph(L):
-    a = np.array([L[0][i],L[1][i]] for i in range(len(L[0])))
+#    a = np.array([L[0][i],L[1][i]] for i in range(len(L[0])))
     fig = plt.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)
     plt.plot(L[0],L[1], "r", label = "BMP_ref")
@@ -80 +81 @@
     plt.show()
 
+def residu(L):
+    fig = plt.figure(1)
+    plt.clf()
+    plt.title("Residu")
+    lr1 = stats.linregress(L[0],L[1]) # return tuple (pente (.slope),ordonnee a l'origine (.intercept), coef de correlation (.rvalue), p-value (.pvalue ?), erreur standard de l'estimation (.stderr))
+    for i in range(len(L[1])):
+        L[1][i] = L[1][i] -  lr1.slope*L[0][i] - lr1.intercept
+    plt.plot(L[0],L[1], "r", label = "BMP_ref")
+    
+    lr2 = stats.linregress(L[0],L[2])
+    for i in range(len(L[2])):
+        L[2][i] = L[2][i] -  lr1.slope*L[0][i] - lr1.intercept
+    plt.plot(L[0],L[2], "b", label = "BMP_impr")
+
+    lr3 = stats.linregress(L[0],L[3])
+    for i in range(len(L[3])):
+        L[3][i] = L[3][i] -  lr1.slope*L[0][i] - lr1.intercept
+    plt.plot(L[0],L[3], "g", label = "BMP_ref")
+ 
+    lr4 = stats.linregress(L[0],L[4])
+    for i in range(len(L[4])):
+        L[4][i] = L[4][i] -  lr1.slope*L[0][i] - lr1.intercept
+    plt.plot(L[0],L[4], "m", label = "BMP_impr")
+    plt.legend()
+    plt.xlabel("length in million motor step")
+    plt.ylabel("tension in mV")
+    plt.show()
     
 """    
@@ -88 +116 @@
     plt.show()
 """
+def print_data(fichier):
+    L = read_data(fichier)
+    print_graph(L)
 
+def print_residu(fichier):
+    L = read_data(fichier)
+    residu(L)
 
-#print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_160153_0.txt")
+print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_160153_0.txt")
+print_residu("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_160153_0.txt")
+
 #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")
+#print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_vertical_acquisition_20180612_0.txt")