PureLightCurve.py 文件源码

def ComputePer(self,beta):
        """
        compute a periodigram for the data and compared it to a theorical line

        """

        datafft=np.fft.fft(self.data)
        per=np.abs(datafft)**2
        py.plot(per)
        x=np.array([100.*i+1 for i in range(self.length*self.res)])
        y=per[0]*(1./x)**(beta/2)
        py.loglog(x,y,c='red')
        py.show()