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