In [7]:

    

from scipy import signal
import numpy as np
import pandas as pd
from matplotlib.pyplot import *

%matplotlib inline


    

In [8]:

    

d = pd.read_csv('data/data_test.csv')
d.head()


    

    
Out[8]:






  

    

      

      
timestamps
      
Unnamed: 1
      
Unnamed: 2
      
Unnamed: 3
      
Unnamed: 4
      
Marker
    
  
  

    

      
0
      
1.489530e+09
      
47.429
      
-265.123
      
7.141
      
13.396
      
0
    
    

      
1
      
1.489530e+09
      
371.733
      
145.523
      
9.855
      
6.313
      
0
    
    

      
2
      
1.489530e+09
      
324.449
      
-208.499
      
24.210
      
30.101
      
0
    
    

      
3
      
1.489530e+09
      
290.334
      
-349.269
      
14.543
      
18.365
      
0
    
    

      
4
      
1.489530e+09
      
409.108
      
31.926
      
-0.342
      
3.121
      
0
    
  

In [9]:

    

X = np.array(d.ix[:, 1:5])
sig = X[:, 1]


    

In [10]:

    

b, a = signal.butter(2, (0.5/100.0, 10/100), btype='bandpass')
# b, a = signal.butter(1, 4/100.0, btype='low')
ecg = signal.filtfilt(b, a, sig)


    

In [6]:

    

plot(ecg)


    

    




---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-6-ee007de21b0e> in <module>()
----> 1 plot(ecg)

NameError: name 'plot' is not defined

In [154]:

    

fftfreq, fft = signal.welch(sig_f, fs=200, nperseg=10000)
plt.plot(fftfreq, fft)
plt.xlim(0,0.2)


    

    
Out[154]:





(0, 0.2)






    

In [168]:

    

np.std(sig_f)*4


    

    
Out[168]:





158.96951852357108

In [124]:

    

_ = plt.specgram(sig, Fs=200.0)


    

In [59]:

    

_ = plt.specgram(sig_f, Fs=200.0, NFFT=256, noverlap=128)


    

In [ ]:

    




    

In [ ]:

    




    

In [ ]: