def add_avg(self): self.D1Ma += self.D1.AvgMag self.D1Pha += self.D1.AvgPhase self.D1vMa += self.D1.vAvgMag self.D1vPha += self.D1.vAvgPh self.D1aPow += self.D1.vAvgPow self.D2Ma += self.D2.AvgMag self.D2Pha += self.D2.AvgPhase self.D2vMa += self.D2.vAvgMag self.D2vPha += self.D2.vAvgPh self.D2aPow += self.D2.vAvgPow if self.doCorrel: self.covAvgMat += getCovMatrix(self.D1.scaledI, self.D1.scaledQ, self.D2.scaledI, self.D2.scaledQ, self.lags)
def add_avg(self): self.D1Ma += self.D1.AvgMag self.D1Pha += self.D1.AvgPhase self.D1vMa += self.D1.vAvgMag self.D1vPha += self.D1.vAvgPh self.D1aPow += self.D1.vAvgPow # Digitizer 2 Values self.D2Ma += self.D2.AvgMag self.D2Pha += self.D2.AvgPhase self.D2vMa += self.D2.vAvgMag self.D2vPha += self.D2.vAvgPh self.D2aPow += self.D2.vAvgPow self.covAvgMat += getCovMatrix(self.D1.scaledI, self.D1.scaledQ, self.D2.scaledI, self.D2.scaledQ, self.lags)
from scipy import signal numtaps = 100 f = 0.2 fir1 = signal.firwin(numtaps, f) plt.ion() I1 = np.random.random(int(1e3)) # dummy test data Q1 = np.random.random(int(1e3)) I2 = np.random.random(int(1e3)) Q2 = np.random.random(int(1e3)) I1new = signal.filtfilt(fir1, 1, I1, padlen=10) plt.figure() plt.plot(fir1) plt.figure() plt.plot(I1) plt.figure() plt.plot(I1new) lags = 20 matrix = getCovMatrix(I1, Q1, I2, Q2) # calculate covariance matrix # plt.figure() # plt.plot(matrix[0]) # plt.figure() # plt.plot(matrix[6])
# Digitizer 2 Values D2M, D2Ph = D2.get_AvgMagPhs() D2Ma = D2Ma + D1M D2Pha = D2Pha + D2Ph D2vM, D2vPh = D2.get_vAvgMagPhs() D2vMa = D2vMa + D2vM D2vPha = D2vPha + D2vPh vdata = vdata + vm.get_val() I1 = D1.scaledI Q1 = D1.scaledQ I2 = D2.scaledI Q2 = D2.scaledQ covAvgMat = covAvgMat + getCovMatrix(I1, Q1, I2, Q2, lags) # D1Lvl = D1Lvl + D1.levelcorr # D2Lvl = D2Lvl + D2.levelcorr # print cz, time()-t0conv # Recording data in Memory DS11.record_data(covAvgMat/np.float(corrAvg),kk,jj,ii) DSP.record_data(vdata/np.float(corrAvg),kk,jj,ii) DSP_LD1.record_data(D1.levelcorr,kk, jj, ii) DS2mD1.record_data(D1Ma/np.float(corrAvg), D1Pha/np.float(corrAvg), kk, jj, ii) DS2vD1.record_data(D1vMa/np.float(corrAvg), D1vPha/np.float(corrAvg) ,kk, jj, ii) DSP_PD1.record_data((D1aPow/np.float(corrAvg)) ,kk, jj, ii)