下面是我的S函数:(帮我看下哪里出了问题)
function [sys,x0,str,ts]=Batch_sim(t,x,u,flag,xx0,ss0,pp0,cc0,um,Km,Ki,Kd,Yx,Yp,m,a,B,K,KLa,c1,Y0,m0)
%xx0,ss0,pp0,cc0分别为x,s,p,c的初始值
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));
end
function [sys,x0,str,ts,simStateCompliance]=mdlInitializeSizes
sizes = simsizes;
sizes.NumContStates = 4;
sizes.NumDiscStates = 0;
sizes.NumOutputs = 4;
sizes.NumInputs = 0;
sizes.DirFeedthrough = 0;
sizes.NumSampleTimes = 1; % at least one sample time is needed
sys = simsizes(sizes);
x0=[xx0 ss0 pp0 cc0];
str = [];
ts = [0 0];
function sys=mdlDerivatives(t,x,u,xx0,ss0,pp0,cc0,um,Km,Ki,Kd,Yx,Yp,m,a,B,K,KLa,c1,Y0,m0)
miu=um*x(2)/(Km+x(2)+x(2)^2/Ki);
rx=miu*x(1);
rd=-Kd*x(1);
rp=a*rx+B*x(1);
rs=(rx/Yx)+(rp/Yp)+m*x(1);
ro2=KLa*(c1-x(4))-(rx/Y0+m0*x(1));
sys = [rx+rd;-rs;rp-K*x(3);ro2];
function sys=mdlOutputs(t,x,u)
sys = [];
![]()
simulik中的模型
![]()
这是显示的错误信息 |