p = tf('p');% декларируем р как опреатор Лапласа
K1=6;% вводим коэффициенты модели
K2=2;
K3=1;
K4=1;
T2=0.025;
T31=0.1;
T32=0.2;
T4=0.01;
W1=tf(K1)% вводим передаточные функции модели
W2=tf(K2/(T2*p+1))
W3=tf(K3*(T31*p+1)/(T32*p+1)/p)
W4=tf(K4/(T4*p+1))
W123=series(series(W1,W2),W3)% находим передаточную функцию прямой ветви системы
Wzam=feedback(W123,W4)% находим передаточную функцию замкнутой системы
pole(Wzam)% вычисляем значения полюсов характеристического уравнения замкнутой системы
pzmap(Wzam)% выводим карту нулей и полюсов передаточной функции замкнутой системы