- 积分
- 21
- 回帖
- 0
- 西莫币
-
- 贡献
-
- 威望
-
- 存款
-
- 阅读权限
- 5
- 最后登录
- 1970-1-1
该用户从未签到
|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?立即注册
×
我编了一个异步电机工作特性的仿真,但是结果有问题,大家能不能帮帮看看。。。
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'); |
|