@@ -0,0 +1,595 @@
+from control import tf, bode_plot%0Afrom matplotlib.pyplot import show%0Afrom scipy import pi%0Afrom scipy.signal import cont2discrete, tf2ss%0Afrom numpy import logspace%0A%0Adef astrom_derivative_filter(Td, N):%0A    num = %5BTd, 0%5D%0A    den = %5BTd/N, 1%5D%0A    return tf(num, den), num, den%0A%0Aomega = logspace(-1, 3, num=100) * 2 * pi%0A%0Anum = %5B1%5D%0Aden = %5B0.01, 1%5D%0A%0Alpf = tf(num, den)%0Amag, phase, f = bode_plot(lpf, omega=omega, dB=True, Hz=True)%0Aprint(mag)%0Aprint(phase)%0Aprint(f)%0A%0Adelay = phase / (360.0 * f)%0A%0Alpf_ss = tf2ss(num, den)%0Aprint(lpf_ss)%0Alpf_ss_d = cont2discrete(lpf_ss, dt=0.05)%0Aprint(lpf_ss_d)%0A%0A%0Ashow()%0A%0A

Add derivative_filter design script.