def frequency_response(filter, n, N ):
s = pysig([ cos((_2pi*k*n)/N) for k in range(N)])
u = filter(s)
return rms(u)/rms(s)
result[k] = sum
return result
if __name__ == "__main__":
import random
import pylab as p
import numpy as n
a = pysig(20, 0, 1)
a.data = numpy.array([ random.random() for x in a.data ])
F = pyfilter( n.ones(20)/20 )
b = F(a)
print a.data
print b.data
#print '\n\n', a.__dict__
p.plot(a.data, 'bo')
p.plot(b.data, 'ro')
p.show()
