新建 文本文档.txt

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

%编写同步发电机有阻尼绕组空在建立电压过程微分方程的M—函数
%将该M—函数定义为sy_ge_damp_noload_ode
function dydt=sy_ge_damp_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;Mafd0=Lmd,MaKd0=Lmd,MakqO=Lmq,Lfd=L1fd+Lmd,Lkd=L1kd+Lmd,Lkq=L1kq+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；]
%下面输入电阻矩阵
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=[0,0,Ufd,0,0]';
%下面列写微分方程
dydt=L\(Udq0-w*G*y-R*y);
%编写求解空ZAI 建立电压微分方程的函数（sy_ge_damp_noload_sol)
function sy_ge_damp_noload_sol
% 设定时间和初始值
tspan=[0 10];
y0=[0;0;0;0;0;];
% 方程求解
[t,y]=ode113('tongbufadianji',tspan,y0);
% 空载建立电压过程中励磁电流的变化规律
ifd=y(:,3);
plot(t,ifd)
%下面指定纵横轴标签
xlabel('Time[s]')
ylabel('Ifd[A]')