def __init__(self, inputsig, outputsig, timevector): self.t = timevector self.input = inputsig self.output = outputsig self.f = makefreqvect(timevector) self.spectra = CalcSpectra(inputsig, outputsig) self.bode = rwkbode.AveBodeFromSpectra(self.spectra) if len(self.bode.mag) < len(self.f): self.f = self.f[0:len(self.bode.mag)]
def __init__(self, signal, time_vector, seedfreq=None, seedphase=None): self.signal = signal self.time = time_vector self.f = makefreqvect(self.time) N = max(shape(signal)) self.comp = squeeze(fft(signal, None, 0)*2/N) self.mag = abs(self.comp) self.phase = arctan2(imag(self.comp),real(self.comp))*180.0/pi self.dBmag = 20.0*log10(self.mag)