tongbufadianji.m

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

% 编写建立两相空载短路微分方程的M—函数
% 将该M—函数定义为tongbufadianji
function dydt=tongbufadianji(t,y)
% 下面输入电机基本数据： 
r=.26;Rfd=0.13;Rkd=0.0224;Rkq=0.02;Ufd=17.8876;w=377;Ll=1.14e-3;Lmd=13.7e-3;Lmq=11.0e-3;
Llfd=2.1e-3;Llkd=1.4e-3;Llkq=1e-3;
Lfd=Llfd+Lmd;Lkd=Llkd+Lmd;Lkq=Llkq+Lmq;Ld=Lmd+Ll;Lq=Lmq+Ll;
%--------------------------------------------------------------------------------------------------------------------------------------   
%电感矩阵 
 L=[  
        Lfd,          Lmd,         0,            -Lmd*sin(w*t);
        Lmd,          Lkd,         0,            -Lmd*sin(w*t);
        0,            0,           Lkq,          -Lmq*cos(w*t);
        Lmd*sin(w*t), Lmd*sin(w*t),Lmq*cos(w*t), -(Ld+Lq)/2+(Ld-Lq)/2*cos(2*w*t)] ;
G=[ 
         0,       0,        0,  -Lmd*cos(w*t);
         0,       0,        0,  -Lmd*cos(w*t);
         0,       0,        0,  Lmq*sin(w*t);
     Lmd*cos(w*t), Lmd*cos(w*t),-Lmq*sin(w*t), -(Ld-Lq)*sin(2*w*t)];
 %--------------------------------------------------------------------------------------------------------------------------------------   
% 电阻矩阵
 R=[ 
       Rfd,  0,   0,  0+0+0;
       0,    Rkd, 0,  0+0+0;
       0,    0,   Rkq,0+0+0;
      0,    0,   0,   -r];
 %--------------------------------------------------------------------------------------------------------------------------------------   
 %电压向量
Ualphabeta=[17.8876,0,0,0]';
 %微分方程
dydt=L\(Ualphabeta-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);