From 0a3d5ea0900ed0f9aaa0b5a5f1b7dcec854fc4f9 Mon Sep 17 00:00:00 2001 From: Carl Hetherington Date: Wed, 16 Sep 2015 15:22:11 +0100 Subject: [PATCH] Add hack. --- hacks/pldec.py | 33 +++++++++++++++++++++++++++++++++ 1 file changed, 33 insertions(+) create mode 100644 hacks/pldec.py 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() + -- 2.30.2