(2.1.1)三相突然短路仿真有阻尼微分方程m-函数程序.m

电机仿真,含阻尼微分方程的M函数程序和仿真图

%编写同步发电机有阻尼绕组空载建立电压过程微分方程的M-函数
%将该M-函数定义为sy_ge_damp_3duanlu_ode
function dydt=sy_ge_damp_3duanlu_ode(t,y)
%下面输入电机基本参数
 r=2.9069,Rfd=5.9013E-01,Rkd=11.900,Rkq=20.081;Ufd=24;w=314,L1=3.0892e-01,Lmd=3.2164,Lmq=9.7153E-01,L1fd=3.0712E-01,
L1kd=4.9076E-01,L1kq=1.0365
Ld=Lmd+L1;Lq=Lmq+L1;MafdO=LMD,Makd0=Lmd,MakqO=Lmq,Lfd=L1fd+Lmd,Lkd=L1kd+Lmd,Lkq=L1kq+Lmp,Mfkd=Lmd;
%下面输入电感系数矩阵
l=[Ld,0,Lmd,lmd,0;
   0,Lq, 0,0,Lmq;
   Lmd, 0, Lfd,Lmd,0;
   Lmd,0,Lmd,Lkd,0;
   0,Lmq,0,0,Lkq];
G=[0, -Lq,0,q,-3/2*Makq0;
   Ld,0, Mafd0,Makd0,0;
   O, 0, 0,0,0;
   0,0,0,0,0];
%下面输入电阻矩阵
R=[r，0，0，0，0；
   0，r,0,0,0;
   0,0,Rfd,0,0;
    0,0,0,Rfd,0;
    0,0,0,0,Rkq];
%下面输入电压向量
Udq0=[-100,0,24,0,0]';
%下面列写微分方程
dydt=L\(Udq0-w*G*y-R*y);