From 3077eae6c9e6bf979e0a5ab56422edb20d7aa551 Mon Sep 17 00:00:00 2001
From: Daniel Friesel <daniel.friesel@uos.de>
Date: Mon, 26 Oct 2020 12:01:32 +0100
Subject: Add --plot=P/U (plot current calculated via power and voltage)

---
 bin/msp430-etv | 20 +++++++++++++++++++-
 1 file changed, 19 insertions(+), 1 deletion(-)

diff --git a/bin/msp430-etv b/bin/msp430-etv
index 18d7622..9cf0fb4 100755
--- a/bin/msp430-etv
+++ b/bin/msp430-etv
@@ -150,7 +150,7 @@ def main():
     parser.add_argument(
         "--plot",
         metavar="UNIT",
-        choices=["U", "I", "P"],
+        choices=["U", "I", "P", "P/U"],
         help="Plot voltage / current / power over time",
     )
     parser.add_argument(
@@ -391,6 +391,24 @@ def main():
             )
             plt.legend(handles=[energyhandle, meanhandle])
             plt.ylabel("Current [A]")
+        elif args.plot == "P/U":
+            (energyhandle,) = plt.plot(
+                data[1:, 0] * 1e-6,
+                power_from_energy / (data[1:, 2] * 1e-3),
+                "b-",
+                label="I=ΔE/(Δt·U)",
+                markersize=1,
+            )
+            (meanhandle,) = plt.plot(
+                data[1:, 0] * 1e-6,
+                smooth_power / (data[1:, 2] * 1e-3),
+                "r-",
+                label="mean(I, 10)",
+                markersize=1,
+            )
+            PELT().get_changepoints((smooth_power / (data[1:, 2] * 1e-3))[:400])
+            plt.legend(handles=[energyhandle, meanhandle])
+            plt.ylabel("Current [A]")
         else:
             (energyhandle,) = plt.plot(
                 data[1:, 0] * 1e-6,
-- 
cgit v1.2.3