% 三相異步電動機的機械特性
%定義為inductor_motor(xiugai)
clc
clear
% 下面輸入電動機參數
U1=220/sqrt(3);
Nphase=3;
P=2;
fN=50;
R1=0.095;
X1=0.680;
X2=0.672;
Xm=18.7;
%下面計算電機同步速度
omegas=2*pi*fN/P;
nS=60*fN/P;
%下面是轉子電阻的循環數值
        R21=0.1;
        hold 
 n=1:20
    s(n)=n/20;
    Tmech(n)=Nphase*P*U1*2*R21./n/omegas/[(R1+R21./n).^2+(X1+X2)^2]
    plot(s(n),Tmech(n));



hold 

%xlabel('轉差率')
%ylabel('電磁轉矩')