Processamento digital de sinal

Índice
Objectivos2
Introdução2
Código de Matlab3
Gráficos4

Objectivo
Este trabalho tem como principal objectivo, desenvolver a componente prática de Processamento Digital de Sinal, aplicando os conhecimentos adquiridos nas aulas teórico-práticas. Com o nosso trabalho pretende-se implementar um filtro IIR, rejeita-banda, no Matlab a um sinal de áudio que apresenta ruído a uma determinada banda de frequência. Posto isto, na conclusão do trabalho, espera-se que estejamos bem familiarizados com as ferramentas de trabalho do Matlab.
Introdução
O sinal consiste numa música com ruído. Por forma a remover o ruído, filtrando o sinal usando um filtro IIR rejeita-banda, chebyshev do tipo 2. A característica utilizada é o rejeita-banda porque queremos eliminar o ruído que se encontra numa banda de frequência distinta. A frequência do ruído, no nosso caso, é de 10515 Hz.
C:\Users\Ze\Desktop\Sem Título.png
Nas nossas especificações, definimos Wp1=10490 Hz e Wp2=10540 Hz, mas não conseguimos visualizar esses parâmetros no freqz devido ao facto de estarmos a utilizar o filtro chebyshev do tipo 2 que é mais difícil de conseguir essas especificações, pois só garante atenuação das frequências na banda de corte.

Código de Matlab y = wavread ('C:\Users\Ze\Desktop\electrotecnica\2ºAno\2ºSemestre\PDS\projecto\MP_4.wav');
[Y, Fs] = wavread ('C:\Users\Ze\Desktop\electrotecnica\2ºAno\2ºSemestre\PDS\projecto\MP_4.wav'); sound(Y,Fs); spectrogram (y,512,400,512,Fs,'yaxis');
Wp=[10490/24000,10540/24000];
Ws=[10500/24000,10530/24000];
[N,Wn]=cheb2ord(Wp,Ws,3,40);
[B,A]=cheby2(N,80,Wn,'stop'); figure(2) freqz(B,A,3000,48000); figure(3) y1=filter(B,A,Y); spectrogram (y1,512,400,512,Fs,'yaxis'); sound(y1,Fs); Gráficos
C:\Users\Ze\Desktop\Sem Título2.png
Fig. 1- Espectrograma do sinal do ruído.
Especificações do filtro:
Wp1=10490/24000 Hz
Wp2=10540/24000 Hz
Ws1=10500/24000 Hz
Ws2=10530/24000 Hz
Rp=3 dB
Rs=40 dB
N=5
C:\Users\Ze\Desktop\Sem Título2.png