MSEI 5

nblocks=4;
n1=1;d1=1;
n2=1;d2=1;
n3=3;d3=[0.01 3];
n4=1;d4=[1 2 3];
blkbuild;
q=[1 0 0; 2 0 0; 3 1 -4; 4 2 3];
iu=[1 2];
iy=[4];
[A,B,C,D]=connect(a,b,c,d,q,iu,iy);
[am,bm,cm,dm]=minreal(A,B,C,D);
[y,x,timp]=step(am,bm,cm,dm,1);
plot(timp,y,'r');

%%timpul primului maxim
n=length(timp);
for j=1:n
    if y(j)==max(y)
    t_primMax=timp(j)
    jmax=j;
    end
end

%%suprareglajul
f=max(y)-y(n)

%%suprareglajul exprimat in procente
fprocent=f*100/y(n)

%% timpul de reglare
for i=1:jmax
    if y(i+1)>0.95*y(n) & y(i)<0.95*y(n)
         t_reglare=timp(i+1)
    end
end

k=0;

for i=1:n-1
    if y(i+1)<=y(n) & y(i)>y(n) & k==0
        jmax2=i;k=1;
    end
end

for i=jmax2:n
    v(i)=y(i);
end

g_amortizare=1-(max(v)-y(n))/(max(y)-y(n))