%?編寫建立自激過程微分方程的M—函數(shù)
%?將該M—函數(shù)定義為sh_ge_se_ex_ode（shunt_generator_self_excited_ode）
function?dydt=tongbufadianji（ty）
%?首先指定全局變量
%global???a1?a2?a3?a4?Rf??
%?下面輸入電機基本數(shù)據(jù)：
r=2.9069E-03Rfd=5.9013E-04Rkd=1.1900E-02Rkq=2.0081E-02;Ufd=24;w=314Ll=3.0892E-04Lmd=3.2164E-03Lmq=9.7153E-04Llfd=3.0712E-04Llkd=4.9076E-04Llkq=1.0365E-03
Laa0=1/3*（Lmd+Lmq）+Ll;Laa2=1/3*（Lmd-Lmq）;Mab0=1/2*Laa0;;Mab2=Laa2;Mafd0=2/3*LmdMakd0=2/3*LmdMakq0=LmqLfd=Llfd+LmdLkd=Llkd+LmdLkq=Llkq+LmqMfkd=Lmd;
L=[?-（Laa0+Laa2*cos（2*w*t））????????????-（-Mab0+Mab2*cos（2*（w*t+2*pi/3）））??-（-Mab0+Mab2*cos（2*（w*t-2*pi/3）））??Mafd0*cos（w*t）	????????Makd0*cos（w*t）??????????-Makq0*sin（w*t）;?
????-（-Mab0+Mab2*cos（2*（w*t+2*pi/3）））?-（Laa0+Laa2*cos（2*（w*t-2*pi/3）））???-（-Mab0+Mab2*cos（2*（w*t）））?????????Mafd0*cos（w*t-2*pi/3）	Makd0*cos（w*t-2*pi/3）??-Makq0*sin（w*t-2*pi/3）;
????-（-Mab0+Mab2*cos（2*（w*t-2*pi/3）））?-（-Mab0+Mab2*cos（2*（w*t）））?????????-（Laa0+Laa2*cos（2*（w*t+2*pi/3）））???Mafd0*cos（w*t+2*pi/3）	Makd0*cos（w*t+2*pi/3）??-Makq0*sin（w*t+2*pi/3）;
????-Mafd0*cos（w*t）???????????????????-Mafd0*cos（w*t-2*pi/3）?????????????-Mafd0*cos（w*t+2*pi/3）?????????????Lfd?????????????????????Mfkd?????????????????????0;
????-Makd0*cos（w*t）???????????????????-Makd0*cos（w*t-2*pi/3）?????????????-Makd0*cos（w*t+2*pi/3）?????????????Mfkd????????????????????Lkd??????????????????????0;
?????Makq0*sin（w*t）????????????????????Makq0*sin（w*t-2*pi/3）??????????????Makq0*sin（w*t+2*pi/3）?????????????0???????????????????????0????????????????????????Lkq]???