% 編寫同步發(fā)電機無阻尼繞組空載建立電壓過程微分方程的M—函數(shù)
% 將該M—函數(shù)定義為sy_ge_nodamp_noload_ode
function dydt=tongbufadianji(t,y)
% 首先指定全局變量
%global   a1 a2 a3 a4 Rf  
% 下面輸入電機基本數(shù)據(jù)：
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,       ]    
% 下面進(jìn)行曲線擬合：

R=[ r,   0,   0,    ;
    0,   r,   0,    ;
    0,   0,   Rfd,  ]

Uabc=[0,0,Ufd]';

dydt=L\(Uabc-w*G*y-R*y);





% 下面列寫標(biāo)準(zhǔn)形式微分方程并將感應(yīng)電勢用勵磁電流的擬合多項式函數(shù)表示：
% 這里的a(1)、a(2)、a(3)、a(4)是指上一步的多項式擬合系數(shù)a(i)
%dydt=[1/Lf*(a(1)*iff^3+a(2)*iff^2+a(3)*iff^1+a(4))-Rf/Lf*iff];
% 給全局變量賦值
%a1=a(1);a2=a(2);a3=a(3);a4=a(4);