tongbufadianji.m

来自「MATLAB电机仿真精华50例,源代码！！！！！！！！」· M 代码 · 共 27 行

% 编写同步发电机无阻尼绕组三相突然短路微分方程的M—函数
% 将该M—函数定义为sy_ge_nodamp_3duanlu_ode
function dydt=sy_ge_nodamp_3duanlu_ode(t,y)
% 下面输入电机基本数据：
r=2.9069,Rfd=5.9013E-01,Rkd=11.900,Rkq=20.081;Ufd=24;w=377,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,	; 
     0,            Lq,          0,	    ;
     Lmd,          0,           Lfd   ];                        
G=[  0,     -Lq,    0,       ; 
     Ld,     0,     Mafd0,   ;
     0,      0,     0,       ];   
% 下面输入电阻矩阵：
R=[ r,   0,   0,    ;
    0,   r,   0,    ;
    0,   0,   Rfd,  ];
% 下面输入电压向量：
Udq0=[-100,-10,24]';
% 下面列写微分方程：
dydt=L\(Udq0-w*G*y-R*y);