sy_ge_damp_3duanlu_ode.m

MATLAB电机仿真精华50例,源代码！！！！！！！！

% 编写同步发电机有阻尼绕组空载建立电压过程微分方程的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;Ll=3.0892E-01;Lmd=3.2164;Lmq=9.7153E-01;Llfd=3.0712E-01;Llkd=4.9076E-01;Llkq=1.0365;
Ld=Lmd+Ll;Lq=Lmq+Ll;Mafd0=Lmd;Makd0=Lmd;Makq0=Lmq;Lfd=Llfd+Lmd;Lkd=Llkd+Lmd;Lkq=Llkq+Lmq;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,       0,	      -3/2*Makq0; 
     Ld,     0,        Mafd0,    Makd0,	   0;
     0,      0,        0,        0,        0;
     0,      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,    Rkd, 0;
    0,   0,   0,    0,   Rkq];
% 下面输入电压向量：
Udq0=[-100,0,24,0,0]';
% 下面列写微分方程：
dydt=L\(Udq0-w*G*y-R*y);