Apéndice

PRÁCTICA #1:
INTRODUCCIÓN A
MATLAB Y
MANIPULACIÓN DE
SINUSOIDES

Integrantes:
 Baque Borbor Kevin Daniel
 Jimmy Rodríguez Gonzales

PARALELO 102
Fecha de Entrega: 24 de junio del 2020
Práctica #1: Introducción a MATLAB y
Manipulación de Sinusoides.
INTRODUCCIÓN
En relación con lo visto previamente en el componente teórico de la materia, que corresponde a la
“Introducción a Señales y Sistemas en tiempo Continuo y Discreto” se realizó la siguiente práctica de
laboratorio con los objetivos a continuación:

 Introducir los aspectos básicos del trabajo en el entorno MATLAB, incluido el uso del sistema
de ayuda, la sintaxis básica de comandos, los números complejos, la indexación de matrices,
el trazado y el uso de sectorización para evitar bucles ineficientes.
 Usar MATLAB para la manipulación de señales y el trazado de su respectiva gráfica.
 Graficar en MATLAB señales sinusoidales a través de amplitud compleja.

Cumpliéndose estos objetivos con éxito, alcanzando destrezas que permitieron culminar la práctica y
posterior realización de este informe con el desempeño deseado.

MARCO TEÓRICO
Los conceptos teóricos que se describen a continuación fueron la base para poder realizar esta
práctica. Para empezar, se toman las definiciones del texto guía[ CITATION 1 \l 12298 ], donde
definimos primero a un SISTEMA, que por lo general su definición se presta para ambigüedades
dependiendo del contexto en el que sea utilizada. Para los propósitos de esta materia se define a un
Sistema como un editor de señales, es decir, en el momento en que le ingrese una señal, la puede
manipular, cambiar, reproducir y finalmente transmitir estas señales previamente ingresadas.
Profundizando en este concepto, una SEÑAL es una variación de patrones que representan o
codifican algún tipo de información. [ CITATION 1 \l 12298 ]

Ahora bien, existen muchos tipos de Señales y Sistemas. Pero por ahora nos interesa definir un tipo
específico de Señales y Sistemas, de los que precisamente hemos utilizado en esta práctica. Las
SEÑALES CONTINUAS, que generalmente hace referencia a las Señales Continuas en el
Tiempo, son señales que pueden ser representadas matemáticamente como una función de una sola
variable real independiente, que justamente esta variable puede ser el tiempo. Con este
conocimiento previo, se hace más sencillo definir a un SISTEMA CONTINUO como el sistema
donde sus entradas y salidas son señales continuas.[ CITATION 1 \l 12298 ]

Con todo esto, gracias a los Sistemas Continuos, que son los encargados de manipular una señal
continua de varias maneras y transmitir ese resultado, podemos definir las OPERACIONES EN
SEÑALES CONTINUAS. Estas operaciones se dividen en dos grupos: las Transformaciones de la
Variable Dependiente y las Transformaciones de la Variable Independiente. Tiene sentido que
solamente sean dos, ya que al ser continuas su representación habitual es en un eje horizontal para
representar al tiempo y un eje vertical. Entonces, estos dos grupos de operaciones se definen cada
una sobre cada eje de la señal. Las Transformaciones de Variable Independiente se definen en el eje
horizontal cuando hablamos que representa el tiempo en la señal.[ CITATION 2 \l 12298 ] En este

1
grupo tenemos la INVERSIÓN O REFLEXIÓN, que como su nombre lo indica, lo que se hace es
cambiar los valores que la señal tomaba en ese instante de tiempo y reemplazarlos por su opuesto
(en caso de que la señal esté definida en un tiempo negativo). Gráficamente podemos observar la
inversión temporal como si la señal fuera “de atrás para adelante” en comparación con la señal
original. En forma de ecuación se define de la siguiente manera:

y ( t ) =x (−t ) (1)

Donde y representa la señal de entrada al sistema y x representa la señal de salida del sistema.

El escalamiento del tiempo, o el FACTOR DE COMPRESIÓN O EXPANSIÓN, es

literalmente comprimir o expandir una señal. Sea a un número real mayor que 0 y diferente de 1,
entonces la operación de Escalamiento de tiempo se define como:

y ( t ) =x ( at ) (2)

Donde, si a> 1, la señal de salida x es una versión comprimida de la señal de entrada y .

El DESPLAZAMIENTO en la variable independiente, gráficamente se denota por mover la señal

a la derecha o a la izquierda, este movimiento depende del signo de la cantidad que se quiera
mover. Definiendo a t 0 como la cantidad de desplazamiento que se le propinará a la señal, esta
operación se define como:

y ( t ) =x ( t−t 0 ) (3)

Donde, si t 0> 0 (positiva), x se desplaza hacia la derecha y si t 0< 0 (negativa), x se desplaza hacia la
izquierda.

Ahora revisaremos las operaciones que se hacen a las señales modificando de alguna manera su
variable dependiente. Estas son: Escalamiento de Amplitud, Suma, Multiplicación, Diferenciación e
Integración.[ CITATION Uni20 \l 12298 ]

Para el ESCALAMIENTO DE AMPLITUD, definimos a la señal x como la señal de entrada y a

la señal y como la señal de salida. Entonces:

y ( t ) =k∗x ( t ) (4)

Donde c es una constante tal que, si c >1 corresponde a una Amplificación y si 0< c<1 corresponde
a una atenuación de la señal.

En la SUMA DE SEÑALES, ahora tendremos dos señales de entrada x 1 y x 2, y a y como la señal

de salida. Entonces:

y ( t ) =x1 ( t )+ x 2 ( t ) (5)

Donde para un tiempo determinado t , la señal y da como resultado la suma de los resultados
independientes de las señales x 1 y x 2.

2
Para la MULTIPLICACIÓN, podemos interpretar a las señales x 1 y x 2como si x 1 amplificara o
atenuara a x 2 por el valor correspondiente a esta señal en un tiempo determinado t , o viceversa.
Este comportamiento se ve reflejado por la siguiente ecuación.

y ( t ) =x1 ( t )∗x 2 ( t ) (6)

Donde, como ya se lo había planteado anteriormente x 1 y x 2 son las señales de entrada, pero ahora
y es la señal de salida correspondiente a la multiplicación de x 1 y x 2.

La DIFERENCIACIÓN E INTEGRACIÓN de señales son procesos donde se emplean

herramientas del Cálculo Diferencial e Integral donde se aprovecha la naturaleza determinística y
continua de las señales. Este cálculo viene denotado por las siguientes ecuaciones:

d
y (t)= x (t) (7)
dt
t
y ( t ) =∫ x ( τ ) dτ (8)
−∞

Donde (7) corresponde a la Diferenciación de la señal x dando como resultado su correspondiente

señal y . Y (8) es la Integración de la señal x resultando en la señal de salida y .

Como último concepto a definir, y con la ayuda de todo el trasfondo construido a partir de los
conceptos ya expuestos, la REPRESENTACIÓN DE SEÑALES POR EXPONENCIALES
COMPLEJOS[ CITATION ESP20 \l 12298 ], que es el resultado de Una multiplicación de una
Sinusoide Compleja por un Exponencial real:

X ( t )=cos ( ω0 t ) + j sen ( ω 0 t ) (9)

Donde X ( t ) corresponde a la representación de una señal por sinusoide compleja, donde cos ( ω0 t )
representa la parte real y j sen ( ω0 t ) representa la parte compleja. Pero esta función X ( t ) puede
ser representada también como:

X ( t )=e j ω t 0
(10)

Entonces, la representación por exponencial compleja quedaría por la multiplicación de la Sinusoide

en (9) por una exponencial real e ∝t . Pero si tomamos la forma de la sinusoide compleja de (10), se
facilita la multiplicación de estos dos términos dando como resultado:
∝+ jω 0 )t
e St =e ∝t e j ω t =e (
0
; S=( ∝+ j ω0 ) (11)

Donde, si S>0 la señal es creciente y si S<0 la señal es decreciente.

Cabe recalcar que esto es una técnica de representación para señales, las que pueden representarse
con otra fórmula o función y seguir describiendo la misma señal, pero como vemos, a veces suele ser
útil para facilitar las operaciones entre estas.

3
RESULTADOS

1. MANIPULACIÓN DE SEÑALES USANDO MATLAB.

a) Genere un vector tiempo (tt) que cubra un rango de t, aproximadamente dos ciclos de las
sinusoides de 4000 Hz que serán definidas en el literal (b). Defina el vector tt. Si usa T
para denotar el período de las sinusoides, defina el tiempo de inicio del vector tt para
que sea igual a -T, y el tiempo de finalización como T. Luego, los dos ciclos incluirán t=0.
Finalmente, asegúrese de tener al menos 25 muestras por período de onda sinusoidal.
En otras palabras, cuando use el operador de dos puntos para definir el vector tiempo,
haga que el incremento sea lo suficientemente pequeño como para generar 25 muestras
por período.

Figura 1. Definición de variable frecuencia, período y salto del vector tiempo.

b) Genere dos sinusoides de 4000 Hz con amplitud arbitraria y con time-delay (t mi)

x 1 ( t )= A1 cos ⁡¿

x 2 ( t )= A 2 cos ⁡¿

Seleccione el valor de las amplitudes y los cambios de tiempo de la siguiente manera: Sea A1
igual a su edad y configure A2=1.6 A 1. Para los desplazamientos de tiempo, configure

( 37.2
t m 1=
M )
∗T y t m2 =− ( 41.3
D )
∗T ; donde D y M son el día y mes de su cumpleaños, y T
es el período.

A1=su edad

A2=1.2 A1

t m 1= ( 37.2
M )
T

t m 2=− ( 41.3
D )
T

4
 Haga un gráfico de ambas señales en el rango de -T < T < T. Use subplot para hacer una
figura con tres gráficos en donde colocará estos dos trazados en la misma ventana de la
figura.

Figura 2. Definición de señales sinusoidales con variables de amplitud y delay de la señal.

Figura 3. Grafica de señales x1(t) y x2(t)

5
 Figura 4. Muestra de las gráficas sinusoidales X1, X2.

c) Cree una tercera sinusoide como la suma: x 3 ( t )=x 1 (t ) + x 2 ( t ). En MATLAB esto equivale
a sumar los vectores que contienen los valores de cada sinusoide. Haga una gráfica de:
x 3 ( t ) en el mismo intervalo de tiempo que el utilizado en los trazados de la parte (b).
Incluya esto como la tercera gráfica en la misma Figura usando subplot.

Figura 5. Definición de la señal X3 como la suma de señales X1 y X2

d) Antes de imprimir las tres gráficas, coloque un título en cada subplot, nombres a los ejes
xlabel('tiempo'); ylabel('x(t)') y muestre la cuadrícula (grid on), tal como se muestra en la
Figura 1.1

6
 Figura 6. Señales de salida graficadas con subplot

7
 2. CÁLCULOS TEÓRICOS
Recuerde que la fase de una sinusoide puede calcularse después de medir el instante de
tiempo de un pico positivo, si conocemos la frecuencia.

a) Realice mediciones de tiempo de un pico positivo y la amplitud de las gráficas de x 1 ( t ) y

x 2 ( t ), y escriba esos valores para Ai y t mi directamente en las gráficas. Luego calcule (a
mano) las fases de las dos señales, x 1 ( t ) y x 2 ( t ), convirtiendo cada time-shift t mi en fase.
Escriba las fases calculadas φ i directamente en las gráficas tal como se muestra en la
Figura 1.2

Figura 7. Valores de picos positivos.

Señal 1 Señal 2
A1=22 A2=26.34
t m 1=0.0001083s t m 1=0.0001s
ω=2 πf ω=2 πf
φ 1=−t m 1 . ω φ 2=−t m 2 . ω
φ 1=−( 0.0001083 ) ( 2 π ) ( 4000 ) φ 2=−(0.0001)(2 π )(4000)
φ 1=−0.8664 π rad φ 2=−0.8 π rad
x 1 ( t )=22 cos ( 2 π ( 4000 ) ( t−0.0001083 ) ) x 2 (t)=26.34 cos ⁡( 2 π ( 4000)(t−0.0001))

8
 b) Mida la amplitud A3 y el desplazamiento temporal t m 3 de x 3 ( t ) directamente en la
gráfica y luego calcule la fase (φ 3 ¿ a mano. Escriba esos valores directamente en el
gráfico para mostrar cómo se midieron la amplitud y el desplazamiento de tiempo (time-
delay) y cómo se calculó su fase.

Señal 3

x 3 ( t )=x 1 (t ) + x 2 ( t )

A3 =48.12

t m 3=0.0001083 s

ω=2 πf
φ 3=−t m 3 . ω

φ 3=−( 0.0001083 ) ( 2 π ) ( 4000 )

φ 3=−0.8664 π rad

c) Ahora use el teorema de sum fasorial. Realice una adición fasorial de amplitudes
complejas para x 1 ( t ) y x 2 ( t ) para determinar la amplitud compleja de x 3 ( t ). Use la
amplitud compleja de x 3 ( t ) para verificar que sus cálculos previos de A3 y (φ 3) fueron
correctos.

Aplicación del Teorema de suma fasorial:

x 1 ( t )=22 cos ( 2 π ( 4000 ) ( t−0.0001083 ) )

x 2 ( t )=26.34 cos ( 2 π ( 4000 )( t−0.0001 ) )

x 3 ( t )=x 1 (t ) + x 2 ( t )=Re {X e j 2 π (4000 ) t (1) Ambas señales tienen la misma frecuencia.

Donde:
2
X =∑ Ak e jφ k

N =1

X =A 1 e j φ + A2 e j φ
1 2

X =22 e j (−0.8664 π ) + 26.34 e j (−0.8 π )

X =−20.09−8.97 j−21.31−15.48 j Transformando coordenadas polares a rectangulares

9
X =−41.4−24.45 j (2)

Transformando (2) a forma polar:

X =48.08 e j (−0.8454 π ) (3)

Sustitución de (3) en (1):

x 3 ( t )=R e {48.08 e j (−0.8454 π ) e j 2 π ( 4000 ) t

x 3 ( t )=48.08 cos ( 2 π ( 4000 ) t−0.8454 π )

d) Calcule el porcentaje de error de 𝐴3 y 𝜑3 con los resultados obtenidos en Matlab y los

obtenidos matemáticamente.

ATeórica =48.08 A Experimental=48.12

A Teórica− A Experimental 48.08−48.12

Error relativo A= | A Teórica
= ||
48.08 |=8.31 x 10
−4

Error porcentual A=( 8.31 x 10−4 ) ( 100 % )=−0.0831

φ Teórica=0.8454 π φ Experimental =−0.8664 π

φTeórica−φ Experimental 0.8454−0.8664
Error relativo φ= | φTeórica
= ||
0.8454 |
=0.024

Error porcentual φ= ( 0.024 ) ( 100 % ) =2.48 %

3. REPRESENTACIÓN DE SINUSOIDES CON EXPONENCIALES COMPLEJOS

Escriba el código en MATLAB que generará valores de la sinusoide 𝑥1(𝑡) que se muestra en la
Fig. 1.2, use la representación de amplitud compleja:
𝑥1(𝑡)=𝑅𝑒{𝑋𝑒𝑗𝜔𝑡}

10
 Figura 1.2 𝒙𝟏(𝒕), generada usando representación de amplitud compleja

Figura 8. Código para generación de gráfica utilizando amplitud compleja.

11
PREGUNTAS:

Para las preguntas del 1 al 4 utilice MATLAB para resolverlas:

1. INDIQUE LA MAGNITUD Y FASE DEL SIGUIENTE NÚMERO DADO EN

COORDENADAS RECTANGULARES: 𝑧1=3+𝑗6.

12
2. GRAFIQUE LA SIGUIENTE SEÑAL S(T):

13
3. GRAFIQUE LA SIGUIENTE SEÑAL PARA LOS 3 VALORES DE FRECUENCIA:

𝑥(𝑡)=7𝑐𝑜𝑠(2𝜋𝑓𝑜𝑡)

a) 𝑓𝑜=200𝐻𝑧

b) 𝑓𝑜=100𝐻𝑧

c) 𝑓𝑜=0𝐻𝑧

14
15
4. GRAFIQUE EN EL MISMO PLOT LAS SIGUIENTES SEÑALES:

a) 𝑥1(𝑡)=3.5𝑐𝑜𝑠(30𝜋𝑡+35𝜋/180)

b) 𝑥2(𝑡)=6.2𝑐𝑜𝑠(30𝜋𝑡+80𝜋/180)

c) 𝑥1(𝑡)+𝑥2(𝑡)

16
5. MATEMÁTICAMENTE REALICE LA SUMA DE 2 SEÑALES Y EXPRESE EL
RESULTADO EN COORDENADAS POLARES:

Nos valemos de la identidad trigonométrica:

cos ( a+b ) =cos ( a ) cos ( b )−sen ( a ) sen ( b )

( 35180π )−4.5 sen ( 30 πt ) sen( 35180π )+7.2 cos ( 30 πt ) cos( 80180π )−7.2 sen ( 30 πt ) sen( 80180π )
4.5 cos ( 30 πt ) cos

35 π 35 π 80 π 80 π
Reemplazando cos ( ) ≈ 0.819 sen ( ) ≈ 0.574 cos ( ) ≈ 0.174 sen (
180 )
, , , ≈ 0.985 .
180 180 180
Sacamos factores comunes:

cos ( 30 πt ) ( 4.5 ∙ 0.819+7.2∙ 0.174 )−sen(30 πt) ¿)

Multiplicando y sumando nos queda:

x s (t) 4.94 cos ( 30 πt )−9.68 sen(30 πt )

Donde x s=x 1 + x 2.

17
Para convertir a polares usamos la fórmula:

r =√t 2 + x 2s

θ=arctan ( xt )
s

2
r =√t 2 +(4.94 cos ( 30 πt ) −9.68 sen ( 30 πt ) )

θ=arctan ( 4.94 cos ( 30 πt )−9.68

t
sen ( 30 πt )
)
La representación (r , θ) sería:

( √ 2
( r , θ )= t 2 + ( 4.94 cos ( 30 πt )−9.68 sen ( 30 πt ) ) , arctan ( 4.94 cos ( 30 πt )−9.68
t
sen ( 30 πt )
))

CONCLUSIONES
En base a lo realizado en la práctica tanto como el trabajo autónomo en parejas para el reporte se
logró cumplir con los siguientes objetivos:

 Se introdujo los aspectos básicos del trabajo en el entorno MATLAB, incluido el uso del
sistema de ayuda, la sintaxis básica de comandos, los números complejos, la indexación de
matrices, el trazado y el uso de sectorización para evitar bucles ineficientes. Por medio de la
lectura y socialización de la documentación del manejo de MATLAB por parte de la profesora
y como parte del trabajo de investigación de los estudiantes.

 Se usó MATLAB para la manipulación de señales y el trazado de su respectiva gráfica, por

medio de comandos donde se especificaba los parámetros necesarios para una señal
sinusoidal tales como su amplitud, frecuencia y fase.

 Se graficó en MATLAB señales sinusoidales a través de amplitud compleja, con el uso de la

teoría socializada en el componente teórico del curso, además de su directa aplicación a la
codificación en el entorno de MATLAB.

18
RECOMENDACIONES:

1. Buscar la documentación necesaria para poder simplificar la escritura del código lo mayor
posible, además de incorporar nuevas funciones y comandos que añadan características a
una gráfica.

2. Estar al día con los temas expuestos en el componente teórico al momento de realizar la
práctica para que se puedan cubrir lo mas rápido posible y poder centrarse en el mero
contenido práctico planificado para el curso.

BIBLIOGRAFÍA

[1] J. H. McClellan, R. W. Schafer y M. A. Yoder, «DPS FIRST Second Edition,» de DPS FIRST Second
Eddition, Harlow, Pearson Education Limited, 2006, pp. 21-26.

[2] I. Bosch Roig y M. Cabedo Fabrés, «Repositorio Institucional de la Universitat Politécnica de

Valencia,» Universitat Politécnica de Valencia, 13 Diciembre 2011. [En línea]. Available:
https://laboratoriosvirtuales.upv.es/eslabon/Ejercicio?do=Funciones_basicas. [Último acceso: 11
Junio 2020].

[3] Universidad Francisco de Paula Santander, «StuDocu,» [En línea]. Available:

https://www.studocu.com/co/document/universidad-francisco-de-paula-santander/teoria-de-
senales-y-sistemas/resumenes/02-tss-lab02-operaciones-basicas-sobre-senales-en-
matlab/5409497/view. [Último acceso: 22 Junio 2020].

[4] ESPOL, «Office 365,» 15 Junio 2020. [En línea]. Available: https://espolec-
my.sharepoint.com/:b:/g/personal/jimjerod_espol_edu_ec/ERumc24CMExCgSnzs8ugzH0BDAut-
NcGLpGmxP_8QGejOQ?e=B81tG7. [Último acceso: 22 Junio 2020].

19
APÉNDICE
%Práctica 1
%Manipulación de señales
clc
clear all
close all
f = 4000;
T = 1/f; %periodo de la señal
tt = -T:T/30:T; %vector de tiempo de las señales, con saltos pequeños de T/30
A1 = 22;
A2=1.2*A1;
M = 5; %Mes de cumpleaños
D=9; %Día de cumpleaños
tm1=((37.2/M)*T); %delay de la señal 1
tm2 =((-41.3/D)*T); %delay de la señal 2
x1 = A1*cos(2*pi*f*(tt-tm1)); %señal 1
x2 = A2*cos(2*pi*f*(tt-tm2)); %señal 2

figure(1)
subplot 311
plot(tt,x1,'r')
title('Senal x1(t)')
xlabel('tiempo(s)')
ylabel('x1(t)')
grid on %activar las cuadrículas

subplot 312
plot(tt,x2,'g')
title('Senal x2(t)')
xlabel('tiempo(s)')
ylabel('x2(t)')
grid on %activar las cuadrículas

%suma de las 2 senales

x3=x1+x2;
subplot 313
plot(tt,x3,'b');
title('Senal x3(t)')
xlabel('tiempo(s)')
ylabel('x3(t)')
grid on %activar las cuadrículas

figure(2)
tt2 = -T:T/100:T; %Vector tiempo que empieza en -T y termina en T, T/100 pasos
x1c = real(20*exp(-i*0.60*pi)*exp(i*2*pi*f*tt2)); % La fn real devuelve la parte
real de la señal compleja
plot(tt2,x1c,'blue'); %Gráfico
xlim([-T,T]);

20
grid on
tittle('x1(t) Generada usando representación de A Compleja')
xlabel('Time(s)'), ylabel('x1(t)')

PREGUNTA 1
clc
clear all
close all
%Pregunta 1
% Encontrar magnitud y fase del siguiente número 𝑧1=3+𝑗6.
z1 = 3 + 6j
r = abs(z1) % Encontrar magnitud
tetha = angle(z1) % Encontrar fase

PREGUNTA 2
clc
clear all
close all

%Pregunta 2

t= -4:8;

a = ((t>=0) & (t<=(0.5))) .* (t*2) + ((t>=0.5) & (t<=2)) .* ((4-

(2*t))*(1/3)) + 0; %DEFINICION DE INTERVALOS DE t

figure(1)
subplot 311
plot(t,a,'r')
title('Senal x1(t)')
xlabel('tiempo(s)')
ylabel('x1(t)')
grid on %activar las cuadrículas

clc
clear all
close all

PREGUNTA 3
%Pregunta 3

%valores de frecuencia
f = 200;
f2= 100;
f3=0;

A1 = 7; %amplitud

T = 1/f; %periodo de la señal

tt = -T:T/30:T; %vector de tiempo de las señales, con saltos pequeños de

T/3

21
x1 = A1*cos(2*pi*f*(tt)); %señal 1
x2 = A1*cos(2*pi*f2*(tt)); %señal 2
x3 = A1*cos(2*pi*f3*(tt)); %señal 2

figure(1)
subplot 311
plot(tt,x1,'r')
title('Senal x1(t)')
xlabel('tiempo(s)')
ylabel('x1(t)')
grid on %activar las cuadrículas

subplot 312
plot(tt,x2,'g')
title('Senal x2(t)')
xlabel('tiempo(s)')
ylabel('x2(t)')
grid on %activar las cuadrículas

subplot 313
plot(tt,x3,'g')
title('Senal x3(t)')
xlabel('tiempo(s)')
ylabel('x3(t)')
grid on %activar las cuadrículas

PREGUNTA 4

clc
clear all
close all

%Pregunta 4

tt1 = 0:1/100:1; %vector de tiempo de las señales

x11 = 3.5*cos((30*pi*tt1) + ((35*pi)/180)); %señal 1

x22 = 6.2*cos((30*pi*tt1)+ ((80*pi)/180)); %señal 2
x33 = x11 + x22;

figure(4)
subplot 311
plot(tt1,x11,'r')
title('Senal x11(t)')
xlabel('tiempo(s)')
ylabel('x11(t)')
grid on %activar las cuadrículas

subplot 312
plot(tt1,x22,'g')
title('Senal x22(t)')
xlabel('tiempo(s)')
ylabel('x22(t)')
grid on %activar las cuadrículas

22
subplot 313
plot(tt1,x33,'g')
title('Senal x33(t)')
xlabel('tiempo(s)')
ylabel('x33(t)')
grid on %activar las cuadrículas

23