异步电机工作特性仿真问题，请教！！！！

幽灵爱你

我编了一个异步电机工作特性的仿真，但是结果有问题，大家能不能帮帮看看。。。
clear;
Pn=160e3;U1=380;p=2;f=50;P2=100:100:Pn;
r1=0.0302;x1=0.123;
r2=0.0317;x2=0.1494;
rm=0.365;xm=5.3;
Pmecs=3.64e3;
x=x1+x2;
a=(3*U1^2*r2-2*r1*r2.*(P2+Pmecs))./(3*U1^2*r2+(r1^2+x^2).*(P2+Pmecs));
b=r2^2.*(P2+Pmecs)./(3*U1^2*r2+(r1^2+x^2).*(P2+Pmecs));
s=(a-sqrt(a.^2-4.*b))./2;
n=(1-s)*60*f/p;
Zm=sqrt((r1+rm)^2+(x1+xm)^2);
A=(r1*Zm^2+(r1+rm)*(r1^2+x^2)).*s.^2+r2*(Zm^2+2*r1*(r1+rm)).*s+(r1+rm)*r2^2;
B=(x*Zm^2+(x1+xm)*(r1^2+x^2)).*s.^2+2*r1*r2*(x1+xm).*s+(x1+xm)*r2^2;
C=Zm^2*((r1^2+x^2).*s.^2+2*r1*r2.*s+r2^2);
I1=U1.*sqrt(A.^2+B.^2)./C;
cosphi=cos(atan(B./A));
nn=P2/(3*U1*I1.*cosphi);
T2=P2*p./(2*pi*f.*(1-s));
plot(P2,I1,'G',P2,T2,'B',P2,nn,'R',P2,n,'Y',P2,cosphi,'k');

yokel

提问要注意点技巧，扔一大堆代码上来好歹加点注释嘛，另外问题的现象也不描述一下