Analisis

Práctica de Laboratorio Nº8
Analisis de la respuesta en el tiempo

1. Objetivos.
Analizar la respuesta de de procesos



2. Procedimiento
Analizar la respuesta de de procesosSistemas de primer orden

K= 1;
tau = 1;
num = K;
den = [tau 1];


t = [0:0.1:10]’;
ye = step(num,den,t);
plot(t,ye);
title (’Respuesta a un escalon unitario’);
xlabel (’tiempo(seg)’);grid;
yRP = ye(length(ye)); % Valor en regimen permanente
n = 1;
while ye(n) < 0.63*yRP
n=n+1;
end
% Constante de tiempo (0.1 es el intervalo transcurrido entre dos medidas,
% se le resta 1, porquelos indices empiezan en 1):
tauEstim = 0.1*(n-1);




ramp = t;
yr = lsim (num,den,ramp,t);
plot (t,yr,t,ramp);
title (’Respuesta a una rampa’);
xlabel (’tiempo(seg)’);
grid;


yi =impulse (num,den,t);
plot (t,yi);
title (’Respuesta a un impulso’);
xlabel (’tiempo(seg)’);



ye2 = 1 - exp(-t/tau);
yr2 = t - tau + tau * exp(-t/tau);
yi2 = (1/tau) * exp(-t/tau);
plot(t,ye,t,ye2,’o’);
title(’Respuesta teorica a escalon unitario’);
grid;
pause;
plot (t,yr,t,ramp,t,yr2,’o’);
title(’Respuesta teorica a rampa unitaria’);
grid;
pause;
plot (t,yi,t,yi2,’o’);
title(’Respuesta teorica a impulso’);
grid;

Sistemas de segundo orden

Primer caso: 2 raíces reales distintas
t = [0:0.2:20]’;
wn = 1;
d = 2;
num = [wn^2];
den = [1,2*d*wn,wn^2];
ye = step(num,den,t);
plot (t,ye);
title (’Respuesta a un escalon unitario’);
xlabel (’tiempo(seg)’);
grid;

Segundo caso: 2 raíces reales iguales (criticamente amortiguado)

d = 1;
den = [1,2*d*wn,wn^2];
ye= step (num,den,t);
plot (t,ye);
title (’Respuesta a un escalon unitario’);
xlabel (’tiempo(seg)’);
grid;

Tercer caso: Dos raíces complejas conjugadas
d = 0.5;
den = [1,2*d*wn,wn^2];
ye =step (num,den,t);
plot (t,ye);
title (’Respuesta a un escalo’n unitario’);
xlabel (’tiempo(seg)’);
grid


Cuarto caso: Sistema con amotiguamiento crítico
d = 0;
den = [1,2*d*wn,wn^2];
ye...