source: branches/numpy_misc/tools/demos/matplotdemo3.py @ 7244

import numpy as num
from pylab import randn, conv, plot, psd, subplot, show


# Create signal
dt = 0.01
t = num.arange(0, 10, dt)    # Time vector
signal = 0.1*num.sin(2*num.pi*t)

# Create noise
noise = conv(randn(len(t)), num.exp(-t/0.05))*dt
noise = noise[:len(t)]

# Create noisy signal and plot
noisy_signal = signal + noise

subplot(211)
plot(t, noisy_signal)
subplot(212)
psd(noisy_signal, 512, 1/dt)

show()