Filtro de voz no Matlab

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

%Rodrigo Borges de Oliveira
%Eliminar os ruídos de áudios
voz = wavread('voz_audio.wav'); %carrega o áudio voz_audio.wav
vozcomsacola = wavread('voz_sacola_audio.wav'); %carrega o áudio voz_sacola_audio.wav
fs = 44100; %frequência de amostragem do áudio
%Voz e sacola separados
subplot(421);
plot(voz);
title('Voz e sacola separados sem filtro no domínio do tempo');
xlabel('Tempo(s)');
ylabel('Amplitude');
subplot(422);
my_fft(voz,fs); %fft para plotar no domínio da frequência
title('Voz e sacola separados sem filtro no domínio da frequência');
fpass = 800; %frequência de passagem
fstop = 1000; %frequência de parada
wp = (fpass/(fs/2))*pi; %frequência angular relacionada à frequência de passagem
ws = (fstop/(fs/2))*pi; %frequência angular relacionada à frequência de parada
wt = ws - wp; %frequência de transição
wc = (ws + wp)/2; %frequência de corte ou frequência de canto
M = ceil((6.2*pi/wt)) + 1; %comprimento do filtro
hd = ideal_lp(wc,M); %Criação da 
w_han = hann(M)'; %janela
h = hd.*w_han; %Hanning 
voz_filtrada = conv(h, voz); %Multiplicação em forma de convolução da janela com o áudio original, retirando as frequências indesejadas
subplot(423);
plot(voz_filtrada);
title('Voz filtrada no domínio do tempo');
xlabel('Tempo(s)');
ylabel('Amplitude');
subplot(424);
my_fft(voz_filtrada,fs);
title('Voz filtrada no domínio da frequência');
%Voz e sacola juntos
subplot(425);
plot(vozcomsacola);
title('Voz e sacola juntos sem filtro no domínio do tempo');
xlabel('Tempo(s)');
ylabel('Amplitude');
subplot(426);
my_fft(vozcomsacola,fs);
title('Voz e sacola juntos sem filtro no domínio da frequência');
vozcomsacola_filtrada = conv(h, vozcomsacola);
subplot(427);
plot(vozcomsacola_filtrada);
title('Voz com sacola filtrada no domínio do tempo');
xlabel('Tempo(s)');
ylabel('Amplitude');
subplot(428);
my_fft(vozcomsacola_filtrada,fs);
title('Voz com sacola filtrada no domínio da frequência');