基于MATLAB的异步电动机仿真

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

2、arm1=3;%电机相数U1=220;%额定定子相电压n1=-1000,1000;%两个相的同步转速,+为规定正相，-为反相p=3;%电机极对数f1=50;%输入电流频率r1=2.08;%定子侧电阻r2=1.53;%转子侧电阻折合值w1=2*pi*f1/p;%电机同步角速度*1=3.12;%定子侧电抗*2=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+(*1+*2)2);%计算转矩 plot(T,n,-); if i=1 te*t(ma*(T),800,strcat(r1=,num2str(r1),Omega),FontSize,9,Color,black); %曲线标注位置设置，每条曲线的标注横轴上对齐其Tm以便于区分，纵坐标无严格限制 end hold on;end*label(电磁转矩itT/(Ncdotm);ylabel(转速itn/(r/min);title(定子串电阻的机械特性曲线)grid on;hold on;%绘制定子串电阻人为机械特性曲线r1

4、p=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+(*1+*2)2); plot(T1,n,k-); hold on; if i=1 str=strcat(r1=,num2str(int16(r1p),Omega); str_y=str_y-150; te*t(ma*(T1),str_y,strcat(r1=,num2str(r1p),Omega),F

5、ontSize,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+r2./s).2+(*1+*2)2); plot(T,n,-); if i=1 te*t(ma*(T),800,strcat(U1=,num2str(int16(U1),V),FontSize,9,Color,black);end hold on;end*label(电磁转矩itT/(Ncdotm);ylabel(转速itn/(r/min

6、);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*U1p2*r2)./s./(w1.*(r1+r2./s).2+(*1+*2)2); plot(T1,n,k-); if i=1 str=strcat(U1=,num2str(int16(U1p),V); str_y=str_y-150; te*t(ma*(T1),str_y,str,FontSi

7、ze,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+(*1+*2)2); plot(T,n,-); if i=1 te*t(ma*(T),800,strcat(*1=,num2str(*1),Omega),FontSize,9,Color,black); end hold on;end*label(电磁转矩itT/(Ncdotm);ylabel(转速itn/(

8、r/min);title(定子串电抗的机械特性曲线)grid on;hold on;%绘制定子串电抗的人为机械特性曲线*1p=*1;str_y=800;for i=1:length(n1); n=n1(i).*(1-s); for coef=1:3 *1p=*1p+1.0; T1=ns(i)*(m1*p*U12*r2)./s./(w1.*(r1+r2./s).2+(*1p+*2)2); plot(T1,n,k-); if i=1 str=strcat(*1=,num2str(*1p),Omega); str_y=str_y-150; te*t(ma*(T1),str_y,strcat(*1=,

9、num2str(*1p),Omega),FontSize,9,Color,black); end hold on; endend%绘制转子串电阻前的固有机械特性曲线figure(4);for i=1:length(n1); n=n1(i)*(1-s); T=ns(i)*(m1*p*U12*r2)./s./(w1.*(r1+r2./s).2+(*1+*2)2); plot(T,n,-); if i=1 te*t(ma*(T),800,strcat(r2=,num2str(r2),Omega),FontSize,9,Color,black); end hold on;end*label(电磁转矩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+(*1