Definition at line 9 of file memcof.py. 00009 : 00010 n=len(data) 00011 xms=Numeric.dot(data,data)/float(n) 00012 wk1=data[:-1] 00013 wk2=data[1:] 00014 d=Numeric.zeros(poles,Numeric.Float) 00015 00016 for k0 in range(poles): 00017 num=Numeric.dot(wk1,wk2) 00018 denom=Numeric.dot(wk1, wk1)+Numeric.dot(wk2, wk2) 00019 d[k0]=2.0*num/denom 00020 xms*=(1.0-d[k0]**2) 00021 if k0!=0: 00022 d[:k0]=wkm-d[k0]*wkm[-1::-1] 00023 if k0!=poles-1: 00024 wkm=d[:k0+1] 00025 wk1, wk2 = wk1[:-1]-wkm[k0]*wk2[:-1], wk2[1:]-wkm[k0]*wk1[1:] 00026 00027 return xms, d 00028 00029 def evlmem(fdt, d, xms):
|