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