syms p K3 T31 T32;% декларируем опреатор Лапласа 'р' и переменные коэффициенты модели как символьные переменные
K1=6;% вводим постоянные коэффициенты модели
K2=2;
K4=1;
T2=0.025;
T4=0.01;
W1=K1;% вводим передаточные функции (ПФ) модели
W2=K2/(T2*p+1);
W3=K3*(T31*p+1)/(T32*p^2+p);
W4=K4/(T4*p+1);
W2=simplify(W2)% упрощаем ПФ модели
W3=simplify(W3)
W4=simplify(W4)
Wzam=(W1*W2*W3)/(1+W1*W2*W3*W4)% находим ПФ замкнутой системы
[Bp,Ap]=numden(Wzam)% выводим числитель и знаменатель ПФ замкнутой системы
Bp=simplify(Bp)% упрощаем числитель и знаменатель ПФ замкнутой системы
Ap=simplify(Ap)
Wzam=Bp/Ap% выражаем через числитель и знаменатель ПФ замкнутой системы
pretty(Wzam)% упрощаем ПФ замкнутой системы
[Bp,Ap]=numden(Wzam)% выводим числитель и знаменатель упрощенной ПФ замкнутой системы
Bp=collect(Bp,p)% группируем полиномы числителя и знаменателя по множителям 'р^x'
Ap=collect(Ap,p)