Motor-Mechanical Features--电机机械特性计算

function Motor_Mechanical_features(Pn,n,n0,K) %Power,Synchronous speed,Rated speed,Overload ratio
s=linspace(0,1,1000);
Tn=9550*Pn/n;
Tm=K*Tn;
Sn=(n0-n)/n0;
Sm=Sn*(K+sqrt(K^2+1));
T1=2.*Tm./(s/Sm+Sm./s);
plot(T1,1-s);
title('Torsion-Speed difference');
xlabel('Tem');
ylabel('1-s');