Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

GuiaDeLaboratorio SeñalesSistemas

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 25

UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS

(Universidad del Perú, DECANA DE AMÈRICA)


FACULTAD DE INGENIERÍA ELECTRÓNICA Y ELÉCTRICA
E.A.P. INGENIERÍA DE TELECOMUNICACIONES
Av. Venezuela s/n - Lima, Perú
Teléfono: 6197000 anexo 4203 Fax: 4209
_____________

LABORATORIO 01 DE SEÑALES Y SISTEMAS


SEÑALES CONTINUAS Y DISCRETAS EN MATLAB

Profesor: Daniel Díaz Ataucuri

I. OBJETIVOS

Analizar las principales características de las señales continua y discretas utilizando


herramienta de software matemático: Matlab y/o Scilab.

.
II. MARCO TEÓRICO

Para el análisis de funciones matemáticas se puede utilizar herramientas de software propietario


como el Matlab (https://la.mathworks.com/) o software libre u open source como Scilab
(https://www.scilab.org). Un artículo donde describe las características de diferentes herramientas
de software matemático lo pueden encontrar en: http://profs.scienze.univr.it/~caliari/pdf/octave.pdf

A continuación, describiremos de manera breve los principales comandos de Matlab para iniciar
los primeros análisis de señales continuas y discretas.

Asignación de un valor a una variable

x = 2; asigna el valor 2 a la variable x pero no lo visualiza (observe que se escribe ; al final).


x = 2 asigna el valor 2 a la variable x pero lo visualiza (observe que se omite el ; al final).

Definición de una función

Si se desea definir y como función de x, se debe ingresar y = @(x) ingrese la expresión

Asignación de un vector a una variable

x = 1:2:10 asigna valores a la variable x desde 1 hasta 10 incrementando en 2.


Como resultado se obtendría x = 1, 3, 5, 7 y 9.

linspace(a,b,n); genera un vector de n puntos desde a hasta b, cuyos elementos están


espaciados linealmente. Por ejemplo, linspace(1,13,7) genera los siguientes valores 1 3 5 7
9 11 13 y linspace(2,27,6) genera los siguientes valores 2 7 12 17 22 27.

logspace(a,b,n); genera un vector de n puntos desde a hasta b, cuyos elementos están


espaciados logarítmicamente.

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


Dibujar una gráfica

plot(‘colocar variable x’, ‘colocar variable x’); crea un gráfico a partir de vectores y/o columnas
de matrices, con escalas lineales sobre ambos ejes.

stem(x,y); crea una gráfica discreta de una secuencia de datos de x versus la columna y.
Las variables x e y deben ser vectores o matrices del mismo tamaño.

title(‘ingrese el título de la gráfica’); añade un título a la gráfica.

xlabel(‘etiqueta en eje x’); añade una etiqueta al eje de abscisas. Con xlabel off desaparece.

ylabel(‘etiqueta en eje y’); añade una etiqueta al eje de ordenadas. Con ylabel off
desaparece.

text(x,y,‘ingrese mensaje de texto’); introduce un mensaje de texto en el lugar especificado


por las coordenadas x e y. Si x e y son: vectores, el texto se repite por cada par de elementos.
Si el menaje de texto ingresado es también un vector de cadenas de texto de la misma
dimensión, cada elemento se escribe en las coordenadas correspondientes.

grid on; habilita una grilla en la gráfica.

grid off; deshabilita la grilla en la gráfica.

clc; limpia la ventana de comando del Matlab (clear command window).

Ejemplo 01:

>> x = -5:5;
>> stem(x,2*x)

Se obtiene:

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


Ejemplo 02:

Si ahora ingresamos:

>> x=linspace(5,50,10);
>> stem(x,5*x)

Se obtiene:

Ejemplo 03:

Si ahora ingresamos:

>> x=linspace(-4,4,9)
>> plot(x,4*x)
>> grid on

Se obtiene

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


Ejemplo 04:

Graficar la función y = 2x – 4 para valores enteros de x = [-10,10]

Se obtiene:

Superposición de gráficas

hold on; permite superponer una gráfica en otra

hold of; deshabilita la superposición de gráficas

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


Generación de dos sub gráficas

subplot(m.n.p); se genera una matriz de m filas y n columnas, y p genera los ejes de cada
gráfica. Por ejemplo, subplot(2,1,1) generará la ventana superior y subplot(2,1,2) generará la
ventana inferior.

subplot(1,2,1) generará la ventana de la izquierda y subplot(1,2,2) generará la ventana de la


derecha.

Colorear las gráficas

plot(x,y,’z’); en este comando, el parámetro z define el color de la gráfica. Si se desea


cambiar color rojo coloque r, otros colores como y(amarillo), g(verde), w(blanco) y k(negro).
Si la curva de color rojo se delinea con + se debe digitar plot(t,m,’r+’); para delinear la curva
de otras formas se puede reemplazar + por *,o,:, etc.

Control de los ejes

axis.- Devuelve los límites del dibujo actual en un vector fila. Para gráficos de dos
dimensiones tiene los elementos [xmin xmax ymin ymax]

axis(v).- Establece la escala de los ejes conforme al vector v que será [xmin xmax ymin ymax].

axis(axis).- Bloquea la escala actual cuando se añaden dibujos posteriores a uno que se ha
mantenido con el comando hold.

axis(cad).- Establece la escala de los ejes con diferentes resultados dependiendo de la


cadena cad utilizada:

 ’auto’ Vuelve a la escala automática.


 ’equal’ Da la misma escala en ambos ejes.

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


 ’ij’ Intercambia la parte positiva y negativa del eje y.
 ’xy’ Deshace el anterior.
 ’image’ Igual que equal pero se ajusta al dibujo.
 ’square’ Modifica la ventana gráfica para hacer que la caja del dibujo sea cuadrada.
 ’normal’ Modifica la ventana gráfica para que la caja del dibujo vuelva al tamaño
habitual.
 ’off’ No se muestran los ejes.
 ’on’ Deshace el comando anterior

Ejemplo 05:

Si ahora ingresamos:

>> r = 0:pi/100:10*pi;
>> plot(r,sin(r),'g');
>> grid on
>> hold on;
>> plot(r,cos(r),'r');

Si ahora ingresamos:

>> r = 0:pi/100:10*pi;
>> plot(r,sin(r),'g');
>> grid on
>> hold on;
>> plot(r,cos(r),'r');
>> axis('equal');

Ejemplo 06:

Analizar el siguiente programa:

clear
x1 = linspace(0,2*pi,50)';
x2 = linspace(pi,3*pi,50)';
X = [x1, x2];
Y = [cos(x1), 0.5*sin(x2)];
stem(X,Y)

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


III. PROCEDIMIENTO

Utilizando Matlab se pide:

1.- Representar los puntos p1(-5,0), p2(-4,1), p3(-3,0), p4(-2,2), p5(-1,2), p6(1,-1), p7(5,-0):

b.- Como se ilustra en la figura siguiente:

c.- Como se ilustra en la figura siguiente, uniendo los puntos:

2.- Mostrar en una misma ventana las siguientes funciones (observe que la GRAFICA 01 tiene
un rango de 0 a 2π en intervalos de π/100 y la GRAFICA 02 tiene un rango de 0 a 4π en
intervalos de π/100):

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


3.- Entre los valores n=0 y n=4π, considerar 50 valores incluyendo los extremos para graficar en
ventanas separadas las siguientes funciones:

a.- exp(0.3n)

b.- sin(3n)

c.- [exp(0.3n)][sin(3n)]

4.- Del problema anterior, graficar las tres funciones en una sola ventana como se ilustra en la
figura siguiente:

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


4.- Graficar en ventanas separadas las siguientes funciones:

a.- [100exp(n)][sin(20n)] , n=0 hasta n=5, para 1000 valores en este rango. Coloque etiquetas
según su criterio.

b.- f(x) = 2x3 – 4x2 – 18x +36 , n=-5 hasta n=5, para 100 valores en este rango. Coloque
etiquetas según su criterio.

c.- En una misma ventana f1(t) = sen(t) , f2(t) = cos(t) y f3(t) = [sen(t)][cos(t)] , en [0,2π]

5.- Representar la función xsen(1/x) en el intervalo [-2,2], Muestre los detalles de esta curva en
las cercanías de (0,0). Considere unos 100 puntos en este intervalo.

IV. RESULTADOS

Imprimir los resultados obtenidos en cada uno de los pasos anteriores y presentar en el informe
final del presente laboratorio.

V. DISCUSIONES DE RESULTADOS

En cada paso anterior explique los comandos de Matlab utilizados.

VI. CONCLUSIONES

Del presente laboratorio de cuatro (04) conclusiones de manera clara.

VII. BIBLIOGRAFÍA

1.- Información oficial del Matlab, https://la.mathworks.com/help/matlab/ref/stem.html

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
(Universidad del Perú, DECANA DE AMÈRICA)
FACULTAD DE INGENIERÍA ELECTRÓNICA Y ELÉCTRICA
E.A.P. INGENIERÍA DE TELECOMUNICACIONES
Av. Venezuela s/n - Lima, Perú
Teléfono: 6197000 anexo 4203 Fax: 4209
_____________

LABORATORIO 02 DE SEÑALES Y SISTEMAS


SEÑALES BÁSICAS EN MATLAB

Profesor: Daniel Díaz Ataucuri

I. OBJETIVOS

Analizar las principales características de las señales básicas continuas y discretas utilizando
herramienta de software matemático: Matlab y/o Scilab.

.
II. MARCO TEÓRICO

Existen un conjunto de funciones básicas que son de gran importancia en el análisis de señales y
sistemas. Entre ellas están:

II.1 Función discretas impulso unitario

Llamado también como muestra unitaria, se define como:

0, n ≠ 0
δ[n] =
1, n = 0

II.2 Función discretas escalón unitario

0, n < 0
u[n] =
1, n ≥ 0

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


II.3 Función continua impulso unitario

0, t≠0
δ(t) =

∫ ∞δ(t)dt = 1

II.4 Función continua escalón unitario

0, t<0
u(t) =
1, t>0

II.5 Función discreta rampa

n, n ≥ 0
r[n] =
0, n < 0

II.6 Función continua rampa

t, t ≥ 0
r[t] =
0, t < 0

Creación de Script

Matlab le permite cread script para ejecutar los comandos que lo constituyen. Para ello ingresar
en New y luego Script. Una vez finalizado el ingreso de los comandos almacenar en una carpeta
adecuada.

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


Ejercicios:

a.- Graficar la siguiente función cuadrada

%Generación de onda Cuadrada de amplitud A, frecuencia


%fundamental w (medida en radianes por segundo) y de
%ciclo útil cu.
clear
A=1;
w=10*pi;
cu=0.5;
t=-1:0.001:1;
sq=A*square(w*t+cu);
plot(t,sq,'r','LineWidth',2); %Dibuja líneas rojas y gruesas
grid on
ylim([-2,2]) %Se define el eje vertical entre -2 y 2

b.- Graficar la siguiente función triangular

%Generación de onda Triangular de amplitud A, frecuencia


%fundamental w (medida en radianes por segundo) y ancho
%wd
clear
A=2;
w=10*pi;
wd=0;
t=-1:0.001:1;
tri=A*sawtooth(w*t+wd);
plot(t,tri,'b','LineWidth',2); %Dibuja línea azul y gruesas
grid on
ylim([-4,4]) %Se define el eje vertical entre -4 y 4

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


c.- Graficar la siguiente función escalón

clear
t=-1:0.01:1;
u=[zeros(1,100),
ones(1,101)];
plot(t,u,'r','LineWidth',2);
grid on
ylim([-2,2])

d.- Graficar la siguiente función discreta escalón

%Grafica el escalón unitario discreto.


clear
n=-10:10;
u=[zeros(1,10), ones(1,11)];
stem(n,u,'g','LineWidth',2);
grid on
ylim([-2,2])

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


e.- Graficar la siguiente función discreta impulso

%Grafica de la función impulso discreto.


clear
n=-10:10;
delta=[zeros(1,10),1 , zeros(1,10)];
stem(n,delta,'g','LineWidth',2);
grid on
ylim([-2,2])

f.- Graficar la siguiente función continua rampa

clear
t1=0:0.1:10;
rampa1=t1;
rampa=[zeros(1,101), rampa1];
t2=-10:0.1:0;
t=[t2,t1];
plot(t,rampa,'r','LineWidth',2);
grid on
ylim([-2,15])

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


III. PROCEDIMIENTO

Utilizando Matlab se pide:

1.- Representar en tres ventanas simultaneas las siguientes funciones, para -12 ≤ n ≤ 12:

δ[n] δ[n-4] δ[n+5]

2.- Representar en cuatro ventanas simultaneas las siguientes funciones, para -10 ≤ n ≤ 10

δ[n]
u[n]
u[n-1]
u[n] - u[n-1]

3.- Considerando la siguiente señal discreta x[n]

Obtenga de manera analítica las siguientes funciones:

i.- En cada caso anterior, en ventanas individuales graficar realizando las operaciones en Matlab
o SCILAB.

IV. RESULTADOS

Imprimir los resultados obtenidos en cada uno de los pasos anteriores y presentar en el informe
final del presente laboratorio.

V. DISCUSIONES DE RESULTADOS

En cada paso anterior explique los comandos de Matlab utilizados.

VI. CONCLUSIONES

Del presente laboratorio de cuatro (04) conclusiones de manera clara.

VII. BIBLIOGRAFÍA

1.- Información oficial del Matlab, https://la.mathworks.com/help/matlab/ref/stem.html

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
(Universidad del Perú, DECANA DE AMÈRICA)
FACULTAD DE INGENIERÍA ELECTRÓNICA Y ELÉCTRICA
E.A.P. INGENIERÍA DE TELECOMUNICACIONES
Av. Venezuela s/n - Lima, Perú
Teléfono: 6197000 anexo 4203 Fax: 4209
_____________

LABORATORIO 03 DE SEÑALES Y SISTEMAS


CONVOLUCIÓN

Profesor: Daniel Díaz Ataucuri

I. OBJETIVO

Analizar la convolución de señales discretas y sus propiedades conmutativas, distributiva y


asociativa utilizando Matlab y/o Scilab.

.
II. MARCO TEÓRICO

(adaptado de la Escuela Politécnica Superior de la Universidad Autónoma de Madrid-UAM)

La función de Matlab conv calcula la suma de convolución



y  n   h  m x  n  m
m 

Para calcular la suma, Matlab requiere que x  n y h  n sean secuencias de duración finita. Si asumimos

que x  n es no nula solamente en el intervalo nx  n  nx  N x  1 (siendo N x su longitud) y que

h  n es no nula solamente en el intervalo nh  n  nh  Nh  1 (siendo N h su longitud), entonces

y  n es no nula únicamente en el intervalo  nx  nh   n   nx  nh   N x  Nh  2 (siendo su

longitud N x  N h  1 ). Esto significa que conv solamente necesita calcular y  n para las N x  N h  1
muestras de este intervalo.

Si x es un vector N x -dimensional que contiene las muestras de x  n en el intervalo

nx  n  nx  N x  1 y h es un vector N h -dimensional que contiene las muestras de h  n en el

intervalo nh  n  nh  Nh  1 , entonces y  conv( x, h) devuelve en y las N x  N h  1 muestras de

y  n en el intervalo  nx  nh   n   nx  nh   N x  Nh  2 .

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


Ha de tener en cuenta que la función conv no devuelve los índices de las muestras de y  n almacenadas
en el vector y . El usuario de la función conv es el responsable de conocer cuáles son dichos índices en
función de los índices de los vectores de entrada.

III. PROCEDIMIENTO

PROPIEDADES CONMUTATIVA, ASOCIATIVA Y DISTRIBUTIVA DE LA CONVOLUCIÓN.


APLICACIÓN A SISTEMAS LINEALES E INVARIANTES (LTI)

En este ejercicio, comprobará las propiedades conmutativa, distributiva y asociativa de la convolución con
un conjunto específico de señales. Además, examinará las implicaciones de estas propiedades en la
conexión serie y paralelo de sistemas lineales e invariantes.

Los problemas de este ejercicio exploran únicamente sistemas de tiempo discreto, dado que en un
computador sólo podemos almacenar y representar un número finito de valores. Sin embargo, estas mismas
propiedades son válidas también para sistemas de tiempo continuo.

A lo largo del ejercicio emplearemos las tres señales siguientes:

 1 n0

x1  n   2 1 n  4

 0 resto

 1 n0
 n 1
 1
h1  n   
 2 n2
 0 resto

 1 n 1
 2 n2

h2  n   3 n3
 1 n4

 0 resto

1. Defina el vector de MATLAB x1 para representar la señal x1  n en el intervalo 0  n  9 y los

vectores h1 y h2 para representar h1  n y h2  n en el intervalo 0  n  4 . Defina también los


vectores nx1 , nh1 y nh2 para que sean vectores de índices (valores de n ) para las señales
correspondientes. Represente las tres señales discretas con los índices correctos empleando la función
stem y dibújelas en el recuadro siguiente.

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


NOTA: pueden ser de utilidad los ejemplos MATLAB siguientes

9 1 2 
 v = [9 1 2; 5 6 8] define la matriz 5 6 8  en la variable v.
 
 n = 2 : 0.5 : 4 define el vector [2.0 2.5 3.0 3.5 4.0] en la variable n

2. La propiedad conmutativa establece que el resultado de la convolución es el mismo


independientemente del orden de los operandos. Esto implica que la salida de un sistema LTI con

respuesta al impulso h  n cuando la entrada es x  n es igual que la salida de un sistema LTI con

respuesta al impulso x  n cuando la entrada es h  n . Use la función conv para verificar esta
propiedad con los vectores h1 y x1 . Represente ambas salidas en el recuadro siguiente.
 ¿Es la salida de la función conv la misma independientemente del orden de los argumentos
de entrada?
 ¿Cuál es el intervalo de valores de n que debe utilizar para la señal de salida?

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


Gráficas y comentarios:

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


3. La convolución también tiene la propiedad distributiva respecto de la suma. Esto quiere decir que:

x  n    h1  n   h2  n   x  n   h1  n   x  n   h2  n 
Esto implica que la salida de dos sistemas LTI conectados en paralelo es la misma que la salida de un
sistema cuya respuesta al impulso es la suma de las respuestas al impulso de los sistemas conectados
en paralelo. Dibuje en el recuadro siguiente dos diagramas de bloques, uno con dos sistemas
conectados en paralelo y otro con el sistema único equivalente.

Compruebe la propiedad distributiva empleando los vectores x1 , h1 y h2 . Para ello calcule la suma
de las salidas de los dos sistemas LTI con respuestas al impulso h1  n y h2  n cuando x1  n es la
señal de entrada. Posteriormente, calcule la salida de un sistema LTI cuya respuesta al impulso es la

suma de h1  n y h2  n cuando x1  n es la señal de entrada. Represente ambas salidas en el


recuadro adjunto y compárelas.
 ¿Producen estos dos métodos la misma salida?

Gráficas, diagramas de bloques y comentarios:

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


4. La convolución también tiene la propiedad asociativa, es decir:

x  n    h1  n   h2  n    x  n   h1  n   h2  n   x  n   h1  n   h2  n 
Esta propiedad implica que el resultado de procesar una señal con una serie de sistemas LTI
conectados en cascada es equivalente a procesar la señal con un único sistema LTI cuya respuesta al
impulso es la convolución de todas las respuestas al impulso de los sistemas LTI conectados en
cascada.
Dibuje en el espacio siguiente dos diagramas de bloques que ejemplifiquen esta propiedad para dos
sistemas conectados en cascada.

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


Diagramas de bloques

Siga los siguientes pasos para comprobar la propiedad asociativa empleando los vectores x1 , h1 y
h2 .
 Calcule w n como la salida de un sistema LTI con respuesta al impulso h1  n cuando la

entrada al sistema es x1  n . Represente esta señal en el recuadro siguiente.

 Calcule la salida y1  n del sistema conectado en cascada como la salida de un sistema LTI con

respuesta al impulso h2  n cuando la entrada al sistema es w n . Represente esta señal en el


recuadro siguiente.
 Calcule ahora la respuesta al impulso del sistema equivalente a la conexión en cascada de los dos

sistemas h1  n y h2  n . Denomine hserie  n a dicha respuesta al impulso y represéntela en el


recuadro siguiente.

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


Grafica de las respuestas w n , y1  n y hserie  n

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


 Por último, calcule la salida y2  n del sistema equivalente hserie  n cuando la entrada al sistema

es x1  n . Represente esta señal en el recuadro siguiente.

Grafica de las respuestas y2  n y hserie  n

Compare las dos salidas calculadas y1  n y y2  n .

 ¿Obtuvo los mismos resultados al procesar x1  n con los dos sistemas individuales y con
el sistema equivalente?.

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM


UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
(Universidad del Perú, DECANA DE AMÈRICA)
FACULTAD DE INGENIERÍA ELECTRÓNICA Y ELÉCTRICA
E.A.P. INGENIERÍA DE TELECOMUNICACIONES
Av. Venezuela s/n - Lima, Perú
Teléfono: 6197000 anexo 4203 Fax: 4209
_____________

LABORATORIO 04 DE SEÑALES Y SISTEMAS


SERIES DE FOURIER

Profesor: Daniel Díaz Ataucuri

I. OBJETIVO

Analizar la convolución de señales discretas y sus propiedades conmutativas, distributiva y


asociativa utilizando Matlab y/o Scilab.

.
II. MARCO TEÓRICO

Elaborado por Daniel Díaz Ataucuri Profesor Titular de Telecomunicaciones UNMSM

También podría gustarte