confreq.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=0;
J=0.3;
%J=10;
TL=0;
ala=0;
sun=0;
TE=0;
 k=0;
 v=8;
 c=0;
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %fc=-10;
for i=0:H:3
    fc=5
    if (i>10), fc=5, end 
    if (i>20), fc=-5, end
    if (i>30), fc=-10, end
    k=k+1;
    ala=WR;
    lam=ala*4/v/4.4857;
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    %if (abs(lam)>50)
     %   cpp=0 
    %else 
     %   cpp=(0.0002694*lam^5-0.02439*lam^4+0.4049*lam^3+8.1192*lam^2-52.7939*lam^1+164.2884)/10000
    %end 
     %    c(k)=cpp;
     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
           lam1=1/(1/(lam-0.02*c)-0.003/(1+c^3));
        m= 0.73*(151/lam1-0.58*c-0.002*c^2.14-13.2)*exp(-18.4/lam1) ;
        if (m<0)
            cp(k)=0;
        else
            cp(k)=m;
        end
            cpp=cp(k);
    TL=cpp*0.5*1.25*3.14*4*4*v^3/WR;
    T(k)=TL;
    TL=0;
    WR=H*(TE+TL)/J+WR;
    sun=sun+(ala+WR)/2*H;
    %SH=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)=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(fc)>=5 )
   U(3,1)=20*abs(fc)*cos(sun+2*pi*fc*i+pi/3);
    U(4,1)=20*abs(fc)*sin(sun+2*pi*fc*i+pi/3);
elseif (fc<0)
     U(3,1)=(7+14*abs(fc))*cos(sun+2*pi*fc*i+pi/3);
     U(4,1)=(7+14*abs(fc))*sin(sun+2*pi*fc*i+pi/3);
 else
     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 (fc<0 and fc>-5) 
 

    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;
    a(k)=i;
end
%i=0:0.0001:3
plot(a,b,'r')