##// END OF EJS Templates
HG_REPOSITORIES » LPP » INSTRUMENTATION » USERS » PELLION » VHD_Lib
RSS

Clone from http://pc-instru.lpp.polytechnique.fr:5000/VHD_Lib

Few fixes....
Few fixes. Whole LFR simulation WIP.
Jeandet Alexis - - Load All Authors

File last commit:

r647:3005781845db default
r682:c53e1b6b3045 default
Show More
plot_results.ipynb
113 lines | 44.3 KiB | text/plain | TextLexer
/ tests / Validation_LFR_Filters_Offset / plot_results.ipynb
History | Annotation | Raw |Copy content |Copy permalink
In [2]:
import numpy as np
import matplotlib.pyplot as plt
import glob 
import pandas as pds
In [2]:
for folder in glob.glob("./2016*"):
    inputsig = pds.read_csv(folder+"/input.txt",sep=" ",header=None)
    inputsig.columns=[ "CH{}".format(i) for i in range(inputsig.shape[1])]
    outputsig = pds.read_csv(folder+"/output_f0.txt",sep=" ",header=None)
    outputsig.columns=["Tstamp"]+[ "CH{}".format(i) for i in range(outputsig.shape[1]-1)]
    f, axarr = plt.subplots(1,2,figsize=(14, 6))
    (outputsig.filter(regex="CH*")- inputsig*0.8912)[150:].plot(ax=axarr[0])
    axarr[0].legend(loc='upper right')
    (outputsig.filter(regex="CH*")[:100]/0.8912).plot(ax=axarr[1])
    axarr[1].legend(loc='upper right')
    plt.show()

---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-2-9719ee11a4bd> in <module>()
----> 1 for folder in glob.glob("./2016*"):
      2     inputsig = pds.read_csv(folder+"/input.txt",sep=" ",header=None)
      3     inputsig.columns=[ "CH{}".format(i) for i in range(inputsig.shape[1])]
      4     outputsig = pds.read_csv(folder+"/output_f0.txt",sep=" ",header=None)
      5     outputsig.columns=["Tstamp"]+[ "CH{}".format(i) for i in range(outputsig.shape[1]-1)]

NameError: name 'glob' is not defined
In [5]:
folder="./simulation/"
inputsig = pds.read_csv(folder+"/input.txt",sep=" ",header=None)
inputsig.columns=[ "CH{}".format(i) for i in range(inputsig.shape[1])]
outputsig = pds.read_csv(folder+"/output_f2.txt",sep=" ",header=None)
outputsig.columns=["Tstamp"]+[ "CH{}".format(i) for i in range(outputsig.shape[1]-1)]
f, axarr = plt.subplots(1,2,figsize=(14, 6))
(outputsig.filter(regex="CH*")- inputsig*0.8912)[150:].plot(ax=axarr[0])
axarr[0].legend(loc='upper right')
(outputsig.filter(regex="CH*")[:100]/0.8912).plot(ax=axarr[1])
axarr[1].legend(loc='upper right')
plt.show()