Идентификация объектов управления. Семенов А.Д - 111 стр.

w0=tf([-1 1],[1 1]);
ws=a(1);nun=a(1);
pp=[1 1];pm=[-1 1];
den=1;d=1;
% Вычисление передаточной функции объекта по (4.18)
for j=2:(n+1)/2
   den=conv(den,pp);
   d=conv(d,pm);
   nun=conv(nun,pp)+a(j)*d;
   ww=tf(nun,den)
end
ww=minreal(ww)% Получение минимальной реализации передаточной функции
step(ww,wo)
pause
bode (ww,wo)
pause
[wb,g]=balreal(ww); )% Получение сбалансированной реализации передаточной функции
wm=modred(wb,[4:(n-1)/2]);
step(wm,ww,wo)
ww=zpk(tf(wm))
      На рис. 4.6 показана вещественная частотная характеристика объекта для
относительных значений частоты ϕ = 2arctg (ω )




                                    Рис. 4.6.