svpwm的两电平s函数

function [sys,x0,str,ts] = VPWM2LNPC(t,x,u,flag) switch flag, case 0, [sys,x0,str,ts]=mdlInitializeSizes(ThetaStart); case {1,4,9}, sys=[]; case 2, sys=mdlUpdate(t,x,u,r,tfn,ThetaDelta); case 3, sys=mdlOutputs(t,x,u); otherwise error(['Unhandled flag = ',num2str(flag)]); end function [sys,x0,str,ts]=mdlInitializeSizes sizes = simsizes; sizes.NumContStates = 0; sizes.NumDiscStates = 0; sizes.NumOutputs = 6; sizes.NumInputs = 4; sizes.DirFeedthrough = 1; sizes.NumSampleTimes = 1; % at least one sample time is needed sys = simsizes(sizes); x0 = []; str = []; ts = [-2 0]; simStateCompliance = 'UnknownSimState'; function sys=mdlDerivatives(t,x,u) sys = []; function sys=mdlUpdate(t,x,u) sys = []; function sys=mdlOutputs(t,x,u) Ua=u(1); Ub=u(2); Udc=u(3); Ts=u(4); A=x(1); B=x(2); C=x(3); X=x(1); Y=x(2); Z=x(3); sector=0; CMPRA1=0; CMPRA2=0; CMPRB1=0; CMPRB2=0; CMPRC1=0; CMPRC2=0; Va=Ub; Vb=(-1)*Ub + sqrt(3)*Ua; Vc=(-1)*Ub - sqrt(3)*Ua; if (Va>=0) A =1; else A=0; end if (Vb>=0) B=1; else B=0; end if (Vc>=0) C = 1; else C=0; end sector = A + 2*B + 4*C; X=1.732051*Ub/Udc; Y=(0.8660*Ub+1.5*Ua)/Udc; Z=(-0.8660*Ub+1.5*Ua)/Udc; switch sector, case 1, t1 =-Z; t2 =X; Ta =((0.5)*(1-t1-t2)*Ts); Tb =(Ta+t1*Ts); Tc =(Tb+t2*Ts); CMPRA1 =Tb; CMPRB1 =Ta; CMPRC1 =Tc; CMPRA2 =Tb; CMPRB2 =Ta; CMPRC2 =Tc; case 2, t1 =Z; t2 =Y; Ta =((0.5)*(1-t1-t2)*Ts); Tb =(Ta+t1*Ts); Tc =(Tb+t2*Ts); CMPRA1 =Ta; CMPRB1 =Tc; CMPRC1 =Tb; CMPRA2 =Ta; CMPRB2 =Tc; CMPRC2 =Tb; case 3, t1 =X; t2 =-Y; Ta =((0.5)*(1-t1-t2)*Ts); Tb =(Ta+t1*Ts); Tc =(Tb+t2*Ts); CMPRA1 =Ta; CMPRB1 =Tb; CMPRC1 =Tc; CMPRA2 =Ta; CMPRB2 =Tb; CMPRC2 =Tc; case 4, t1 =-X; t2 =Z; Ta =((0.5)*(1-t1-t2)*Ts); Tb =(Ta+t1*Ts); Tc =(Tb+t2*Ts); CMPRA1 =Tc; CMPRB1 =Tb; CMPRC1 =Ta; CMPRA2 =Tc; CMPRB2 =Tb; CMPRC2 =Ta; case 5, t1 =-Y; t2 =-Z; Ta =((0.5)*(1-t1-t2)*Ts); Tb =(Ta+t1*Ts); Tc =(Tb+t2*Ts); CMPRA1 =Tc; CMPRB1 =Ta; CMPRC1 =Tb; CMPRA2 =Tc; CMPRB2 =Ta; CMPRC2 =Tb; case 6, t1 =Y; t2 =-X; Ta =((0.5)*(1-t1-t2)*Ts); Tb =(Ta+t1*Ts); Tc =(Tb+t2*Ts); CMPRA1 =Tb; CMPRB1 =Tc; CMPRC1 =Ta; CMPRA2 =Tb; CMPRB2 =Tc; CMPRC2 =Ta; end sys =[CMPRA1,CMPRA2,CMPRB1,CMPRB2,CMPRC1,CMPRC2,num]; function sys=mdlGetTimeOfNextVarHit(t,x,u) sampleTime = 78.125e-6; % Example, set the next hit to be one second later. sys = t + sampleTime; function sys=mdlTerminate(t,x,u) sys = [];