Procesamiento Digital de Señales Con Matlab

PROCESAMIENTO DIGITAL DE SEALES CON MATLAB
http://lonely113.blogspot.com
Generacin y Grfica de Seales

2
Sea una seal:

Dnde: A: Amplitud Fo: Frecuencia phi: Fase La representacin de esta seal en tiempo discreto est dada por:
Donde: Fs=1/Ts: Frecuencia de muestreo. La grfica de seales discretas se realiza con el comando:
stem(t,xt)
Ejemplo: Generacin y grfica de una seal muestreada a 8 KHz.

3
Se genera una seal sinusoidal de 400 Hz y amplitud 2.

2 1.5 1
>> F0=400; >> A=2; >> phi=pi/4; >> Fs=8000; >> Ts=1/Fs; >> t=-0.002:Ts:0.002; >> xt=A*sin(2*pi*F0*t+phi); >> stem(t,xt)
0.5 0 -0.5 -1 -1.5 -2 -2
-1.5
-1
-0.5
0.5
1.5 x 10
2
-3
Ejemplo: Seales no sinusoidales peridicas.

4
Generacin de una seal diente de sierra.

>> F0=400; >> A=2; >> Fs=8000; >> Ts=1/Fs; >> t=-0.002:Ts:0.002; >> xt=sawtooth(2*pi*F0*t); >> stem(t,xt) >> hold on >> plot(t,xt) >> xlabel('time (s)'); >> ylabel('x(t)');
1 0.8 0.6 0.4 0.2
x(t)
0 -0.2 -0.4 -0.6 -0.8 -1 -2
-1.5
-1
-0.5
0 time (s)
0.5
1.5 x 10
2
-3
Ejemplo: Seal aperidica.

5
Muestreo de la seal Sinc.

2
>> F0=400; >> A=2; >> Fs=8000; >> Ts=1/Fs; >> t=-0.003:Ts:0.003; >> xt=A*sinc(2*F0*t); >> stem(t,xt)
1.5
0.5
-0.5 -3
-2
-1
2 x 10
3
-3
Generacin de Ruido
6
Se utilizan los comandos:

Y=rand(rows,cols): Genera ruido con distribucin uniforme. Y=randn(rows,cols):Genera ruido Gaussiano. Dnde: rows, cols: Indica la dimensin de la matriz de ruido aleatorio a generar. Para observar el histograma se utiliza el comando: hist(y,m): Representa el ruido "y" mediante m "contenedores".
Ejemplo: Generacin y Representacin de Ruido Gaussiano

7
4 3.5 3
x 10
>> Gnoise=randn(1,1e6); >> hist(Gnoise,100);
2.5 2 1.5 1 0.5 0 -6
-4
-2
Submuestreo y Sobremuestro
8
Submuestreo
xtDown=downsample(xt,N)
La seal xDown tendr una frecuencia de muestro Fs/N. Submuestrear la seal significa conservar cada N-sima muestra y eliminar las muestras restantes.
Sobremuestreo
xtUp=upsample(xt,N)
La seal xUp tendr una frecuencia de muestreo NFs.
Sobremuestrear la seal significa introducir N-1 ceros entre muestras consecutivas.
Ejemplo: Submuestreo.
9
Se submuestrea la seal sinusoidal generada anteriormente (Fs=8 KHz) para obtener una seal xtDown muestreada a 2 KHz.
>> F0=400; >> A=2; >> phi=pi/4; >> Fs=8000; >> Ts=1/Fs; >> t=-0.002:Ts:0.002; >> xt=A*sin(2*pi*F0*t+phi); >> stem(t,xt) >>xtDown=downsample(xt,4); >> tDown=downsample(t,4); >> hold on >> stem(tDown,xtDown,'r');
2 1.5 1 0.5 0 -0.5 -1 -1.5 -2 -2
-1.5
-1
-0.5
0.5
1.5 x 10
2
-3
Procesamiento de Audio
10
Para escuchar un tono de seal en Matlab se utiliza el comando: soundsc(xt,Fs) Donde:
xt: Tono a escuchar. Fs: Frecuencia de muestreo.

Por ejemplo: Para escuchar la seal sinusoidal generada anteriormente:
>> soundsc(xt,Fs)
Grabacin de Audio
11
Para grabar una seal audible mediante la tarjeta de sonido y un micrfono se utilizan los siguientes comandos:
r=audiorecorder: Crea un objeto de grabacin. record(r): Inicio de grabacin. pause(r) ,stop(r): Pausa y finalizacin. play(r): Escuchar la grabacin y=getaudiodata(r): Para obtener la matriz que contiene las muestras de la seal audible. sta es la seal que se puede procesar. Fs=r.SampleRate: Para obtener la frecuencia de muestreo.
Ejemplo:
12
>> r=audiorecorder; >> record(r) >> pause(r) >> record(r) >> stop(r) >> play(r) >> y=getaudiodata(r); >> Fs=r.SampleRate Fs = 8000
Guardar como Archivo de Audio

13
Para guardar la seal de audio que se acaba de grabar se utiliza el comando:

wavwrite(y,Fs,grabacion)
Se guarda en formato wav en la carpeta de trabajo de Matlab.

>> wavwrite(y,Fs,'grabacion'); >> which grabacion.wav C:\Users\Casa\Documents\MATLAB\ grabacion.wav
Leer un Archivo de Audio

14
Para leer un archivo de audio en formato wav que se encuentra en la carpeta de trabajo de Matlab se utiliza la funcin:
[xt,Fs]=wavread(nombre_de_archivo)
Se guardan las muestras y la frecuencia de muestro en xt y Fs respectivamente.
>> [xt,Fs]=wavread('tuner1');
Anlisis en Frecuencia
15
Transformada Rpida de Fourier: xF=fft(xt,nFFT) Donde: nFFT=nmero de puntos a representar.
Comandos adicionales:
magxF=abs(xF): Magnitud de la fft de xt. phasexF=angle(xF): Fase de la fft de xt. phasexF=unwrap(phasexF): Para obtener la apropiada informacin de fase.
Ejemplo: Obtencin de La FFT de la Funcin Compuerta. 16

Magnitud de xF 8
>> Fs=2000; >> Ts=1/Fs; >> t=-0.005:Ts:0.005; >> xt=rectpuls(t,0.004); >> nFFT=64; >> xF=fft(xt,nFFT); >> magxF=abs(xF); >> phasexF=angle(xF); >> phasexF=unwrap(phasexF); >> figure,stem(magxF); >>title(Magnitud de xF) >> figure,stem(phasexF) >>title(Fase de xF)
7 6 5 4 3 2 1 0
10
20
30
40
50
60
70
Fase de xF 0 -5 -10 -15 -20 -25 -30 -35 -40
10
20
30
40
50
60
70

Magnitud de xF 8
El comando fft obtiene la transformada de Fourier en la regin comprendida entre 0 y la frecuencia de muestro Fs. Se requiere desplazar una parte de la grfica para obtener la regin negativa de la transformada de Fourier.
>> magxF=fftshift(magxF); >> phasexF=fftshift(phasexF); >> stem(magxF); >>title(Magnitud de xF) >> figure,stem(phasexF); >>title(Fase de xF)
7 6 5 4 3 2 1 0
10
20
30
40
50
60
70
Fase de xF 0 -5 -10 -15 -20 -25 -30 -35 -40
10
20
30
40
50
60
70

Magnitud de xF 8
Los grficos anteriores no presentan el eje de frecuencia correcto. Se debe crear el vector de eje de frecuencias para graficar correctamente.
>> f_esp=Fs/nFFT; >> fNyquist=Fs/2; >> f_inicio=-fNyquist; >> f_fin=fNyquist-f_esp; >> f_eje=f_inicio:f_esp:f_fin; >> stem(f_eje,magxF) >> figure,stem(f_eje,phasexF)
7 6 5 4 3 2 1 0 -1000
0 -5 -10 -15 -20 -25 -30 -35 -40 -1000
-800
-600
-400
-200
0 200 Fase de xF
400
600
800
1000
-800
-600
-400
-200
200
400
600
800
1000
Sistemas LTI
19
Todos los sistemas LTI pueden ser analizados como Filtros LTI discretos.
Ecuacin de Diferencias Para Filtros Causales

20
Dnde: x[n]: Entrada y[n]: Salida
En el dominio Z
Funcin de transferencia en el dominio Z
Representacin de Filtros en Matlab

21
Para representar la funcin de transferencia del filtro:
en Matlab:
>> num=[1 1.1]; >> den=[1 0 -0.1];
Respuesta al Impulso y Convolucin

22
Para obtener la respuesta al impulso se utilizan los comandos:
Impz(num,den): Grafica la respuesta al impulso del filtro. Por defecto se considera Fs=1 Hz y
grafica las 10 primeras muestras.
Impz(num,den,n,Fs): Grafica la respuesta al impulso. Se especifica el nmero de muestras a

graficar y la frecuencia de muestro utilizada.
[h,t]= Impz(num,den,n,Fs): Guarda los vectores de magnitud y eje de tiempo en h y t

respectivamente.
Ejemplo: Graficando la respuesta al impulso del filtro creado anteriormente.

23
Impulse Response 1.4 1.2
Amplitude
0.8
0.6
>> num=[1 1.1]; >> den=[1 0 -0.1]; >> impz(num,den); >> Fs=8000; >> [h,t]=impz(num,den,32,Fs); >> figure,stem(t,h)
0.4
0.2
3 4 n (samples)
1.4
1.2
0.8
0.6
0.4
0.2
0.5
1.5
2.5
3.5 x 10
4
-3
Seal de Salida del Filtro Para una Seal de Entrada

24
Se realiza mediante los comandos:

y=conv(h,x) Donde: h es la respuesta al impulso del filtro. x es la seal de entrada.
y=filter(num,den,x)
Donde: x es la seal de entrada.
Ejemplo:
25
Espectro de xt
Se obtendr la seal de salida del filtro para una seal de entrada audible:
>> [xt,Fs]=wavread('tuner1'); >> [h,t]=impz(num,den,10,Fs); >> y1=convn(h,xt); >> soundsc(xt,Fs); >> soundsc(y1,Fs); >> y2=filter(num,den,xt); >> soundsc(xt,Fs); >> soundsc(y2,Fs);
35 30 25 20 15 10 5 0 0 500 1000 1500 2000 2500 3000 3500 4000
Espectro de y2
80 70 60 50 40 30 20 10 0 500 1000 1500 2000 2500 3000 3500 4000
Anlisis de Filtros en el Dominio de la Frecuencia

26
Relacin Entrada Salida de un filtro en el dominio de la frecuencia.
Respuesta de magnitud: Respuesta de fase:
Retardo de grupo:
Respuesta de Magnitud y Fase de Filtros

27
Se utilizan los comandos: [H,F]=freqz(num,den): Por defecto asume 512 puntos equidistantemente espaciados y el eje de frecuencias es normalizada. [H,F]=freqz(num,den,n,Fs): Donde se especifica nmero de puntos n y la frecuencia de muestreo Fs. [H,F]=freqz(num,den,f,Fs): Donde f es un vector en el que se especifica las frecuencias particulares en las que se obtendr la respuesta de magnitud y fase.
Ejemplo: Respuesta en frecuencia del filtro creado anteriormente.

28
10
Magnitude (dB)
-10
-20
0.1
0.2
0.3 0.4 0.5 0.6 0.7 0.8 Normalized Frequency ( rad/sample)
0.9
>> num=[1 1.1]; >> den=[1 0 -0.1]; >> Fs=8000; >> freqz(num,den) >>figure,freqz(num,den,1024,Fs); >> f=[0,200,400]; >> [H,F]=freqz(num,den,f,Fs);
Phase (degrees)
-50 -100 -150 -200
0.1
0.2
0.9
10
Magnitude (dB)
-10
-20
500
1000
1500 2000 2500 Frequency (Hz)
3000
3500
4000
Phase (degrees)
-50 -100 -150 -200
500
1000
1500 2000 2500 Frequency (Hz)
3000
3500
4000
Grfico de Retardo de Grupo

29
Se utilizan los comandos:

[G,F]=grpdelay(num,den) [G,F]=grpdelay(num,den,n,Fs) [G,F]=grpdelay(num,den,f,Fs)
Los argumentos son iguales a los del caso anterior.
Ejemplo: Continuando con el ejemplo anterior

30
12 10
Group delay (samples)
>> [G,F]=grpdelay(num,den); >> grpdelay(num,den); >> figure,grpdelay(num,den,1024,Fs);

Group delay (samples)
0.1
0.2
0.9
12
10
500
1000
1500 2000 2500 Frequency (Hz)
3000
3500
4000
Anlisis de Polos y Ceros de Funciones de Transferencia

31
Para estudiar la estabilidad de los sistemas LTI.

Funcin de transferencia:
Ganancia polo cero:
Obtencin de Polos y Ceros

32
Se realiza una transformacin del dominio tf al dominio zp.
[z,p,k]=tf2zp(num,den)
Dnde los vectores: z: ceros p: polos k: ganancia de cada numerador de tf.
De manera inversa:
[num,den]=zp2tf(z,p,k)
Ejemplo: Con el filtro implementado anteriormente.

33
>> num=[1 1.1]; >> den=[1 0 -0.1]; >> [z,p,k]=tf2zp(num,den) z= -1.1000 p= 0.3162 -0.3162 k= 1 >> [num,den]=zp2tf(z,p,k) num = 0 1.0000 1.1000 den = 1.0000 0 -0.1000
Grfico de Polos y Ceros

34
Se utiliza el comando:
zplane(z,p)
Si hay polos fuera del crculo unitario significa que el sistema es inestable. Continuando con el ejemplo anterior:
1 0.8 0.6 0.4
Imaginary Part
0.2 0 -0.2 -0.4 -0.6 -0.8 -1 -1 -0.5 0 Real Part 0.5 1
>> zplane(z,p)
Uso del FVTool

35
FVTool (Filter Visualization Tool) Es una herramienta que simplifica el trabajo de estudiar el comportamiento de los filtros.
Se trata de una interfaz GUI para el anlisis de filtros.
Llamada al FVtool
36
Se llama a la herramienta con el comando:
fvtool(num,den)
Donde: num y den corresponden al filtro en estudio. Continuando con el ejemplo anterior:
>>fvtool(num,den)
Configuracin de FVtool
37
El eje de frecuencias por defecto no se corresponde con los valores correctos. Para configurarlo correctamente:
Men Analysis>Sampling Frecuency. En el cuadro Fs se introduce el valor de la frecuencia de muestreo.
La herramienta Fvtool permite visualizar diferentes respuestas del filtro en estudio

38
Respuesta de Magnitud
Respuesta de fase Respuestas de magnitud y fase
Respuesta de retardo de grupo

Respuesta de retardo de fase Respuesta al impulso
Respuesta al escaln
Grfico de polos y ceros Etc
Ejemplo:
39
Respuesta al impulso
Respuestas de Magnitud y fase
Visualizacin de Mltiples Filtros

40
Es posible visualizar la respuesta de varios filtros en simultneo. Ejemplo:
>> num1=[1 1]; >> den1=[1]; >> num2=[1 -1]; >> den2=[1]; >> fvtool(num1,den1,num2,den2)
Diseo de Filtros Digitales en Matlab

41
Los filtros elementales son los siguientes:
Diseo de Filtros IIR

42
IIR: Infinite Impulse Response.
Respuesta en frecuencia:
Tipos y Caractersticas de Filtros IIR

43
Especificaciones Para un Filtro Pasabanda

44
Para el diseo de un filtro pasabanda se especifica:

fstop1: Frecuencia stop 1. fpass1: Frecuencia pass 1. fpass2: Frecuencia pass 2. fstop2: Frecuencia stop 2. Rp: Rizado en la banda de paso en dB. Rs: Atenuacin en la banda de supresin en dB.
Diseo de un Filtro Pasabanda

45
El diseo de filtros en Matlab se realiza con valores normalizados. Para esto se realiza:

Especificar Fs. Fnyquist=Fs/2 Ws(1)=fstop1/Fnyquist Wp(1)=fpass1/Fnyquist Wp(2)=fpass2/Fnyquist Ws(2)=fstop2/Fnyquist
Dnde Ws y Wp son vectores de frecuencias normalizadas.

Diseo de un Filtro Butterworth

46
1. Determinar el orden del filtro para las especificaciones:
[nB,WnB]=buttord(Wp,Ws,Rp,Rs)
Wp, Ws, Rp y Rs corresponden a las especificaciones. WnB: Frecuencia a -3 dB.
2. Obtener los coeficientes del filtro.
[numB,denB]=butter[nB,WnB]
En el diseo de filtros Pasabanda o rechazo de banda el orden del filtro resulta ser en realidad el doble del especificado por el comando buttord.
Ejemplo: Diseo de un filtro Butterworth

47
>> fstop1=1000; >> fpass1=1200; >> fpass2=2200; >> fstop2=2400; >> Rp=1; >> Rs=60; >> Fs=8000; >> Fnyquist=Fs/2; >> Ws(1)=fstop1/Fnyquist; >> Wp(1)=fpass1/Fnyquist; >> Wp(2)=fpass2/Fnyquist; >> Ws(2)=fstop2/Fnyquist; >> [nB,WnB]=buttord(Wp,Ws,Rp,Rs) nB = 22 WnB = 0.2958 0.5551 >>[numB,denB]=butter(nB,WnB);
Diseo de un Filtro Elptico

48
Se realiza un procedimiento similar al caso anterior.
>> [nE,WnE]=ellipord(Wp,Ws,Rp,Rs) nE = 6 WnE = 0.3000 0.5500 >> [numE,denE]=ellip(nE,Rp,Rs,Wp); Se puede observar que la cantidad de coeficientes es mucho menor que en el caso de un filtro Butterworth.
Verificacin del Diseo de Filtros

49
Con la herramienta Fvtool.

>> fvtool(numB,denB,numE,denE)
Se observa que ambos filtros son estables.
Ejemplo: Filtrando el archivo tuner1.wav

50
Espectro de xt 35 30 25 20 15 10
>> [xt,fs]=wavread('tuner1'); >> yt=filter(numE,denE,xt); >> soundsc(xt,fs) >> soundsc(yt,fs)
5 0 0 500 1000 1500 2000 2500 3000 3500 4000
Espectro de yt
10
0 500 1000 1500 2000 2500 3000 3500 4000
Llamada a los Comandos de Diseo de Filtros

51
Los ms usuales son los siguientes:

Butterworth Pasabajas
[num,den]=butter(nB,Wn)
Elptico Pasabajas
[num,den]=ellip(nE,Rp,Rs,Wp)
Butterworth Pasaaltas
[num,den]=butter(nB,Wn,high)
Elptic0 Pasaaltas
[num,den]=ellip(nE,Rp,Rs,Wp,high)
Diseo de Filtros FIR

52
FIR: Finite Impulse Response.
Respuesta en frecuencia:
Diseo Mediante Mtodo de Ventanas

53
Se multiplica la respuesta en el tiempo de un filtro ideal con una funcin ventana. Truncando la respuesta ideal a un nmero finito de elementos.
1 0.8 0.6 0.4 0.2 0 -0.2 -0.4 -50
Se realiza mediante el comando:
-40
-30
-20
-10
10
20
30
40
50
num=fir1(n,Wp)
Dnde:
n: Orden del filtro Wp: Vector que contiene las frecuencias de corte. Si slo tiene un elemento se obtiene un filtro pasabajas y si contiene dos elementos se obtiene un filtro pasabanda.
Ejemplo:
54
Se disea un filtro mediante el mtodo de ventanas para las siguientes especificaciones: Fs=8000 Hz fpass1=2200 Hz fpass2=3200 Hz Rp=1 dB Rs=60 dB Orden=127
Ejemplo:
55
>> Fs=8000; >> Fnyquist=Fs/2; >> fpass1=2200; >> fpass2=3200; >> Rp=1; >> Rs=60; >> Wp(1)=fpass1/Fnyquist; >> Wp(2)=fpass2/Fnyquist; >> n=127; >> num=fir1(n,Wp); >>fvtool(num)
Wintool Para diseo de Funciones Ventana

56
>> wintool Initializing Window Design & Analysis Tool ..... done.
Herramienta con interfaz GUI que permite disear y visualizar la respuesta en el tiempo y en frecuencia de varios tipos de funciones ventana.
Propiedades de Algunas Funciones Ventana de Longitud M

57
Tipo
Anchura de Lbulo Principal
Mnima Atenuacin en la Banda de Rechazo 20.9
Anchura de Banda de Transicin
Rectangular
2Fs/(2M+1)
0.46Fs/M
Hann
Hamming Blackman
4Fs/(2M+1)
2Fs/(2M+1) 6Fs/(2M+1)
43.9
54.5 75.3
1.55Fs/M
1.55Fs/M 2.78Fs/M
Diseo de un Filtro FIR Utilizando una Ventana Diseada Previamente

58
Se disea la funcin ventana Blackman de longitud 128 mediante la herramienta wintool y se exporta a Matlab con Save to workspace. El nombre por defecto es window_1.
Continuacin
59
>> numBlackman=fir1(n,Wp,window_1); >>fvtool(numBlackman)
Ejemplo: Filtrando el archivo tuner1.wav con el filtro diseado.

60
Espectro de xt 35 30 25 20 15 10
de
audio
>> [xt,Fs]=wavread('tuner1'); >> yt=filter(numBlackman,1,xt); >> soundsc(xt,Fs) >> soundsc(yt,Fs)
5 0 0 500 1000 1500 2000 2500 3000 3500 4000
Espectro de yt
0 500 1000 1500 2000 2500 3000 3500 4000
Diseo de Filtros Mediante la Herramienta FDATool

61
Fdatool (Filter Design and Analysis tool) es una herramienta con interfaz GUI para el diseo de filtros digitales.
>> fdatool
Diseo de un Filtro usando FDATool

62
Se disear un filtro pasaaltas con las siguientes especificaciones:

Pasaaltas Tipo IIR Elptico Mnimo orden posible Fs=8000 Hz Fstop=2800 Hz Fpass=3200 Hz Rs=60 dB Rp=1 dB Se introducen todas estas Especificaciones en la herramienta y un clic en Design Filter.
Exportacin al rea de Trabajo

63
Se realiza el siguiente procedimiento: 1. Clic derecho en Structure en el campo Current filter information y en seguida se escoge convert to single section.
2. Men File>Export Se le da nombres al numerador y denominador del filtro. Luego clic en Export.
Ejemplo: Filtrando el archivo de audio tuner1.Wav

64
Espectro de xt 35 30 25 20 15 10
>> [xt,Fs]=wavread('tuner1'); >> yt=filter(numHp,denHp,xt); >> soundsc(xt) >> soundsc(yt)
5 0 0 500 1000 1500 2000 2500 3000 3500 4000
Espectro de yt
1.2
0.8
0.6
0.4
0.2
0 500 1000 1500 2000 2500 3000 3500 4000
65

Procesamiento Digital de Señales Con Matlab

Cargado por

Copyright:

Formatos disponibles

Procesamiento Digital de Señales Con Matlab

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Procesamiento Digital de Señales Con Matlab

Cargado por

Copyright:

Formatos disponibles

PROCESAMIENTO DIGITAL DE SEALES CON MATLAB

Generacin y Grfica de Seales

Sea una seal:

Ejemplo: Generacin y grfica de una seal muestreada a 8 KHz.

Se genera una seal sinusoidal de 400 Hz y amplitud 2.

0.5 0 -0.5 -1 -1.5 -2 -2

Ejemplo: Seales no sinusoidales peridicas.

Generacin de una seal diente de sierra.

0 -0.2 -0.4 -0.6 -0.8 -1 -2

Ejemplo: Seal aperidica.

Muestreo de la seal Sinc.

Se utilizan los comandos:

Ejemplo: Generacin y Representacin de Ruido Gaussiano

>> Gnoise=randn(1,1e6); >> hist(Gnoise,100);

2.5 2 1.5 1 0.5 0 -6

Sobremuestrear la seal significa introducir N-1 ceros entre muestras consecutivas.

Para escuchar un tono de seal en Matlab se utiliza el comando: soundsc(xt,Fs) Donde:

xt: Tono a escuchar. Fs: Frecuencia de muestreo.

Guardar como Archivo de Audio

Para guardar la seal de audio que se acaba de grabar se utiliza el comando:

Se guarda en formato wav en la carpeta de trabajo de Matlab.

Leer un Archivo de Audio

Transformada Rpida de Fourier: xF=fft(xt,nFFT) Donde: nFFT=nmero de puntos a representar.

Ejemplo: Obtencin de La FFT de la Funcin Compuerta. 16

Fase de xF 0 -5 -10 -15 -20 -25 -30 -35 -40

Ejemplo: Obtencin de La FFT de la Funcin Compuerta. 17

Fase de xF 0 -5 -10 -15 -20 -25 -30 -35 -40

Ejemplo: Obtencin de La FFT de la Funcin Compuerta. 18

Ecuacin de Diferencias Para Filtros Causales

Dnde: x[n]: Entrada y[n]: Salida

Representacin de Filtros en Matlab

Para representar la funcin de transferencia del filtro:

Respuesta al Impulso y Convolucin

Para obtener la respuesta al impulso se utilizan los comandos:

Impz(num,den,n,Fs): Grafica la respuesta al impulso. Se especifica el nmero de muestras a

[h,t]= Impz(num,den,n,Fs): Guarda los vectores de magnitud y eje de tiempo en h y t

Ejemplo: Graficando la respuesta al impulso del filtro creado anteriormente.

Seal de Salida del Filtro Para una Seal de Entrada

Se realiza mediante los comandos:

35 30 25 20 15 10 5 0 0 500 1000 1500 2000 2500 3000 3500 4000

80 70 60 50 40 30 20 10 0 500 1000 1500 2000 2500 3000 3500 4000

Anlisis de Filtros en el Dominio de la Frecuencia

Relacin Entrada Salida de un filtro en el dominio de la frecuencia.

Respuesta de magnitud: Respuesta de fase:

Respuesta de Magnitud y Fase de Filtros

Ejemplo: Respuesta en frecuencia del filtro creado anteriormente.

0.3 0.4 0.5 0.6 0.7 0.8 Normalized Frequency ( rad/sample)

-50 -100 -150 -200

0.3 0.4 0.5 0.6 0.7 0.8 Normalized Frequency ( rad/sample)

1500 2000 2500 Frequency (Hz)

-50 -100 -150 -200

1500 2000 2500 Frequency (Hz)

Grfico de Retardo de Grupo

Se utilizan los comandos:

Los argumentos son iguales a los del caso anterior.

Ejemplo: Continuando con el ejemplo anterior

Group delay (samples)

>> [G,F]=grpdelay(num,den); >> grpdelay(num,den); >> figure,grpdelay(num,den,1024,Fs);

0.3 0.4 0.5 0.6 0.7 0.8 Normalized Frequency ( rad/sample)