Лаб3 / LR3

k=1.38;
n1=[1 2.07];
n2=[1 0.05 0.066];
d1=[1 0];
d2=[1 0.38 1.813];
d3=[1 2.09];
d4=[1 -0.004];
num=k*conv(n1,n2);
den=conv(conv(d1,d2),conv(d3,d4));
F=tf(num,den)
F1=zpk(F)


da=[1 1];
dena=da;
numa=10;
Fa=tf(numa,dena);
C=feedback(F,Fa)
C1=zpk(C)
Css=ss(C);


% B=eye(6)+2*inv(Css.a-eye(6));
% my=[];
% K=1000;
% for (k=1:1:K)
% my(k)=norm(B^k);
% end
% my(k)-my(k-1)

Smax=[]; Ka=[];
i=1;
for (ka=100:-0.1:0)
numa=ka;
Fa=tf(numa,dena);
C=feedback(F,Fa);
Css=ss(C);
smax=max(real(eig(Css.a)));
Ka(i)=ka;
Smax(i)=smax;
i=i+1;
end
plot(Ka,Smax); grid on; ylabel('Smax'); xlabel('Ka');