In [29]:

    

import numpy as np
from scipy import fftpack
import matplotlib.pyplot as plt


    

In [67]:

    

data = np.loadtxt("TEST01.csv")
Ts = 0.200
N = len(data)
T_end = (N * Ts) / 60
print T_end
time = np.linspace(0, T_end, N)


    

    




132.583333333

In [41]:

    

F = fftpack.fft(data)
xf = np.linspace(0.0, 1.0/(2.0*Ts), N/2)


    

    




2.5

In [81]:

    

# THRESHOLD DETECTOR
idx = np.where(data > 137)
print len(idx[0])


    

    




120

In [45]:

    

%matplotlib inline
plt.rcParams['figure.figsize'] = (10, 6)

fig, (a, b) = plt.subplots(2,1)

a.plot(time, data)

a.set_xlim([0, T_end])
a.grid(True)
a.set_xlabel('Time [min.]')
a.set_ylabel('ADC val.')

b.plot(xf, 2.0/N * np.abs(F[:N//2]))
b.set_ylim([0, 5])
b.set_xlim([0, 0.03])
b.set_ylabel('Amplitude')
b.set_xlabel('Freq. [Hz]')
b.grid(True)