INFORME 4

NOMBRE: CARLOS EDUARDO VASQUEZ BOLAÑOS

INFORME PRACTICA 4


UCSM-PPIE
PROCESAMIENTO DIGITAL DE SEÑALES
CUARTA PRÁCTICA DE LABORATORIO
PROCESAMIENTO DE SEÑALES DE AUDIO SPTOOLlo

I. OBJETIVOS
Procesar señales de audio
Implementar efectos sobre dichas señales
Utilizar sptool



fs=11025;
[y fs]=wavread('ejemplo.wav');
fs
sptool
sig2
g=sig2.data
fsg=sig2.Fs
sound(g,fsg)
filt1
filt1.tf
spect1plot(spect1.signal)















II. ACTIVIDADES

1. Diseño de un filtro pasabanda
Iniciar sptool
En la lista de Filters click en New. Esto crea un nuevo diseño de filtro
En la caja de sampling frequency ingrese 11025.
En la opción algorithm, seleccione Butterworth IIR filter.
Seleccione bandpass en el filter type
Ingrese 200 y 1200 para low y high stopband frequencies
Ingrese 250 y 1000 para lowy high passband frequencies
Click Apply en el panel general.
Cierre filter design tool.


2. Filtrado de una señal de audio
Cargue en MATLAB un archivo de al menos 15 segundos de duración de voz en formato .wav’.
El sonido muestreado está a 11025Hz.
>> sptool
>> [y fs]=wavread('ejemplo.wav');
>> fs

fs =

11025

>>

En sptool, en File seleccione Import
Seleccione desde el workspace elnombre del vector que contiene el sonido e importarlo usando la flecha
Seleccione sampling frequency desde el workspace e impórtela.
Ingrese un nombre, por defecto es sig1. Click OK


Desde la lista de Filters seleccione el diseño hecho en el paso 1 y click apply. Click OK.
En la lista Signals, podrá ver el resultado del filtrado, selecciónela, presionando control key y seleccione la señal deentrada. Click View.


Si no puede ver ambas señales, click en el selector de trace en el toolbar y seleccione la señal de salida. Click play y escuche la salida.




Cambie a la señal de entrada, escúchela y evalúe la diferencia.
Se escucha el sonido con menos calidad, esto es porque que se le aplico el filtro de banda.



Regresar a la pantalla de diseño de filtros y compare diferentes métodos dediseño. Verifique que R p esté puesto a 3 dB y R s a 20 dB.


Observe la curva de la respuesta del filtro Butterworth y especifique la frecuencia de corte.



Ahora cambie a un filtro Chebychev type I. Tome nota de la frecuencia de corte y la respuesta del filtro.





Repetir para el filtro Chebychev type II y el filtro elliptic. Interprete sus observaciones



Chebychev type II:Elliptic:



Para ver la impulse response de cada filtro, regresar a la ventana principal de sptool seleccione su filtro y presione en View. Una nueva ventana aparece donde se puede ver la amplitud y la fase. Desde el panel izquierdo seleccione Impulse Response. Use la herramienta zoom para ver la primera porción de la respuesta, repita en otras áreas de la respuesta.

FILTRO BUTTERWORTH:CHEVYCHEV TYPE I:





CHEVYSHEV TYPE II:





ELLIPTIC:





Desmarque Magnitude and Phase y seleccione Zeros and Poles, en la ventana de diseño cambie el tipo de filtro, ¿qué sucede con la posición de los ceros y polos?

FILTRO BUTTERWORTH








CHEVEYSHEV TYPE I

CHEVYSHEV TYPE II:

ELLIPTIC:


La posición de los ceros y polos cambia dependiendo del filtro

3. Filtrado de ruido
En esta parte seusará un filtro lowpass para eliminar ruido a altas frecuencias en una señal
>> fs=11025;
>> [y fs]=wavread('ejemplo.wav');
>> sound(y,fs)

En el prompt de MATLAB use el comando size para determinar el número de muestras del archivo anterior.
>> [y fs]=wavread('ejemplo.wav');
>>sound(y)
>>ruido=0.6*randn(size(y));
>>yruido=y+ruido
>>sound(yruido)

Genere un vector aleatorio de ruido n con el mismonúmero de muestras. Use el commando randn para generar el ruido que tiene una distribución de probabilidad Gausiana. Escálelo multiplicando el vector por 0.1.
>> [y fs]=wavread('ejemplo.wav');
>> sound(y)
>> ruido=0.1*randn(size(y));
>> yruido=y+ruido;
>> sound(yruido)


Adicione el ruido al vector original en un nuevo vector, no destruya el sonido original.
>> [y fs]=wavread('ejemplo.wav');
>>...