ch2example23sfun.m

通信程序源代码和模型通信程序源代码和模型通信程序源代码和模型

function [sys,x0,str,ts] =ch2example23Sfun(t,x,u,flag,v0, theta0,L,m,g,k) 
%  单摆的实现
switch flag, 
    case 0                 % flag=0 初始化 
       sizes = simsizes;           % 获取SIMULINK仿真变量结构
       sizes.NumContStates  = 2;   % 连续系统的状态数是2
       sizes.NumDiscStates  = 0;   % 离散系统的状态数是0
       sizes.NumOutputs     = 1;   % 输出信号数目是1
       sizes.NumInputs      = 0;   % 输入信号数目是0
       sizes.DirFeedthrough = 0;   % 该系统不是直通的 
       sizes.NumSampleTimes = 1;   % 这里必须为1 
       sys = simsizes(sizes); 
       str = [];           % 通常为空矩阵 
       x0  = [v0, theta0]; % 初始状态矩阵x0
       ts  = [0 0];        % 连续系统
   case 1                  % flag=1 计算连续系统状态方程
       vdot=g*sin(x(2))-k/m*x(1);
       thetadot=-1/L*x(1);
       sys=[vdot; thetadot]; %状态方程
   case 3                    % flag=3 计算输出 
       sys=x(2) ;    % 输出方程 
   case {2,4,9}      % 其他作不处理的flag 
       sys=[];       % 无用的flag时返回sys为空矩阵
otherwise % 异常处理 
       error(['Unhandled flag = ',num2str(flag)]); 
end