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