subplot(312)
plot(t)
title("trends")
xlabel("time")
ylabel("Quantized current")
subplot(313)
plot(n)
title("noise")
xlabel("time")
ylabel("Quantized current")
tight_layout(pad=1.0)
savefig('ssa2.pdf')
show()
# remove noise
tm = time.time()
ovssa = SSA.OVSSA(traces[1], None, 201, 100)
print(time.time()-tm)
title("SSA vs OVSSA trace")
xlabel('time')
ylabel("Quantized current")
plot((s+t)[50:], alpha=0.7)
plot(ovssa, alpha=0.5)
legend(['ssa', 'ovssa'])
tight_layout(pad=1.1)
savefig('ssavsovssa.pdf')
show()
# TODO: check various modes / trends oscilations without noise
# change threshholds, use at least 200 files to process
# modes include n, s, t, s+t, 1/2s+t, 1/2t+s, try different groupings
