statistics.py 文件源码

def get_rv(self,signal):

        #extract parameters
        N=signal.S.shape[0]
        r,p=self.detector.A.shape

        if self.type=="PFA":

            if self.detector.estimate_sigma2==True:
                rv=f(r,N-p)
            else:
                rv=chi2(r)

        if self.type=="PD":

            #construct matrices
            A=np.matrix(self.detector.A)
            b=np.matrix(self.detector.b).T
            H=np.matrix(signal.H)
            x=np.matrix(signal.X)
            sigma2=signal.sigma2

            #compute the non central parameter
            term1=A*x-b
            nc=term1.T*lg.inv(A*lg.inv(H.T*H)*A.T)*term1/sigma2
            nc=min(np.asscalar(nc),self.nc_max)

            if self.detector.estimate_sigma2==True:
                rv=ncf(r,N-p,nc)
            else:
                rv=ncx2(r,nc)

        return rv


###### LINEAR MODEL WITH INTERFERENCE ##################################