Download - Practica 1 a PDS Terminado
-
8/12/2019 Practica 1 a PDS Terminado
1/21
Alumno: Luis Brian Alvarez Ponce
Curso: Procesamiento Digital de Seales
Practicas 01 A
III. ACTIVIDADES
Para los ejercicios siguientes:
Indicar los comandos de MatLab a utilizar Consignar las tabulaciones y las grficas obtenidas para las funciones
construidas
Grficos Bsicos 2D
1. Construir la grfica siguientex = [1 3 0.5 2.5 2];y = [0 1 2 3 4];
plot(x,y)
2. Cambiar el comandoplotpor el comandostem, que sucedi con la grafica?Cuando usamos plot (x,y) la figura es:
0.5 1 1.5 2 2.5 30
0.5
1
1.5
2
2.5
3
3.5
4
-
8/12/2019 Practica 1 a PDS Terminado
2/21
Al usar plot es porque estamos indicando una seal continua.
Al usar stem estamos indicado un muestreo es por esto que la figura cambia el resultado
es este:
3. Colquele ttulo y nombre a los ejes de la grfica, usando title, ylabel, xlabel
En Plot:
title('Usando Plot')xlabel('t imagen 1')ylabel('f(t)')
0.5 1 1.5 2 2.5 30
0.5
1
1.5
2
2.5
3
3.5
4
0.5 1 1.5 2 2.5 30
0.5
1
1.5
2
2.5
3
3.5
4
-
8/12/2019 Practica 1 a PDS Terminado
3/21
En Stem:
title('Usando Stem')xlabel('t imagen 2')ylabel('f(n)')
Impulso unitario
4. Usando el editor de MatLab crear la funcin, comente cada lneafunction[x,n]=impseq(n0,n1,n2)%Generacion x(n) = delta(n-n0); n1
-
8/12/2019 Practica 1 a PDS Terminado
4/21
-
8/12/2019 Practica 1 a PDS Terminado
5/21
7. Modifique la funcin para darle al pulso una amplitud A, donde A sera el cuartoparmetro de entrada
function[x,n]=impseq(n0,n1,n2,A) % tipo de funcin y cantidad de
parametros
%Generacion x(n) = delta(n-n0); n1 n=[-5:15];
>> x=[(n-5)==0];
>> x=(x*-4);
>> stem(n,x)
0 1 2 3 4 5 6 7 80
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
t
imagen 4
f(n)
impseq
-
8/12/2019 Practica 1 a PDS Terminado
6/21
9. Modifique las propiedades del grfico en trminos de colores y smbolo (genere almenos tres variantes.
-5 0 5 10 15-4
-3.5
-3
-2.5
-2
-1.5
-1
-0.5
0
t
imagen 5
f(n)
Agregando A
-
8/12/2019 Practica 1 a PDS Terminado
7/21
>>stem(n,x,'+g')
>> stem(n,x,'oy')
>> stem(n,x,':r')
-5 0 5 10 15-4
-3.5
-3
-2.5
-2
-1.5
-1
-0.5
0
t
imagen 6
f(n)
colores y smbolos
-5 0 5 10 15-4
-3.5
-3
-2.5
-2
-1.5
-1
-0.5
0
t
imagen 7
f(n)
colores y simbolos
-
8/12/2019 Practica 1 a PDS Terminado
8/21
10.Agregar leyenda al grfico usando el comnado legend
>> stem(n,x,'+g')
>> legent
>> stem(n,x,'oy')
>> legend
-5 0 5 10 15-4
-3.5
-3
-2.5
-2
-1.5
-1
-0.5
0
t
imagen 8
f(n)
colores y simbolos
-5 0 5 10 15-4
-3.5
-3
-2.5
-2
-1.5
-1
-0.5
0
t
imagen 6
f(n)
colores y simbolos
data1
-
8/12/2019 Practica 1 a PDS Terminado
9/21
>> stem(n,x,':r')
>> legend
-5 0 5 10 15-4
-3.5
-3
-2.5
-2
-1.5
-1
-0.5
0
t
imagen 7
f(n)
colores y simbolos
data1
-5 0 5 10 15-4
-3.5
-3
-2.5
-2
-1.5
-1
-0.5
0
t
imagen 8
f(n)
colores y simbolos
data1
-
8/12/2019 Practica 1 a PDS Terminado
10/21
Escaln unitario11.Usando el editor de MatLab crear la funcin, comente cada lnea
function[x,n]=stepseq(n0,n1,n2)%Generacion x(n) = u(n-n0); n1> n=[-2:10];
>> x=[(n-5)>=0];
>> stem(n,x,':*r');
>> title('Escaln unitario');
>> xlabel('t imagen 9');
>> ylabel('f(n)')
>> legend muestreo
-2 0 2 4 6 8 100
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1Escaln unitario
t imagen 9
f(n)
muestreo
-
8/12/2019 Practica 1 a PDS Terminado
11/21
13.Modifique la funcin para darle al escaln una amplitud A, donde A sera el cuartoparmetro de entrada
function[x,n]=stepseq(n0,n1,n2,A) % tipo de funcin y cantidad de
parametros
%Generacion x(n) = u(n-n0); n1=0 en caso contrario n< 0x=(x*A); %agregando un cuarto parametro
stem(n,x); %graficar en muestreo
14.Repita 12 para generar un escaln de amplitud -2 a partir de n=4 visualizandomuestras entre -2 y 15
>> n=[-2:15];
>> x=[(n-4)>=0];>> x=(x*-2);
>> stem(n,x,':*r')
>> title('Escaln unitario add A');>> xlabel('t imagen 9');
ylabel('f(n)')
legend muestreo
ExponencialUna secuencia exponencial responde a la expresin:
x[n] = C z^n Donde C y z son, en general, nmeros complejos.
-2 0 2 4 6 8 10 12 14 16-2
-1.8
-1.6
-1.4
-1.2
-1
-0.8
-0.6
-0.4
-0.2
0Escaln unitario add A
t imagen 9
f(n)
muestreo
-
8/12/2019 Practica 1 a PDS Terminado
12/21
15.Construya las siguientes exponenciales% Exponencial
n1=0;n2=20;n=[n1:n2];x=0.5.^n;
stem(n,x);
>> n=[0:20];
>> x=0.5.^n;
>> stem(n,x,':*r')
title('Exponencial');
xlabel('t imagen 10');
ylabel('f(n)')
legend muestreo
% Exponencialn = 0:0.1:10;
x=0.5.^n;stem(n,x);
>> n=[0:0.1:10];
>> x=0.5.^n;
>> stem(n,x,':*r')
>> title('Exponencial add A');
>> xlabel('t imagen 11');
ylabel('f(n)')
legend muestreo
0 2 4 6 8 10 12 14 16 18 200
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1Exponencial
t imagen 10
f(n)
muestreo
-
8/12/2019 Practica 1 a PDS Terminado
13/21
.
El comando n = 0:0.1:10; define el vector con componentes de rango de 0 a 10 en pasos de 0.1;La funcin x= 0.9.^n; define un vector con componentes 0.9.^(0.1), 0.9.^(0.2), 0.9.^(0.3), etc.
>> n=[0:0.1:10];>> x = 0.9.^n;>> stem(n,x,':*r')title('Exponencial add A');xlabel('t imagen 12');
>> ylabel('f(n)')legend muestreo
0 1 2 3 4 5 6 7 8 9 100
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1Exponencial add A
t
imagen 11
f(n)
muestreo
0 1 2 3 4 5 6 7 8 9 100
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1Exponencial add A
t imagen 12
f(n)
muestreo
-
8/12/2019 Practica 1 a PDS Terminado
14/21
% Exponencialn = 0:20;
x=5.^n;stem(n,x);
>> n = 0:20;>> x=5.^n>> stem(n,x,':*r')title('Exponencial');xlabel('t imagen 13');>> ylabel('f(n)')legend muestreo>> grid on
%Exponencial Compleja
n1=0;n2=10;
n=[n1:n2];x=exp((2+3j)*n);stem(n,x);
>> n=[0:10];
x=exp((2+3j)*n);
stem(n,x);
stem(n,x,':*r')
title('Exponencial');xlabel('t imagen 14');
ylabel('f(n)')
legend muestreo
grid on
0 2 4 6 8 10 12 14 16 18 200
1
2
3
4
5
6
7
8
9
10x 10
13 Exponencial
t imagen 13
f(n)
muestreo
-
8/12/2019 Practica 1 a PDS Terminado
15/21
-
8/12/2019 Practica 1 a PDS Terminado
16/21
Seal senoidal17.Graficar la siguiente seal senoidal
0 2 4 6 8 10-2
0
2
4
6
8x 10
7
Muestreo
Amplitud
Parte real
0 2 4 6 8 10-5
-4
-3
-2
-1
0
1x 10
8
Muestreo
Amplitud
Parte imaginaria
%Seal Senoidal
n1=0;n2=60;
n=[n1:n2];
x=2*sin(0.2*pi*n+pi/5);stem(n,x);
-
8/12/2019 Practica 1 a PDS Terminado
17/21
n=[0:60];x=2*sin(0.2*pi*n+pi/5);stem(n,x);stem(n,x,':*r')title('Senoidal');xlabel('t imagen 17');ylabel('f(n)')
legend muestreogrid on
18.Graficar la siguiente seal senoidal, observe que la suma tambin es peridica
n=[0:120];x=2*cos(0.2*pi*n+pi/5)+sin(0.5*pi*n);stem(n,x);stem(n,x,':*r')title('Senoidal');xlabel('t imagen 1');ylabel('f(n)')legend muestreogrid on
0 10 20 30 40 50 60-2
-1.5
-1
-0.5
0
0.5
1
1.5
2Senoidal
t imagen 17
f(n
)
muestreo
%Seal Senoidal
n1=0;n2=120;
n=[n1:n2];
x=2*cos(0.2*pi*n+pi/5)+sin(0.5*pi*n);stem(n,x);
-
8/12/2019 Practica 1 a PDS Terminado
18/21
CUESTIONARIO FINAL
1. Crear un programa en matlab donde se grafiquen las siguientes ondas: (use subplot)cos(2n/16)
cos(2n2/16)cos(2n4/16)
cos(2n8/16)
cos(2n12/16)
n1=0;n2=60;n=[n1:n2];x=cos(2*pi*n/16);subplot(5,1,1)stem(n,x,':*r')xlabel('Muestreo')ylabel('Amplitud')title('Funcion coseno')legend muestreogrid on
x=cos(2*pi*2*n/16);subplot(5,1,2)
stem(n,x,':*r')xlabel('Muestreo')ylabel('Amplitud')title('Funcion coseno')legend muestreogrid on
x=cos(2*pi*4*n/16);subplot(5,1,3)stem(n,x,':*r')xlabel('Muestreo')
ylabel('Amplitud')title('Funcion coseno')legend muestreo
0 20 40 60 80 100 120-3
-2
-1
0
1
2
3Senoidal
t imagen 1
f(n)
muestreo
-
8/12/2019 Practica 1 a PDS Terminado
19/21
grid on
x=cos(2*pi*8*n/16);subplot(5,1,4)stem(n,x,':*r')xlabel('Muestreo')
ylabel('Amplitud')title('Funcion coseno')legend muestreogrid on
x=cos(2*pi*12*n/16);subplot(5,1,5)stem(n,x,':*r')xlabel('Muestreo')ylabel('Amplitud')title('Funcion coseno')legend muestreogrid on
2. Forme y grafique el muestreo de las siguientes seales usando una razn de muestreo de l0Hz (10 muestras por segundo). Incluya las graficas en su reporte:
a. y1 = cos(6t);
0 10 20 30 40 50 60-1
0
1
Muestreo
Amplitud
Funcion coseno
0 10 20 30 40 50 60-1
0
1
Muestreo
Amplitud
Funcion coseno
0 10 20 30 40 50 60-1
0
1
Muestreo
Amplitud
Funcion coseno
0 10 20 30 40 50 60-1
0
1
Muestreo
Amplitud
Funcion coseno
0 10 20 30 40 50 60-1
0
1
Muestreo
Amplitud
Funcion coseno
muestreo
muestreo
muestreo
muestreo
muestreo
-
8/12/2019 Practica 1 a PDS Terminado
20/21
b. y2 = 3exp(-4t)*cos(5t);
c. sume un ruido escalado a 0.2 en y1. (use la funcin randn, y escriba help randn sino sabe cmo usarla).
t=[0:10];r = cos(6*t)+ 0.2.*randn(size(t));
stem(t,r,'k');
0 1 2 3 4 5 6 7 8 9 10
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
-
8/12/2019 Practica 1 a PDS Terminado
21/21
CONCLUSIONES
Emita al menos cinco conclusiones en torno al manejo de seales en MatLab