clear
m1=3;
U1=220*sqrt(3);
R1=0.28;
R2=0.35;
p=2;
f=50;
omega=2*pi*f/p;
x1=1.19;
x2=1.19;
s=0.005:0.005:1
Te=(m1*p*U1^2*R2)./s./(omega.*((R1+R2./ s).^2.+(x1+x2)^2));
figure(1)
plot(s,Te,'k-');
xlabel('电磁转矩');
ylabel('转速');
strx=0.02;
text(strx,max(Te)+10,strcat('U1=',num2str(U 1),'V'),'color','black');
hold on;
for coef=0.75:-0.25:0.25
U1p=U1*coef;
Te1=(m1*p*U1p^2*R2)./s./(omega.*((R1+R 2./s).^2.+(x1+x2)^2));
plot(s,Te1,'k-');
str=strcat('U1=',num2str(U1p),'V');
stry=max(Te1)+10;
text(strx,stry,str,'color','black');
end
figure(2)
plot(s,Te,'k-');
xlabel('电磁转矩');
ylabel('转速');
strx=0.75;
stry=Te(length(Te))+10;
text(strx,stry,strcat('R2=',num2str(R2),'\omeg a'),'color','black');
hold on;
for coef=2:1:5
R2p=R2*coef;
Te1=(m1*p*U1^2*R2p)./s./(omega.*((R1+R 2p./s).^2.+(x1+x2).^2));
plot(s,Te1,'k-');
str=strcat('R2=',num2str(R2p),'\omega');
stry=Te1(length(Te1))+10;
text(strx,stry,str,'color','black');
end