同步电机的三相突然短路的 MATLAB 计算 姓名：林启星 专业:电力电子与电力传动 学号：1001202591、新建一个 M 文件存放主程序如下：％＊＊*******＊*＊****＊＊＊**＊**＊*＊*＊＊*＊＊＊＊%定义变量及变量赋值%＊*＊**＊***＊***＊*＊*＊＊*＊***＊＊＊**＊*＊***global A B U0R=0.005；Rf=0。000656；RD=0。00151;RQ=0。00159；xd=1.0;xq=0。60；xf=1。03;xD=0。95;xQ=0。70;xad=0.85;xaq=0。45；u0=1;i0=1；w0=0。5548；A=［—xd，0,xad,xad，0;0，-xq,0,0，xaq；-xad，0，xf，xad,0;—xad,0,xad，xD,0;0,—xaq，0,0,xQ］；B=[-R ， xq,0,0 ， — xaq;—xd ， — R,xad ， xad ， 0;0 ， 0,Rf ， 0 ， 0 ； 0 ， 0,0,RD,0 ；0,0，0，0,RQ］；％＊**＊*＊*****＊＊＊*＊**＊*****＊＊*＊＊＊*＊**＊％求出功率角 a0 及 ud0,uq0，id0,iq0,if0%＊＊*＊＊*＊*＊＊＊＊****＊＊*＊＊＊**＊**＊*＊＊＊**＊a0=atan(（u0*sin(w0）+xq*i0)/（u0＊cos(w0)+R＊i0））-w0;ud0=u0*sin(a0）；uq0=u0*cos（a0）;id0=i0＊sin(a0+w0);iq0=cos(a0+w0);Eom=uq0+R＊iq0+xd*id0;If0=Eom/xad;U0=［-ud0;-uq0;0；0;0］；%***＊**＊**＊*＊*＊＊＊＊**＊**＊***＊*＊＊＊**＊＊％通过使用 ode45 来计算各绕组电流％*＊*＊**＊＊＊***＊****＊**＊**＊＊**＊***＊***I0=[0;0；0;0;0］;[t，i］=ode45（'IdqfDQ'，[0，2000]，I0）；id=i(:，1)+id0;iq=i（:,2)+iq0;If=i(:,3)+If0；iD=i(：，4)；iQ=i（：，5);ia=id.＊cos（t+pi）-iq。＊sin(t+pi）;%＊＊＊*＊＊＊*****＊＊***＊＊＊**＊＊＊**＊＊＊**＊**%三相短路时的电磁转矩 Te%＊*****＊**＊**＊＊*＊*＊*＊＊*＊****＊＊**＊＊**Wd=—xd*id+xad＊If+xad*iD;Wq=—xq*iq+xaq＊iQ；Te=Wd.*iq-Wq.*id；％＊＊＊*****＊***＊＊*＊**＊**＊*＊***＊＊*＊*＊**%计算 ia 及 Te％＊**＊***＊＊*****＊＊＊＊＊*＊＊********＊*＊*＊Ten=Te(end，1）ian=ia(end,1)％＊*＊＊*＊＊＊**＊＊*＊＊＊**＊＊***＊＊＊**＊******％各变量图形输出%**＊＊＊＊＊＊＊*＊*＊*＊*＊**＊＊＊*＊**＊＊*＊＊＊＊*＊figure(1）;plot(t，ia）title（'ia=f（t）'）；xlabel(’t’);ylabel(’ia');figure（2）；plot(t,If)title（’If=f(t)’);xlabel('t')；ylabel('If'）;figure（3）;plot(t，iD)title（'iD=f（t）’)；xlabel（’t’）；ylabel（’iD’）；figure(4）；plot（t，iQ）title（'iQ=f（t）’）;xlabel（’t’）；ylabel('iQ')；figure(5）；plot（t,Te)title（'Te=f（t）'）;xlabel（’t'）;ylabel('Te’)；2、新建一个M文件（IdqfDQ。m）存放子程序程序如下：function dI=IdqfDQ(t，i)global A B U0dI=A\U0-A\B＊i;3、 运行结果如下：Ten =: 0。0160ian =0.6706
