statistics.py 文件源码

def compute(self,signal):

        check_MD(signal.X)

        N=signal.N
        M,L=signal.X.shape
        sigma2=signal.sigma2
        H=np.matrix(signal.H)
        X=np.matrix(signal.X)
        Rx=X*X.H/L
        PH=np.eye(N)-H*lg.pinv(H)
        D=np.matrix(np.diag(1j*np.arange(N)))
        M=H.H*D.H*PH*D*H

        CRB=(sigma2/2)*lg.inv(np.real(np.multiply(M,Rx)))
        self.w=np.diag(CRB)

        return self