scalar.m

无刷双馈电机建模仿真的matlab、simulink代码。控制绕组的频率一定时

R=[0.81,0,0,0,0,0;0,0.81,0,0,0,0;0,0,0.81,0,0,0;0,0,0,0.81,0,0;0,0,0,0,0.00157,0;0,0,0,0,0,0.00157];
G=[0,0.24,0,0,0,0.00267;-0.24,0,0,0,-0.00267,0;0,0,0,0.63,0,0.0043;0,0,-0.63,0,0.0043,0;0,0,0,0,0,0;0,0,0,0,0,0];
L1=[177.62,0,-101.31,0,-14842,0;0,177.62,0,101.31,0,-14842;-101.31,0,63.74,0,9106.1,0;0,101.31,0,63.74,0,-9106.1;-14842,0,9106.1,0,1.3341e+006,0;0,-14842,0,-9106.1,0,1.3341e+006];
X=[0;0;0;0;0;0];
X1=[0;0;0;0;0;0];
U=[380;0;0;0;0;0];
H=0.0001;%步长
T=10;
WR=47.1;
J=10;
TL=0;
ala=0;
sun=0;
sun1=0;
TE=0;
 k=0;
 wc=20*2*3.14;
 a0=0.5;
 %fc=-10;
 %R radius=4
for i=0:H:10
      v=8
    %if (i>10), v=9,end 
    k=k+1;
    ala=WR;
  lam=ala*4/v;
    cpp=(0.0002694*lam^5-0.02439*lam^4+0.4049*lam^3+8.1192*lam^2-52.7939*lam^1+164.2884)/10000
    TL=-cpp*0.5*1.25*3.14*4*4*v^3/WR;
    wc=wc+(-a0*wc-29*a0*v+a0*314)*H
    WR=H*(TE-TL)/J+WR;
    sun=sun+(ala+WR)/2*H;
    sun1=sun1+wc*H;
   U(1,1)=380*cos(314*i-3*sun);
    %U(1,1)=380*cos(314*i-3*sun);
    U(2,1)=-380*sin(314*i-3*sun);
    U(3,1)=20*abs(wc/2/3.14)*cos(sun+sun1+pi/3);
    U(4,1)=20*abs(wc/2/3.14)*sin(sun+sun1+pi/3);
   % U(3,1)=7*cos(sun+2*pi*fc*i+pi/3);
   % U(4,1)=7*sin(sun+2*pi*fc*i+pi/3); this is the direct case
   % U(3,1)=0*cos(sun+2*pi*fc*i+pi/3);% this is the short circuit case
   % U(4,1)=0*sin(sun+2*pi*fc*i+pi/3);  above case function
  % if (fc>0 and fc<5 ) 
   %    U(3,1)=(7+20*abs(fc))*cos(sun+2*pi*fc*i+pi/3);
   % U(4,1)=(7+20*abs(fc))*sin(sun+2*pi*fc*i+pi/3);
   %end 
% {           if (abs(wc)>=5*2*3.14)
 %          U(3,1)=20*abs(wc/2/3.14)*cos(sun+wc*i+pi/3);
  %        U(4,1)=20*abs(wc/2/3.14)*sin(sun+wc*i+pi/3);
%  elseif (wc<0)
 %    U(3,1)=(7+14*abs(wc/2/3.14))*cos(sun+wc*i+pi/3);
  %   U(4,1)=(7+14*abs(wc/2/3.14))*sin(sun+wc*i+pi/3);
 %else
  %   U(3,1)=(7+20*abs(wc/2/3.14))*cos(sun+wc*i+pi/3);
   %  U(4,1)=(7+20*abs(wc/2/3.14))*sin(sun+wc*i+pi/3);
%end
%{ if (fc<0 and fc>-5) 
 %     
  %end }

    X1=L1*(U-WR*G*X-R*X);
    X=X+H*X1;
    TE=3*0.00089*(X(1,1)*X(6,1)-X(2,1)*X(5,1))+0.0043*(X(3,1)*X(6,1)+X(4,1)*X(5,1));
    
    N=WR*30/3.14;
    b(k)=N;
    c(k)=cpp;
    a(k)=i;
end
%i=0:0.0001:3
plot(a,b,'r')