% 三相異步電動機的機械特性
%定義為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;
%下面是轉子電阻的循環數值
        R2=0.1;
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)
hold on
        R2=0.2;
        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)
hold on

        R2=0.5;
        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)
hold on

        R2=1.0;
        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)
hold on
        R2=1.5;
        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)
hold on


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