Capitulo 3 Ejercicios

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 54

CARATULA UNIVERSIDAD NACIONAL DE LOJA

AREA DE LA ENERGÍA Y LOS RECURSOS NATURALES NO RENOVABLES

INGENIERIA EN ELECTRÓNICA Y TELECOMUNICACIONES

Procesamiento digital de
señales

Ejercicios cap. #3

ING. Diego Orellana

Estudiante: Andrea Stefania Enriquez


Gunsha
Ciclo: 6to “A”

Fecha de Entrega Planificada: 17-12-2020.

Loja – Ecuador
G 2020-2021
1. Utilizando el método de multiplicación matriz-vector analizado en este capítulo,
escriba una función MATLAB para calcular la DTFT de una secuencia de duración finita.
El formato de la función debe ser

Utilice esta función para calcular la DTFT X (ejω) de las siguientes secuencias de
duración finita sobre −π ≤ ω ≤ π. Trace gráficas de magnitud y ángulo DTFT en una
ventana de figura.
1. Utilice esta función para calcular la DTFT X (ejω) de las siguientes
secuencias de duración finita sobre −π ≤ ω ≤ π. Trace gráficas de magnitud
y ángulo DTFT en una ventana de figura.
2. x (n) = n (0,9) n [u (n) - u (n - 21)].

3. x (n) = [cos (0.5πn) + j sin (0.5πn)] [u (n) - u (n - 51)]. Comente sobre la


gráfica de magnitud.
4. x (n) = {4, 3, 2, 1, 1, 2, 3, 4}. Comente sobre la gráfica de ángulos.

5. x (n) = {4, 3, 2, 1, −1, −2, −3, −4}. Comente sobre la gráfica de ángulos.
2. Sea

. Se forma una nueva secuencia 𝑥2 (𝑛) usando

1. Exprese X2 (ejω) en términos de X1 (ejω) sin calcular explícitamente X1


(ejω).

2. Verifique su resultado usando MATLAB calculando y trazando magnitudes


de la respectivas DTFT.
3. Determine analíticamente la DTFT de cada una de las siguientes secuencias. Grafique la
magnitud y el ángulo de X (ejω) sobre 0 ≤ ω ≤ π.
1. x(n) = 2(0.5)n u(n + 2)
2. x(n) = (0.6)|n| [u(n + 10) − u(n − 11)]
3. x(n) = n (0.9)n u(n + 3)
4. x(n) = (n + 3) (0.8)n−1 u(n − 2)
5. x(n) = 4(−0.7)n cos(0.25πn)u(n)
Las siguientes secuencias de duración finita se denominan ventanas y son muy útiles en DSP:

4. Para cada una de estas ventanas, determine sus DTFT para M = 10, 25, 50, 101. Escala
transformar valores para que el valor máximo sea igual a 1. Grafique la magnitud de la
DTFT normalizada sobre −π ≤ ω ≤ π. Estudie estas tramas y comente su
comportamiento como una función de M.
Ventana rectangular: RM (n) = u (n) - u (n - M)
|𝑀−1−2𝑛|
Triangular Window: 𝑇𝑀 (𝑛) = [1 − ] 𝑅𝑀 (𝑛)
𝑀−1

2𝜋𝑛
Ventana de Hann: 𝐶𝑀 (𝑛) = 0.5 [1 − 𝑐𝑜𝑠 ] 𝑅𝑀 (𝑛)
𝑀−1
2𝜋𝑛
Ventana de Hamming: 𝐻𝑀 (𝑛) = [0.54 − 0.46𝑐𝑜𝑠 ] 𝑅𝑀 (𝑛)
𝑀−1
5. Utilizando la definición de DTFT en (3.1), determine las secuencias correspondientes a
las siguientes DTFT.
1. X (ejω) = 3 + 2 cos (ω) + 4 cos (2ω)

2. X (ejω) = [1 - 6 cos (3ω) + 8 cos (5ω)] e−j3ω


3. X (ejω) = 2 + j4 sin (2ω) - 5 cos (4ω)

4. X (ejω) = [1 + 2 cos (ω) + 3 cos (2ω)] cos (ω / 2) e−j5ω / 2

5. X (ejω) = j [3 + 2 cos (ω) + 4 cos (2ω)] sin (ω) e −j3ω

6. Usando la definición de DTFT inversa en (3.2), determine las secuencias


correspondientes a las siguientes DTFT.
Recuerde que las transformadas anteriores son periódicas en ω con un período igual a
2π. Por lo tanto, las funciones se dan solo durante el período primario de −𝜋 ≤ 𝜔 ≤
𝜋.
7. Una secuencia de valores complejos x (n) se puede descomponer en una parte simétrica
conjugada xe (n) y una parte antisimétrica conjugada xo (n) como se discutió en el
Capítulo 2. Demuestre que F [xe (n)] = XR (ejω) y F [xo (n)] = jXI (ejω) donde XR (ejω) y
XI (ejω) son las partes real e imaginaria de la DTFT X (ejω), respectivamente. Verifique
esta propiedad en x (n) = 2 (0.9) −n [cos (0.1πn) + j sin (0.9πn)] [u (n) - u (n - 10)] usando
las funciones de MATLAB desarrolladas en el Capítulo 2.
8. Una DTFT X de valor complejo (𝒆𝒋𝝎 ) también se puede descomponer en su simétrica
conjugada parte 𝑿 𝒆 (𝒆𝒋𝝎 )) y la parte antisimétrica conjugada Xo (𝒆𝒋𝝎 ) es decir,

dónde

Muestra esa

donde 𝒙𝑹 (𝒏) y 𝒙𝑰 (𝒏) son las partes real e imaginaria de 𝒙(𝒏). Verifique esta propiedad en

utilizando las funciones de MATLAB desarrolladas en el Capítulo 2.


9. Usando la propiedad de desplazamiento de frecuencia de la DTFT, demuestre que la
parte real de X (ejω) de un pulso sinusoidal

donde RM (n) es el pulso rectangular dado en el problema P3 .4 viene dado por

Calcule y grafique XR (ejω) para ωo = π / 2 y M = 5, 15, 25, 100. Utilice el intervalo de


trazado [−π, π]. Comente sus resultados.
10. Sea 𝒙 (𝒏) = 𝓣𝟏𝟎 (−𝒏) un pulso triangular dado en el problema P3.4. Utilizando las
propiedades de la DTFT, determine y grafique la DTFT de las siguientes secuencias.
11. Para cada uno de los sistemas lineales invariantes de desplazamiento descritos por la
respuesta al impulso, determine la función de respuesta de frecuencia H (ejω). Grafique
la respuesta de magnitud | H (ejω) | y la respuesta de fase H (ejω) en el intervalo [−π,
π].
1. h(n) = (0.9)|n|
2. h(n) = sinc(0.2n)[u(n + 20) − u(n − 20)], Donde sinc(0) = 1
3. h(n) = sinc(0.2n)[u(n) − u(n − 40)]

4. h(n) = [(0.5)n + (0.4)n]u(n)


5. h(n) = (0.5)|n| cos(0.1πn)
12. Sea 𝒙(𝒏) = 𝑨 𝒄𝒐𝒔 (𝝎𝟎𝒏 + 𝜽𝟎 ) una secuencia de entrada a un sistema LTI descrito por
la respuesta al impulso 𝒉(𝒏). Demuestre que la secuencia de salida 𝒚(𝒏) está dada por

13. Sea x (n) = 3 cos (0.5πn + 60 °) + 2 sin (0.3πn) la entrada a cada uno de los sistemas
descritos en el problema P3.11. En cada caso, determine la secuencia de salida y (n).
1. La entrada al sistema h (n) = (0.9) | n | es x (n) = 3 cos (0.5πn + 60◦) +2 sin (0.3πn).
2. La respuesta de estado estable y (n) se calcula usando Matlab.
14. Sea 𝒙(𝒏) = 𝑨 𝒄𝒐𝒔 (𝝎𝟎𝒏 + 𝜽𝟎 ) una secuencia de entrada a un sistema LTI descrito por
la respuesta al impulso 𝒉 (𝒏). Entonces la secuencia de salida 𝒚 (𝒏) viene dada por

donde 𝝎𝒄 se llama frecuencia de corte y 𝜶 se llama retardo de fase.

1. Determine la respuesta al impulso ideal 𝒉𝒅 (𝒏) utilizando la relación IDTFT (3.2).

2. Determine y grafique la respuesta al impulso truncada


3. Determine y grafique la función de respuesta de frecuencia 𝑯(𝒆𝒋𝝎 ), y compárela con
la respuesta ideal del filtro de paso bajo 𝑯𝒅 (𝒆𝒋𝝎 ). Comente sus observaciones.
15. Un filtro de paso alto ideal se describe en el dominio de frecuencia por

donde ωc se llama frecuencia de corte y α se llama retardo de fase.

1. Determine la respuesta al impulso ideal hd (n) utilizando la relación IDTFT (3.2).

2. Determine y grafique la respuesta al impulso truncada

para N = 31, α = 15 y ωc = 0.5π.


3. Determine y grafique la función de respuesta de frecuencia H (e jω), y compárela con la
respuesta ideal del filtro de paso alto Hd (ejω). Comente sus observaciones.
Observamos que la respuesta truncada es una versión manchada o borrosa de la respuesta
ideal.

16. Para un sistema lineal invariante de desplazamiento descrito por la ecuación en


diferencias

la función de respuesta en frecuencia está dada por

Escriba una función de MATLAB freqresp para implementar esta relación. El formato de esta
función debe ser

Solución:
17. Determine H (ejω) y grafique su magnitud y fase para cada uno de los siguientes
sistemas.
1.

2. y(n) = x(n) − x(n − 2) + 0.95y(n − 1) − 0.9025y(n − 2)


3. y(n) = x(n) − x(n − 1) + x(n − 2) + 0.95y(n − 1) − 0.9025y(n − 2)
4. y(n) = x(n) − 1.7678x(n − 1) + 1.5625x(n − 2) + 1.1314y(n − 1) − 0.64y(n − 2)

5.
18. Un sistema lineal invariante de desplazamiento se describe mediante la ecuación en
diferencias

En cada caso, genere x (n), 0 ≤ n ≤ 200 y procese a través de la función de filtro para obtener
y (n). Compare su y (n) con las respuestas de estado estable en cada caso.

Determine la respuesta de estado estable del sistema a las siguientes entradas


19. Se muestrea una señal analógica xa (t) = sin (1000πt) utilizando los siguientes intervalos
de muestreo. En cada caso, trace el espectro de la señal de tiempo discreto resultante.

20. Implementamos el siguiente filtro analógico usando un filtro discreto:


La frecuencia de muestreo en A/D y D/A es 8000 sam/seg, y la respuesta al impulso es

1. ¿Cuál es la frecuencia digital en 𝒙(𝒏) si 𝒙𝒂 (𝒕) = 𝟏𝟎𝐜𝐨 𝐬(𝟏𝟎, 𝟎𝟎𝟎𝝅𝒕) ?

2. Determine la salida en estado estable 𝒚𝒂 (𝒕) si 𝒙𝒂 (𝒕) = 𝟏𝟎 𝒄𝒐𝒔 (𝟏𝟎, 𝟎𝟎𝟎𝝅𝒕)

3. Determine la salida de estado estable 𝒚𝒂 (𝒕) si 𝒙𝒂 (𝒕) = 𝟓 𝒔𝒊𝒏 (𝟖, 𝟎𝟎𝟎𝝅𝒕).

4. Para evitar el aliasing, se necesitaría un prefiltros para procesar 𝒙𝒂 (𝒕)antes de que


pase al convertidor A / D. ¿Qué tipo de filtro debería usarse y cuál debería ser la
frecuencia de corte más grande que funcionaría para la configuración dada?
El prefiltro debe ser un filtro de paso bajo con una frecuencia de corte de 4 KHz.

21. Considere una señal analógica xa (t) = cos (20πt), 0 ≤ t ≤ 1. Se muestrea en intervalos Ts
= 0.01, 0.05 y 0.1 seg para obtener x (n).

1. Para cada Ts, grafique x (n).


2. Reconstruya la señal analógica ya (t) a partir de las muestras x (n) utilizando la
interpolación sinc (use Δt = 0.001) y determine la frecuencia en ya (t) a partir de su
gráfico. (Ignore los efectos finales).
3. Reconstruya la señal analógica ya (t) a partir de las muestras x (n) utilizando la
interpolación spline cúbica y determine la frecuencia en ya (t) a partir de su gráfico.
(Nuevamente, ignore los efectos finales).
4. Comente sus resultados.

A partir de las gráficas de las figuras, está claro que las reconstrucciones de muestras en
Ts = 0.01 y 0.05 representan la frecuencia original (excluyendo los efectos finales) pero
las reconstrucciones para Ts = 0.1 muestran la frecuencia original con un alias a cero.
Además, la interpolación spline cúbica es una mejor reconstrucción que la interpolación
sinc, es decir, la interpolación sinc es más susceptible al efecto de contorno.

22. Considere la señal analógica 𝑥𝑎 (𝑡) = 𝑐𝑜𝑠 (20𝜋𝑡 + 𝜃), 0 ≤ t ≤ 1. Se muestrea a intervalos
𝑇𝑠 = 0.05 seg para obtener x (n). Deje θ = 0, π / 6, π / 4, π / 3, π / 2. Para cada uno de estos
valores de θ, realice lo siguiente.
1. Grafique 𝑥𝑎 (𝑡) y superponga 𝑥 (𝑛) usando la función plot (n, x, ’o’).
2. Reconstruya la señal analógica 𝑦𝑎 (𝑡) a partir de las muestras 𝑥(𝑛) usando la interpolación
sinc (Use ∆t = 0.001) y superponga 𝑥(𝑛) sobre ella.
3. Reconstruya la señal analógica 𝑦𝑎 (𝑡) a partir de las muestras 𝑥(𝑛) utilizando la interpolación
spline cúbica y superponga 𝑥(𝑛) sobre ella.
4. Debe observar que la reconstrucción resultante en cada caso tiene la frecuencia correcta pero
una amplitud diferente. Explique esta observación. Comente sobre el papel de la fase de 𝑥𝑎 (𝑡)
en el muestreo y reconstrucción de señales.

También podría gustarte