%編寫同步發(fā)電機(jī)無阻尼空在建立電壓過程微分方程的M-函數(shù)
%將該M-函數(shù)定義為sy_ge_nodamp_noload_ode
function dydt=sy_ge_nodamp_noload_ode(t,y)
%下面輸入電機(jī)基本參數(shù)
 r=2.9069,Rfd=5.9013E-01,Rkd=11.900,Rkq=20.081;Ufd=24;w=377,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,MAKQ0=LMD,MAKQO=LMQ,LFD=L1FD+LMD,LKD=L1KD+LMD,LKQ=L1KQ+LMP,MFKD=LMD;
%下面輸入電感系數(shù)矩陣
l=[ld,0,lmd;
   0,lq, 0;
   lmd, 0, lfd]
G=[0, -LQ,0;
   LD,0, MAFD0;
   O, 0, RFD]
%下面輸入電壓向量
Uabc=[0,0,Ufd]';
% 下面列寫微分方程
dydt=L\(Uabe-w*g*y-r*y);