In [53]:
signal = np.sin(np.linspace(0,10 * pi, 1000))
signal2 = np.sin(np.linspace(10,110 * pi, 1000))
signal3 = np.sin(np.linspace(0, 1000* pi, 1000))
#signal= signal2
#signal2 = signal
sig = pr.Signal(signal + signal2 + signal3, 200)
sig
Out[53]:
In [54]:
coefs = pywt.wavedec(sig, 'db1')
In [55]:
for c in coefs:
pl.figure()
print len(c), sum(c **2)
pl.plot(c)
In [ ]: