UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS

FACULTAD DE INGENIERÍA ELECTRÓNICA Y ELÉCTRICA

INFORME N°1:

TEMA: PROCESAMIENTO DE SEÑALES Y FUNCIONES

PERIODICAS USANDO MATLAB

ALUMNO CÓDIGO

 Lucas Huaman Jonathan 08190017

 LABORATORIO DE INTRODUCCION A LAS TELECOMUNICACIONES

PRACTICA Nº 1

TEMA : PROCESAMIENTO DE SEÑALES Y FUNCIONES PERIODICAS USANDO MATLAB

I. OBJETIVO:

Simular e investigar en forma experimental el procesamiento de señales usando los comandos e

instrucciones del software MATLAB, desarrollando los ejercicios propuestos en función de los
ejercicios planteados.

INTRODUCCIÓN TEÓRICA:

Clasificación de las señales:

o Continuas: toman valor en cualquier instante del tiempo.
o Discretas: toman valor en algunos instantes concretos.
o Analógicas: en un instante de tiempo dado toman un valor real
o Digitales: toman un valor de un conjunto de valores discretos

 Una señal continua en el tiempo es periódica si los valores se repiten cada T0 unidades: f (t+T0)=f
(t).
 Período fundamental o periodo T0
 Frecuencia N = 1/ T0 se mide en ciclos por unidad de
 Tiempo (Hz= ciclos por segundo) frecuencia angular, pulsación o frecuencia w0 = 2π / T0.

La aplicación de las técnicas de procesamiento digital de señales (DSP) al tratamiento de las

señales de audio se ha convertido en el conjunto de herramientas más poderosas actualmente
disponibles para la transformación y análisis del sonido, y sus aplicaciones en la música. Gracias
al desarrollo y expansión de la tecnología informática, estas herramientas son cada vez más
accesibles a compositores, músicos e ingenieros de sonido, pudiendo ser implementadas hasta
en sencillas computadoras domésticas. Hoy en día estas técnicas son ampliamente utilizadas en
todo tipo de aplicaciones, que van desde la composición electroacústica hasta la post-producción
fonográfica.
Sin embargo, la mayoría de los usuarios desconoce los fundamentos teóricos y las bases
matemáticas de los procesos involucrados, los que quedan reducidos a una caja negra, limitando
enormemente la capacidad de manejarlos de manera óptima. Los programas que implementan
estas técnicas suelen presentar una interfaz que al usuario resulta opaca, y los parámetros son
ajustados en el mejor de los casos a base de prueba y error, cuando no son simplemente
aceptados los valores que el programa adquiere por defecto. El desconocimiento de las
formulaciones matemáticas subyacentes convierte además a la "caja negra" en una especie de
 "caja mágica", de la que se ignora cuáles son sus verdaderas posibilidades reales, y cuáles sus
limitaciones intrínsecas.
El problema es que los fundamentos teóricos de las técnicas de DSP son presentados
generalmente en términos de estudios relativamente avanzados de ingeniería, y resultan por
tanto inaccesibles a la mayoría de los músicos, que no disponen de los conocimientos
matemáticos necesarios para abordar la literatura sobre el tema.

II. EQUIPOS Y MATERIALES:

1. MATLAB
2. PC Pentium II-actual.
3. Manual de MATLAB.

III. PROCEDIMIENTO:

1. Desarrolle en el programa principal de MATLAB los siguientes ejemplos y anotar sus resultados.

Funciones matemáticas: y2 = 2

>> x=sqrt(2)/2 y3 = 3

x = 0.7071 Números complejos:

>> y=asin(x) >> a=1;b=5;c=6

y = 0.7854 c= 6

>> y_deg=y*180/pi %convierte grados a >> x1=(-b+sqrt(b^2-4*a*c))/(2*a)

radianes
x1 = -2
y_deg = 45.0000
>> x2=(-b-sqrt(b^2-4*a*c))/(2*a)
>> x=sqrt(3^2+4^2)
x2 = -3
x= 5
>> a=1,b=4,c=13;
>> y=rem(23,4)
a=1
y= 3
b=4
>>
x=2.6,y1=fix(x),y2=floor(x),y3=ceil(x),y4=roun >> x1=(-b+sqrt(b^2-4*a*c))/(2*a)
d(x); x1 = -2.0000 + 3.0000i
x = 2.6000 >> x2=(-b-sqrt(b^2-4*a*c))/(2*a)
y1 = 2 x2 = -2.0000 - 3.0000i
>> a*x1^2+b*x1+c >> imag_c1=imag(c1)

ans = 0 imag_c1 = -2

>> a*x2^2+b*x2+c construccion de arrays

ans = 0 >> x=[0 .1*pi .2*pi .3*pi .4*pi .5*pi .6*pi

.7*pi .8*pi .9*pi pi]
>> C1=1-2i
x = Columns 1 through 4
C1 = 1.0000 - 2.0000i
0 0.3142 0.6283 0.9425
>> C1=1-2j
Columns 5 through 8
C1 = 1.0000 - 2.0000i
1.2566 1.5708 1.8850 2.1991
>> C2=3*(2-sqrt(-1)*3)
Columns 9 through 11
C2 = 6.0000 - 9.0000i
2.5133 2.8274 3.1416
>> c3=sqrt(-2)
>> y=sin(x)
c3 = 0 + 1.4142i
y = Columns 1 through 4
>> c4=6+sin(.5)*i
0 0.3090 0.5878 0.8090
c4 = 6.0000 + 0.4794i
Columns 5 through 8
>> c5=6+sin(.5)*j
0.9511 1.0000 0.9511 0.8090
c5 = 6.0000 + 0.4794i
Columns 9 through 11
>> C1=1-2i
0.5878 0.3090 0.0000
C1 = 1.0000 - 2.0000i
>> x(1:5),
>> mag_c1=abs(c1)
ans = Columns 1 through 4
mag_c1 = 2.2361
0 0.3142 0.6283 0.9425
>> angle_c1=angle(c1)
Column 5
angle_c1 = -1.1071
1.2566
>> deg_c1=angle_c1*180/pi
>> x=(0:0.1:1)*pi
deg_c1 = -63.4349
x = Columns 1 through 4
>> real_c1=real(c1)
0 0.3142 0.6283 0.9425
real_c1 = 1
 Columns 5 through 8 z=cos(x);

1.2566 1.5708 1.8850 2.1991 figure(2)

Columns 9 through 11 plot(x,z,x,y)

2.5133 2.8274 3.1416 figure(3)

>> x=linspace(0,pi,11) plot(x,y,x,y,'*')

x = Columns 1 through 4 figure(4)

0 0.3142 0.6283 0.9425 plot(y,z)

Columns 5 through 8 figure(5)

1.2566 1.5708 1.8850 2.1991 plot(x,y,x,2*y.*z,'--')

Columns 9 through 11 grid

2.5133 2.8274 3.1416 xlabel('variable independiente x')

Grafica de funciones trigonométricas ylabel('variables dependientes')

x=linspace(0,4*pi,60); title('2sen(x)cos(x)sen(2x)')

y=sin(x); figure(6)

figure(1) plot3(y,z,x),grid

plot(x,y)
Generación de señales discontinuas

x=linspace(0,16,100);

y=sin(x);

z=(y>=0).*y;

z=z+0.5*(y<0);

z=(x<=14).*z;

plot(x,z)

xlabel('x'),ylabel('Z=f(x)')

title('señal discontinua')

Manipulación de gráficos

x=linspace(0,4*pi,60);

y=sin(x);

z=cos(x);

plot(x,y)

hold on

plot(x,z,'m')

hold off
Subdivisión de ventanas de graficas

x=linspace(0,4*pi,60);

y=sin(x);

z=cos(x);

a=2*sin(x).*cos(x);

b=sin(x)./(cos(x)+eps);

subplot(2,2,1)

plot(x,y),axis([0 4*pi -1 1]),title('sin(x)')

subplot(2,2,2)

plot(x,z),axis([0 4*pi -1 1]),title('cos(x)')

subplot(2,2,3)

plot(x,a),axis([0 4*pi -1 1]),title('2sin(x)cos(x)')

subplot(2,2,4)

plot(x,b),axis([0 4*pi -40 40]),title('tg=sin(x)/cos(x)')

Cuestionario

1- Cambie el valor de cada uno de los ejercicios y desarrolle nuevas aplicaciones

Funciones matemáticas
>> y=rem(23,4)
>> x=sqrt(2)/3
y= 3
x = 0.4714
>>
>> y=acos(x) x=2.6,y1=fix(x),y2=floor(x),y3=ceil(x
),y4=round(x);
y = 1.0799
x = 2.6000
>> y_deg=y*180/pi%convierte
grados a radianes y1 = 2

y_deg = 61.8745 y2 = 2

>> x=sqrt(3^2+4^2)
y3 = 3
x= 5
números complejos
>> c3=sqrt(-2)
>> a=1;b=3;c=5
c3 =
c= 5
0 + 1.4142i
>> x1=(-b+sqrt(b^2-4*a*c))/(2*a)
>> c4=6+sin(.5)*i
x1 = -1.5000 + 1.6583i
c4 =
>> x2=(-b-sqrt(b^2-4*a*c))/(2*a)
6.0000 + 0.4794i
x2 = -1.5000 - 1.6583i
>> c5=6+sin(.5)*j
>> a=1,b=3,c=12;
c5 =
>> x1=(-b+sqrt(b^2-4*a*c))/(2*a)
6.0000 + 0.4794i
x1 = -1.5000 + 3.1225i
>> C1=1-2i
>> x2=(-b-sqrt(b^2-4*a*c))/(2*a)
C1 =
x2 = -1.5000 - 3.1225i
1.0000 - 2.0000i
>> a*x1^2+b*x1+c,
>> c1=1-2j
ans = 0
c1 =
>> a*x2^2+b*x2+c
1.0000 - 2.0000i
ans = 0
>> mag_c1=abs(c1)
>> C1=1-2i mag_c1 = 2.2361

C1 = 1.0000 - 2.0000i >> angle_c1=angle(c1)

>> C1=1-2j angle_c1 = -1.1071

>> deg_c1=angle_c1*180/pi
C1 = 1.0000 - 2.0000i
deg_c1 = -63.4349
>> C2=3*(2-sqrt(-1)*3)
>> real_c1=real(c1)
C2 =
real_c1 = 1
6.0000 - 9.0000i
>> imag_c1=imag(c1)
imag_c1 = -2

Construccion de arrays

>> x=[0 .1*90 .2*90 .3*90 .4*90 .5*90 .6*90 .7*90 .8*90 .9*90 90]
x=
Columns 1 through 8
0 9.0000 18.0000 27.0000 36.0000 45.0000 54.0000 63.0000
Columns 9 through 11
72.0000 81.0000 90.0000
>> y=cos(x*pi/180)
y=
Columns 1 through 8
1.0000 0.9877 0.9511 0.8910 0.8090 0.7071 0.5878 0.4540
Columns 9 through 11
0.3090 0.1564 0.0000
>> x(2)
ans =
9
>> x(4)
ans =
27
>> x(6)
ans =
45
>> x=(0:0.1:1)*90
x=
Columns 1 through 8
0 9.0000 18.0000 27.0000 36.0000 45.0000 54.0000 63.0000
Columns 9 through 11
72.0000 81.0000 90.0000
>> x=linspace(0,90,11) % Ambos arrays deben de ser similares
x=
 0 9 18 27 36 45 54 63 72 81 90

i.
x=[0 .1*10 .2*10 .3*10 .4*10 .5*10 .6*10 .7*10 .8*10 .9*10 10]
x=
0 1 2 3 4 5 6 7 8 9 10
>> y=cos(x)
y=
Columns 1 through 8
1.0000 0.5403 -0.4161 -0.9900 -0.6536 0.2837 0.9602 0.7539
Columns 9 through 11
-0.1455 -0.9111 -0.8391
>> x(2)
ans =
1
>> x(8)
ans =
7
>> x=(0:0.1:1)*10
x=
Columns 1 through 8
0 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 7.0000
Columns 9 through 11
8.0000 9.0000 10.0000
>> x=linspace(0,10,11) % Ambos arrays deben de ser similares
x=
0 1 2 3 4 5 6 7 8 9 10

1) Escribir en el block de notas o depurador del MATLAB los siguientes ejercicios y copiar en el
programa principal, para graficar las funciones periódicas:

a. Grafica de funciones trigonométricas:

x=linspace(0,8*pi,100); y=sin(x);
figure(1) figure(4)

plot(x,y) plot(y,z)

z=cos(x); figure(5)

figure(2) plot(x,y,x,2*y.*z,'--')

plot(x,z,x,y) grid

figure(3) xlabel('variable independiente x')

plot(x,y,x,y,'*') ylabel('variables dependientes')

title('2sen(x)cos(x)sen(2x)')

figure(6)

plot3(y,z,x),grid
Generaciopn de señales discontinuas

x=linspace(0,20,200);

y=sin(x);

z=(y>=0).*y;

z=z+0.5*(y<0);

z=(x<=14).*z;

plot(x,z)

xlabel('x'),ylabel('Z=f(x)')
title('señal discontinua')

Manipulacion de graficos

x=linspace(0,8*pi,100);

y=sin(x);

z=cos(x);

plot(x,y)

hold on

plot(x,z,'m')

hold off
Subdivision de ventana de graficas

x=linspace(0,8*pi,100);

y=sin(x);

z=cos(x);

a=2*sin(x).*cos(x);

b=sin(x)./(cos(x)+eps);

subplot(2,2,1)

plot(x,y),axis([0 4*pi -1 1]),title('sin(x)')

subplot(2,2,2)
plot(x,z),axis([0 4*pi -1 1]),title('cos(x)')

subplot(2,2,3)

plot(x,a),axis([0 4*pi -1 1]),title('2sin(x)cos(x)')

subplot(2,2,4)

plot(x,b),axis([0 4*pi -40 40]),title('tg=sin(x)/cos(x)')

2- Explique las funciones de cada una de la principales instrucciones y comandos

utilizados en los ejercicios anteriores.

 sqrt(x): obtiene la raíz cuadrada de la variable x.

 asin(x): obtiene el arcoseno de la variable x.
 rem(x,y): obtiene el residuo de la división x/y.
 fix(x): redondea a x hacia cero.
 floor(x): redondea a x hacia menos infinito.
 ceil(x): redondea a x hacia infinito.
 round(x): redondea a x hacia el entero más próximo.
 abs(c): obtiene la magnitud del número complejo c.
 angle(c): obtiene el ángulo de fase del número complejo c.
 real(c): obtiene la parte real del número complejo c.
 imag(c): obtiene la parte imaginaria del número complejo c.
 x= [0 .1*pi .2*pi .3*pi .4*pi .5*pi .6*pi .7*pi .8*pi .9*pi pi]: crea un arreglo
con los valores dentro de los corchetes.
 y=sin(x): obtiene el seno de cada uno de los valores de x en el arreglo
anterior.
 x(3): accede al tercer elemento del arreglo x.
 y(5): accede al quinto elemento del arreglo y.
 x= (0:0.1:1)*pi: crea un arreglo con intervalo de 0.1, con primer elemento
0, último elemento 1, y los multiplica por π.
 x=linspace (0, pi, 11): crea un arreglo con primer elemento 0, último
elemento pi, y contenga en total 11 elementos.
 plot(x,y):crea una grafica bidimensional que coloca en el eje de las
abscisas a los valores de x, y en el eje de las ordenadas a los valores
de y.
 plot(x,z,x,y): crea los plot(x,z) y plot(x,y) simultáneamente.
 plot(x,y,x,y,’*’): crea los plot(x,y) y plot(x,y) simultáneamente, este ultimo
formado de “*”.
 plot(x,2*y.*z,’--’): coloca en el eje de las ordenadas el doble producto de
los valores de y x z, y forma la grafica con “-“.
 grid: traza una cuadricula en la grafica .
 xlabel(nombre): coloca un nombre al eje x de la grafica.
 ylabel(nombre): coloca un nombre al eje y de la grafica.
 title(titulo): coloca un titulo al grafico existente.
 plot3(y,z,x): realiza un grafica en tres dimensiones con los valores de
y,z,x.
 z=(y>=0).*y: fija a cero los valores negativos de y.
 z=z+0.5*(y<0): suma 0.5 a los valores negativos de y.
 z=(x<=14).*z: fija a cero los valores de x mayores a 14.
 hold on: matlab conserva las curvas o graficasn existentes.
 hold off: libera la ventana de la figura actual para nuevas graficas.
 subplot(2,2,1): selecciona la subgrafica superior izquierda.
 subplot(2,2,2): selecciona la subgrafica superior derecha.
 subplot(2,2,3): selecciona la subgrafica inferior izquierda.
 subplot(2,2,4): selecciona la subgrafica inferior derecha.
  axis(0 4*pi -1 1): localiza la grafica para los valores del eje x entre 0 y 4π,
y los del eje y entre -1 y 1.

3 Presentar su informe final con una breve introducción teórica, con el resultado
de cada ejercicio plantado y sus modificaciones experimentales.
Desarrollado al principio.