ty.m

来自「这是多电平svpwm空间电压矢量控制的算法」· M 代码 · 共 107 行

function [sys,x0,str,ts]=s_function(t,x,u,flag)
switch flag,
   case 0,
    [sys,x0,str,ts]=mdlInitializeSizes;
   case 3,
    sys=mdlOutputs(t,x,u);
   case {2,4,9}
    sys=[];
otherwise
    error(['Unhandled flag= ',num2str(flag)]);
end
 
function [sys,x0,str,ts]=mdlInitializeSizes
sizes=simsizes;
sizes.NumContStates=0;
sizes.NumDiscStates=0;
sizes.NumOutputs=1;
sizes.NumInputs=10;
sizes.DirFeedthrough=1;
sizes.NumSampleTimes=0;

sys=simsizes(sizes);
x0=[];
str=[];
ts=[];

function sys=mdlOutputs(t,x,u)
if u(6)<=0
    sys=0;
end

if  u(6)==1
    if (u(3)+u(4))<=u(1)
        sys=u(10)*u(9);
    elseif (u(3)-u(4))>=u(1)
        sys=u(10)*u(9)-u(10);
    elseif u(5)>u(2)
        sys=u(10)-u(10)*u(9);
    else
        sys=u(10)*u(8);
    end
end
    
if  u(6)==2
    if (u(3)+u(4))<=u(1)
        sys=u(10)*u(8);
    elseif (u(3)-u(4))>=u(1)
        sys=u(10)*u(8)-u(10);
    elseif u(5)>u(2)
        sys=u(10)-u(10)*u(8);
    else
        sys=u(10)*u(8);
    end
end

if  u(6)==3
    if (u(3)+u(4))<=u(1)
        sys=u(10)*u(7);
    elseif (u(3)-u(4))>=u(1)
        sys=u(10)*u(7)-u(10);
    elseif u(5)>u(2)
        sys=u(10)-u(10)*u(7);
    else
        sys=u(10)*u(7);
    end
end

if  u(6)==4
    if (u(3)+u(4))<=u(1)
        sys=2*u(10)+u(10)*u(9);
    elseif (u(3)-u(4))>=u(1)
        sys=-u(10)-u(10)*u(9);
    elseif u(5)>u(2)
        sys=u(10)+u(10)*u(9);
    else
        sys=-u(9)*u(10);
    end
end

if  u(6)==5
    if (u(3)+u(4))<=u(1)
        sys=-u(10)*u(8);
    elseif (u(3)-u(4))>=u(1)
        sys=-u(10)*u(8)-u(10);
    elseif u(5)>u(2)
        sys=u(10)*u(8)+u(10);
    else
        sys=-u(10)*u(8);
    end
end

if  u(6)==6
    if (u(3)+u(4))<=u(1)
        sys=-u(10)*u(7);
    elseif (u(3)-u(4))>=u(1)
        sys=-u(10)-u(10)*u(7);
    elseif u(5)>u(2)
        sys=-u(10)+u(10)*u(9);
    else
        sys=-u(10)*u(7);
    end
end

if u(6)>6
    sys=0;
end