clear all; %Declaracion de los parametros: wn=1;%la frecuencia sera fija zeta=[0 0.5 1 2];% Variaremos zeta t=0:15/200:(15-15/200); %Por cada parametro zeta almacenaremos la respuesta %del sistema en una matriz y: for n=1:length(zeta) num=wn^2; den=[1 2*zeta(n)*wn wn^2]; W=tf(num,den); y(:,n)=step(W,t); end;
%Representacion en 2D de la respuesta figure(2) plot(t,y);xlabel('t');ylabel('y(t)'); title('Respuesta de un sistema de segundo orden variando \zeta'); legend('\zeta=0 (sin amortiguamiento)',... '\zeta=0.5 (Subamortiguado)',... '\zeta=1 (Criticamente Amortiguado)','\zeta=2 (Sobreamortiguado)');