datAc.saveStates("../datasets/MPU/"+dayString+'_'+sstring+"cylA12.txt", estAngc_A) print "time c_A12: ", c_A ''' 14 ''' startT = time.time() estAngd_A = modDA.estimateSeries(s1_fit, fingList, sensList, jointList, bnds=True, met=1) d_A = time.time()-startT print "time d_A14: ", d_A datAc.saveStates("../datasets/MPU/"+dayString+'_'+sstring+"dipA44.txt", estAngd_A) startT = time.time() estAngc_A = modCA.estimateSeries(s1_fit, fingList, sensList, jointList, bnds=True, met=1) c_A = time.time()-startT datAc.saveStates("../datasets/MPU/"+dayString+'_'+sstring+"cylA44.txt", estAngc_A) print "time c_A14: ", c_A plo.plotLeapVsMag((tim,estAngd_A),(tim,s1,s1_fit[:,:3]),head="estAngd_A vs B-field "+sstring) (timLeap, angInd, angMid, angRin, angPin) = datAc.readLeap("../datasets/160210/"+dayString+'_'+sstring+"_leap" ) #plo.plotLeapVsMag((timLeap,angInd),(tim,estAngd_A),head="leap state vs estAngd_A "+sstring,dif=False) #plo.plotLeapVsMag((tim,estAngd_A),(tim,estAngc_A),head="dip vs cyl ad-ab "+sstring,dif=False) #plo.plotLeapVsMag((tim,estAngd),(tim,estAngc),head="dip vs cyl NO ad-ab "+sstring,dif=False) plo.plotLeapVsMag((timLeap,angInd,angMid,angRin,angPin), (tim,estAngd_A[:,:4],estAngd_A[:,4:8],estAngd_A[:,8:12],estAngd_A[:,12:]),head="estAngd") plo.plotLeapVsMag((timLeap,angInd,angMid,angRin,angPin), (tim,estAngc_A[:,:4],estAngc_A[:,4:8],estAngc_A[:,8:12],estAngc_A[:,12:]),head="estAngc")
import numpy as np import matplotlib.pyplot as plt import plotting as plo import dataAcquisitionMulti as datAc #def readMagState(mStateF): sstring = 'set4' leapFile = "../datasets/160217/160217_"+'set1'+"_leap" magFile = "../datasets/160217/160217_"+sstring+"_mag" #caliFile = "../datasets/160217/160217_"+'cali3' #(tC,s1c,s2c,s3c,s4c) = datAc.readMag(caliFile) (tLeap,ind,mid,rin,pin) = datAc.readLeap(leapFile) #(tMag,s1,s2,s3,s4) = datAc.readMag(magFile1) (tMag,s1,s2,s3,s4) = datAc.readMag(magFile) #plo.plotAngles(tLeap,(ind,mid,rin,pin),head='leapStates set5') #plo.plotLeapVsMag((np.zeros((3,1)),np.zeros((3,3))),(tMag,s1,s2,s3,s4),head=sstring) #plo.plotLeapVsMag((tLeap,ind),(tC,s1c,s2c,s3c,s4c),head="cli") plo.plotLeapVsMag((tLeap,ind),(tMag,s1,s2,s3,s4),head=sstring)
cnt += 1 # add the DIP states estAng_I = estAng[:,:3] dip_I = estAng_I[:,1]*(2./3.) estAng_I = np.insert(estAng_I,2,dip_I,1) estAng_M = estAng[:,3:6] dip_M = estAng_M[:,1]*(2./3) estAng_M = np.insert(estAng_M,2,dip_M,1) estAng_R = estAng[:,6:9] dip_R = estAng_R[:,1]*(2./3) estAng_R = np.insert(estAng_R,2,dip_R,1) estAng_P = estAng[:,9:] dip_P = estAng_P[:,1]*(2./3) estAng_P = np.insert(estAng_P,2,dip_P,1) plo.plotLeapVsMag((tim,estAng_I[1:]),(tim,s1,s1_fit[:,:3]),head="estAng vs B-field "+sString) (timLeap, angInd, angMid, angRin, angPin) = datAc.readLeap("../datasets/160210/160210_"+sString+"_leap") plo.plotLeapVsMag((timLeap,angInd,angMid,angRin,angPin), (tim,estAng_I[1:],estAng_M[1:],estAng_R[1:],estAng_P[1:]),head="leap state vs est "+sString,dif=False)