扰动观察法.docx

function sys,x0,str,ts,simStateCompliance = PO_MPPT_Boost(t,x,u,flag)switch flag, case 0, sys,x0,str,ts,simStateCompliance=mdlInitializeSizes; case 1, sys=mdlDerivatives(t,x,u); case 2, sys=mdlUpdate(t,x,u); case 3, sys=mdlOutputs(t,x,u); case 4, sys=mdlGetTimeOfNextVarHit(t,x,u); case 9, sys=mdlTerminate(t,x,u); otherwise DAStudio.error(Simulink:blocks:unhandledFlag, num2str(flag);endfunction sys,x0,str,ts,simStateCompliance=mdlInitializeSizessizes = simsizes;sizes.NumContStates = 0;sizes.NumDiscStates = 3;sizes.NumOutputs = 1;sizes.NumInputs = 2;sizes.DirFeedthrough = 1;sizes.NumSampleTimes = 1; % at least one sample time is neededsys = simsizes(sizes);x0 = 0,0,0.5;%U I Urefstr = ;ts = 1e-6;simStateCompliance = UnknownSimState;% end mdlInitializeSizesfunction sys=mdlDerivatives(t,x,u)sys = ;% end mdlDerivativesfunction sys=mdlUpdate(t,x,u)DU=0.001; %步长dU=u(1)-x(1);dP=u(1)*u(2)-x(1)*x(2);if dP0 if dU0 Uref=x(3)+DU; else Uref=x(3)-DU; endelseif dP0 Uref=x(3)-DU; else Uref=x(3)+DU; endelse Uref=x(3);endif Uref=1 Uref=1;endif Uref=0 Uref=0;endx(1)=u(1);x(2)=u(2);x(3)=Uref;sys = x(1),x(2),x(3);% end mdlUpdatefunction sys=mdlOutputs(t,x,u)sys = x(3);% end mdlOutputsfunction sys=mdlGetTimeOfNextVarHit(t,x,u)sampleTime = 1; % Example, set the next hit to be one second later.sys = t +