(23:57:19) NiK: Михаил Гутник (01 июня, 00:00): k1=130;
k2=120;
k=0.7;
b1=4;
b2=5;
w01=13;
w02=8;
num1=[k1 0 (k1*w02^2)-(k2*w01^2)];
num2=[k2 0 (k2*w01^2)-(k1*w02^2)];
den1=[1 0 (w01^2)+(w02^2) 0 0];
den2=[1 0 (w01^2)+(w02^2) 0 0];
[num,den]=parallel(num1,den1,num2,den2);

w = 0:.01:50;
[RE,IM] = nyquist(num, den, w);
[A,F]=bode(num, den,w);

t = 0:0.01:10;
u = sin(20*t);
Y=lsim(num, den, u, t);

subplot(3,2,1);
plot(t,u);

title('Внешнее воздействие');
subplot(3,2,2);
plot(t,Y);

title('Выходной сигнал');
subplot(3,2,3);
plot(w,RE);

title('ВЧХ');
subplot(3,2,4);
plot(w,IM);

title('МЧХ');
subplot(3,2,5);
plot(w,A);

title('АЧХ');
subplot(3,2,6);
plot(w,F);

title('ФЧХ');