Add hack.

diff --git a/hacks/pldec.py b/hacks/pldec.py
new file mode 100644 (file)
index 0000000..165fd6b
--- /dev/null
+++ b/hacks/pldec.py
@@ -0,0 +1,33 @@
+#!/usr/bin/python
+
+import numpy as np
+import scipy.signal
+import matplotlib.pylab as plt
+import math
+
+radians = np.pi * 16
+points = 1024
+cut = math.pow(10, -3 / 20.0)
+print cut
+
+t = np.linspace(0, radians, points)
+inL = np.zeros((points,))
+inC = np.sin(t)
+inR = np.zeros((points,))
+inS = np.zeros((points,))
+
+# Encode
+Lt = inL + inC * cut + np.imag(scipy.signal.hilbert(inS * cut))
+Rt = inR + inC * cut - np.imag(scipy.signal.hilbert(inS * cut))
+
+# Decode
+outL = Lt
+outR = Rt
+outS = Lt - Rt
+
+plt.plot(t, outL, label='L')
+plt.plot(t, outR, label='R')
+plt.plot(t, outS, label='S')
+plt.legend()
+plt.show()
+