#!/usr/bin/env python3

# Analysis script motor measurements
#
# 2024-07-09 ch
#
import numpy as np
import matplotlib.pyplot as plt

from numpy import genfromtxt

data = genfromtxt('data.txt', delimiter=' ')

#print(data)

dt= data[1:,0]        
iField= data[1:,1]
iRotor= data[1:,2]
phi= data[1:,3]

t=np.cumsum(dt/1000);

plt.subplot(211)
plt.title( "current sensors ")
plt.plot(t, iField)
plt.plot(t, iRotor)
plt.xlabel("t[ms]");
plt.legend(['field', 'rotor'])
plt.grid()
plt.tight_layout()

plt.subplot(212)
plt.title( "encoder ")
plt.plot(t,phi)
plt.xlabel("t[ms]");
plt.grid()
plt.tight_layout()

plt.show()

print("mean iField ")
print(np.mean(iField))

print("mean iRotor ")
print(np.mean(iRotor))