w=logspace(-1,2,1000);	% Vector de frecuencias
% Análisis de la influencia de la ganancia en la respuesta en frecuencia de un sistema de 2º orden
psi = 0.1; wn=2;
for k=[2 1 0.5];
    num=[k*wn^2]; den=[1 2*psi*wn wn^2];
    
    figure(1), hold on
    bode(num,den,w), grid minor
    title('Influencia de la ganancia en la respuesta en fecuencia de un sistema de 2º orden');
    legend('K=2','K=1','K=0.5');
    hold off
    
    figure(2), hold on
    step(num,den);grid minor
    title('Influencia de la ganancia en la respuesta en el tiempo de un sistema de 2º orden');
    legend('K=2','K=1','K=0.5');
    hold off
    
    disp('Pulsa una tecla para continuar'), pause
end;
close all

% Análisis de la influencia del coeficiente de amortiguamiento en la respuesta en frecuencia de un sistema de 2º orden
wn=2; k=10;
for psi=[0.1 0.3 0.5 0.707 0.9 1]
    num=[k*wn^2]; den=[1 2*psi*wn wn^2];
    
    figure(1), hold on
    bode(num,den,w), grid minor
    title('Influencia de \zeta en la respuesta en fecuencia de un sistema de 2º orden');
    legend('\zeta=0.1','\zeta=0.3','\zeta=0.5','\zeta=0.707','\zeta=0.9','\zeta=1');
    hold off
    
    figure(2);, hold on
    step(num,den);grid minor
    title('Influencia de \zeta en la respuesta en el tiempo de un sistema de 2º orden');
    legend('\zeta=0.1','\zeta=0.3','\zeta=0.5','\zeta=0.707','\zeta=0.9','\zeta=1');
    hold off
    
    disp('Pulsa una tecla para continuar'), pause
end;
close all


% Análisis de la influencia de la frecuencia antural en la respuesta en frecuencia de un sistema de 2º orden
psi=0.1; k=10;
for wn=[1 3 5 10]
   num=[k*((wn)^2)]; den=[1 2*psi*wn wn^2];
   
    figure(1), hold on
    bode(num,den,w), grid minor
    title('Influencia de \omega_n en la respuesta en fecuencia de un sistema de 2º orden');
    legend('\omega_n=1','\omega_n=3','\omega_n=5','\omega_n=10');
    hold off
    
    figure(2);, hold on
    step(num,den);grid minor
    title('Influencia de \omega_n en la respuesta en el tiempo de un sistema de 2º orden');
    legend('\omega_n=1','\omega_n=3','\omega_n=5','\omega_n=10');
    hold off
    
    disp('Pulsa una tecla para continuar'), pause
end;
close all
% Adición de polos 1/(1+ps) a un sistema de 2º orden
wn=2;psi=0.1; k=10;			

for p=[0 1 2 4]
   num=[k*wn^2]; den=[p (1+p*2*psi*wn) wn*(wn*p+2*psi) wn^2];
   
    figure(1), hold on
    bode(num,den,w), grid minor
    title('Influencia de un polo 1/(1+pS) en la respuesta en frecuencia de un sistema de 2º orden');
    legend('p=0','p=1','p=2','p=4');
    hold off
    
    figure(2);, hold on
    step(num,den);grid minor
    title('Influencia de un polo 1/(1+pS) en la respuesta en el tiempo de un sistema de 2º orden');
    legend('p=0','p=1','p=2','p=4');
    hold off
    
    disp('Pulsa una tecla para continuar'), pause
end;
close all
% Adición de cero (1+as) a un sistema de 2º orden
wn=2;psi=0.1;k=10;			
for a=[0 1 2 4]

    num=[a*wn^2*k k*wn^2]; den=[1 2*psi*wn wn^2];
    
    figure(1), hold on
    bode(num,den,w), grid minor
    title('Influencia de un cero (1+aS) en la respuesta en frecuencia de un sistema de 2º orden');
    legend('a=0','a=1','a=2','a=4');
    hold off
    
    figure(2);, hold on
    step(num,den);grid minor
    title('Influencia de un cero (1+aS) en la respuesta en el tiempo de un sistema de 2º orden');
    legend('a=0','a=1','a=2','a=4');
    hold off
    
    disp('Pulsa una tecla para continuar'),pause
end;
close all