pac.py 文件源码

def estimate(self, phases, phases_lohi):
        num_ts, ts_len = np.shape(phases)

        self.pairs = [(r1, r2) for r1 in xrange(0, num_ts) for r2 in xrange(r1, num_ts)]

        pacs_ts = np.zeros((num_ts, num_ts, ts_len), dtype=np.complex)
        pacs_avg = np.zeros((num_ts, num_ts))

        for pair in self.pairs:
            p1, p2 = pair
            phase1 = phases[p1, ]
            phase1_lohi = phases_lohi[p2, ]

            ts, avg = self.estimator.estimate_pair(phase1, phase1_lohi)

            pacs_ts[pair] = ts
            pacs_avg[pair] = avg

        return pacs_ts, pacs_avg