% 編寫同步發電機有阻尼繞組空載建立電壓過程微分方程的M—函數
% 將該M—函數定義為sy_ge_damp_noload_ode
function dydt=sy_ge_damp_noload(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,	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=[0,0,Ufd,0,0]';
% 下面列寫微分方程：
dydt=L\(Udq0-w*G*y-R*y);