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)
