Simulink壁咚(二)——What and How

最新推荐文章于 2024-08-31 12:14:30 发布

注释远方

最新推荐文章于 2024-08-31 12:14:30 发布

阅读量2.1k

点赞数

文章标签： MBD

本文链接：https://blog.csdn.net/qinze5857/article/details/129439428

版权

一、Simulink 知多少

1.1 记录每次运行结果

1.2 子系统封装

右击子系统，点击Mask，Create Mask 创建并编辑

数字信号处理X——EasySimulink封装和库_simulink函数封装库_注释远方的博客-CSDN博客

找实例， demo----->输入关键字

1.3 代码生成

打开设置齿轮

二、S-Function

在命令行中输入demo并搜索s function examples，打开level1中的csfunc获得模板。

function [sys,x0,str,ts] = sfunc(t,x,u,flag,para1)
%CSFUNC An example MATLAB file S-function for defining a continuous system.  
%   Example MATLAB file S-function implementing continuous equations: 
%      x' = Ax + Bu
%      y  = Cx + Du
%   
%   See sfuntmpl.m for a general S-function template.
%
%   See also SFUNTMPL.
  
%   Copyright 1990-2009 The MathWorks, Inc.
% t表示仿真时间，x当前的状态值，u当前系统的输入，flag是阶段标志位


A=[-0.09   -0.01
    1          0];

B=[ 1   -7
    0   -2];

C=[ 0    2
    1   -5];

D=[-3    0
    1    0];

switch flag,

  %%%%%%%%%%%%%%%%%%
  % Initialization %
  %%%%%%%%%%%%%%%%%%
  case 0,
    [sys,x0,str,ts]=mdlInitializeSizes(A,B,C,D);

  %%%%%%%%%%%%%%%
  % Derivatives %
  %%%%%%%%%%%%%%%
  case 1,
    sys=mdlDerivatives(t,x,u,A,B,C,D);

  %%%%%%%%%%%
  % Outputs %
  %%%%%%%%%%%
  case 3,
    sys=mdlOutputs(t,x,u,A,B,C,D,para1);

  %%%%%%%%%%%%%%%%%%%
  % Unhandled flags %
  %%%%%%%%%%%%%%%%%%%
  case { 2, 4, 9 },
    sys = [];

  %%%%%%%%%%%%%%%%%%%%
  % Unexpected flags %
  %%%%%%%%%%%%%%%%%%%%
  otherwise
    DAStudio.error('Simulink:blocks:unhandledFlag', num2str(flag));

end
% end csfunc

%
%=============================================================================
% mdlInitializeSizes
% Return the sizes, initial conditions, and sample times for the S-function.
%=============================================================================
%
function [sys,x0,str,ts]=mdlInitializeSizes(A,B,C,D)

sizes = simsizes;
sizes.NumContStates  = 2;%连续系统状态的数量,A
sizes.NumDiscStates  = 0;%离散系统
sizes.NumOutputs     = 2;% C D
sizes.NumInputs      = 2;% B
sizes.DirFeedthrough = 1;%直馈0/1，如果输入直接影响输出为1（y=Cx+Du），否则为0(y=Cx)
sizes.NumSampleTimes = 1;%采样时间的个数

sys = simsizes(sizes);
x0  = zeros(2,1);%初始状态的定义
str = [];
ts  = [0 0];%[采样间隔 采样截止时间]，连续系统为[0 0],离散系统[0.1 0]

% end mdlInitializeSizes
%
%=============================================================================
% mdlDerivatives
% Return the derivatives for the continuous states.
%=============================================================================
%
function sys=mdlDerivatives(t,x,u,A,B,C,D)

sys = A*x + B*u;%实际上为x' = Ax + Bu

% end mdlDerivatives
%
%=============================================================================
% mdlOutputs
% Return the block outputs.
%=============================================================================
%
function sys=mdlOutputs(t,x,u,A,B,C,D,para1)

sys = C*x + D*u; % 实际为 y= Cx +Du
disp(para1);
% end mdlOutputs