From 3b17facf2e63aa6b3f550add70d18d7a312a2c03 Mon Sep 17 00:00:00 2001 From: efish Date: Thu, 8 Nov 2018 15:32:58 +0100 Subject: [PATCH 1/2] add eod_cl.py --- code/eod_cl.py | 0 1 file changed, 0 insertions(+), 0 deletions(-) create mode 100644 code/eod_cl.py diff --git a/code/eod_cl.py b/code/eod_cl.py new file mode 100644 index 0000000..e69de29 From 8d39a21777636dbb1e07000bf34b7a280f84644e Mon Sep 17 00:00:00 2001 From: efish Date: Thu, 8 Nov 2018 16:28:03 +0100 Subject: [PATCH 2/2] -a --- code/eod_sim.py | 31 +++++++++++++++++++++++++++++++ 1 file changed, 31 insertions(+) diff --git a/code/eod_sim.py b/code/eod_sim.py index e69de29..861f62e 100644 --- a/code/eod_sim.py +++ b/code/eod_sim.py @@ -0,0 +1,31 @@ +import numpy as np +import matplotlib.pyplot as plt + +freq = 800 +freq2 = 820 +dt = 0.00001 + +x = np.arange(0.0,1.0,dt) +eod = np.sin(x * 2 *np.pi *freq) + np.sin(x * 2* np.pi *freq * 2)* 0.5 +eod2 = np.sin(x * 2 *np.pi *freq2) + np.sin(x * 2* np.pi *freq2 * 2)* 0.5 + +fig = plt.figure(figsize = (5., 7.5)) + +ax = fig.add_subplot(311) +ax.plot(x, eod, color = 'darkgreen', linewidth = 1.0) +ax.set_xlabel('time[s]') +ax.set_ylabel('voltage [mV]') + +ax = fig.add_subplot(312) +ax.plot(x, eod2, color = 'crimson', linewidth = 1.0) +ax.set_xlabel('time[s]') +ax.set_ylabel('voltage [mV]') + +ax = fig.add_subplot(313) +ax.plot(x, eod+ eod2 *0.05, color = 'lightblue', linewidth = 1.0) +ax.set_xlabel('time[s]') +ax.set_ylabel('voltage [mV]') + +plt.tight_layout() +plt.show() +plt.savefig('eodsTest.pdf')