[adsense:responsibe:9545213979]
s=%s/(2*%pi); g=4/(s*(s^2+1)); gs=syslin('c',g); clf(); bode(gs,0.1,10);
s=%s/(2*%pi); g=4/(s*(s^2+1)); gcad=(5*s+1); gs=syslin('c',g); gs2=syslin('c',g*gcad); clf(); bode(gs2,'compensado');
s=%s/(2*%pi); g=4/(s*(s^2+1)); gcad=(5*s+1); gcad2=(0.25*s+1); gs=syslin('c',g); gs2=syslin('c',g*gcad*gcad2); clf(); bode(gs2,0.01,100,'compensado');
s=%s; g=4/(s*(s^2+1)); gcad=(5*s+1); gcad2=(0.25*s+1); gt1=g*gcad*gcad2; glc=g /. 1; glc2=gt1 /. 1; gs=syslin('c',glc); gs2=syslin('c',glc2); t=0:0.1:14; y=csim('step',t,gs); y2=csim('step',t,gs2); clf(); subplot(2,1,1); plot(t,y,'g'); legends('no compensado',3,opt=1) xgrid; xtitle('Respuesta ante una entrada escalon','','y(t)') subplot(2,1,2) plot(t,y2,'b'); legends('compensado',2,opt=4); xtitle('','t(seg)','y(t)') xgrid;
s=%s; g=4/(s*(s^2+1)); gcad=(5*s+1); gcad2=(0.25*s+1); gt1=g*gcad*gcad2; glc=g /. 1; glc2=gt1 /. 1; gs=syslin('c',glc); gs2=syslin('c',glc2); t=0:0.1:20; y=csim(t,t,gs); y2=csim(t,t,gs2); clf(); subplot(2,1,1); plot(t,y,'g'); legends('no compensado',3,opt=1) xgrid; xtitle('Respuesta ante una rampa','','y(t)') subplot(2,1,2) plot(t,y2,'b'); plot(t,t,'r'); legends(['compensado';'rampa'],[2;5],opt=4); xtitle('','t(seg)','y(t)') xgrid;
Añadir nuevo comentario