Ingmate03t4 Lab

Asignatura Datos del alumno Fecha
Apellidos:
Procesamiento de
Señales, Sonido e
Imágenes Digitales Nombre:
Actividades
Laboratorio: Representación de señales y convolución
Objetivo
El objetivo de este laboratorio es familiarizarnos con la representación de señales en

Octave, estudiar las propiedades de periodicidad de señales continuas y discretas, y
familiarizarnos con la suma de convolución. Antes de acudir al laboratorio deberás haber
instalado Octave en tu ordenador. No está permitido el uso de estructuras de control (if,
for, etc).
Tarea 1: representación de señales discretas
Para representar señales discretas en Octave debemos definir un vector (habitualmente

llamado n) con los valores de la variable independiente espaciados una unidad y otro
vector (habitualmente llamado x) con los valores de la variable dependiente. Para dibujar
la señal discreta se usa stem. Por ejemplo:
n = [-3:3];
x = e.^n;
stem(n,x);
En Octave los operadores producto, división y exponenciación (*, /, ^) hacen operaciones

a nivel de matriz. Si queremos que la operación sea elemento a elemento debemos usar
.*, ./ y .^.
Use Octave para representar la siguiente señal discreta en el intervalo n = [-4,4]:
2, 𝑛 = −2
−2, 𝑛=1
𝑥[𝑛] = 3, 𝑛=0
1, 𝑛=1
{0, resto
TEMA 4 – Actividades © Universidad Internacional de La Rioja. (UNIR)

Apellidos:
Procesamiento de
Señales, Sonido e
Entrega
Entrega la solución en un fichero tarea1.m y la gráfica resultante en tarea1.png.
Tarea 2: representación de señales continuas y complejas
Una señal continua se representa en Octave mediante vectores que contengan valores
muestreados de la señal en intervalos cortos. Para ello, habitualmente se define primero
el intervalo de muestreo ts (inct), habitualmente menor que la unidad, y después sus
variables dependientes e independientes.
Para dibujar la señal continua se usa el comando plot. Por ejemplo:
inct = 0.1;
t = [-2:inct:3];
x = 0.5 * (sign(t) + 1);
plot(t,x);
En ocasiones el rango de la coordenada horizontal (abscisas) y vertical (coordenadas) no

coincide con el que deseamos representar. En este caso podemos usar el comando axis
para modificar estos rangos. Por ejemplo, para que el eje de abscisas sea [-2..3] y el de
coordenadas [-1..1] hacemos:
axis([-2 3 -1 1])
También podemos etiquetar la gráfica función usando title, xlabel e ylabel. Por ejemplo:
title('x = escalon(t)');
xlabel('t (seg)');
ylabel('x(t)');
En Octave las muestras del eje de coordenadas (vertical) pueden ser números complejos,
en cuyo caso las operaciones stem y plot por defecto muestran la parte real. Los

Apellidos:
Procesamiento de
Señales, Sonido e
comandos real, imag, abs y angle permiten calcular la parte real, imaginaria, módulo y
argumento de la función compleja discreta o continua que estemos representando.
Otro aspecto importante que se empleará en esta tarea es la representación de varias

funciones en una misma gráfica. Una vez que dibujamos una primera señal en una
gráfica, para que al dibujar la segunda no se borre la primera debemos de usar el
comando hold on. Además, el comando legend permite etiquetar cada una de las señales
dibujando una pequeña leyenda dentro de la gráfica. A continuación, se presenta un
ejemplo de cómo se usan los comandos hold on y legend representando una exponencial
y una sinusoidal solapadas en la misma gráfica:
n = -4:4;
x1 = exp(2.0/16.0*pi*n);
x2 = sin(2.0/16.0*pi*n);
hold off;
stem(n,x1);
hold on;
stem(n,x2,'Color','red', 'MarkerEdgeColor', 'red'); % El primer parámetro de color es el
color del segmento de la muestra y el segundo es el color del círculo de la muestra.
legend('exp(2*pi*n/16)', 'sin(2*pi*n/16)')
Dada la siguiente señal de tiempo continuo en el intervalo t=[-8..8]:
𝜋
𝑥(𝑡) = 𝑒 𝑗 4 𝑡
Usa subplot para dividir la figura en 1 fila x 2 columnas y dibujar:
En la primera gráfica (1,1) la parte real e imaginaria de la señal solapadas.

En la segunda gráfica (1,2) el valor absoluto y la fase de la señal solapadas.
Para ello podemos usar las funciones Octave: real, imag, abs y angle.

Apellidos:
Procesamiento de
Señales, Sonido e
Entrega
Entrega la solución en un fichero tarea2.m y la gráfica resultante en tarea2.png.
Tarea 3: representación de señales periódicas
En la tarea anterior hemos visto que se podía parametrizar el color de las señales
representadas. Por ejemplo, para representar x(t)=sen(πt) en rojo, haríamos:
t= [-2*pi:0.1:2*pi];
x = sin(pi*t);
plot(t,x,'Color','red');
Utiliza Octave para representar, en 1 columna y 2 filas con subplot, las siguientes
señales en el rango de la variable independiente -10 a 10 con la señal discreta en color
por defecto y la señal continua en rojo:
14𝜋 14𝜋
𝑥1 (𝑡) = 𝑅𝑒 {𝑒 𝑗 16 𝑡 } 𝑥1 [𝑛] = 𝑅𝑒 {𝑒 𝑗 16 𝑛 }
1 1
𝑥2 (𝑡) = cos ( 𝑡) 𝑥2 [𝑛] = cos ( 𝑛)
4 4
En vista de los resultados obtenidos indica qué señales de las anteriores son periódicas
o aperiódicas. En las periódicas indique su periodo fundamental.
Entrega
Entrega la solución en un fichero tarea3.m, la gráfica resultante en tarea3.png y las

respuestas a las preguntas en Actividad1_blanco.docx (no copiar el enunciado, solo
poner las respuestas).

Apellidos:
Procesamiento de
Señales, Sonido e
Tarea 4: convolución de señales
Dadas la siguiente señal y respuesta al impulso:
𝑛 + 1, 0 ≤ 𝑛 ≤ 4 1, 0 ≤ 𝑛 ≤ 2
𝑥[𝑛] = { ℎ[𝑛] = {
0, resto 0, resto
Representa estas señales con subplot en n [0,10] usando dos filas y una columna.
Vamos a calcular la convolución y[n]=x[n]*h[n]. Para ello, observando que x[n] tiene
valor solo un en un rango finito [0,4], podemos escribir la suma de convolución como
suma de cinco términos:
𝑦[𝑛] = ∑ 𝑥[𝑚]ℎ[𝑛 − 𝑚] = 𝑥[0]ℎ[𝑛 − 0] + 𝑥[1]ℎ[𝑛 − 1]+. . . +𝑥[4]ℎ[𝑛 − 4]

𝑚=0
Donde:
𝑦𝑚 [𝑛] = 𝑥[𝑚]ℎ[𝑛 − 𝑚]
Representa lo que ocurre con el valor en la posición x[m] (que corresponde a un escalar)
al pasarlo por un sistema con respuesta al impulso desplazada h[n-m] (que corresponde
a un vector). Lógicamente ym[n] también será un vector.
Para calcular h[n-m] deberás implementar una función en Octave que posteriormente
llamarás en el programa principal. Dicha función se almacenará como un fichero .m.
Para ello el fichero debe comenzar con la palabra reservada function y el nombre del
fichero coincidir con el de la función. Por ejemplo, podemos crear un fichero
real_imag.m con este contenido:
function [Xr,Xi] = real_imag(X)

Xr = real(X);
Xi = imag(X);
end

Apellidos:
Procesamiento de
Señales, Sonido e
Las funciones en Octave pueden retornar más de un valor (en este ejemplo, dos valores)
y debemos asignar un valor a cada valor retornado antes de que acabe la función.
Posteriormente, para ejecutar la función, suponiendo que en X1 tengamos representada
una señal, obtendríamos su parte real (X1r) e imaginaria (X1i) haciendo:
[X1r,X1i] = real_imag(X1);
Calcula h[n-m] implementando la siguiente función en Octave que posteriormente

llamarás en el programa principal:
function [y] = desplaza(x,n)

y = [zeros(1,n) x(1:end-n)];
end
Después calcula las sumas parciales ym[n] en los vectores y0, y1, ... y4 y represéntalos
con subplot. Por ejemplo, para calcular y2 haríamos:
y2 = x(2+1).*desplaza(h,2);
Calcula la convolución y[n] en el vector y, sumando los vectores anteriores que

contienen las sumas parciales:
𝑦[𝑛] = ∑ 𝑦𝑚 [𝑛]
𝑚
Representa con subplot en 6 filas y una columna (o en 6 columnas y una fila,

escogiendo la distribución que mejor se vea) la suma total y junto con las sumas
parciales anteriores. Indica en Actividad1_blanco.docx cuál es en general la posición
de comienzo y duración de una convolución en función de la posición de inicio y
duración de cada una de las señales convolucionadas. Si esta no es la posición de inicio
o duración de y[n] justifica por qué es así.
Calcula la convolución z[n]=x[n]*h[n] usando el comando conv de Octave y

represéntela con stem. Deberías obtener el mismo valor que obtuviste en y[n], en caso

Apellidos:
Procesamiento de
Señales, Sonido e
contrario revisa los apartados anteriores. Ten en cuenta que conv no define los índices
de la variable independiente, luego es nuestra responsabilidad determinar los índices
de los valores obtenidos.
Entrega
Entrega la solución en un fichero tarea4.m, la gráfica resultante en tarea4.png y las

respuestas a las preguntas en Actividad1_blanco.docx (no copiar el enunciado, solo
poner las respuestas).
Rúbrica
Tarea 1 (2.5 puntos)

Gráfica correcta: 1.25 puntos.
Código correcto: 1.25 puntos.


Cálculo de los periodos de las señales: 0.5 puntos (0.5/4 cada uno).
Pregunta de periodicidad: 0.5 puntos.

Gráfica correcta de representación de x[n] y h[n]: 0.25 puntos.
Código correcto de representación de x[n] y h[n]: 0.25 puntos.
Gráfica correcta de representación de sumas parciales: 0.5 puntos.
Código correcto de representación de sumas parciales: 0.5 puntos.

Apellidos:
Procesamiento de
Señales, Sonido e
Gráfica correcta de representación de y[n] con el comando conv: 0.25 puntos.

Código correcto de representación de y[n] con el comando conv: 0.25 puntos.
Pregunta de posición de comienzo de la convolución: 0.25 puntos.
Pregunta de duración de la convolución: 0.25 puntos.

Ingmate03t4 Lab

Cargado por

Copyright:

Formatos disponibles

Ingmate03t4 Lab

Cargado por

Información del documento

Descripción original:

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

Ingmate03t4 Lab

Cargado por

Copyright:

Formatos disponibles

Asignatura Datos del alumno Fecha

Laboratorio: Representación de señales y convolución

El objetivo de este laboratorio es familiarizarnos con la representación de señales en

Tarea 1: representación de señales discretas

Para representar señales discretas en Octave debemos definir un vector (habitualmente

En Octave los operadores producto, división y exponenciación (*, /, ^) hacen operaciones

Use Octave para representar la siguiente señal discreta en el intervalo n = [-4,4]:

TEMA 4 – Actividades © Universidad Internacional de La Rioja. (UNIR)

Entrega la solución en un fichero tarea1.m y la gráfica resultante en tarea1.png.

Tarea 2: representación de señales continuas y complejas

Para dibujar la señal continua se usa el comando plot. Por ejemplo:

En ocasiones el rango de la coordenada horizontal (abscisas) y vertical (coordenadas) no

TEMA 4 – Actividades © Universidad Internacional de La Rioja. (UNIR)

Otro aspecto importante que se empleará en esta tarea es la representación de varias

Dada la siguiente señal de tiempo continuo en el intervalo t=[-8..8]:

Usa subplot para dividir la figura en 1 fila x 2 columnas y dibujar:

En la primera gráfica (1,1) la parte real e imaginaria de la señal solapadas.

TEMA 4 – Actividades © Universidad Internacional de La Rioja. (UNIR)

Entrega la solución en un fichero tarea2.m y la gráfica resultante en tarea2.png.

Tarea 3: representación de señales periódicas

Entrega la solución en un fichero tarea3.m, la gráfica resultante en tarea3.png y las

TEMA 4 – Actividades © Universidad Internacional de La Rioja. (UNIR)

Tarea 4: convolución de señales

Dadas la siguiente señal y respuesta al impulso:

𝑦[𝑛] = ∑ 𝑥[𝑚]ℎ[𝑛 − 𝑚] = 𝑥[0]ℎ[𝑛 − 0] + 𝑥[1]ℎ[𝑛 − 1]+. . . +𝑥[4]ℎ[𝑛 − 4]

function [Xr,Xi] = real_imag(X)

TEMA 4 – Actividades © Universidad Internacional de La Rioja. (UNIR)

Calcula h[n-m] implementando la siguiente función en Octave que posteriormente

function [y] = desplaza(x,n)

Calcula la convolución y[n] en el vector y, sumando los vectores anteriores que

Representa con subplot en 6 filas y una columna (o en 6 columnas y una fila,

Calcula la convolución z[n]=x[n]*h[n] usando el comando conv de Octave y

TEMA 4 – Actividades © Universidad Internacional de La Rioja. (UNIR)

Entrega la solución en un fichero tarea4.m, la gráfica resultante en tarea4.png y las

Tarea 1 (2.5 puntos)

Tarea 2 (2.5 puntos)

Tarea 3 (2.5 puntos)

Tarea 4 (2.5 puntos)

TEMA 4 – Actividades © Universidad Internacional de La Rioja. (UNIR)

Gráfica correcta de representación de y[n] con el comando conv: 0.25 puntos.

TEMA 4 – Actividades © Universidad Internacional de La Rioja. (UNIR)

También podría gustarte