s函数问题

zzz.dan · 发表于 2020-9-3 10:35

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

这是dynamic_.m文件的内容
function [sys,x0,str,ts,simStateCompliance] = dynamic_(t,x,u,flag)
switch flag,
  case 0,
[sys,x0,str,ts,simStateCompliance]=mdlInitializeSizes;
  case 3,
sys=mdlOutputs(t,x,u);
  case {1, 2, 4, 9 }
sys=[];
  otherwise
DAStudio.error('Simulink:blocks:unhandledFlag', num2str(flag));
end
function [sys,x0,str,ts,simStateCompliance]=mdlInitializeSizes
sizes = simsizes;
sizes.NumContStates  = 0;
sizes.NumDiscStates  = 0;
sizes.NumOutputs    = 3;
sizes.NumInputs    = 3;
sizes.DiRFeedthrough = 1;
sizes.NumSampleTimes = 1; % at least one sample time is needed
sys = simsizes(sizes);
x0  = [];
str = [];
ts  = [0 0];
simStateCompliance = 'UnknownSimState';
% end mdlInitializeSizes
function sys=mdlOutputs(t,x,u)
z1 = u(1);
z2 = u(2);
ut = u(3);
thet = 0.95;
a1 = 0.5;
a2 = 1.9;
d1_z1 = z2-1/8*(z1^4*thet^3);
d1_z2 = ut+1/10*(z2*exp(z1*thet))-(a1*sin(z2)*z1*log*(1+a2^2*z1^2));
y1 = thet*z1;
sys = [d1_z1;d1_z2;y1];

问题如下：