%начальные значения
k=1;w=1;
om=[];B=[]
%параметры передаточной функции
k1=2;k2=0.5;T=k2/k1;
while w<=25,
 s=i*w;
 %вид передаточной функции
 f=k1/s+k2;
 %массив значений для фазовой характеристики
 B(k)=-pi/2+atan(T*w);
 om(k)=w;
 w=w+0.01;
 k=k+1;
end;
%построение графика
plot(om,B,'-k');
grid off