файлы матлаб / lab33

clear,clc;
a11=-0.089;
a12=-2.19;
a13=0.328;
a14=0.319;
a15=0.0327;
a16=0.089;
b11=0.076;
b12=-0.217;
b13=-0.166;
b14=0.0264;
b15=-0.151;
b16=-0.076;
d11=-0.602;
d12=0.327;
d13=-0.975;
d14=0.227;
d15=0.0636;
d16=0.602;
i11=0.15;
f=2.21;

A=[a11, a12, 0, a13, a14, 0; b11, b12, 0, b13, 0, 0; 0, 1, 0, 0, 0, 0;...
d11, d12, 0, d13, 0, 0; 0, i11, 0, 1, 0, 0; 1, 0, f, 0, 0, 0];
B=[0, a15, a16; b14, b15, b16; 0, 0, 0; d14, d15, d16; 0, 0, 0; 0, 0, 0];
C=eye(6);
D=[0, 0, 0;0, 0, 0;0, 0, 0;0, 0, 0;0, 0, 0;0, 0, 0];
S1=ss(A,B,C,D);
DD=[-1.947, -3.59, -1.421, -1.672, -7.29, -0.859; 1.263, 6.42, 0.799, 1.424, 6.08, 0.487];
S2=ss(DD);
set(S1,'B',[B(:,3),B(:,1),B(:,2)]);
S=lft(S1,S2,2,6);
set(S,'C', [0, 0, 0, 0, 1, 0], 'D',[0]);
P=lyap(S.A,eye(size(S.A)));
q=[];
for i=1:6 %цикл составляет ряд из миноров матрицы P
q=[q det([P(1:i,1:i)])];
i=i+1;
end
rooots=roots(poly(S.A));
plot(real(rooots),imag(rooots),'*') %построение полюсов
grid on
hold on
t=[-2:2];
x=0*t;
plot(t,x,'black');
plot(x,t,'black'); % строит оси
hold off
k=[];
N=[];
for i=-10:0.1:100 % цикл создает ряд наибольших вещественных корней объекта REZ
DD2=[-1.947, -3.59, -1.421, -1.672, -7.29, -0.859; 1.263, 6.42, i, 1.424, 6.08, 0.487];
S3=ss(DD2);
REZ=lft(S1,S3,2,6);
k=[k i];
N=[N max(real(eig(REZ)))];
end
figure(2);
plot(k,N);
grid on