From 3b47758d5f3b51827bca2afe8b3e72a0d526c35b Mon Sep 17 00:00:00 2001
From: su530201 <olivier.kaufmann@umons.ac.be>
Date: Sun, 23 Apr 2023 11:47:24 +0200
Subject: [PATCH] Adds an sp property
---
test_measure_with_ohmpi_card_3_15.py | 16 +++-------------
1 file changed, 3 insertions(+), 13 deletions(-)
diff --git a/test_measure_with_ohmpi_card_3_15.py b/test_measure_with_ohmpi_card_3_15.py
index eca45fbb..b03669f0 100644
--- a/test_measure_with_ohmpi_card_3_15.py
+++ b/test_measure_with_ohmpi_card_3_15.py
@@ -36,19 +36,9 @@ ax2.plot(k.readings[:,0], k.readings[:,2]*k.readings[:,4], '-b', marker='.', lab
ax2.set_ylabel('Vmn [mV]')
fig.legend()
plt.show()
-# compote mean voltages inside pulses
-pulses = k.readings[:,1]
-mean_vmn=[]
-mean_iab=[]
-for i in range(int(np.min(pulses)), int(np.max(pulses))+1):
- mean_vmn.append(np.mean(k.readings[k.readings[:,1]==i, 4]))
- mean_iab.append(np.mean(k.readings[k.readings[:,1]==i, 3]))
-mean_vmn = np.array(mean_vmn)
-mean_iab = np.array(mean_iab)
-print(f'Vmn: {mean_vmn}, Iab: {mean_iab}')
-sp = np.mean(mean_vmn[np.ix_([0,2,4])]-mean_vmn[np.ix_([1,3,5])])/2
-print(f'SP: {sp} mV, sp property: {k.sp}')
-r = ((k.readings[:,4]-k.readings[:,2]*sp)/k.readings[:,3])
+# compute resistances corrected for SP
+print(f'SP: {k.sp} mV')
+r = ((k.readings[:,4]-k.sp)/k.readings[:,3])
print(f'Mean resistance with sp correction : {np.mean(r):.3f} Ohms, Dev. {100*np.std(r)/np.mean(r):.1f} %')
change_config('config_default.py', verbose=False)
--
GitLab