（1.2.1）无阻尼微分方程的m-函数程序.m

电机仿真,含阻尼微分方程的M函数程序和仿真图

%编写同步发电机无阻尼空在建立电压过程微分方程的M-函数
%将该M-函数定义为sy_ge_nodamp_noload_ode
function dydt=sy_ge_nodamp_noload_ode(t,y)
%下面输入电机基本参数
 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;
%下面输入电感系数矩阵
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);