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-')