06.08.2013 Views

内容简介作者简介 - 科学与工程计算系

内容简介作者简介 - 科学与工程计算系

内容简介作者简介 - 科学与工程计算系

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

%<br />

%<br />

switch flag,<br />

case 0,<br />

[sys,x0,str,ts]=mdlInitializeSizes(angzzy);<br />

case 1,<br />

sys=mdlDerivatives(t,x,u,dampzzy,gravzzy);<br />

case 2,<br />

sys=mdlUpdate(t,x,u);<br />

case 3,<br />

sys=mdlOutputs(t,x,u);<br />

case 9,<br />

sys=mdlTerminate(t,x,u);<br />

otherwise<br />

error(['Unhandled flag = ',num2str(flag)]);<br />

end<br />

% ------ mdlInitializeSizes ------------------------<br />

function [sys,x0,str,ts]=mdlInitializeSizes(angzzy)<br />

sizes = simsizes;<br />

sizes.NumContStates = 2;<br />

sizes.NumDiscStates = 0;<br />

sizes.NumOutputs = 1;<br />

sizes.NumInputs = 1;<br />

sizes.DirFeedthrough = 0;<br />

sizes.NumSampleTimes = 1;<br />

sys = simsizes(sizes);<br />

x0 = angzzy;<br />

str = [ ];<br />

ts = [0, 0];<br />

% -----mdlDerivatives -----------------------------<br />

function sys=mdlDerivatives(t,x,u,dampzzy,gravzzy)<br />

dx(1)=-dampzzy*x(1)-gravzzy*sin(x(2))+u;<br />

dx(2)=x(1);<br />

sys =dx ;<br />

% -----mdlUpdate ----------------------------------<br />

function sys=mdlUpdate(t,x,u)<br />

sys = [ ];<br />

% -----mdlOutputs ---------------------------------<br />

function sys=mdlOutputs(t,x,u)<br />

sys = x(2);<br />

% -----mdlTerminate -------------------------------<br />

function sys=mdlTerminate(t,x,u)<br />

sys = [ ];<br />

(4)<br />

34

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!