哈工大 电机学 MATLAB 仿真 实验报告.docx

《哈工大 电机学 MATLAB 仿真 实验报告.docx》由会员分享，可在线阅读，更多相关《哈工大 电机学 MATLAB 仿真 实验报告.docx（11页珍藏版）》请在第一文库网上搜索。

1、基于MAT1AB的电机学计算机辅助分析与仿真实验报告班级：学号：姓名：完成时间：HARBININSTITUTEOFTECHNO1OGY一、实验内容1.1 单相变压器不同负载性质的相量图通过MAT1AB画出单相变压器带感性，阻性，容性三种不同性质负载的变压器向量图1.2 感应电机的T-S曲线通过MAT1AB画出三相感应电动机的转矩转差率曲线二、实验要求2.1 单相变压器不同负载性质的相量图根据给定的仿真实例画出负载相位角上=-30,0,30三种情况下得向量图，观察电压大小与相位的关系，了解总结负载性质不同对向量图的影响2.2 感应电机的T-S曲线根据给定的实例，画出s=T.31.3的T-S曲线，

2、了解感应电机临界转差率的大小和稳定工作区间的大小，给出定性分析三、实验方法3.1 单相变压器不同负载性质的相量图1 .单相变压器不同负载性质的相量图(1)先画出负载电压的相量；(2)根据负载的性质和阻抗角画出二次电流(规算值)的相量(3)在巴上加上一个与电流方向相同的压降，其大小为二次电流规算值4与二次漏电阻规算值R2之积；再加上一个超前电流方向90。的压降，其大小为二次电流I2规算值与二次漏电抗规算值Z2之积；(4)根据上一步结果连线，得出耳；超前反方向90。画出加；(6)根据励磁电阻与电抗的大小得出励磁阻抗角，并超前中加一个励磁阻抗角的大小得出的方向；(7)根据平行四边形法则，做出与，的和

3、，即为4；根据片=瓦得出用，并得出一与。(9)在-4上加上一个与电流方向相同的压降，其大小为一次电流人与一次漏电阻与之积；再加上一个超前电流方向90。的压降，其大小为一次电流人与一次漏电抗小之积；(10)根据上一步结果连线，得出G；3.2 感应电机的7-S曲线实验采用matIab对转矩转差率曲线进行仿真。由转矩转差率关系公式知，CS（+c&）2+（玉+扇）2S只有S为自变量，其他参数均为已知。编程时，先取S在0.01T.3正区间的r-S,进行绘图；再取相应负区间对T-S绘图；最后加入（0,0）四、实验源程序(1分)4.1 单相变压器不同负载性质的相量图见附录4.2 感应电机的TS曲线%T-S曲

4、线绘制%定义常量R2=0.04；R1=0.06;M1=3;U1=380;W=2*pi*1485/60;X1=0.27;X2=0.56;C=1+X1/16.4；%画出s=0.011.3的T-S曲线s=0.01:0.01:1.3;T=ones(1,Iength(s);fori=1:Iength(s)T2(i)=R2s(i);T1(i)=1(R1+C*T2(i)2+(X1+C*X2)2)；T3(i)=1/W；T(i)=M1*U2*T1(i)*T2(i)*T3(i)；end%画出s=-0.01.3的T-S曲线s1=-1.3:0.01:-0.01;T1=ones(1,1ength(s1);fori=1:

5、130T21(i)=R2s1(i);T11(i)=1(R1+C*T21(i)2+(X1+C*X2)2)；T31(i)=1/W；T1(i)=M1*2*T11(i)*T21(i)*T31(i);end%将(0,0)加入曲线中，合并成完整的向量ss=s10s;TT=T10T;pIot(ss,TT);%改变X轴方向set(gca,XDir,rev,);titIe(感应电机转矩转差率曲线)；XIabeI(,rotatingspeeds1ip(%),)yIabeI(,Torque(Nm),)五、实验结果5.1 单相变压器不同负载性质的相量图不同性质负载对相量图带来的影响：(1)二次电流I2在负载性质为感性

6、、容性时分别超前、滞后负载电压U2一个阻抗角的大小；纯阻性负载时I2与U2同方向；(2)在容性负载时，负载电压U2有可能大于感应电动势E2。感应电机的TS曲线实验结果如下：通过MatIab仿真得到感应电机的7-S曲线如下图；可以看出实验的曲线和理论预测基本相符，实验的仿真结果较为成功。曲线的特点：(1)曲线不是单调的而是像一个山峰一样，存在最大转矩点和对应的临界转差率。(2)从特性上看，电机的最大转矩很大，但是对应S=I的起动点，转矩却很小；(3)曲线的正负半区基本对称六、心得体会做负载性质的的实验时，最先遇到的问题就是如何画出向量图，我们查找了许多资料,例如使用COmPaSS函数，但CO叩a

7、ss函数是从原点开始画向量图，且有背景，又如quiver函数，也不能满足要求，最后决定自己在线段终点处画出两个微小的与该线段成一比较小的角度的线段作为箭头，这也是我本次一个很大的收获，学会灵活的运用知识，而不是一味想使用程序内置的函数，做实验，搞科研的时候要灵活。在绘制不同负载的曲线时，我也进一步的理解了曲线的绘制步骤，对单相变压器有了更深的理解。在绘制T-S时，我们同样也遇到了问题，当程序第一次运行时，会出现错误的提示，就是向量的长度不匹配，但是第二次运行却能够正常运行。调试了多次后还是没有结果，最后在老师的帮助下，我们终于找到了原因，所使用的矩阵应该先声明一个长度与其等式另一侧的矩阵维数相

8、匹配，才能够正常运行。在绘制T-S时，我也更加理解了转矩与转差率的关系，也了解了电机的起动转矩和最大转矩，以及稳定运行的区间。通过这次实验，我更加意识到了MAT1AB的重要性，虽然上学期我们学习了数学实验这门课，并且初步学习了Mat1ab的使用方法，但是数学实验上所学的毕竟只是基本的操作和一些简单的程序，没有经过大的检验和锻炼。当这次作业再一次重新使用MAT1AB时，起初感觉十分的不熟练，之后抽出一部分时间练习Mat1abo正所谓磨刀不误砍柴工，通过训练，我的Mat1ab越来越熟练，后面的工作也越来越快。未来我们的工作离不开计算机，未来熟练的掌握计算机基础知识对于我们来说是十分必要的。这次实验

9、，不仅提高了我的电机学知识水平，也帮助我更加熟练地掌握的MAT1AB的使用方法，对以后的工作学习有着巨大的帮助。最重要的是，我们的成功离不开老师的耐心指导，我们的几个难点都是在老师的耐心帮助下才得以解决的。可以说，老师的帮助与指导与我们的成功密不可分。七、附录%单相变压器负载运行时不同负载阻抗下的向量图c1ear;c1f;c1c;f=50;N1=876;N2=260;u2=6000;i2=180;k=N1N2;%调整阻、抗的比例r1n=5.5;x1n=12.4;rm=850;xm=8600;r2n=k2*0.45x2n=k2*0.964;%变压器二次侧负载阻抗角，取-90-30O3090the

10、ta=-30;%为了使图形更加直观，放大原副端阻抗三角形beta=8;r1=beta*r1n;r2=beta*r2n;x1=beta*x1n;x2=beta*x2n;Z1=r1+j*x1;Z2=r2+j*x2;u22=k*u2;i22=i2/k;%U2的大小U2=u22;%U2与x轴的夹角theta_U2=atan2(imag(U2),rea1(U2);%I2落后与U2I2=i22*(cos(theta*pi180)+j*sin(theta*pi180);theta_12=atan2(imag(12),rea1(12);,求出E2E2=U2+I2*Z2;theta_E2=atan2(imag(

11、E2),rea1(E2);E1=E2?theta_E1=atan2(imag(E1),rea1(E1);%为了在图形更加直观，放大kim=1.5;告求出主磁通fai=kim*E1/(-j*4.44*f*N1);theta_fai=atan2(imag(fai),rea1(fai);theta_fai1=theta_fai-pi/2;需求出Im=beta*(-E1)/(rm+j*xm);theta_Im=atan2(imag(Im),rea1(Im);括求出I1I1=1m-I2;theta_I1=atan2(imag(I1),rea1(I1);行求出U1U1=-E1+I1*r1+j*I1*x1;

12、theta_U1=atan2(imag(U1),rea1(U1);ur1=I1*r1;theta_ur1=atan2(imag(ur1),rea1(ur1);ux1=j*I1*x1;theta_ux1=atan2(imag(ux1),rea1(ux1);ur2=I2*r2;theta_ur2=atan2(imag(ur2),rea1(ur2);ux2=j*I2*x2;theta_ux2=atan2(imag(u2),rea1(ux2);rot_fai=abs(fai)*i;fai_amp=rot_fai*1e5;fai_ampr=rea1(fai_amp);fai_ampi=imag(fai_

13、amp);rot_U2=abs(U2)*(cos(theta_U2-theta_fai1)+j*sin(theta_U2-theta_fai1);U2r=rea1(rot_U2);U2i=imag(rot_U2);rot_E1=abs(E1)*(cos(theta_E1-theta_fai1)+j*sin(theta_E1-theta_fai1);rE1r=rea1(-rot_E1);rE1i=imag(-rot_E1);rot_U1=abs(U1)*(cos(theta_U1-theta_fai1)+j*sin(theta_U1-theta_fai1);U1r=rea1(rot_U1);U1

14、i=imag(rot_U1);rot_E2=abs(E2)*(cos(theta_E2-theta_fai1)+j*sin(theta_E2-theta_fai1);E2r=rea1(rot_E2);E2i=imag(rot_E2);rot_ur1=abs(ur1)*(cos(theta_ur1-theta_fai1)+j*sin(theta_ur1-theta_fai1);ur1r=rea1(rot_ur1);ur1i=imag(rot_ur1);rot_ux1=abs(ux1)*(cos(theta_ux1-theta_fai1)+j*sin(theta_ux1-theta_fai1);u

15、x1r=rea1(rot_ux1);ux1i=imag(rot_ux1);rot_ur2=abs(ur2)*(cos(theta_ur2-theta_fai1)+j*sin(theta_ur2-thetafai1);ur2r=rea1(rot_ur2);ur2i=imag(rot_ur2);rot_ux2=abs(ux2)*(cos(theta_ux2-theta_fai1)+j*sin(theta_ux2-theta_fai1);ux2r=rea1(rot_ux2);ux2i=imag(rot_ux2);rot_I1=abs(I1)*(cos(theta_I1-theta_fai1)+j*sin(theta_I1-theta_fai1);I1r=rea1(rot_I1);I1i=imag(rot_I1);rot_I2=abs(12)*(cos(theta_I2-theta_fai1)+j*sin(theta_I2