def fitExp(r11,maxlag=None): if maxlag==None: maxlag=len(r11) xdata=np.arange(len(r11[:maxlag])) ydata=np.array(r11[:maxlag]) try: popt, pcov = tt.fit_exp_correlation(xdata,ydata) T=popt except RuntimeError: print("Error - curve_fit failed") T=0 return T
def fitExp(r11, maxlag=None): if maxlag == None: maxlag = len(r11) xdata = np.arange(len(r11[:maxlag])) ydata = np.array(r11[:maxlag]) try: popt, pcov = tt.fit_exp_correlation(xdata, ydata) T = popt except RuntimeError: print("Error - curve_fit failed") T = 0 return T
def doAutoCorr(xkey,ykey,Tkey,Lkey=None,threshold=0.1): xdata=self.data[xkey] ydata=self.data[ykey] if checkAutoCorr(xdata,ydata,threshold=threshold): try: popt, pcov = tt.fit_exp_correlation(xdata,ydata) #self.data['Txx']=abs(popt[0])*np.sqrt(np.pi)*0.5*self.data['dt'] self.data[Tkey]=popt[0]*self.data['dt'] if Lkey: self.data[Lkey]=self.data[Tkey]*self.Umean() except RuntimeError: print("Error - curve_fit failed") self.data[Tkey]=0 if Lkey: self.data[Lkey]=0 else: self.data[Tkey]=0 if Lkey: self.data[Lkey]=0