def analysis.memcof.memcof (   data,
  poles 
)

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):


Generated on Wed Nov 21 10:18:33 2007 for analysis by  doxygen 1.5.4