% Исходные данные
J=0.35;
La=0.01;
%La=0.1;
Ra=0.3;
R_v=145;
U_c=220 ;
w=4000;
Ce=205;
Cm=200;
%Номинальные данные
F_nom=U_c*w/R_v;
Potok_nom=0.007;
w_nom=100;
i_nom=50;
M_nom=70;
a=[w*Potok_nom U_c R_v*F_nom/w La*i_nom Ce*Potok_nom*w_nom i_nom*Ra U_c J*w_nom M_nom Cm*Potok_nom*i_nom];
u=[1 1 1];
u4=[0 0 0];
x0=[0.983 1.0173 2.16];
% x0=[0 0 0];
c=[0.7479 -0.1896 0.5022];
A4=[-27 0 0; 287*x0(3) -30 287*x0(1); -2.001376*x0(2) -2*x0(1) 0];
B4=[U_c/(w*Potok_nom) -U_c/(La*i_nom) M_nom/(J*w_nom)];
r_A=eig(A4);
sim('laba3.slx');
sim('laba4.slx');
%Графики
figure('Color','w');
subplot(3,1,1);
hPlot=plot(x(:,1),x(:,4),x(:,1),x1(:,2),':');
set( hPlot, 'LineWidth', 1 );
xlabel('t');
ylabel('Ф')
legend ('СНДУ','СЛДУ' )
grid on;
subplot(3,1,2);
h2Plot=plot(x(:,1),x(:,3),x(:,1),x2(:,2),':');
set( h2Plot, 'LineWidth', 1 );
xlabel('t');
ylabel('i');
grid on;
subplot(3,1,3);
h3Plot=plot(x(:,1),x(:,2),x(:,1),x3(:,2),':');
set( h3Plot, 'LineWidth', 1 );
xlabel('t');
ylabel('w')
grid on;
%Графики
figure(1);
plot(real(r_A), imag(r_A),'*');
grid on;