c6exsf1.m

基于MATLAB_Simulink的系统仿真技术与应用大家可以学习下

function [sys,x0,str,ts] = c6exsf1(t,x,u,flag,A,B,C,D)
switch flag,
case 0 % 初始化设置
   [sys,x0,str,ts]=mdlInitializeSizes(A,D); 
case 1 % 连续状态变量计算
   sys = mdlDerivatives(t,x,u,A,B); 
case 3 % 输出量计算
   sys = mdlOutputs(t,x,u,C,D); 
case { 2, 4, 9 } % 未定义标志
   sys = [];
otherwise % 处理错误
   error(['Unhandled flag = ',num2str(flag)]); 
end
%==============================================================
% mdlInitializeSizes 进行初始化，设置系统变量的大小
%==============================================================
function [sys,x0,str,ts] = mdlInitializeSizes(A,D)
sizes = simsizes;  % 取系统默认设置
sizes.NumContStates = size(A,1); % 设置连续变量个数
sizes.NumDiscStates = 0; % 设置离散状态个数，因为无离散状态，故设其为0
sizes.NumOutputs = size(A,1)+size(D,1); % 设置输出变量个数，为D的行数加系统的阶次
sizes.NumInputs = size(D,2); % 设置输入变量的个数，为D的列数
sizes.DirFeedthrough = 1; % 输出量的计算取决于输入量D
sizes.NumSampleTimes = 1; % 采样周期的个数，
sys = simsizes(sizes); % 设置系统的大小参数
x0 = zeros(size(A,1),1); % 设置为零初始状态
str = []; % 设置字符串矩阵
ts = [-1 0]; % 采样周期设置，前面的-1表示继承输入信号的采样周期
%==============================================================
% mdlDerivatives 计算系统的状态变量
%==============================================================
function sys = mdlDerivatives(t,x,u,A,B)
sys = A*x + B*u;
%==============================================================
% mdlOutputs 计算系统输出
%==============================================================
function sys = mdlOutputs(t,x,u,C,D)
sys = [C*x+D*u; x]; % 系统的增广输出