def plot_pink_autocorr(beta, label):
signal = PinkNoise(beta=beta)
wave = signal.make_wave(duration=1.0, framerate=10000)
lags, corrs = autocorr(wave)
plt.plot(lags, corrs, label=label)
import sys
sys.path.insert(1, 'dsp-modulo')
import numpy as np
import thinkplot as plt
from thinkdsp import decorate
from thinkdsp import PinkNoise
pendientes = [0.0, 1.5, 2.0]
for pendiente in pendientes:
senal = PinkNoise(beta=pendiente)
wave = senal.make_wave(duration=0.5, framerate=22050)
#wave.write("ruido_rosa.wav")
#wave.plot()
#decorate(xlabel="Tiempo", ylabel= "Amplitud")
#plt.show()
etiqueta = f'pendiente= {pendiente}'
espectro = wave.make_spectrum()
espectro.plot_power(label=etiqueta, alpha=0.5)
loglog = dict(xscale="log", yscale="log")
decorate(xlabel="Frecuencia", ylabel="Poder", **loglog)
plt.show()