def plot_sample(ref, samp):
N = min(len(samp), 256)
fs = 48000
pl.subplot(211)
pl.grid()
pl.plot(ref)
pl.plot(samp)
pl.subplot(212)
pl.grid()
#pl.ylim((-50, 50))
pl.psd(ref, N, window=pl.blackman(N), Fs=fs)
pl.psd(samp, N, window=pl.blackman(N), Fs=fs)
def getWindowedData(self,windowlength_time=1e-12):
#returns the blackmanwindowed tdData
N=int(windowlength_time/self.dt)
w=py.blackman(N*2)
w=py.hstack((w[0:N],py.ones((self.getLength()-N*2),),w[N:]))
windowedData=self.tdData
#this could be written more elegantly?!
windowedData[:,1]*=w
windowedData[:,2]*=w
windowedData[:,3]*=w
windowedData[:,4]*=w
return windowedData