From: Carl Hetherington <cth@carlh.net>
Date: Wed, 16 Sep 2015 14:22:11 +0000 (+0100)
Subject: Add hack.
X-Git-Tag: v2.3.3~4
X-Git-Url: https://main.carlh.net/gitweb/?p=dcpomatic.git;a=commitdiff_plain;h=0a3d5ea0900ed0f9aaa0b5a5f1b7dcec854fc4f9

Add hack.
---

diff --git a/hacks/pldec.py b/hacks/pldec.py
new file mode 100644
index 000000000..165fd6bce
--- /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()
+