"""
pyule
=======================


"""
###########################################################
# Context
# ----------

##############################################
# Example 
# --------
#
# In the following example, we use pyule to analyse 1000 noisy signals.
# We plot them in gray and the average in red.
#
from spectrum import pyule
from scipy.signal import lfilter
from pylab import plot, mean, log10, randn
from numpy import array
allpsds = []
a = array([1,-2.2137,2.9403,-2.1697,0.9606])

for this in range(1000):
    x = lfilter([1], a, randn(1,256))
    p = pyule(x[0],4, scale_by_freq=False, NFFT=256)
    p.plot(color='gray', alpha=0.1)
    allpsds.append(p.psd)

plot(p.frequencies(), 10*log10(mean(allpsds, axis=0)), "red")