matlab(Simulink)中S-function函数编写规则.pdfVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
matlab(Simulink)中S-function 函数编写规则 s 函数是system Function 的简称,用它来写自己的simulink 模块。(够简单吧,^_^ ,详细的概念介绍大伙看帮助吧)可以用matlab 、C 、C++ 、Fortran 、Ada 等语 言来写,这儿我只介绍怎样用matlab 语言来写吧(主要是它比较简单) 先讲讲为什么要用s 函数,我觉得用s 函数可以利用matlab 的丰富资源,而不仅仅局限于simulink 提供的模块,而用c 或c++等语言写的s 函数还可以实现对 硬件端口的操作,还可以操作windows API 等的 先介绍一下simulink 的仿真过程(以便理解s 函数),simulink 的仿真有两个阶段:一个为初始化,这个阶段主要是设置一些参数,像系统的输入输出个数、 状态初值、采样时间等;第二个阶段就是运行阶段,这个阶段里要进行计算输出、更新离散状态、计算连续状态等等,这个阶段需要反复运行,直至结束。 在matlab 的workspace 里打edit sfuntmpl(这是matlab 自己提供的s 函数模板),我们看它来具体分析s 函数的结构。它的第一行是这样的:function [sys,x0,str,ts]=sfuntmpl(t,x,u,flag) 先讲输入与输出变量的含义:t 是采样时间,x 是状态变量,u 是输入(是做成simulink 模块的输入),flag 是仿真过程中的状态标志(以它来判断当前是初始化还 是运行等);sys 输出根据flag 的不同而不同(下面将结合flag 来讲sys 的含义),x0 是状态变量的初始值,str 是保留参数(mathworks 公司还没想好该怎么用它, 嘻嘻,一般在初始化中将它置空就可以了,str=[]),ts 是一个 1×2 的向量,ts(1)是采样周期,ts(2)是偏移量。 下面结合sfuntmpl.m 中的代码来讲具体的结构: switch flag, %判断flag,看当前处于哪个状态 case 0, [sys,x0,str,ts]=mdlInitializeSizes; flag=0 表示处于初始化状态,此时用函数mdlInitializeSizes 进行初始化,此函数在sfuntmpl.m 的149 行 我们找到他,在初始化状态下,sys 是一个结构体,用它来设置模块的一些参数,各个参数详细说明如下 size = simsizes;%用于设置模块参数的结构体用simsizes 来生成 sizes.NumContStates = 0;%模块连续状态变量的个数 sizes.NumDiscStates = 0;%模块离散状态变量的个数 sizes.NumOutputs = 0;%模块输出变量的个数 sizes.NumInputs = 0;%模块输入变量的个数 sizes.DirFeedthrough = 1;%模块是否存在直接贯通(直接贯通我的理解是输入能%直接控制输出) sizes.NumSampleTimes = 1;%模块的采样时间个数,至少是一个 sys = simsizes(sizes); %设置完后赋给sys 输出 举个例子,考虑如下模型: dx/dt=fc(t,x,u) 也可以用连续状态方程描述:dx/dt=A*x+B*u x(k+1)=fd(t,x,u) 也可以用离散状态方程描述:x(k+1)=H*x(k)+G*u(k) y=fo(t,x,u) 也可以用输出状态方程描述:y=C*x+D*u 设上述模型连续状态变量、离散状态变量、输入变量、输出变量均为 1 个,我们就只需改上面那一段代码为: (一般连续状态与离散状态不会一块用,我这儿是为了方便说明) sizes.NumContStates=1;sizes.NumDiscStates=1;sizes.NumOutputs=1;sizes.NumInpu ts=1; 其他的可以不变。继续在mdlInitializeSizes 函数中往下看: x0 = []; %状态变量设置为空,表示没有状态变量,以我们上面的假设,可改 %为x0=[0,0](离散和连续的状态变量我们都设它初值为0) str = []; %这个就不用说了,保留参数嘛,置[]就可以了,反正没什么用,可%能7.0 会给它一些意义 ts = [0 0]; %采样周期设为0 表示是连续系统,如果是离散系

文档评论(0)

zsmfjy + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档