Filtro Fir - Bandstop

[Gleiciane Silva 0]

Pessoal,

 

Alguém sabe  me ajudar sobre filtro fir (Bandstop) em python?

[_Isis_ 202]

Sem você explicar o que você quer, não.

[Gleiciane Silva 0]

O código é esse aqui:

 

Preciso fazer um filtro de rejeita banda de um arquivo de áudio em wave

 

import pyaudio  
import wave
import struct
import matplotlib.pyplot as plt
import numpy as np
import scipy.signal as sig
from matplotlib.mlab import window_hanning

# define stream chunk (Define o pedaco do fluxo)
chunk = 20 * 4096  

# open a wav format music  (Abrir uma musica no formato wav)
f = wave.open("Daft_Punk_Get_Lucky_Cover.wav", "rb")  

# instantiate PyAudio  (Cria uma instancia de PyAudio)
p = pyaudio.PyAudio()  
formato = p.get_format_from_width(f.getsampwidth())
frame_rate = f.getframerate()
                                                                
print('formato ', formato, 'width ', f.getsampwidth(), 'channels ', f.getnchannels())

# open stream  (Abre um fluxo de audio)
stream = p.open(format=formato,
                channels=f.getnchannels(),
                rate=frame_rate,
                output=True)  
# read data  (Le os dados)
f.readframes(chunk)  
data = f.readframes(chunk)  
# Tamanho
size = 2 * chunk - 1
fmt = str(size) + 'HH'
print('len(data)', len(data), "calcsize ", struct.calcsize(fmt))
data_ = np.reshape(np.fromstring(data, 'Int16'), [chunk, 2])
x = np.arange(len(data_[:, 0]))
plt.ion()

# play stream  (Toca o audio)
while data:    
 
    
    data_ = np.reshape(np.fromstring(data, 'Int16'), [chunk, 2])
    
    plt.plot(data_[:, 0])
    
    plt.ylim((-40000, 40000))
    
    plt.draw()
    plt.pause(0.0000001)
    plt.gcf().clear()
    
    stream.write(data)
    data = f.readframes(chunk)  
    
# stop stream (Parar o fluxo)
stream.stop_stream()

# Feche o fluxo quando termina
stream.close()  

# close PyAudio (Fecha o PyAudio)
p.terminate()