《基于MATLAB的异步电动机仿真.doc》由会员分享,可在线阅读,更多相关《基于MATLAB的异步电动机仿真.doc(8页珍藏版)》请在金锄头文库上搜索。
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);ylabel(转速itn/(r/min);title(定子串电阻的机械特性曲线)grid on;hold on;%绘制定子串电阻人为机械特性曲线r1p
4、=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=strcat(r1=,num2str(int16(r1p),Omega); str_y=str_y-150; text(max(T1),str_y,strcat(r1=,num2str(r1p),Omega),Fo
5、ntSize,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+(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
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+(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,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+(x1+x2)2); plot(T,n,-); if i=1 text(max(T),800,strcat(x1=,num2str(x1),Omega),FontSize,9,Color,black); end hold on;endxlabel(电磁转矩itT/(Ncdotm);ylabel(转速itn/(
8、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+(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=,
9、num2str(x1p),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+(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(电磁转矩i
10、tT/(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-); if i=1 str=strcat(r2=,num2str(r2p),Omega); str_y=str_y-150; text(max(T),str_y,str,FontSize,9,Color,black); end hold on; endend3. 仿真波形