基于MATLAB的异步电动机仿真

课程报告COURSE REPORT课程名称: ——课程号: ——授课教师: ——学号：——姓名：西木小卒所属：上大桂树林基于MATLAB的异步电动机仿真1.电机仿真模型一台三相六机鼠笼式异步电动机，定子绕组为Y型连接，额定电压为U N=380V，额定转速n N=975r/min，电源频率f1=50Hz，定子电阻R1=2.08Ω，定子漏电抗X1=3.12Ω，转子电阻折合值=1.53Ω，转子漏电抗折合值=4.25Ω。

要求：绘制以上参数电动机的固有机械特性曲线、定子串电阻人为特性曲线、电子串电抗人为特性曲线、转子串电阻人为特性曲线、降电压人为特性曲线；给出仿真源代码。

2. 仿真代码实现clcclearm1=3;%%电机相数U1=220;%%额定定子相电压n1=[-1000,1000];%%两个相的同步转速,+为规定正相，-为反相p=3;%%电机极对数f1=50;%%输入电流频率r1=2.08;%%定子侧电阻r2=1.53;%%转子侧电阻折合值w1=2*pi*f1/p;%%电机同步角速度x1=3.12;%%定子侧电抗x2=4.25;%%转子侧折合电抗s=-1:0.001:1;%%设定转差率，是画图的关键自变量ns=[-1,1];%用来标定转矩的方向，规定ns=1为正，-1为反%%绘制定子串电阻前的固有机械特性曲线figure(1);for i=1:length(n1);n=n1(i)*(1-s);%%计算转速T=ns(i)*(m1*p*U1^2*r2)./s./(w1.*((r1+r2./s).^2+(x1+x2)^2));%%计算转矩plot(T,n,'--');if i==1text(max(T),800,strcat('r1=',num2str(r1),'\Omega'),'FontSize',9,'Color','black');%%曲线标注位置设置，每条曲线的标注横轴上对齐其Tm以便于区分，纵坐标无严格限制endhold on;endxlabel('电磁转矩{\itT}/(N\cdotm)');ylabel('转速{\itn}/(r/min)');title('定子串电阻的机械特性曲线')grid on;hold on;%%绘制定子串电阻人为机械特性曲线r1p=r1;str_y=800;for i=1:length(n1);%%绘制正反相序的波形图，n=n1(i).*(1-s);for coef=1:3%%三条人为曲线r1p=r1p+1;%%串电阻的步进值T1=ns(i)*(m1*p*U1^2*r2)./s./(w1.*((r1p+r2./s).^2+(x1+x2)^2));plot(T1,n,'k-');hold on;if i==1str=strcat('r1=',num2str(int16(r1p)),'\Omega');str_y=str_y-150;text(max(T1),str_y,strcat('r1=',num2str(r1p),'\Omega'),'FontSize',9,'Color','black' );endendend%%绘制降电压前的固有机械特性曲线figure(2);for i=1:length(n1);n=n1(i)*(1-s);T=ns(i)*(m1*p*U1^2*r2)./s./(w1.*((r1+r2./s).^2+(x1+x2)^2));plot(T,n,'--');if i==1text(max(T),800,strcat('U1=',num2str(int16(U1)),'V'),'FontSize',9,'Color','black');endhold on;endxlabel('电磁转矩{\itT}/(N\cdotm)');ylabel('转速{\itn}/(r/min)');title('降低定子电压的机械特性曲线')grid on;hold on;%%绘制降电压人为机械特性曲线str_y=800;for i=1:length(n1);n=n1(i).*(1-s);for coef=0.75:-0.25:0.25;U1p=U1*coef;T1=ns(i)*(m1*p*U1p^2*r2)./s./(w1.*((r1+r2./s).^2+(x1+x2)^2));plot(T1,n,'k-');if i==1str=strcat('U1=',num2str(int16(U1p)),'V');str_y=str_y-150;text(max(T1),str_y,str,'FontSize',9,'Color','black');endhold on;endend%%绘制定子串电抗前的固有机械特性曲线figure(3);for i=1:length(n1);n=n1(i)*(1-s);T=ns(i)*(m1*p*U1^2*r2)./s./(w1.*((r1+r2./s).^2+(x1+x2)^2));plot(T,n,'--');if i==1text(max(T),800,strcat('x1=',num2str(x1),'\Omega'),'FontSize',9,'Color','black');endhold on;endxlabel('电磁转矩{\itT}/(N\cdotm)');ylabel('转速{\itn}/(r/min)');title('定子串电抗的机械特性曲线')grid on;hold on;%%绘制定子串电抗的人为机械特性曲线x1p=x1;str_y=800;for i=1:length(n1);n=n1(i).*(1-s);for coef=1:3x1p=x1p+1.0;T1=ns(i)*(m1*p*U1^2*r2)./s./(w1.*((r1+r2./s).^2+(x1p+x2)^2));plot(T1,n,'k-');if i==1str=strcat('x1=',num2str(x1p),'\Omega');str_y=str_y-150;text(max(T1),str_y,strcat('x1=',num2str(x1p),'\Omega'),'FontSize',9,'Color','black' );endhold on;endend%%绘制转子串电阻前的固有机械特性曲线figure(4);for i=1:length(n1);n=n1(i)*(1-s);T=ns(i)*(m1*p*U1^2*r2)./s./(w1.*((r1+r2./s).^2+(x1+x2)^2));plot(T,n,'--');if i==1text(max(T),800,strcat('r2=',num2str(r2),'\Omega'),'FontSize',9,'Color','black');endhold on;endxlabel('电磁转矩{\itT}/(N\cdotm)');ylabel('转速{\itn}/(r/min)');title('转子串电阻的机械特性曲线')grid on;hold on;%%绘制转子串电阻的人为机械特性曲线r2p=r2;str_y=800;for i=1:length(n1);n=n1(i).*(1-s);for coef=1:3r2p=r2p+0.75;T1=ns(i)*(m1*p*U1^2*r2p)./s./(w1.*((r1+r2p./s).^2+(x1+x2)^2));plot(T1,n,'k-');if i==1str=strcat('r2=',num2str(r2p),'\Omega');str_y=str_y-150;text(max(T),str_y,str,'FontSize',9,'Color','black');endhold on;endend3.仿真波形.. . .. ...。