tongbufadianji.m

MATLAB电机仿真精华50例,源代码！！！！！！！！

% 同步发电机有阻尼绕组转子绕组短路,定子绕组突然加上对称电压
% 将该M—函数定义为sy_ge_damp_rduanlu_ode
function dydt=sy_ge_damp_rduanlu_ode(t,y)
% 下面输入电机基本数据：
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,	Lmd,      0; 
     0,            Lq,          0,	    0,          Lmq;
     Lmd,    0,           Lfd,      Lmd,       0;
     Lmd,    0,           Lmd,     Lkd,        0;
     0,            Lmq,   0,        0,          Lkq]                        
 G=[  0,     -Lq,    0,       0,	      -3/2*Makq0; 
     Ld,     0,     Mafd0,   Makd0,	   0;
     0,      0,     0,       0,        0;
     0,      0,     0,       0,        0;
     0,      0,     0,       0,        0]
% 下面输入电阻矩阵：
R=[ r,   0,   0,    0,   0;
    0,   r,   0,    0,   0;
    0,   0,   Rfd,  0,   0;
    0,   0,   0,    Rkd, 0;
    0,   0,   0,    0,   Rkq]

% 下面输入电压向量：
Udq0=[100,0,0,0,0]';
% 下面列写微分方程：
dydt=L\(Udq0-w*G*y-R*y);





% 下面列写标准形式微分方程并将感应电势用励磁电流的拟合多项式函数表示：
% 这里的a(1)、a(2)、a(3)、a(4)是指上一步的多项式拟合系数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);