def goDTW2(n1,n2,f0,dt,rmsNoise,maxStrain,g,f): #for uo fr = 2 kmin = -2 kmax = 2 rmin = -1.5 rmax = 1.5 nl = 1+fr*(n2-n1) dst = dt/fr sf = Sampling(n1,dt,0.0) sl = Sampling(nl,dst,0.0) dw = DynamicWarpingWT(nl,fr,kmin,kmax) uoi,udmin = dw.finduo(f,g) uo = uoi[0] #uo = 350 print 'uo= ',uo,'\n' e = dw.computeErrors(f,g) d = dw.accumulate(e) for ir in range(4): dr = 0.001*pow(10,ir) nr = int(floor((rmax-rmin)/dr)) + 1 sr = Sampling(nr,dr,rmin) u = dw.findShiftsX(e,sr,uo) u = mul(u,dt) dmin = dw.findShiftsX(e,sr,uo,True) print 'dr= ',dr print 'nr= ',nr print 'esum= ',dmin[0],'\n' #plotMatrix(etran(d),sf,sl) plotMatrix(etran(d),sf,sl,u=u,title='dr= '+str(dr))