Convolución discreta
Publicado por David (1 intervención) el 17/03/2013 18:12:51
Hola tengo un programa y tengo unas dudas a ver si me podeis orientar un poco.
Entiendo que es una convolución entre un sin(t) y una función impulso.
Tengo que modificar h para que el sistema funcione como un sistema de primer orden con 3 segundos de tiempo de establecimiento.
Después tengo que ejecutar el programa modificado y analizar el resultado obtenido cuando la señal de entrada es un escalón que comienza en tiempo t=2 y tiene una amplitud de 5 unidades.
Estoy muy verde en Matlab, gracias de antemano.
t_sim=20; % Tiempo total simulación
T=0.001; % Periodo de muestreo
t=(0.0:T:t_sim); % Instantes en los que conoce la señal muestreada
nk=length(t); % Número de muestras de la señal de entarda
u=sin(t); % Señal de entrada muestreada
h=[100 -100]; % Núcleo de convolución (respuesta ante un impulso)
n=length(h); % Longitud del núcleo convolución
% Calculo de la respuesta para cada punto de la señal de entrada
for k=(1+n:nk) % En cada iteración procesamos un punto de la entrada
% Calculo de la convolución para cada punto de la señal
y(k)=0;
for i=(1:n) % Para cada punto del núcleo de convolución
y(k)=y(k)+h(i)*u(k-i+1);
end
end
plot(t,u,t,y);
T=0.001; % Periodo de muestreo
t=(0.0:T:t_sim); % Instantes en los que conoce la señal muestreada
nk=length(t); % Número de muestras de la señal de entarda
u=sin(t); % Señal de entrada muestreada
h=[100 -100]; % Núcleo de convolución (respuesta ante un impulso)
n=length(h); % Longitud del núcleo convolución
% Calculo de la respuesta para cada punto de la señal de entrada
for k=(1+n:nk) % En cada iteración procesamos un punto de la entrada
% Calculo de la convolución para cada punto de la señal
y(k)=0;
for i=(1:n) % Para cada punto del núcleo de convolución
y(k)=y(k)+h(i)*u(k-i+1);
end
end
plot(t,u,t,y);
Entiendo que es una convolución entre un sin(t) y una función impulso.
Tengo que modificar h para que el sistema funcione como un sistema de primer orden con 3 segundos de tiempo de establecimiento.
Después tengo que ejecutar el programa modificado y analizar el resultado obtenido cuando la señal de entrada es un escalón que comienza en tiempo t=2 y tiene una amplitud de 5 unidades.
Estoy muy verde en Matlab, gracias de antemano.
Valora esta pregunta
0