Convolucion

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

Practica N°7

Respuesta al impulso y convolución


Maicol Nava, Brayan Ramírez
Fundación Universitaria Los Libertadores

Abstract Se introduce en Matlab por medio de


polinomios.
En este laboratorio el objetivo es
obtener la respuesta de un sistema en el >>y= [1];
tiempo ante una entrada estándar, >>x= [1 0.25 1];
para esto usaremos Matlab donde se Para obtener la respuesta del sistema en el
define el sistema mediante una función tiempo para una entrada de escalón se usa
de transferencia propia del sistema o el comando step
de ecuaciones de estado. >>step(y,x)
MatLab despliega la respuesta en el
Palabras claves
tiempo en la ventana de figuras: Figura
Matlab, convolución, Respuesta, Función N° 1.
de transferencia Para obtener la respuesta en el tiempo
para una entrada impulso unitario se usa
I. Introducción el comando impulse, con sintaxis idéntica
En este documento se desarrolla la a la utilizada con el comando step:
respuesta de un sistema frente a diferentes Si se define el sistema en MatLab por los
entradas y luego se hace la convolución polinomios del numerador y denominador
de dos funciones de la función de transferencia tenemos:
» y=[1 5 4];
De esta manera se ayuda a comprender el » x=[1 6 11 6];
funcionamiento de la respuesta y la » impulse(y,x)
convolución.
MatLab despliega la respuesta en el
II. Desarrollo tiempo en la ventana de figuras (Fig. N°
2).
La función de transferencia es la relación
formada por un numerador y un
denominador, en Matlab debe definirse el
numerador y el denominador como
vectores, cuyos elementos son los
coeficientes de los polinomios y el
denominador en potencias decrecientes de
S.
Tomando como ejemplo la función de
transferencia:
1
H ( s )= 2
s + 0.25 s+ 1
Fig. 1. Respuesta al escalón

Fig. 3. Respuesta a la rampa

Fig. 2. Respuesta al impulso Por otro lado, se puede hacer una


aproximación de Dirac mediante la
MatLab permite, además de obtener la función de Matlab Sinc lo cual se puede
respuesta en el tiempo para una entrada observar en la figura 3.
escalón o impulso, también obtener Escalar y aumentar la amplitud del seno
respuesta para otras entradas tal como cardinal por 100:
rampas o sinusoides. El comando lsim
>> t=-10:0.1:10;
permite obtener la respuesta en el tiempo
para un sistema con una entrada x, donde >> diracsinc=100*sinc(100*t);
x se define como una función del tiempo. >> plot(t,diracsinc)
La sintaxis de este comando es:
𝑙𝑠i𝑚(𝑁𝑈𝑀, 𝐷𝐸𝑁, X, 𝑇) usando la función >> title('Dirac aproximado con un seno cardinal')
de transferencia. >> xlabel('t')
Para obtener la respuesta en el tiempo >> ylabel('100*sinc(1oo*t)')
para una función rampa, se define X de
la siguiente forma:
>>T=0:0.1:10
>>X=T;
>>NUM=[1];
>>DEN=[1 0.25 1];
>>[Y,X]=lsim(NUM,DEN,X,T);
>>PLOT(T,Y,T,X)
Al hacer X=T se está definiendo la
función rampa. T es el vector de tiempo
variando desde 0 hasta 10 seg. En la
variable Y se almacena la salida del
sistema en función del tiempo T. El
comando plot permite desplegar en la
ventana de figuras la variable Y (salida) y
la entrada X (rampa) en función del
tiempo. Fig. 4. Aproximación de un delta de Dirac
Convolución:

Realizar la convolución de 2 pulsos entradas, por otro lado, también


rectangulares. El resultado es un pulso podemos hacer un acercamiento al
triangular. delta de Dirac que como podemos
𝑦 = 𝑟𝑒𝑐𝑡(𝑡) * 𝑟𝑒𝑐𝑡(𝑡) observar que esta correcto ya que
tiende a infinito en x=0 y por ultimo
Equivalentemente: tenemos la convolución obtenida de
las gráficas la cual demostró el pico

que forman.
𝑦(𝑡) = ∫ rect (r )rect (t−r ) dr
−∞ Conclusiones
 MatLab permite el análisis de
t=-10:0.01:10; funciones de transferencia y el
x1=3*exp(-t).*rectpuls(t-5,10); tratamiento de datos por medio
subplot(3,1,1)
plot(t,x1) de convolución.
ylabel('x1(t)')  Se comprobó que los resultados
x2=(t/2).*rectpuls(t-1,2) obtenidos son correctos por
subplot(3,1,2)
plot(t,x2) medio de las gráficas obtenidas
ylabel('x2(t)')  Se obtuvo la convolución de
y=conv(x1,x2) más de una función lo cual
t=-10:0.01/2:10;
subplot(3,1,3) facilita el análisis de la gráfica.
plot(t,y)
title('Convolucion de dos funciones')
xlabel('t')

Fig. N° 4. Convolución de dos funciones

III. Análisis de resultados


Se puede observar en la figura 1,2 y 3
que por medio de Matlab se puede
observar como responder un sistema a
diferentes tipos de entradas, además
de ser capaz de graficar las diferentes

También podría gustarte