[adsense:responsibe:9545213979]
k | 0 | 1 | 2 | 3 | 4 | 5 | 6 |
y(k) | 0 | 2 | 2.2 | 2.34 | 2.4106 | 2.4188 | 2.42 |
clear; z=%z; g=2*(z-0.2)/((z+0.1)*(z-0.4)); r=z/(z-1); gt=g*r gs=syslin('d',g); gts=syslin('d',gt); //Calculo y comprobacion de descomposicion fraciones simples a=horner((z+0.1)*gts,-0.1) b=horner((z-0.4)*gts,0.4) c=horner((z-1)*gts,1) gt2=a/(z+0.1)+b/(z-0.4)+c/(z-1) gf=tf2ss(gs); se=pfss(gf) //Calculo y comprobacion de la salida for i=1:1:10 y(i)=a*0.1^(i-1)*cos((i-1)*%pi)+b*0.4^(i-1)+c*1^(i-1); end; aux=[0 y']; k=0:1:40; x=ones(1,41); y2=flts(x,gs); [aux;y2] //Calculo del valor final y sobreoscilacion aux2=horner(gs,1) aux3=(max(y2)-aux2)/aux2 //Representacion grafica clf; plot(k,y2,'o'); plot(k,aux); mtlb_axis([0 10 0 2.1]); xgrid; xtitle('Respuesta a un escalon','k','y(k)')
Añadir nuevo comentario