Fix some bugs and add systheo2functions.py to notebooks/V01... and notebooks/V02...

notebooks/systheo2functions.py

@@ -2,6 +2,7 @@ import numpy as np
 import control
 from scipy import signal
 from scipy import linalg
+from itertools import zip_longest
 import matplotlib.pyplot as plt
 inv = np.linalg.inv
 pi = np.pi
@@ -14,6 +15,18 @@ log10 = np.log10
 
 s = control.TransferFunction.s
 
+def deconvolve(N,D,i):
+    f,R = signal.deconvolve(N,D)
+    print("deconvolve (im ersten Schritt N/D sonst R/D):")
+    print("f"+i+":"+str(f))
+    print("R:"+str(R))
+    return(f,R)
+
+def convolve(R):
+    R = signal.convolve(R,[1, 0])
+    print("convolve with z^(-1): R=R*z^(-1):"+str(R)+"\n")
+    return(R)
+
 def print_det(M):
     print("Determinante: "+str(np.around(np.linalg.det(M))))