From fef26c3cd47c1e6d3a9dc7592db5721b996ead6b Mon Sep 17 00:00:00 2001
From: Julian T <julian@jtle.dk>
Date: Wed, 23 Sep 2020 17:53:34 +0200
Subject: Added signal processing assignment

---
 sem5/sig/mm6/default.nix |  5 +++++
 sem5/sig/mm6/opg2.py     | 30 ++++++++++++++++++++++++++++++
 sem5/sig/mm6/opg3.py     | 29 +++++++++++++++++++++++++++++
 3 files changed, 64 insertions(+)
 create mode 100644 sem5/sig/mm6/default.nix
 create mode 100644 sem5/sig/mm6/opg2.py
 create mode 100644 sem5/sig/mm6/opg3.py

diff --git a/sem5/sig/mm6/default.nix b/sem5/sig/mm6/default.nix
new file mode 100644
index 0000000..5ba7198
--- /dev/null
+++ b/sem5/sig/mm6/default.nix
@@ -0,0 +1,5 @@
+with import <nixpkgs> {};
+stdenv.mkDerivation {
+    name = "plotting";
+    buildInputs = [ python38Packages.matplotlib python38Packages.numpy ];
+}
diff --git a/sem5/sig/mm6/opg2.py b/sem5/sig/mm6/opg2.py
new file mode 100644
index 0000000..38d65fe
--- /dev/null
+++ b/sem5/sig/mm6/opg2.py
@@ -0,0 +1,30 @@
+#!/usr/bin/env python3
+# 2020-09-23:s5/sig/mm6/opg/2 
+
+import matplotlib.pyplot as plt
+import sys
+import numpy as np
+
+lookahead = 100
+values = 50
+
+res = np.empty(values)
+
+for n in range(values):
+    s = 0
+    for k in range(lookahead):
+        h = np.exp(-k/8)
+
+        x = 1
+        if (n - k) < 0:
+            x = 0
+
+        s += x * h
+    res[n] = s
+
+x = np.arange(values)
+
+plt.plot(x, res)
+plt.show()
+
+
diff --git a/sem5/sig/mm6/opg3.py b/sem5/sig/mm6/opg3.py
new file mode 100644
index 0000000..c338c49
--- /dev/null
+++ b/sem5/sig/mm6/opg3.py
@@ -0,0 +1,29 @@
+#!/usr/bin/env python3
+# 2020-09-23:s5/sig/mm6/opg/3 #cc7f
+
+import matplotlib.pyplot as plt
+import numpy as np
+import sys
+
+# Load data
+data = []
+with open(sys.argv[1], "r") as f:
+    for line in f:
+        data.append(float(line))
+
+data = np.array(data)
+datalen = len(data)
+
+x = np.arange(datalen)
+before = plt.subplot(2, 1, 1)
+after  = plt.subplot(2, 1, 2)
+before.plot(x, data)
+
+res = np.empty(datalen)
+
+for i in range(datalen):
+    last = res[i-1] if i > 0 else 0
+    res[i] = np.exp(-1/8) * last + 1/8 * data[i]
+
+after.plot(x, res)
+plt.show()
-- 
cgit v1.2.3