% 三相異步電動機的機械特性

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;
%下面是轉子電阻的循環數值
for m=1:5
    if m==1
        R2=0.1;
    elseif m==2
        R2=0.2;
        elseif m==3
        R2=0.5;
        elseif m==4
        R2=1.0;
    else 
        R2=1.5;
    end
%下面是轉差率計算
for n=1:2000
    s(n)=n/2000;
    Tmech=Nphase*P*U1*2*R2/s(n)/omegas/[(R1+R2/s(n))^2+(X1+X2)^2];
%繪圖
plot(s(n),Tmech)
end
if m==1
    hold
end
end
hold
xlabel('轉差率')
ylabel('電磁轉矩')