基于MATLAB的异步电动机仿真

1、课程报告COURSE REPORT课程名称: 课 程 号: 授课教师: 学 号： 姓 名：西木小卒 所 属：上大桂树林基于MATLAB的异步电动机仿真1. 电机仿真模型一台三相六机鼠笼式异步电动机，定子绕组为Y型连接，额定电压为UN=380V，额定转速nN=975r/min，电源频率f1=50Hz，定子电阻R1=2.08，定子漏电抗X1=3.12，转子电阻折合值R2、=1.53，转子漏电抗折合值X2、=4.25。要求：绘制以上参数电动机的固有机械特性曲线、定子串电阻人为特性曲线、电子串电抗人为特性曲线、转子串电阻人为特性曲线、降电压人为特性曲线；给出仿真源代码。2. 仿真代码实现clcclea

2、rm1=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(

3、i)*(m1*p*U12*r2)./s./(w1.*(r1+r2./s).2+(x1+x2)2);%计算转矩 plot(T,n,'-'); if i=1 text(max(T),800,strcat('r1=',num2str(r1),'Omega'),'FontSize',9,'Color','black'); %曲线标注位置设置，每条曲线的标注横轴上对齐其Tm以便于区分，纵坐标无严格限制 end hold on;endxlabel('电磁转矩itT/(Ncdotm)');ylab

4、el('转速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*U12*r2)./s./(w1.*(r1p+r2./s).2+(x1+x2)2); plot(T1,n,'k-'); hold on; if i=1 str=strca

5、t('r1=',num2str(int16(r1p),'Omega'); str_y=str_y-150; text(max(T1),str_y,strcat('r1=',num2str(r1p),'Omega'),'FontSize',9,'Color','black'); end endend%绘制降电压前的固有机械特性曲线figure(2);for i=1:length(n1); n=n1(i)*(1-s); T=ns(i)*(m1*p*U12*r2)./s./(w1.*(r1

6、+r2./s).2+(x1+x2)2); plot(T,n,'-'); if i=1 text(max(T),800,strcat('U1=',num2str(int16(U1),'V'),'FontSize',9,'Color','black'); end hold on;endxlabel('电磁转矩itT/(Ncdotm)');ylabel('转速itn/(r/min)');title('降低定子电压的机械特性曲线')grid on;hold o

7、n;%绘制降电压人为机械特性曲线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*U1p2*r2)./s./(w1.*(r1+r2./s).2+(x1+x2)2); plot(T1,n,'k-'); if i=1 str=strcat('U1=',num2str(int16(U1p),'V'); str_y=str_y-150; text(max(T1),str_y,str,'FontSiz

8、e',9,'Color','black'); end hold on; endend%绘制定子串电抗前的固有机械特性曲线figure(3);for i=1:length(n1); n=n1(i)*(1-s); T=ns(i)*(m1*p*U12*r2)./s./(w1.*(r1+r2./s).2+(x1+x2)2); plot(T,n,'-'); if i=1 text(max(T),800,strcat('x1=',num2str(x1),'Omega'),'FontSize',9,

9、9;Color','black'); end hold on;endxlabel('电磁转矩itT/(Ncdotm)');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:3 x1p=x1p+1.0; T1=ns(i)*(m1*p*U12*r2)./s./(w1.*(r1+r2./s).2+(

10、x1p+x2)2); plot(T1,n,'k-'); if i=1 str=strcat('x1=',num2str(x1p),'Omega'); str_y=str_y-150; text(max(T1),str_y,strcat('x1=',num2str(x1p),'Omega'),'FontSize',9,'Color','black'); end hold on; endend%绘制转子串电阻前的固有机械特性曲线figure(4);for i=1:leng

11、th(n1); n=n1(i)*(1-s); T=ns(i)*(m1*p*U12*r2)./s./(w1.*(r1+r2./s).2+(x1+x2)2); plot(T,n,'-'); if i=1 text(max(T),800,strcat('r2=',num2str(r2),'Omega'),'FontSize',9,'Color','black'); end hold on;endxlabel('电磁转矩itT/(Ncdotm)');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:3 r2p=r2p+0.75; T1=ns(i)*(m1*p*U12*r2p)./s./(w1.*(r1+r2p./s).2+(x1+x2)2); plot(T1,n,'k-')