syms p w % задание переменных
Wlc=6/(0.5*p^2+p)*exp(-0.5*p) % ввод операторной передаточной функции линейной части системы
pretty(Wlc) % представление операторной передаточной функции линейной части в дробном виде
Wjw=subs(Wlc,p,j*w)% подмена оператора Лапласа на частотный комплексный оператор jw
pretty(Wjw) % представление частотной передаточной функции линейной части в дробном виде
imWjw0=subs(imag(Wjw),w,0) % проверка доп. условия Попова (только если есть нулевой корень)
Wjwmod=real(Wjw)+j.*w.*imag(Wjw) % вычисление модифицированной частотной передаточной функции линейной части
w1=[0.01:0.05:100]' % задание значений частоты сигнала - w
Wjw1=subs(Wjw,w,w1) % вычисление значений частотной передаточной функции линейной части
Wjwmod1=subs(Wjwmod,w,w1) % вычисление значений модифицированной частотной передаточной функции линейной части
plot(real(Wjw1),imag(Wjw1)) % построение частотного годографа линейной части
hold all % удержание построенного графика
plot(real(Wjwmod1),imag(Wjwmod1)) % построение модифицированного частотного годографа линейной части