Приложение 1. Программа построения АФЧХ и ЛФХ.
% lax.m
function lax
K0 = 100;
T1 = 1;
T2 = 4;
Kp = 2;
Tp = 0.2;
num = [K0*Kp*Tp Kp*K0];
den = [T1*T2 T1+T1*T2+T2 T1+T2+1 1];
sys = tf(num, den);
figure(2)
nyquist(sys)
title('АФЧХ для схемы 1')
grid on
figure(3)
asymp(sys)
title('ЛАХ и ЛФХ для схемы 1')
Приложение 2. Программа кривой Михайлова.
% кривая Михайлова
K0 = 100;
T1 = 1;
T2 = 4;
Kp = 2;
Tp = 0.2;
% разомкнутая система
num=[T1 T1+T1*T2+T2 T1+T2+1 1];
den=[1];
w=0:0.01:3;
apk=freqs(num,den,w);
u=real(apk);
v=imag(apk);
plot(u,v)
xlabel('Re')
ylabel('Im')
grid on
hold on
% замкнутая система
figure(2)
num=[T1 2*T1+T2 T1+T2+1+K0*Kp*Tp 1+K0*Kp];
den=[1];
w=0:0.01:3;
apk=freqs(num,den,w);
u=real(apk);
v=imag(apk);
plot(u,v)
xlabel('Re')
ylabel('Im')
grid on
hold on