2018-11-13 19:01

4 files changed, 42 insertions, 0 deletions

diff --git a/data/2018-11-08/ohm.png b/data/2018-11-08/ohm.png
new file mode 100644
index 0000000..8142011
--- /dev/null
+++ b/data/2018-11-08/ohm.png
diff --git a/data/2018-11-08/ohm.txt b/data/2018-11-08/ohm.txt
new file mode 100644
index 0000000..7d8a582
--- /dev/null
+++ b/data/2018-11-08/ohm.txt
@@ -0,0 +1,9 @@
+1	0.013
+2	0.023
+5	0.054
+10	0.102
+20	0.200
+50	0.497
+100	0.993
+200	1.985
+500	4.996
diff --git a/data/2018-11-08/ohm_symlog.png b/data/2018-11-08/ohm_symlog.png
new file mode 100644
index 0000000..6847b86
--- /dev/null
+++ b/data/2018-11-08/ohm_symlog.png
diff --git a/data/2018-11-08/workup_ohm.py b/data/2018-11-08/workup_ohm.py
new file mode 100644
index 0000000..95ebc42
--- /dev/null
+++ b/data/2018-11-08/workup_ohm.py
@@ -0,0 +1,33 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+# get data
+xi, yi = np.genfromtxt('ohm.txt').T
+
+# prepare figure
+plt.figure()
+ax = plt.gca()
+
+# raw
+ax.plot(xi, yi)
+ax.grid()
+ax.set_xlim(0, 500)
+ax.set_ylim(0, 5 * 1.1)
+ax.set_xlabel('set point')
+ax.set_ylabel('measured', color='C0')
+ax.set_xscale('symlog')
+
+# diff
+ax = ax.twinx()
+error = yi - (xi*0.01)
+error /= xi
+error *= 100  # to percent
+error = np.abs(error)  # absolute value
+print(error)
+ax.plot(xi, error, c="C1")
+ax.set_ylabel('percent error', color='C1')
+ax.grid()
+ax.set_ylim(0, 0.3 * 1.1)
+ax.set_yticks(np.linspace(0, 0.3, 6))
+
+plt.savefig('ohm_symlog.png')