#!/usr/bin/env python
from matplotlib.pyplot import *
from numpy import sin,exp,absolute,pi,arange
from numpy.random import normal
def f(t):
s1 = sin(2*pi*t)
e1 = exp(-t)
return absolute((s1*e1))+.05
t = arange(0.0, 5.0, 0.1)
s = f(t)
nse = normal(0.0, 0.3, t.shape) * s
plot(t, s+nse, 'b^')
vlines(t, [0], s)
xlabel('time (s)')
title('Comparison of model with data')
show()
|