matlab电机仿真同步发电机

% 编写建立自激过程微分方程的M—函数
% 将该M—函数定义为sh_ge_se_ex_ode（shunt_generator_self_excited_ode）
function dydt=tongbufadianji（ty）
% 首先指定全局变量
%global   a1 a2 a3 a4 Rf  
% 下面输入电机基本数据：
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]