Taller Final

Señales y Sistemas Lineales – Facultad de Minas – Universidad Nacional de Colombia

Taller Final Transformada z y Análisis de Fourier

Señales y Sistemas Lineales – Eliana I. Arango Z.Taller Final Transformada z y Análisis de Fourier

Objetivos: Obtener la transformada Z de funciones en el dominio de n utilizando Matlab. Obtener la transforma Z inversa de funciones en el dominio de Z utilizando Matlab. Analizar los sistemas descritos en el dominio Z Analizar la estabilidad de sistemas en el dominio Z. Aplicar los teoremas del análisis de Fourier a las señales representadas. Identificar las principales cuestiones implicadas en el muestreo.

Orientación para la realización del taller:

A continuación se describe el trabajo que debe realizar. Encontrará algunos cálculos teóricos solicitados que puede hacerlos a mano y escanearlos para facilitar el trabajo. Debe presentar un informe detallado en un archivo en Word que incluya todos los cálculos teóricos, los análisis solicitados, todas las respuestas obtenidas en matlab, todas las conclusiones solicitadas y los pantallazos que usted considere necesarios para darle solución a las preguntas del taller.

Parte 1 - Transformada Z de funciones en el dominio discreto

La transformada Z desempeña el mismo papel para las señales y sistemas en tiempo discreto que el que tiene la transformada de Laplace para señales y sistemas de tiempo continuo. El análisis de un sistema discreto en el dominio de Z permite obtener algunas características de la respuesta del sistema y así mismo inferir sobre su estabilidad.La transformada Z se puede calcular en Matlab haciendo uso de las variables simbólicas y de la función ztrans, para hacer uso de este comando se recomienda que la función a la cual se desea calcular la transformada Z sea función de n.

Ejemplo 1

syms n; % Se define la variable simbólica nf = n^4; % Se define una función f que depende de nztrans(f) % Se calcula la transforma Z

Actividad 1

Consulte la ayuda del comando ztrans. Implemente el Ejemplo 1, en un archivo .m y verifique la respuesta. Calcule la transformada Z para las siguientes funciones básicas usando el comando ztrans de Ma-

tlab. Realice esta actividad en un archivo .m

a) nu [n ]b) nα nu [n ]c) cos (αn)u [n ]

Compare la respuesta con los datos que se encuentran en la siguiente tabla.



1 nu [n ]z

( z−1 )2|z|>1

2nα nu [n]

zα

( z−α )2|z|>|α|

3cos (nΩ)u [n]

z2−zcos(Ω)z2−2 zcos (Ω)+1

|z|>1

Transforma Z inversa

Análogamente a la transformada inversa de Laplace, la transformada Z inversa permite obtener la respuesta en el dominio del tiempo discreto de un sistema en el dominio de Z.

La transformada Z inversa se puede calcular en Matlab haciendo uso de las variables simbólicas y del comando iztrans, para el uso de dicho comando se recomienda que la función se encuentre descrita en términos de Z. Veamos el siguiente ejemplo obtenido de la ayuda del Matlab.

Ejemplo 2

syms z % Se define la variable simbólica nf = 2*z/(z-2)^2; % Se define una función f que depende de ziztrans(f) % Se calcula la transformada inversa de la función f

Actividad 2

Consulte la ayuda del comando iztrans. Implemente el Ejemplo 2 en un archivo .m y analice el resultado. Calcule la transformada Z inversa de las siguientes funciones haciendo uso del comando iztrans de

Matlab. Esta actividad se debe entregar en un archivo .m.

a. X [ z ]= 2Z3+Z(Z−2 )2 (Z−1 )

b. H [z ]= 2ZZ−0.8

c.F [ z ]= 2 Z(Z+1)

(Z−13 )(Z2+ 14 ) ( z2+4 Z+5 )

Funciones de Transferencia

El poder real de la transformada de Laplace se halla en el análisis del comportamiento dinámico de sistemas en tiempo continuo. De una manera análoga, el de la transformada Z está en el análisis del comportamiento dinámico de sistemas en tiempo discreto



La respuesta y [n ] de un sistema con una respuesta al impulso h [n ], ante una entrada arbitraría x [n], está

dada por la convolución y [n ]=x [n]∗h[n]. Dado que la operación de convolución se convierte en un

producto con la transformada Z. Se tiene:

H ( z )=Y ( z )X (z )

Los sistemas en tiempo discreto se pueden describir de una manera muy conveniente mediante diagramas de bloques, y es posible escribir las funciones de transferencia de manera directa a partir del diagrama de bloques.

Ejemplo 3.

Considere el siguiente diagrama de bloques, a partir del diagrama obtenga:

La ecuación en diferencias que representa el sistema. La función de transferencia del sistema. Represente la función de transferencia en Matlab. Realice en Simulink la representación del diagrama de bloques y de la función de transferencia;

aplique un escalón discreto con un Ts=0.01s.

Solución

La ecuación en diferencias que describe el sistema que está representando en el diagrama de blo-ques es:

y [n ]=2 x [n ]−x [n−1 ]−12y [n−1 ]

Y es equivalente a:

Y (z )=2 X (z )−z−1 X (z )−12z−1Y (z )

La función de transferencia corresponde a:

H ( z )=Y ( z )X (z )

=2 z−1

2+12

La función de transferencia de un sistema en tiempo discreto se representa de manera similar a las correspondientes en sistemas en tiempo continuo, la única diferencia es que se debe especificar el tiempo de muestreo (o el tiempo entre cada muestra) del sistema.

fs=100; % Se define la frecuencia de muestreo del sistemaTs=1/fs; % Se define el tiempo de muestreo del sistemaH=tf([2 -1],[2 1/2], Ts) % Se define los componentes de la función de transferencia



La implementación en Simulink del diagrama de bloques y de la función de transferencias se obser-va en la siguiente figura.

Actividad 3

Considere el siguiente diagrama de bloques, a partir del diagrama obtenga:

La ecuación en diferencias que representa el sistema. La función de transferencia del sistema. Represente la función de transferencia en Matlab. (considere un tiempo de muestreo igual 0.01s) Realice en Simulink la representación del diagrama de bloques y de la función de transferencia;

aplique un escalón discreto con un Ts=0.01s. Compare las respuestas.

Ejemplo 4

Sea y [n ]=0.6 y [n−1 ]+4 x [n ], se pide:

Hallar la función de transferencia Hallar los polos y los ceros de la función de transferencia, y determinar la estabilidad del sistema. Realizar un análisis gráfico de estabilidad. Obtener el diagrama de bode de la respuesta del sistema, use el comando dbode, el tiempo de

muestro es de 0.01ms Hallar la respuesta del sistema ante una entrada tipo escalón unitario. Utilizando el comando step.



Solución

La función de transferencia del sistema se halla de la siguiente manera:

Y ( z )=0.6 z−1Y ( z )+4 X (z )

Y ( z )X ( z )

= 4 zz−0.6

→Función de transferencia

Los polos y los ceros de un sistema discreto se hallan de manera análoga para un sistema en tiem-po continuo, esta operación se puede realizar en Matlab mediante el siguiente comando:

[ceros,polos,ganancia]=tf2zpk([4 0],[1 -0.6])

De acuerdo a la respuesta este sistema no tiene ceros, tiene una ganancia de 4 y tiene un polo en z=0.6.

Como presenta un único polo en z=0.6<1 entonces se puede decir que el sistema es estable, esto se corrobora usando el comando isstable.

F=tf([4 0],[1 -0.6], 0.01)isstable(F)

Usando el comando rlocus se puede observar gráficamente la posición de los polos y ceros dentro del circulo unitario, por lo que se podría determinar la estabilidad de un sistema de una manera sencilla

F=tf([4 0],[1 -0.6], 0.01)rlocus(F)

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

-0.8

-0.6

-0.4

-0.2

0

0.2

0.4

0.6

0.8

1Root Locus

Real Axis

Imag

inar

y A

xis

En la gráfica anterior se observa la posición del único polo dentro del circulo unitario, por lo que se confirma la estabilidad del sistema.

El diagrama de bode se obtienen utilizando la función dbode.



6

8

10

12

14

16

18

20

Ma

gn

itud

(d

B)

10-1

100

101

102

-40-35

-30-25

-20-15

-10-5

05

Fa

se (

de

g)

Diagrama de Bode

Frequencia (Hz)

La respuesta del sistema al escalón unitario se realizó de la siguiente manera:

F=tf([4 0],[1 -0.6], 0.01)step(F)

Respuesta al escalon

Tiempo (sec)

Am

plit

ud

0 0.05 0.1 0.15 0.2 0.254

5

6

7

8

9

10

Actividad 4

Sea F (z)=

2Z (Z+1)

(Z−13 )(Z2+ 14 ) ( z2+4 Z+5 ) , se pide

Hallar los polos y los ceros de la función de transferencia, y determinar la estabilidad del sistema. Realizar un análisis gráfico de estabilidad. Consulte la ayuda del comando rlocus. Obtener el diagrama de bode de la respuesta del sistema, use el comando dbode, el tiempo de

muestro es de 0.01ms. Cómo se interpreta el diagrama de bode. Qué información se puede extraer a partir del diagrama de boe



Hallar la respuesta del sistema ante una entrada tipo escalón unitario. Utilizando el comando step. Halle los valores del tiempo de establecimiento y el tiempo de levantamientos, y determine si estos valores se ven afectador por el tiempo de muestreo. Realice una explicación coherente de todas sus conclusiones.

Actividad 5

Evalúe la estabilidad de los siguientes sistemas de manera analítica, posteriormente evalúe la estabilidad con el comando isstable de Matlab. Finalmente compare los resultados.

Nota: Antes de realizar esta actividad consulte la ayuda del comando isstable

a) H ( z )= zz−2

b)H ( z )= z

z2−78

c)H ( z )= z

z2−32z+98

d) H ( z )= z2−1z3−2 z2+3.75 z−0.5625

Parte 2 - Análisis de Fourier

Esta parte se subdivide en una parte teórica que consistirá en responder una serie de preguntas, luego tiene una parte práctica para la cual deberá utilizar unos “applet” que se le especifican al final.

1. Identifique la señal a la que pertenece cada uno de los espectros dados a continuación:

(a) (b) (c)

2. Calcule la densidad espectral de potencia de una señal senoidal x(t) = Asen(2 fo t)3. Haga una figura en la que se explique claramente el proceso de muestreo ideal de una señal cualquiera en el dominio

del tiempo. (La figura debe incluir la señal analógica, la función de muestreo y la señal muestreada)4. Se tiene una señal cualquiera con un espectro dado por la siguiente figura, donde B=1000 Hz:



Dibuje el espectro de la señal muestreada idealmente (usando un tren de impulsos) para los siguientes casos:a. Se usa una frecuencia de muestreo S=2000 Hz. ¿Qué caso de muestreo se observa?b. Se usa una frecuencia de muestreo S=2500 Hz. ¿Qué caso de muestreo se observa?c. Se usa una frecuencia de muestreo S=1000 Hz. ¿Qué caso de muestreo se observa?

5. Para los ejercicios del tema de Series de Fourier vamos a utilizar un applet que se encuentra en la siguiente dirección web:

http://www.jhu.edu/signals/fourier2/index.html

He probado el applet en mi computador y me ha funcionado correctamente. Deben usar el internet explorer e instalar el pluggin de java, de ser necesario, para poder obtener el funcionamiento del applet.

Les adiciono una figura de cómo se ve el applet:

Inicialmente vamos a verificar el fenómeno de Gibbs para cada una de las señales que se pueden calcular en el applet.

Seleccione “Rectangular Pulse” Escriba el número 1 en la casilla “Fourier Series Coefficientes” Click en “Calculate” Click en “Table”, abre una pequeña ventana con los valores numéricos de los coeficientes calculados. Observe la gráfica roja que acaba de aparecer sobre el pulso azul y verifique que es muy poco aproximado al

pulso rectangular azul Ahora haga Click en el símbolo “+” para aumentar el número de coeficientes Debe observar que la gráfica roja se aproxima cada vez más a la gráfica azul al aumentar el número de coefi-

cientes En la parte inferior también puede observar los espectros de magnitud y fase correspondientes.

6. Complete la siguiente TABLA 3

http://www.jhu.edu/signals/fourier2/index.html



Función a aproximar en la sumatoria de armónicos de la Serie de Fourier

Número de Coeficientes de la serie de Fourier que usted considera suficientes para lograr una buena aproximación

Rectangular Pulse

Bipolar Pulse

Sawtooth

Triangle

Exponential

Noise

El espectro de Magnitud de cada aproximación de Fourier que se observa debajo de la Figura principal y la tabla de coeficientes de la Serie que aparece al hacer Click sobre “Table” nos puede orientar acerca de cuántos coeficientes de la serie de Fourier correspondiente serán necesarios para obtener un alto porcentaje de la potencia de la señal, utilizando el Teorema de Parseval.

7. Complete la TABLA 4 utilizando el Teorema de Parseval:Función aproximada por la Serie de Fourier

Número de Coeficientes de Mayor Magnitud en el espectro

Potencia contenida hasta el armónico representado en el coeficiente de la columna anterior

Rectangular Pulse

Bipolar Pulse

Sawtooth

Triangle

Exponential

Noise

8. Para ejemplificar el teorema del muestreo de Shannon y Nyquist utilizaremos un applet que se encuentra en la siguiente dirección web:

http://www.jhu.edu/signals/sampling/index.html

He probado el applet en mi computador y me ha funcionado correctamente. Deben usar el internet explorer e instalar el pluggin de java, de ser necesario, para poder obtener el funcionamiento del applet.Les adiciono una figura de cómo se ve el applet:

http://www.jhu.edu/signals/sampling/index.html



Inicialmente vamos a seleccionar la señal SINC.

Como podemos observar la gráfica azul es la señal en el dominio del tiempo y la señal roja es el espectro de magnitud de la transformada de Fourier de la señal, es decir está graficada en el dominio de la frecuencia.

Podemos seleccionar entre las señales Pulse, Sen, Sinc. También podemos dibujar una señal usando el cursor.

Podemos introducir la frecuencia de muestreo: “Sampling frecuency” para la señal seleccionada y al hacer click en “Sample”, nos aparece la gráfica del espectro de Fourier de la señal muestreada.

Luego con el objetivo de recuperar la señal muestreada podemos introducir una frecuencia de corte para un filtro pasabajas ideal, (ideal low-pass cutoff frequency), y al hacer click en “filter” nos aparece el espectro de Fourier filtrado (en rojo) y la señal en el dominio del tiempo (en azul), que se debe aproximar bastante a las señal original si hemos utilizado una frecuencia de muestreo adecuada.

Recordemos el Teorema del Muestreo: Una señal limitada en banda a una frecuencia B (es decir que tiene un ancho de banda igual a B= Máxima frecuencia contenida en la señal y que se puede observar en el espectro de magnitud de su transformada de Fourier), puede recuperarse completamente de sus muestras tomadas a una frecuencia de dos veces su ancho de banda, es decir a una frecuencia 2B que se denomina Frecuencia de Nyquist ó frecuencia crítica de muestreo.

9. Vamos a comprobar los resultados para los tres casos del teorema de muestreo para la señal Sinc. Utilizaremos copy-paste para copiar las figuras obtenidas en el applet para cada una de las frecuencias de muestreo utilizadas.

Seleccione la señal Seno Cardenal (click en Sinc) En la gráfica roja puede observar que el ancho de banda de la señal Sinc es aproximadamente B=15 rad/sec.

Este dato nos informa sobre la frecuencia del filtro pasabajas que se utilizará para recuperar la señal y que se mantiene constante aunque se modifique la frecuencia de muestreo. También nos informa sobre la Frecuencia de Nyquist para esta señal.

a. Frecuencia de muestreo crítica F=2B. Copie la gráfica obtenidab. Submuestreo F<2B. Copie la gráfica obtenidac. Sobremuestreo F>2B. Copie la gráfica obtenidaObserve los resultados obtenidos y explique con sus propias palabras que ha sucedido en cada uno de los casos de muestreo:



10. Ahora observe el caso del pulso rectangular. Click en Pulse En la gráfica roja puede observar que el ancho de banda de esta señal es muy amplio, ya que su espectro se ex-

tiende hasta frecuencias muy altas. Por tanto para poderla recuperar con buena exactitud tendremos que utilizar una frecuencia de muestreo bastan-te alta.

Comience utilizando una frecuencia de muestreo que considere aproximada y defina para el filtro pasabajas una frecuencia de corte igual a la mitad de la frecuencia de muestreo. Observe las gráficas obtenidas y si es necesa-rio aumente las frecuencias de muestreo y de corte del filtro, hasta obtener que la señal recuperada después del filtro le parece lo suficientemente aproximada al pulso rectangular inicial.

a. Utilice Copy-paste para copiar la figura obtenida en el applet:

b. Qué Frecuencias de muestreo y de corte utilizó?

Con este ejercicio del pulso rectangular se puede entender por qué el ancho de banda de los sistemas digitales debe ser tan amplio, ya que las señales utilizadas en estos sistemas son formas diferentes de pulsos rectangulares.

Taller Final

Documents

Transcript of Taller Final