In [33]:
# Scientific modules
%pylab inline
# Goes into current dir
%cd "C:\Users\JH218595\Documents\Notebooks\TP Master Fusion"
In [34]:
mes1 =np.loadtxt('Data_2014/tp4-s1.csv', skiprows=3)
f = mes1[:,0]/1e9
S00_dB = mes1[:,1]
S10_dB = mes1[:,5]
mes2 =np.loadtxt('Data_2014/tp4-s2.csv', skiprows=3)
S20_dB = mes2[:,5]
mes3 =np.loadtxt('Data_2014/tp4-s3.csv', skiprows=3)
S30_dB = mes3[:,5]
In [35]:
plot(f, S10_dB, f, S20_dB, f, S30_dB)
Out[35]:
In [37]:
power_conservation = 10**(S00_dB/10)+10**(S10_dB/10)+10**(S20_dB/10)+10**(S30_dB/10)
plot(f, power_conservation, lw=2)
Out[37]:
In [24]: