added convolutional theorem implementation in python

Author: pillarxyz

contents/convolutions/convolutional_theorem/code/python/convolutional_theorem.py

@@ -0,0 +1,16 @@
+from scipy.fft import fft, ifft
+import numpy as np
+
+def convolve_fft(signal1, signal2):
+    return ifft(np.multiply(fft(signal1),fft(signal2)))
+
+
+x = np.random.rand(100)
+y = [np.exp(-((i-50)/100)**2/.01) for i in range(1,101)]
+
+x /= np.linalg.norm(x)
+y /= np.linalg.norm(y)
+
+fft_output = convolve_fft(x, y)
+np.savetxt("fft.dat", fft_output)
+

contents/convolutions/convolutional_theorem/convolutional_theorem.md

@@ -44,6 +44,13 @@ Also note that the Fourier Transform is a periodic or cyclic operation, so there
 [import, lang:"julia"](code/julia/convolutional_theorem.jl)
 {% endmethod %}
 
+
+{% method %}
+{% sample lang="py" %}
+[import, lang:"python"](code/python/convolutional_theorem.py)
+{% endmethod %}
+
+
 <script>
 MathJax.Hub.Queue(["Typeset",MathJax.Hub]);
 </script>