0% encontró este documento útil (0 votos)
22 vistas16 páginas

Informe Lab 2 Señales - Manuel Vergara - Daniel Torres

Este informe describe algoritmos para resolver la convolución de señales continuas y discretas manualmente y con MATLAB. Se propone calcular la convolución de dos señales pulso a mano y verificar los resultados con MATLAB. Adicionalmente, se analiza un sistema LTI usando ecuaciones diferenciales para encontrar la respuesta al impulso.

Cargado por

Andres Genoy
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
22 vistas16 páginas

Informe Lab 2 Señales - Manuel Vergara - Daniel Torres

Este informe describe algoritmos para resolver la convolución de señales continuas y discretas manualmente y con MATLAB. Se propone calcular la convolución de dos señales pulso a mano y verificar los resultados con MATLAB. Adicionalmente, se analiza un sistema LTI usando ecuaciones diferenciales para encontrar la respuesta al impulso.

Cargado por

Andres Genoy
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

Informe de laboratorio 2

Daniel Torres González


Manuel Vergara
(13 de abril 2023)

Resumen— III. EJERCICIO 1


El propósito de este informe es desarrollar un algoritmo para
resolver la operación de convolución entre señales continuas y
discretas y verificar los resultados con el software De las siguientes señales obtenga la convolución h ( t )=f 1 ( t )
MATLAB/SIMULINK. Sugiere resolver manualmente la * f 2 ( t ) a mano y luego compruebe sus resultados usando
convolución de la señal y verificar los resultados con Matlab
MATLAB, definir matemáticamente la operación de
convolución de la señal discreta y escribir los algoritmos del Descripción y definición.
método gráfico y del método de tira deslizante en MATLAB
sin utilizar la función conv(). Además, usando el sistema LTI La convolución es una operación matemática utilizada para
dado por las ecuaciones diferenciales, encontramos el combinar dos señales y obtener información sobre cómo una
polinomio característico, construimos la solución de la afecta a la otra. En este caso, se recomienda calcular
ecuación diferencial en estado cero y entrada cero, trazamos la manualmente la convolución entre las dos señales f 1 ( t )y
respuesta homogénea con condiciones iniciales específicas y
f 2 ( t )y luego usar el software MATLAB para verificar los
determinamos la respuesta de impulso del sistema.
resultados. La convolución se denota con el símbolo *, que se
Palabras clave—MATLAB, Señal, Variable, Convolución expresa matemáticamente de la siguiente manera:

h ( t )=f 1 ( t )* f 2 ( t )
I. INTRODUCCIÓN
El resultado de la convolución es una tercera señal ℎ(𝑡) que
E STE documento se encuentra basado y regido bajo el
proporciona información sobre cómo interactúan las dos
señales originales. Lograr la convolución a mano requiere una
formato Institute of Electrical and Electronics Engineers, serie de operaciones matemáticas, incluida la inversión de una
(IEEE), realizado en el procesador de texto, Microsoft Word, señal y la integración de la multiplicación de dos señales.
con fuente ´´Time New Romance´´. Después de la adquisición manual de la convolución, los
resultados se pueden verificar con el software MATLAB, que
La convolución es una operación matemática esencial en el tiene una función integrada para calcular la convolución de la
procesamiento de señales que permite combinar dos señales señal. De esta manera, puede verificar manualmente los
para obtener información sobre cómo una afecta a la otra. Este resultados y asegurarse de que sean correctos. Lograr la
informe describe la construcción de algoritmos para tratar convolución de la señal es la base del procesamiento de la
sistemáticamente las operaciones de convolución entre señales señal y se puede utilizar en varios campos, como las
continuas y discretas. Además, se exploran conceptos como comunicaciones, la ingeniería eléctrica y la acústica.
respuesta al impulso y convolución de señales continuas y
discretas, y se verifican resultados teóricos y prácticos con el a.)
software MATLAB/SIMULINK. La revisión se divide en
objetivos específicos, que incluyen soluciones manuales y de
software de convoluciones de señales, definición matemática
de operaciones de convolución de señales discretas y
programación de algoritmos apropiados, y análisis de un
sistema LTI dado usando ecuaciones diferenciales. Se requiere
una computadora con el software MATLAB/SIMULINK Fig. 1 Señales pulso, punto a.
instalado para completar este proyecto.
Para este punto tenemos dos señales pulso. Primero debemos
II. OBJETIVOS expresarlas de forma matemática para poderles hacer la
1. Construir algoritmos para resolver de forma sistemática la convolución. Recordando que las señales pulso son la
operación de convolución entre dos señales. diferencia de 2 señales escalón unitario. Respectivamente
2. Estudiar conceptos como la respuesta impulsiva y quedarían:
convolución para señales continuas y discretas
3. Verificar los resultados teóricos y prácticos a través de la f 1 ( t )=u ( t−4 )−u ( t−6 )
herramienta de software. f 2 ( t )=u ( t+5 )−u ( t+ 4 )
2

t +5

Ahora el paso siguiente es cambiar de f 1 ( t ), la variable t por ∫ u ( t+1 ) dτ


4
τ y de nuestra f 2 ( t ), la variable t por t−τ . Entonces
tendríamos:
Queda u ( t+1 ) ya que cuando despejamos t nos dio –1 por lo
f 1 ( τ )=u ( τ−4 )−u ( τ−6 ) que gráficamente estaría en el lado negativo por lo que en
f 2 ( t−τ )=u ( t+5−τ ) −u ( t+ 4−τ ) estos casos se representa con una suma para indicar esto. Si
nos fijamos u ( t+1 ) se toma como constante ya que la
Ahora lo que haremos será reemplazar en la fórmula de constante de integración es, luego se resuelve:
convolución la cual nos dice que:
t +5
u ( t+1 ) ∫ dτ

∫ f ( τ ) ⋅g (t−τ ) dτ 4
−∞
Reemplazando. t +5
¿ u ( t+1 ) [ τ ]4

∫ u ( τ−4 )−u ( τ−6 ) ⋅u ( t+ 5−τ )−u ( t+ 4−τ ) dτ ¿ u ( t+1 ) [ t+5−4 ]


−∞

¿ u ( t+1 ) [ t+1 ]
Ahora lo que haremos será resolver la integral por partes y con
Este será nuestro primer trozo de la señal de convolución,
el método analítico. Por lo tanto, nuestra primera integral
Haremos lo mismo con las demás integrales:
tendríamos:

∞ Ahora nuestra segunda integral tendríamos:


∫ u ( τ−4 ) ⋅u ( t+5−τ ) dτ ∞
−∞
∫ u ( τ−4 ) ⋅u ( t+ 4−τ ) dτ
−∞
El siguiente paso es hacer que la función que está dentro de la
integral la llevemos a 0, esto con el objetivo de establecer los
nuevos límites de integración. Si observamos solo hay 2
El siguiente paso es hacer que la función que está dentro de la
formas de hacerlo:
integral la llevemos a 0, esto con el objetivo de establecer los
nuevos límites de integración. Si observamos solo hay 2
u ( τ −4 )=0 formas de hacerlo:
τ =4 u ( τ −4 )=0
τ =4
u ( t+5−τ )=0
τ =t+5 u ( t+ 4−τ ) =0
τ =t+ 4
Ahora se crearán unos nuevos intervalos en donde por
naturaleza de ser escalones unitarios tendrán los intervalos 1 y
0:
Ahora se crearán unos nuevos intervalos en donde por
naturaleza de ser escalones unitarios tendrán los intervalos 1 y
{ 1 ; 4 <τ <t +5 } 0:
{ 0 ; Otro caso }
{ 1 ; 4 <τ <t +4 }
El siguiente paso es despejar t del intervalo de 1, despejaremos { 0 ; Otro caso }
t:
El siguiente paso es despejar t del intervalo de 1, despejaremos
4 < τ< t+5 t:
4−5< τ−5 <t
−1< τ−5<t 4 < τ< t+ 4
t >−1 4−4< τ−4< t
0< τ−4<t
Con esto ya podemos plantear nuestra nueva integral:
t >0
3

Con esto ya podemos plantear nuestra nueva integral: Con esto ya podemos plantear nuestra nueva integral:

t +4 t +5

∫ u ( t ) dτ ∫ u ( t−1 ) dτ
4 6

Queda u(t) ya que cuando despejamos t nos dio 0 por lo que Queda u ( t−1 ) ya que cuando despejamos t nos dio 1 por lo
gráficamente estaría en 0 por lo que en estos casos se que gráficamente estaría en el lado positivo por lo que en estos
representa de esta forma. Si nos fijamos u(t) se toma como casos se representa con una resta para indicar esto. Si nos
constante ya que la constante de integración es, luego se fijamos u ( t−1 ) se toma como constante ya que la constante
resuelve:
de integración es, luego se resuelve:

t +4
u ( t ) ∫ dτ t +5
u ( t−1 ) ∫ dτ
4

t +4 6
¿ u ( t ) [ τ ]4
t+5
¿ u ( t−1 ) [ τ ] 6
¿ u ( t ) [ t+ 4−4 ]
¿ u ( t−1 ) [ t +5−6 ]
¿ u (t )[ t ]
¿ u ( t−1 ) [ t−1 ]
Ahora nuestra tercera integral tendríamos:

∫ −u ( τ−6 ) ⋅u ( t+ 5−τ ) dτ Ahora nuestra cuarta integral tendríamos:


−∞

∫ −u ( τ−6 ) ⋅−u ( t+ 4−τ ) dτ


−∞
El siguiente paso es hacer que la función que está dentro de la
integral la llevemos a 0, esto con el objetivo de establecer los
nuevos límites de integración. Si observamos solo hay 2
El siguiente paso es hacer que la función que está dentro de la
formas de hacerlo:
integral la llevemos a 0, esto con el objetivo de establecer los
u ( τ −6 )=0 nuevos límites de integración. Si observamos solo hay 2
τ =6 formas de hacerlo:
u ( τ −6 )=0
u ( t+5−τ )=0 τ =6
τ =t+5
u ( t+ 4−τ ) =0
Ahora se crearán unos nuevos intervalos en donde por
τ =t+ 4
naturaleza de ser escalones unitarios tendrán los intervalos 1 y
0:
Ahora se crearán unos nuevos intervalos en donde por
naturaleza de ser escalones unitarios tendrán los intervalos 1 y
{ 1 ; 6< τ< t+5 } 0:
{ 0 ; Otro caso }
{ 1 ; 6< τ< t+ 4 }
El siguiente paso es despejar t del intervalo de 1, despejaremos
t:
{ 0 ; Otro caso }
El siguiente paso es despejar t del intervalo de 1, despejaremos
6< τ <t+ 5 t:
6−5<τ −5<t
1< τ−5<t 6< τ <t+ 4
t >1 6−4 <τ −4 <t
4

2< τ−4 <t


t >2
Con esto ya podemos plantear nuestra nueva integral:

t +4

∫ u ( t−2 ) dτ
6

Queda u ( t−2 ) ya que cuando despejamos t nos dio 1 por lo


que gráficamente estaría en el lado positivo por lo que en estos
casos se representa con una resta para indicar esto. Si nos
fijamos u ( t−2 ) se toma como constante ya que la constante Fig. 3 Gráfica señales.
de integración es, luego se resuelve:
Como podemos observar nuestra señal f1(t) es la de color
naranja y la de azul es nuestra f2(t). La señal se amarillo es
nuestra convolución de las otras dos señales.
t+ 4
u ( t−2 ) ∫ dτ Para verificar que la señal amarilla corresponde a los
6 resultados obtenidos de las integrales se creó un código en el
t+ 4
¿ u ( t−2 ) [ τ ] 6 cuál se graficó los resultados obtenidos por trozos:

¿ u ( t−2 ) [ t +4−6 ]

¿ u ( t−2 ) [ t−2 ]

Demostración en Matlab.

Primero se creó un código en el cuál podamos graficar ya


hacer la convolución de nuestras dos señales f1(t) y f2(t), la
cual nos da:
Fig. 4 Código Verificación.
El anterior código nos da como resultado:

Fig. 2 Código ejercicio 1

Como resultados nos arrojó lo siguiente:

Fig. 5 Gráfica convolución


5

Como podemos observar la anterior gráfica es la misma que


obtuvimos anteriormente por lo cual podemos afirmar que el u ( t+5−τ )=0
procedimiento y los resultados son correctos. τ =t+5
b.) Ahora se crearán unos nuevos intervalos en donde por
naturaleza de ser escalones unitarios tendrán los intervalos 1 y
0:

{ 1 ; 3<τ <t +5 }
{ 0 ; Otro caso }
Fig. 6 Señales pulso punto b. El siguiente paso es despejar t del intervalo de 1, despejaremos
t:

Para este punto tenemos dos señales pulso. Primero debemos 3< τ <t+5
expresarlas de forma matemática para poderles hacer la
convolución. Recordando que las señales pulso son la
3−5<τ −5<t
diferencia de 2 señales escalón unitario. Respectivamente −2< τ−5<t
quedarían: t >−2

f 1 ( t )=u ( t−3 )−u ( t−5 ) Con esto ya podemos plantear nuestra nueva integral:
f 2 ( t )=u ( t+5 )−u ( t+ 3 ) t +5

Ahora el paso siguiente es cambiar de f 1 ( t ), la variable t por


∫ u ( t−2 ) dτ
3
τ y de nuestra f 2 ( t ), la variable t por t−τ . Entonces
tendríamos: Queda u ( t−2 ) ya que cuando despejamos t nos dio –1 por lo
f 1 ( τ )=u ( τ−3 )−u ( τ −5 ) que gráficamente estaría en el lado positivo por lo que en estos
casos se representa con una resta para indicar esto. Si nos
f 2 ( t−τ )=u ( t+5−τ ) −u ( t+3−τ )
fijamos u ( t−2 ) se toma como constante ya que la constante
Ahora lo que haremos será reemplazar en la fórmula de de integración es τ , luego se resuelve:
convolución la cual nos dice que:
t+5
∞ u ( t−2 ) ∫ dτ
∫ f ( τ ) ⋅g (t−τ ) dτ 3
−∞
t+ 5
Reemplazando. ¿ u ( t−2 ) [ τ ] 3

¿ u ( t−2 ) [ t +5−3 ]
∫ u ( τ−3 )−u ( τ−5 ) ⋅u ( t+5−τ )−u ( t+3−τ ) dτ
−∞
¿ u ( t+2 ) [ t+ 2 ]

Ahora lo que haremos será resolver la integral por partes y con Este será nuestro primer trozo de la señal de convolución,
el método analítico. Por lo tanto, nuestra primera integral Haremos lo mismo con las demás integrales:
tendríamos:
Ahora nuestra segunda integral tendríamos:

∫ u ( τ−3 ) ⋅u ( t+5−τ ) dτ ∞
−∞ ∫ u ( τ−3 ) ⋅−u (t +3−τ ) dτ
−∞
El siguiente paso es hacer que la función que está dentro de la
integral la llevemos a 0, esto con el objetivo de establecer los
nuevos límites de integración. Si observamos solo hay 2 El siguiente paso es hacer que la función que está dentro de la
formas de hacerlo: integral la llevemos a 0, esto con el objetivo de establecer los
nuevos límites de integración. Si observamos solo hay 2
u ( τ −3 )=0 formas de hacerlo:
τ =3 u ( τ −3 )=0
6

τ =3 τ =5

u ( t+ 3−τ )=0 u ( t+5−τ )=0


τ =t+3 τ =t+5

Ahora se crearán unos nuevos intervalos en donde por Ahora se crearán unos nuevos intervalos en donde por
naturaleza de ser escalones unitarios tendrán los intervalos 1 y naturaleza de ser escalones unitarios tendrán los intervalos 1 y
0: 0:

{ 1 ; 3<τ <t +3 } { 1 ; 5<τ <t +5 }


{ 0 ; Otro caso } { 0 ; Otro caso }
El siguiente paso es despejar t del intervalo de 1, despejaremos El siguiente paso es despejar t del intervalo de 1, despejaremos
t: t:

3< τ <t+3 5< τ <t+5


3−3<τ −3<t 5−5<τ −5<t
0< τ−3< t 0< τ−5< t
t >0 t >0
Con esto ya podemos plantear nuestra nueva integral: Con esto ya podemos plantear nuestra nueva integral:

t +3 t +5

∫ u ( t ) dτ ∫ u ( t ) dτ
3 5

Queda u(t) ya que cuando despejamos t nos dio 0 por lo que Queda u(t) ya que cuando despejamos t nos dio 0 por lo que
gráficamente estaría en 0 por lo que en estos casos se gráficamente estaría en 0 por lo que en estos casos se
representa de esta forma. Si nos fijamos u(t) se toma como representa de esta forma. Si nos fijamos u(t) se toma como
constante ya que la constante de integración es τ , luego se constante ya que la constante de integración es τ , luego se
resuelve: resuelve:

t +5
t +3 u ( t ) ∫ dτ
u ( t ) ∫ dτ 5
3
t +5
t +3
¿ u ( t ) [ τ ]5
¿ u (t )[ τ ] 3
¿ u ( t ) [ t+5−5 ]
¿ u ( t ) [ t+ 3−3 ]
¿ u (t )[ t ]
¿ u (t )[ t ]

Ahora nuestra tercera integral tendríamos: Ahora nuestra cuarta integral tendríamos:

∞ ∞

∫ −u ( τ−5 ) ⋅−u (t +5−τ ) dτ ∫ −u ( τ−5 ) ⋅−u (t +3−τ ) dτ


−∞ −∞

El siguiente paso es hacer que la función que está dentro de la El siguiente paso es hacer que la función que está dentro de la
integral la llevemos a 0, esto con el objetivo de establecer los integral la llevemos a 0, esto con el objetivo de establecer los
nuevos límites de integración. Si observamos solo hay 2 nuevos límites de integración. Si observamos solo hay 2
formas de hacerlo: formas de hacerlo:
u ( τ −5 )=0 u ( τ −5 )=0
7

τ =5

u ( t+ 3−τ )=0
τ =t+3

Ahora se crearán unos nuevos intervalos en donde por


naturaleza de ser escalones unitarios tendrán los intervalos 1 y
0:

{ 1 ; 5<τ <t +3 }
{ 0 ; Otro caso }
El siguiente paso es despejar t del intervalo de 1, despejaremos
t:

5< τ <t+3 Fig. 6 Código señales.


5−3<τ −3<t
2< τ−3<t Como resultados nos arrojó lo siguiente:
t >2
Con esto ya podemos plantear nuestra nueva integral:

t +3

∫ u ( t−2 ) dτ
5

Queda u ( t−2 ) ya que cuando despejamos t nos dio 1 por lo


que gráficamente estaría en el lado positivo por lo que en estos
casos se representa con una resta para indicar esto. Si nos
fijamos u ( t−2 ) se toma como constante ya que la constante
de integración es τ , luego se resuelve:

t+3
u ( t−2 ) ∫ dτ
5
Fig. 7 Gráficas señales.
t+ 3
¿ u ( t−2 ) [ τ ] 5 Como podemos observar nuestra señal f1(t) y f2(t) coinciden
con la del ejercicio propuesto. Tener en cuenta la 3 gráfica la
¿ u ( t−2 ) [ t +3−5 ] cual corresponde a la de la convolución.

¿ u ( t−2 ) [ t−2 ] Para verificar que la tercera señal corresponde a los resultados
obtenidos de las integrales se creó un código en el cual se
graficó los resultados obtenidos por trozos:
Demostración en Matlab.

Primero se creó un código en el cuál podamos graficar ya


hacer la convolución de nuestras dos señales f1(t) y f2(t), la
cual nos da:

Fig. 8 Código señal verificación.

El anterior código nos da como resultado:


8

∫ u ( τ+ 2 ) ⋅u ( t +2−τ )−u ( t−1−τ ) dτ


−∞

Ahora lo que haremos será resolver la integral por partes y con


el método analítico. Por lo tanto, nuestra primera integral
tendríamos:

∫ u ( τ+ 2 ) ⋅u ( t +2−τ ) dτ
−∞

El siguiente paso es hacer que la función que está dentro de la


integral la llevemos a 0, esto con el objetivo de establecer los
nuevos límites de integración. Si observamos solo hay 2
formas de hacerlo:
Fig. 9 Código señal convolución.
u ( τ +2 ) =0
τ =−2
Como podemos observar la anterior gráfica es la misma que
obtuvimos anteriormente por lo cual podemos afirmar que el u ( t+ 2−τ )=0
procedimiento y los resultados son correctos.
τ =t+ 2

c.) Ahora se crearán unos nuevos intervalos en donde por


naturaleza de ser escalones unitarios tendrán los intervalos 1 y
0:

{ 1 ;−2< τ <t+ 2 }
{ 0 ; Otro caso }
El siguiente paso es despejar t del intervalo de 1, despejaremos
Fig. 10 Señales punto c. t:
Para este punto tenemos una señal escalón unitario y otra señal
pulso. Primero debemos expresarlas de forma matemática para −2< τ< t+2
poderles hacer la convolución. Recordando que las señales −2−2< τ−2< t
pulso son la diferencia de 2 señales escalón unitario. −4< τ−2<t
Respectivamente quedarían: t >−4

f 1 ( t )=u ( t+2 ) Con esto ya podemos plantear nuestra nueva integral:


f 2 ( t )=u ( t+ 2 )−u ( t−1 ) t +2

Ahora el paso siguiente es cambiar de f 1 ( t ), la variable t por


∫ u ( t+ 4 ) dτ
−2
τ y de nuestra f 2 ( t ), la variable t por t−τ . Entonces
tendríamos: Queda u ( t+ 4 ) ya que cuando despejamos t nos dio –4 por lo
f 1 ( τ )=u ( τ +2 ) que gráficamente estaría en el lado negativo por lo que en
f 2 ( t−τ )=u ( t+2−τ ) −u (t−1−τ ) estos casos se representa con una resta para indicar esto. Si
nos fijamos u ( t+ 4 ) se toma como constante ya que la
Ahora lo que haremos será reemplazar en la fórmula de constante de integración es τ , luego se resuelve:
convolución la cual nos dice que:
t +2
∞ u ( t+ 4 ) ∫ dτ
∫ f ( τ ) ⋅g (t−τ ) dτ −2
−∞
Reemplazando. t +2
¿ u ( t+ 4 ) [ τ ]−2
9

t +1

¿ u ( t+ 4 ) [ t+ 2+ 2 ] u ( t+1 ) ∫ dτ
−2

¿ u ( t+ 4 ) [ t+ 4 ] t +1
¿ u ( t+1 ) [ τ ]−2
Este será nuestro primer trozo de la señal de convolución,
Haremos lo mismo con las demás integrales: ¿ u ( t+1 ) [ t+1−2 ]
Ahora nuestra segunda integral tendríamos: ¿ u ( t+1 ) [ t−1 ]

∫ u ( τ+ 2 ) ⋅−u ( t−1−τ ) dτ
−∞

Demostración en Matlab.
El siguiente paso es hacer que la función que está dentro de la
integral la llevemos a 0, esto con el objetivo de establecer los
Primero se creó un código en el cuál podamos graficar ya
nuevos límites de integración. Si observamos solo hay 2
hacer la convolución de nuestras dos señales f1(t) y f2(t), la
formas de hacerlo:
cual nos da:
u ( τ +2 ) =0
τ =−2

u ( t−1−τ )=0
τ =t−1

Ahora se crearán unos nuevos intervalos en donde por


naturaleza de ser escalones unitarios tendrán los intervalos 1 y
0:

{ 1 ;−2< τ <t−1 }
{ 0 ; Otro caso }
El siguiente paso es despejar t del intervalo de 1, despejaremos
t:

−2< τ< t−1


1−2< τ +1<t Fig. 11 Código señales
−1< τ+ 1< t
t >−1 .
Como resultados nos arrojó lo siguiente:
Con esto ya podemos plantear nuestra nueva integral:

t −1

∫ u ( t +1 ) dτ
−2

Queda u ( t+1 ) ya que cuando despejamos t nos dio –1 por lo


que gráficamente estaría en el lado negativo por lo que en
estos casos se representa con una resta para indicar esto. Si
nos fijamos u ( t+1 ) se toma como constante ya que la
constante de integración es τ , luego se resuelve:
10

poderles hacer la convolución. Recordando que las señales


Fig. 12 Gráfica señales pulso son la diferencia de 2 señales escalón unitario.
Respectivamente quedarían:
Como podemos observar nuestra señal f1(t) y f2(t) coinciden
con la del ejercicio propuesto. Tener en cuenta la 3 gráfica la −t
f 1 ( t )=e u ( t )
cual corresponde a la de la convolución.
f 2 ( t )=u ( t+3 )−u ( t )
Para verificar que la tercera señal corresponde a los resultados
obtenidos de las integrales se creó un código en el cual se Ahora el paso siguiente es cambiar de f 1 ( t ), la variable t por
graficó los resultados obtenidos por trozos: τ y de nuestra f 2 ( t ), la variable t por t−τ . Entonces
tendríamos:

−τ
f 1 ( τ )=e u ( τ )
f 2 ( t−τ )=u ( t+3−τ ) −u ( t−τ )

Ahora lo que haremos será reemplazar en la fórmula de


convolución la cual nos dice que:
Fig. 13 Código señal verificación. ∞

∫ f ( τ ) ⋅g (t−τ ) dτ
−∞
El anterior código nos da como resultado:
Reemplazando.

∫ e−t u ( τ ) ⋅u ( t +3−τ )−u ( t−τ ) dτ


−∞

Ahora lo que haremos será resolver la integral por partes y con


el método analítico. Por lo tanto, nuestra primera integral
tendríamos:

∫ e−t u ( τ ) ⋅u ( t +3−τ ) dτ
−∞

El siguiente paso es hacer que la función que está dentro de la


integral la llevemos a 0, esto con el objetivo de establecer los
nuevos límites de integración. Si observamos solo hay 2
formas de hacerlo:

Fig. 14 Señal verificación. u ( τ ) =0


τ =0
Como podemos observar la anterior gráfica es la misma que
obtuvimos anteriormente por lo cual podemos afirmar que el u ( t+ 3−τ )=0
procedimiento y los resultados son correctos. τ =t+3
d.) Con esto ya podemos separar por integrales quedándonos así:

t +3 t

∫e −τ
u ( t+ 3 ) dτ−∫ e u ( t ) dτ
−τ

0 0

Ahora resolvemos la integral:


Fig. 15 Señales punto d.
u ( t+ 3 ) ¿
Para este punto tenemos una señal exponencial y otra señal
pulso. Primero debemos expresarlas de forma matemática para ¿ u ( t+3 ) [ 1−e−(t +3) ] −u ( t ) 1−[e−t ]
11

Demostración en Matlab.

Primero se creó un código en el cuál podamos graficar ya


hacer la convolución de nuestras dos señales f1(t) y f2(t), la
cual nos da: Fig. 18 Código señal verificación
El anterior código nos da como resultado:

Fig. 16 Código señales.

Como resultados nos arrojó lo siguiente: Fig. 19 Gráfica señal verificación

Como podemos observar la anterior gráfica es la misma que


obtuvimos anteriormente por lo cual podemos afirmar que el
procedimiento y los resultados son correctos.

f.)

Fig. 20 Señales punto f.

. Para este punto tenemos una señal exponencial y otra señal


pulso. Primero debemos expresarlas de forma matemática para
poderles hacer la convolución. Recordando que las señales
pulso son la diferencia de 2 señales escalón unitario.
Respectivamente quedarían:

−t
f 1 ( t )=e u ( t )
Fig. 17 Gráfica señales.
f 2 ( t )=u ( t )−u ( t−3 )

Como podemos observar nuestra señal f1(t) y f2(t) coinciden Ahora el paso siguiente es cambiar de f 1 ( t ), la variable t por
con la del ejercicio propuesto. Tener en cuenta la 3 gráfica la τ y de nuestra f 2 ( t ), la variable t por t−τ . Entonces
cual corresponde a la de la convolución. tendríamos:

Para verificar que la tercera señal corresponde a los resultados −τ


obtenidos de las integrales se creó un código en el cual se
f 1 ( τ )=e u ( τ )
graficó los resultados obtenidos por trozos: f 2 ( t−τ )=u ( t−τ )−u ( t−τ−3 )
12

Ahora lo que haremos será reemplazar en la fórmula de


convolución la cual nos dice que:

∫ f ( τ ) ⋅g (t−τ ) dτ
−∞
Reemplazando.

∫ e−t u ( τ ) ⋅u ( t−τ )−u(t−τ−3)dτ


−∞

Ahora lo que haremos será resolver la integral por partes y con


el método analítico. Por lo tanto, nuestra primera integral
tendríamos:

∞ ∞

∫ e u ( τ )∗u(t−τ )dτ −∫ e u ( τ )∗u (t−τ −3)dτ


−τ −τ
Fig. 21 Código señales.
−∞ −∞ Como resultados nos arrojó lo siguiente:
El siguiente paso es hacer que la función que está dentro de la
integral la llevemos a 0, esto con el objetivo de establecer los
nuevos límites de integración. Si observamos solo hay 2
formas de hacerlo:

u ( τ ) =0
τ =0

u ( t−3−τ )=0
τ =t−3
Con esto ya podemos separar por integrales quedándonos así:

t t −3

∫e −τ
u ( t ) dτ− ∫ e u ( t−3 ) dτ
−τ

0 0

Ahora resolvemos la integral: Fig. 22 Gráficas señales.

u (t )¿ Como podemos observar la anterior gráfica es la misma que


obtuvimos anteriormente por lo cual podemos afirmar que el
¿ u ( t ) [ 1−e−t ] −u ( t −3 ) [1−e−(t−3 )] procedimiento y los resultados son correctos.

Demostración en Matlab.
IV. EJERCICIO 2
Primero se creó un código en el cuál podamos graficar ya
hacer la convolución de nuestras dos señales f1(t) y f2(t), la Definir matemáticamente la operación de convolución para
cual nos da: señales discretas. Consulte el método gráfico y el método de
cinta deslizante (Sliding Tape Method), programe ambos
algoritmos en Matlab sin usar la función conv().

Para definir una convolución de manera matemática usamos la


siguiente formula:
13

∞ Ahora hablemos del método de la cinta, este es un algoritmo


¿
X 1 [n] X 2 [n]=h[n] ∑ X 1 [k ] X 2 [n−k ] para convolucionar dos señales discretas. Este enfoque es
k=−∞ especialmente útil cuando se trabaja con señales grandes, ya
que reduce la cantidad de pasos necesarios para realizar la
Donde tomamos X1 y X2 como cada una de las funciones las convolución. La idea básica del método de cinta es dividir la
cuales en este caso serian: señal en segmentos más pequeños y apilarlos para calcular la
convolución. En lugar de hacer toda la convolución a la vez, la
X 1 =[ A 1 , A2 , A3 , A 4 , ...... An , ] señal de entrada se divide en segmentos más pequeños y cada
X 2 =[B 1 , B2 , ] segmento se convolución con la señal de respuesta. Los
resultados de cada convolución luego se suman para obtener la
Las convoluciones las podemos realizar de varios métodos por respuesta completa.
ahora usaremos el método de cinta para eso se hace la El proceso de superposición se realiza utilizando una técnica
multiplicación de las funciones y luego se suman para obtener llamada "desplazamiento de cinta", donde una cinta virtual se
C la cual será la respuesta de la convalidación de la siguiente mueve sobre la señal de entrada y se calcula una convolución
manera en cada posición. A medida que avanza la cinta, las piezas se
superponen y se suman para obtener la respuesta completa.
El enfoque de cinta es eficiente porque la mayoría de las
operaciones de convolución se realizan solo en pequeños
segmentos y no en toda la señal. Esto reduce el número total
A1 A2 A3 A4 de pasos necesarios y, por lo tanto, acelera el proceso de
B1 B2 plegado.

[B3 . A 1 ]=C 1
Demostración en Matlab del método de cinta.
A1 A2 A3 A4 Ahora para graficar mediante el método de cinta se usó el
B1 B2 siguiente código:

[B2 . A 1 +B 1 . A 2 ]=C 2

A1 A2 A3 A4
B1 B2

[B2 . A 2 +B 1 . A 3 ]=C3

A1 A2 A3 A4
B1 B2

[B2 . A 3 +B 1 . A 4 ]=C 4

A1 A2 A3 A4
B1 B2

[B1 . A 4 ]=C5

De esta manera quedaría armada la función en convolución.

h=(C 1 ,C 2 ,C 3 ,C 4 , C4 )
Fig.23 Código método de cinta.
Método de cinta.
El cual después de correr el código con valores la gráfica
daría:
14

Fig.24 Grafica del código de cinta

Método de grafica

El método gráfico consiste en graficar las dos señales en el


plano cartesiano, mover una de las señales y superponer la
otra, multiplicar y sumar los valores en cada punto de
superposición para obtener la respuesta completa. Esta es una
forma intuitiva de comprender la convolución que es útil para
cálculos rápidos en señales pequeñas, pero puede volverse
engorrosa para señales grandes o complejas. Por lo tanto, el
método de la banda divide la señal en pequeños segmentos y
superpone cada uno de ellos para calcular la convolución en
cada ubicación de la banda. Los resultados de cada
convolución luego se suman para obtener la respuesta
completa. Es muy eficiente y funciona bien para señales
grandes, ya que reduce el número de pasos necesarios.

Demostración en Matlab

Fig.25 Código método de gráfico.

El cual después de correr el código con los mismos valores


que se le dinero a método de cinta para corroborar que ambos
métodos funcionan y daría lo mismo:
15

h(t)=C 1 e(−2 t ) +C 2 e(−3 t )


la solución completa de la ecuación diferencial es igual a la
solución homogénea asociada:

(−2 t ) (−3 t)
y ( t ) =h(t )=C 1 e +C 2 e

Supongamos que en t=0, la posición y la velocidad del sistema


son cero, es decir:
y ( 0 )=0 y ( 0 )=0

Sustituyendo estas condiciones iniciales en la solución


general, obtenemos:

0 = 𝐶1 + 𝐶2 0 = −2𝐶1 − 3𝐶2

Fig.26 Grafica del código de grafico Resolviendo este sistema de ecuaciones lineales, encontramos
que 𝐶1 = 3 y 𝐶2 = −3.

Por lo tanto, la solución de la ecuación diferencial a estado


cero y entrada cero es:

V. EJERCICIO 3 y ( t ) =3 e
(−2 t)
−3 e
(−3 t )

Un sistema LTI viene dado por el siguiente modelo en


operador diferencial (𝐷^2 + 5𝐷 + 6)𝑦(𝑡) = (𝐷 + 1)𝑢(𝑡)
VI. REFERENCIAS
a.) Encuentre el polinomio característico del sistema y [1] A. V. Oppenheim and A. S. Willsky, Signals and Systems, Prentice
determine sus raíces Hall, 1996.J. U.
[2] B. P. Lathi, Signal Processing and Linear Systems, 2nd
El sistema LTI dado por el modelo en operador ed., Berkeley-Cambridge Press, 2005.
diferencial: [3] S. T. Karris, "Signals and Systems: With MATLAB
Applications", Orchard Publications, 2010.
( D2 +5 D+6 ) y ( t )=( D+1 ) u(t) [4] R. Figueroa and A. Cornejo, Análisis complejo: teoría,
métodos y aplicaciones, 1ra edición, Pearson Educación,
Se puede representar en el dominio de la frecuencia como: 2010.
[5] J. R. Pineda, Señales y Sistemas, 3ra edición, Ecoe
Y(s) = H(s)U(s) Ediciones, 2018.

Para encontrar el polinomio característico del sistema, se


sustituye la variable D por s en la ecuación del modelo:

( s2 +5 s+6 ) y ( t )=( s+1 ) u(t)


Es necesario resolver la ecuación diferencial asociada
igualándola a 0 y así conocemos los coeficientes dónde
fácilmente podemos hallar sus raíces:

2
s +5 s +6=0
Por lo tanto, las raíces del polinomio
característico y del sistema LTI son R1 = -2 y R2 = -3.

b.) Construya la solución de la ecuación diferencial a


estado cero y entrada cero.
Se resuelve la ecuación diferencial homogénea asociada:
16

También podría gustarte