import LCD
from math import sin, cos, exp, pi
from time import sleep

FileName = 'Current_Halogen.txt'

f = open(FileName, "rt")
Data = f.readlines()
f.close()

n = len(Data)
N = 384

t = [0]*N
X = [0]*N
for i in range(0,N):
    X[i] = float(Data[i])
    t[i] = i
    
Navy = LCD.RGB(0,0,5)
White = LCD.RGB(200,200,200)
LtBlue = LCD.RGB(100,100,250)

LCD.Init()
LCD.Clear(Navy)
LCD.Plot(t,X)

sleep(3)

c = [0] * 25
w0 = 2*pi/N
for n in range(1,25):
    a = b = 0
    for i in range(0,N):
        a += X[i] * cos(n*i*w0) / N
        b += X[i] * sin(n*i*w0) / N
    c[n] = (a**2 + b**2) ** 0.5
    
LCD.Bar(c, White, LtBlue)
LCD.Title(FileName	, White, Navy)