s函数是system Function的简称，用它来写自己的simulink模块。（够简单吧，^_^，详细的概念介绍大伙看帮助吧）可以用matlab、C、C++、Fortran、Ada等语言来写，这儿我只介绍怎样用matlab语言来写吧（主要是它比较简单）先讲讲为什么要用s函数，我觉得用s函数可以利用matlab的丰富资源，而不仅仅局限于simulink提供的模块，而用c 或c++等语言写的s函数还可以实现对硬件端口的操作，还可以操作windows API等的先介绍一下simulink的仿真过程（以便理解s函数），simulink的仿真有两个阶段：
1.初始化：这个阶段主要是设置一些参数，像系统的输入输出个数、状态初值、采样时间等；
2.运行阶段：这个阶段里要进行计算输出、更新离散状态、计算连续状态等等。这个阶段需要反复运行，直至结束。
在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.NumInp
u ts=1; 其他的可以不变。
继续在mdlInitializeSizes函数中往下看：
x0 = []; %状态变量设置为空，表示没有状态变量，以我们上面的假设，可改%为x0=[0,0](离散和连续的状态变量我们都设它初值为0) str = []; %这个就不用说了，保留参数嘛，置[]就可以了，反正没什么用，可能7.0会给它一些意义
ts = [0 0]; %采样周期设为0表示是连续系统，如果是离散系统在下面的mdlGet %TimeOfNextVarHit函数中具体介绍嘻嘻，总算讲完了初始化，后面的应该快了。
在sfuntmpl的106行继续往下看：
case 1, sys=mdlDerivatives(t,x,u); flag=1表示此时要计算连续状态的微分，即上面提到的dx/dt=fc(t,x,u)中的dx/dt,找到mdlDerivatives函数（在193行）如果设置连续状态变量个数为0，此处只需sys=[]; 就可以了（如sfuntmpl中一样），按我们上述讨论的那个模型，此处改成sys=fc(t,x(1),u)或sys=A*x(1)+B*u %我们这儿x(1)是连续状态变量，而x(2)是离散的，这儿只用到连续的，此时的输出sys就是微分继续，
在sfuntmpl的112行：
case 2, sys=mdlUpdate(t,x,u); flag=2表示此时要计算下一个离散状态，即上面提到的x(k+1)=fd(t,x,u)，找到mdlUpdate函数（在206行）它这儿sys=[];表示没有离散状态，我们这而可以改成sys=fd(t,x(2),u)或sys=H*x(2)+G*u;%sys即为x(k+1) 看来后面几个一两句话就可了，呵呵，
在sfuntmpl的118行：
case 3, sys=mdlOutputs(t,x,u); flag=3表示此时要计算输出，即y=fo(t,x,u),找到mdlOutputs函数（在218行），如上，如果sys=[]表示没有输出，我们改成sys=fo(t,x,u)或sys=C*x+D*u %sys此时为输出y 好像快完了，嘻嘻，
在sfuntmpl的124行
case 4, sys=mdlGetTimeOfNextVarHit(t,x,u); flag=4表示此时要计算下一次采样的时间，只在离散采样系统中有用(即上文的mdlInit ializeSizes中提到的ts设置ts(1)不为0) 连续系统中只需在mdlGetTimeOfNextVarHit函数中写上sys=[];这个函数主要用于变步长的设置，具体实现大家可以用edit vsfunc看vsfunc.m这个例子最后一个，
在sfuntmpl的130行
case 9, sys=mdlTerminate(t,x,u); flag=9表示此时系统要结束，一般来说写上在mdlT erminate函数中写上sys=[]就可，如果你在结束时还要设置什么，就在此函数中写关于sfuntmpl这个s函数的模板讲完了。s函数还可以带用户参数，下面给个例子，和simulink 下的gain模块功能一样，大伙自己看吧，我睡觉去了，累了
function [sys,x0,str,ts] = sfungain(t,x,u,flag,gain)
switch flag,
case 0,
sizes = simsizes;
sizes.NumContStates = 0;
sizes.NumDiscStates = 0;
sizes.NumOutputs = 1;
sizes.NumInputs = 1;
sizes.DirFeedthrough = 1;
sizes.NumSampleTimes = 1;
sys = simsizes(sizes);
x0=[]; str=[]; ts=[0,0];
case 3, sys=gain*u;
case {1,2,4,9}, sys = [];
end
SIMULINK s-function的设计Simulink为用户提供了许多内置的基本库模块，通过这些模块进行连接而构成系统的模型。对于那些经常使用的模块进行组合并封装可以构建出重复使用的新模块，但它依然是基于Simulink原来提供的内置模块。而Simulink s-function 是一种强大的对模块库进行扩展的新工具。
(一)、s-function的概念
s-function是一个动态系统的计算机语言描述，在MATLAB里，用户可以选择用m文件编写，也可以用c或mex文件编写，在这里只给大家介绍如何用m文件编写s-function。
S-function提供了扩展Simulink模块库的有力工具，它采用一种特定的调用语法，使函数和Simulink解法器进行交互。S-function最广泛的用途是定制用户自己的Simulink模块。它的形式十分通用，能够支持连续系统、离散系统和混合系统。
(二)、建立m文件s-function
1、使用模板文件：sfuntmp1. m 格式：[sys,x0]=function(t,x,u,flag) 该模板文件位于MATLAB根目录下toolbox/simulink/blocks目录下。模板文件里s-function的结构十分简单，它只为不同的flag的值指定要相应调用的m文件子函数。
比如当flag=3时，即模块处于计算输出这个仿真阶段时，相应调用的子函数为sys=mdloutputs(t,x,u)。模板文件使用switch语句来完成这种指定，当然这种结构并不唯一，用户也可以使用if语句来完成同样的功能。
而且在实际运用时，可以根据实际需要来去掉某些值，因为并不是每个模块都需要经过所有的子函数调用。模板文件只是Simulink为方便用户而提供的一种参考格式，并不是编写s-function的语法要求，用户完全可以改变子函数的名称，或者直接把代码写在主函数里，但使用模板文件的好处是，比较方便，而且条理清晰。
使用模板编写s-function，用户只需把s-函数名换成期望的函数名称，如果需要额外的输入参量，还需在输入参数列表的后面增加这些参数，因为前面的4个参数是simulink调用s-function时自动传入的。对于输出参数，最好不做修改。接下去的工作就是根据所编s-function要完成的任务，用相应的代码去替代模板里各个子函数的代码即可。Simulink 在每个仿真阶段都会对s-function进行调用，在调用时，Simulink会根据所处的仿真阶段为flag传入不同的值，而且还会为sys这个返回参数指定不同的角色，也就是说尽管是相同的sys变量，但在不同的仿真阶段其意义却不相同，这种变化由simulink自动完成。
m文件s-function可用的子函数说明如下：
mdlInitializeSizes(flag=0)：定义s-function模块的基本特性，包括采样时间、连续或者离散状态的初始条件和sizes数组。
mdlDerivatives(flag=1)：计算连续状态变量的微分方程。
mdlUpdate(flag=2)：更新离散状态、采样时间和主时间步的要求。
mdlOutputs(flag=3)：计算s-function的输出。
mdlGetTimeOfNextVarHit(flag=4)：计算下一个采样点的绝对时间，这个方法仅仅是在用户在mdlInitializeSizes 里说明了一个可变的离散采样时间。
概括说来，建立s-function可以分成两个分离的任务：初始化模块特性包括输入输出信号的宽度，离散连续状态的初始条件和采样时间。将算法放到合适的s-function子函数中去。
2、定义s-function的初始信息为了让Simulink识别出一个m文件s-function，用户必须在s-函数里提供有关s-函数的说明信息，包括采样时间、连续或者离散状态个数等初始条件。这一部分主要是在mdlInitializeSizes子函数里完成。Sizes数组是s-function函数
信息的载体，它内部的字段意义为：NumContStates(sys(1))：连续状态的个数（状态向量连续部分的宽度）NumDiscStates(sys(2))：离散状态的个数（状态向量离散部分的宽度）NumOutputs(sys(3))：输出变量的个数（输出向量的宽度）NumInputs(sys(4))：输入变量的个数（输入向量的宽度）DirFeedthrough(sys(5))：有不连续根的数量NumSampleTimes(sys(6))：采样时间的个数，有无代数循环标志如果字段代表的向量宽度为动态可变，则可以将它们赋值为－1。注意DirFeedthrough是一个布尔变量，它的取值只有0和1两种，0表示没有直接馈入，此时用户在编写mdlOutputs子函数时就要确保子函数的代码里不出现输入变量u；1表示有直接馈入。NumSampleTimes表示采样时间的个数，也就是ts变量的行数，与用户对ts的定义有关。需要指出的是，由于s-function会忽略端口，所以当有多个输入变量或多个输出变量时，必须用mux模块或demux模块将多个单一输入合成一个复合输入向量或将一个复合输出向量分解为多个单一输出。3、输入和输出参量说明S-function默认的4个输入参数为t、x、u和flag，它们的次序不能变动，代表的意义分别为：t：代表当前的仿真时间，这个输入参数通常用于决定下一个采样时刻，或者在多采样速率系统中，用来区分不同的采样时刻点，并据此进行不同的处理。x：表示状态向量，这个参数是必须的，甚至在系统中不存在状态时也是如此。它具有很灵活的运用。u：表示输入向量。flag：是一个控制在每一个仿真阶段调用哪一个子函数的参数，由Simulink在调用时自动取值。S-function默认的4个返回参数为sys、x0、它们的次序不能变动，代表的意义分别为：sys：是一个通用的返回参数，它所返回值的意义取决于flag 的值。x0：是初始的状态值（没有状态时是一个空矩阵[]），这个返回参数只在flag值为0时才有效，其他时候都会被忽略。一、有一系统如下：dx1=x2 dx2=9.81*sin(x(1))-2*x(2)+u 求出系统在单位阶跃输入下的x1的状态变化曲线，假设x1,x2初值为0。function [sys,x0]=dong(t,x,u,flag) if flag==0 sys=[2;0;2;1;0;0]; x0=[0;0]; elseif flag==1 sys=[x(2);9.81*sin(x(1))-2*x(2)+u]; elseif flag==3 sys=[x(1);x(2)]; else sys=[]; end
s函数一下是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个，我们就只需改上面那一段代码为：
MATLAB 学习指南 第六章．编写与调用函数 在这一章中，我们讨论如何用多源代码文件来构造一个程序。首先，解释代码文件在MATLAB中如何工作。在编译语言中，例如FORTRAN，C ，或C++，代码被存储在一个或多个源文件中，在进行编译的时候，这些源文件组合在一起 形成了一个单独的可执行文件。作为一种解释型语言，MATLAB以一种更广泛的方式来处理多个源文件。MATLAB代码被放入带有扩展名.m的ASCII文件（或称m-文件）中。MATLAB 6 有一个集成字处理与调试应用程序，尽管会用到其它编辑程序如vi或emacs，集成字处理与调试应用程序仍是编译m-文件的首选程序。 有两种不同的m-文件。一种是脚本文件，它是一种最简单的文件，仅仅将MATLAB中的指令收集在一起。当在交互提示符处输入文件名执行脚本文件时，MATLAB在m-文件内读取并执行指令，就好像指令是我们输入的。而且，似乎我们能够削减m-文件的内容并将削减过的内容传到MATLAB指令窗口中。这种m-文件的用法将在6.1节中给予概述。 在6.2节中要讨论的第二种m-文件包含一个单一函数，此函数名与此m-文件名相同。这种m-文件包含一段独立的代码，这段代码具有一个明确规定的输入/输出界面；那就是说，传给这段代码一列空变量arg1，arg2，…，这段独立代码就能够被调用，然后返回输出值out1，out2，…。一个函数m-文件的第一个非注释行包含函数标头，其形式如下： 此m-文件以返回指令结束，将执行程序返回到函数被调用的位置。或者在交互指令提示符处或者在另一个m-文件内，无论何时用下列指令调用函数代码，函数代码都将被执行。 输入映射到空变量：arg1=var1，arg2=var2，等等。在函数主体内，输出值被分配给了变量out1，out2，等等。当遇到返回值时，当前值out1，out2，…在函数被调用处被映射到变量outvar1，outvar2，…。在用可变长度自变量和输出变量列表编写函数时，MATLAB允许更多的自由。例如，也可以使用下列指令来调用函数。 在此情况下，仅返回一个单一输出变量，这个变量在出口处包含函数变量out1的值。输入和输出自变量可能是字符串，数值，向量，矩阵，或者更高级的数据结构。 为什么使用函数呢？因为从每门计算机科学课程中可知，把一个大的程序分割 成多个可以单独执行一个被明确规定的和被注释过的任务的小程序会使大程序 易读，易于修改，不易于出错。在MATLAB中，先为程序编写一个主文件，或者是一个脚本文件或者更好的话，是一个能够返回一个单一整数的函数m-文件（返回1表示程序执行成功，0表示不完全程序执行，负值表示出现运行误差），这个主文件是程序的进入点。通过把m-文件当作函数来调用，此程序文件可以
S函数编辑流程及相关 1. S函数模板编辑环境进入： 在MA TLAB主界面中直接输入：edit sfuntmpl 即可弹出S函数模板编辑的M文件环境，修改即可。 在MA TLAB主界面中直接输入：sfundemos，即可调出S 函数的许多编程例子。 2. S函数模板的相关基础： 1)M文件S函数的引导语句为： t x f p u sys Function [p x str ts flag , ,1 2 ,...) , , ,0 , ] ,( , S函数默认的四个输入参数：t ,x ,u ,flag S函数默认的四个输出函数：sys ,x0 ,str ,ts 各个参数的含义如下： T ：代表当前的仿真时间，该输入决定了下一个采样时间； X ：表示状态向量，行向量，引用格式：X(1),X(2) U ：表示输入向量; Flag ：控制在每一个仿真阶段调用哪一个子函数的参数，由SIMULINK在调用时自动取值； Sys ：通用的返回变量，返回的数值决定Flag值，mdlUpdates 里：列向量，引用格式：Sys(1,1),Sys(2,1)；mdlOutputs里：行向量，引用格式：Sys =x. X0 ：初始的状态值；列向量，引用格式：X0=[ 0；0；0 ]
Str ：空矩阵，无具体含义； Ts ：包含模块采样时间和偏差的矩阵。[period, offset] 当Ts为-1时，表示与输入信号同采样周期。 2)S函数工作方式： Flag = 0时，调用mdlInitializeSizes函数，定义S函数的基本特性，包括采样时间，连续或者离散状态的初始条件和Sizes 数组； Flag = 1时，调用mdlDerivatives函数，计算连续状态变量的微分方程；求所给表达式的等号左边状态变量的积分值的过程。 Flag = 2时，调用mdlUpdate函数，用于更新离散状态，采样时间和主时间步的要求； Flag = 3时，调用mdlOutputs函数，计算S函数的输出； Flag = 4时，调用mdlGetTimeOfNextV arHit函数，计算下一个采样点的绝对时间，这个方法仅仅是使用户在mdlInitializeSize 里说明一个可变的离散采样时间； Flag = 9时，调用mdlTerminate函数，实现仿真任务的结束。 3)S函数仿真过程： ①初始化：mdlInitializeSizes，初始化S函数 ●初始化SimStruct，包含了S函数的所有信息； ●设置输入、输出端口数； ●设置采样时间； ●分配存储空间。
利用matlab编写S函数求解微分方程
自动化专业综合设计报告 设计题目：利用matlab编写S函数求解微分方程 所在实验室：自动化系统仿真实验室 指导教师：郭卫平 学生姓名律迪迪 班级文自0921 学号200990519114 成绩评定：
一、设计目的 了解使用simulink的扩展工具——S-函数，s函数可以利用matlab的丰富资源，而不仅仅局限于simulink提供的模块，而用c或c++等语言写的s函数还可以实现对硬件端口的操作，还可以操作windows API等的，它的魅力在于完美结合了simulink框图简洁明快的特点和编程灵活方便的优点，提供了增强和扩展sinulink能力的强大机制，同时也是使用RTW实现实时仿真的关键。 二、设计要求 求解解微分方程 y’=y-2x/y
y(0)=1 要求利用matlab编写S函数求解 三、设计内容（可加附页） 【步骤1】获取状态空间表达式。 在matlab中输入 dsolve(‘Dy=y-2*x/y’,’y(0)=1’，’x’) 得到 y=(2*x+1).^(1/2); 【步骤2】建立s函数的m文件。 利用21·用S函数模板文件。 以下是修改之后的模板文件sfuntmpl.m的内容。function [sys,x0,str,ts] = sfuntmpl(t,x,u,flag) %SFUNTMPL General M-file S-function template % With M-file S-functions, you can define you own ordinary differential % equations (ODEs), discrete system equations, and/or just about % any type of algorithm to be used within a Simulink block diagram. %
matlab中s函数编写心得(转) Part I: 所谓s函数是system Function的简称, 用它来写自己的simulink模块. s函数可以用matlab、C、C++、Fortran、Ada等语言来写，这儿我只介绍怎样用matlab语言来写吧（主要是它比较简单）< xmlnamespace prefix ="o" ns ="urn:schemas- microsoft-com:office:office" /> 先讲讲为什么要用s函数，我觉得用s函数可以利用matlab的丰富资源，而不仅仅局限于simulink提供的模块，而用c或c++等语言写的s函数还可以实现对硬件端口的操作，还可以操作windows API等 先介绍一下simulink的仿真过程（以便理解s函数），simulink 的仿真有两个阶段：一个为初始化，这个阶段主要是设置一些参数，像系统的输入输出个数、状态初值、采样时间等；第二个阶段就是运行阶段，这个阶段里要进行计算输出、更新离散状态、计算连续状态等等，这个阶段需要反复运行，直至结束. 在matlab的workspace里输入edit sfuntmpl(这是matlab自己提供的s函数模板)，我们看它来具体分析s函数的结构.
1. 函数的函数头 函数的第一行：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)是偏移量 2. 函数分析 下面结合sfuntmpl.m中的代码来讲具体的结构： switch flag, %判断flag，看当前处于哪个状态 case 0, [sys,x0,str,ts]=mdlInitializeSizes; // 解释说明 flag=0表示当前处于初始化状态，此时调用函数mdlInitializeSizes进行初始化，此函数在该文件的第149行定义. 其中的参数sys是一个结构体，它用来设置模块的一些参数，各个参 数详细说明如下 size = simsizes;%用于设置模块参数的结构体用simsizes来生
MATLAB 程序设计方法及若干程序实例 樊双喜 (河南大学数学与 信息科学学院开封475004) 摘要本文通过对 MATLAB 程序设计中的若干典型问题做简要的分析和总结,并在此基础上着重讨论了有关算法设计、程序的调试与测试、算法与程序的优化以及循环控制等方面的问题.还通过对一些程序实例做具体解析,来方便读者进行编程训练并掌握一些有关MATLAB 程序设计方面的基本概念、基本方法以及某些问题的处理技巧等.此外,在文章的最后还给出了几个常用数学方法的算法程序, 供读者参考使用.希望能对初学者进行 MATLAB 编程训练提供一些可供参考的材料,并起到一定的指导和激励作用,进而为MATLAB 编程入门打下好的基础. 关键字算法设计；程序调试与测试；程序优化；循环控制 1 算法与程序 1.1 算法与程序的关系算法被称为程序的灵魂,因此在介绍程序之前应先了 解什么是算法.所谓算 法就是对特定问题求解步骤的一种描述.对于一个较复杂的计算或是数据处理的问题,通常是先设计出在理论上可行的算法,即程序的操作步骤,然后再按照算法逐步翻译成相应的程序语言,即计算机可识别的语言. 所谓程序设计,就是使用在计算机上可执行的程序代码来有效的描述用于解决特定问题算法的过程.简单来说,程序就是指令的集合.结构化程序设计由于采用了模块分化与功能分解,自顶向下,即分而治之的方法,因而可将一个较复杂的问题分解为若干子问题,逐步求精.算法是操作的过程,而程序结构和程序流程则是算法的具体体现. 1.2MATLAB 语言的特点 MATLAB 语言简洁紧凑,使用方便灵活,库函数极其丰富,其语法规则与科技人员的思维和书写习惯相近,便于操作.MATLAB 程序书写形式自由,利用其丰富
S函数的简介及编写 摘自恒润科技S-function的编写 1. S函数模板编辑环境进入： 在MATLAB主界面中直接输入：edit sfuntmpl 即可弹出S函数模板编辑的M文件环境，修改即可。 在MATLAB主界面中直接输入：sfundemos，即可调出S 函数的许多编程例子。 2. S函数模板的相关基础： 1)M文件S函数的引导语句为： x t flag f u Function [p sys x str ts p , ,1 2 ,...) , , ,0 , ] ,( , S函数默认的四个输入参数：t ,x ,u ,flag S函数默认的四个输出函数：sys ,x0 ,str ,ts 各个参数的含义如下： T ：代表当前的仿真时间，该输入决定了下一个采样时间； X ：表示状态向量，行向量，引用格式：X(1),X(2) U ：表示输入向量; Flag ：控制在每一个仿真阶段调用哪一个子函数的参数，由SIMULINK在调用时自动取值； Sys ：通用的返回变量，返回的数值决定Flag值，mdlUpdates里：列向量，引用格式：Sys(1,1),Sys(2,1)；mdlOutputs里：行向量，引用格式：Sys =x. X0 ：初始的状态值；列向量，引用格式：X0=[ 0；0；0 ] Str ：空矩阵，无具体含义；
Ts ：包含模块采样时间和偏差的矩阵。[period, offset] 当Ts为-1时，表示与输入信号同采样周期。 2)S函数工作方式： Flag = 0时，调用mdlInitializeSizes函数，定义S函数的基本特性，包括采样时间，连续或者离散状态的初始条件和Sizes数组； Flag = 1时，调用mdlDerivatives函数，计算连续状态变量的微分方程；求所给表达式的等号左边状态变量的积分值的过程。 Flag = 2时，调用mdlUpdate函数，用于更新离散状态，采样时间和主时间步的要求； Flag = 3时，调用mdlOutputs函数，计算S函数的输出； Flag = 4时，调用mdlGetTimeOfNextVarHit函数，计算下一个采样点的绝对时间，这个方法仅仅是使用户在mdlInitializeSize 里说明一个可变的离散采样时间； Flag = 9时，调用mdlTerminate函数，实现仿真任务的结束。 3)S函数仿真过程： ①初始化：mdlInitializeSizes，初始化S函数 ●初始化SimStruct，包含了S函数的所有信息； ●设置输入、输出端口数； ●设置采样时间； ●分配存储空间。 ②数值积分：mdlDerivatives ●用于连续状态的求解和非采样过零点；
、最小二乘拟合原理 x= xl x2 ... xn y= yl y2 ... yn 求m 次拟合 ?力* y 卅…I ZA ； A T A = ZX 茁 X x i - X x i +1 ,- ? ? ? ［函Oi …备F =⑷矿丄? A T y 所以m 次拟合曲线为y = a 0 +勿?怎+吐■审+???? +如■牙皿 二、 Matlab 实现程序 function p=funLSM (x, y, m) %x z y 为序列长度相等的数据向量，m 为拟合多项式次数 format short; A=zeros(m+l,m+l); for i=0:m for j=0:m A(i + 1, j + 1)=sum(x.A (i+j)); end b(i+1)=sum(x.A i.*y); end a=A\b 1; p=fliplr (a'); 三、 作业 题1：给出如下数据，使用最小二乘法球一次和二次拟合多项式(取小数点后3位) X 1.36 1.49 1.73 1.81 1.95 2.16 2.28 2.48 Y 14.094 15.069 16.844 17.378 18.435 19.949 20.963 22.495 解:
? x=[1.36 1.49 1.73 1. 81 1. 95 2. 16 2. 28 2. 48]: ? y=[14.094 15.069 16.844 17. 378 18.435 19.949 20.963 22.495]; >> p=funLSM(x, y? 1) P = 7.4639 3.9161 >> p=funLSM(x, y? 2) P = 0.3004 6.3145 4.9763 一次拟合曲线为： y = 7.464x+ 3.91S 二次拟合曲线为： y = +6.315^4-4.976 一次拟合仿真图
1.三角波产生器 t=-3:0.01:3; f1=tripuls(t); subplot(3,1,1); plot(t,f1); axis([-3,3,-0.2,1.2]) set(gcf,'color','w'); f2=tripuls(t,4); subplot(3,1,2); plot(t,f2); axis([-3,3,-0.2,1.2]) %set(gcf,'color','w'); f3=tripuls(t,4,-1); subplot(3,1,3); plot(t,f3); axis([-3,3,-0.2,1.2]) 2.离散序列的相加与相乘 function[x,n]=jxl(x1,x2,n1,n2) n=min(min(n1),min(n2)):max(max(n1),max(n2)); s1=zeros(1,length(n));s2=s1; s1(find((n>=min(n1))&(n<=max(n1))==1))=x1; s2(find((n>=min(n2))&(n<=max(n2))==1))=x2; x=s1+s2;//x=s1.*s2:%序列乘 axis([(min(min(n1),min(n2))-1),(max(max(n1),max(n2))+1),(min(x)-0.5), (max(x)+0.5)]) 3.序列的反摺 function[x,n]=xlfz(x1,n1) x=fliplr(x1);n=fliplr(n1); stem(n,x,'filled') axis([min(n)-1,max(n)+1,min(x)-0.5,max(x)+0.5]) 4.序列的卷积 function[x,n]=gghconv(x1,x2,n1,n2) x=conv(x1,x2) ns=n1(1)+n2(1); leg=length(x1)+length(x2)-2; n=ns:(ns+leg) subplot(2,2,1) stem(n1,x1,'filled') title('x1(n)') xlabel('n') subplot(2,2,2)
MATLAB作业一1、试求出如下极限。 （1） 23 25 (2)(3) lim (5) x x x x x x x ++ + →∞ ++ + ，（2） 23 3 1 2 lim () x y x y xy x y →- → + + ，（3） 22 22 22 1cos() lim ()x y x y x y x y e+ → → -+ + 解：（1）syms x； f=((x+2)^(x+2))*((x+3)^(x+3))/((x+5)^(2*x+5)) limit（f，x，inf） =exp（-5） （2）syms x y； f=(x^2*y+x*y^3)/(x+y)^3； limit(limit(f,x,-1),y,2) =-6； （3）syms x y； f=(1-cos(x^2+y^2))/(x^2+y^2)*exp(x^2+y^2)； limit(limit(f,x,0),y,0) =0 2、试求出下面函数的导数。 （1 ）() y x=, (2)22 atan ln() y x y x =+ 解; （1）syms x; f=sqrt(x*sin(x)*sqrt(1-exp(x))); g= diff(f,x); g== (sin(x)*(1 - exp(x))^(1/2) + x*cos(x)*(1 - exp(x))^(1/2) - (x*exp(x)*sin(x))/(2*(1 - exp(x))^(1/2)))/(2*(x*sin(x)*(1 - exp(x))^(1/2))^(1/2)) pretty(g)= (2)syms x y; f=atan(y/x)-log(x^2+y^2) pretty(-simple(diff(f,x)/diff(f,y)))= 2 x + y =------- x - 2 y (3) 假设1 cos u- =，试验证 22 u u x y y x ?? = ???? 。 解：syms x y; u=1/cos(sqrt(x/y)); diff(diff(u,x),y)-diff(diff(u,y),x)=0; 所以： 22 u u x y y x ?? = ????
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 sf untmpl(这是matlab自己提供的s函数模板)，我们看它来具体分析s函数的结构。它的第一行是这样的：f unction [sys,x0,str,ts]=sf untmpl(t,x,u,f lag) 先讲输入与输出变量的含义：t是采样时间，x是状态变量，u是输入（是做成simulink模块的输入）,f lag是仿真过程中的状态标志（以它来判断当前是初始化还是运行等）；sys输出根据f lag的不同而不同（下面将结合f l ag来讲sys的含义），x0是状态变量的初始值，str是保留参数（mathworks公司还没想好该怎么用它，嘻嘻，一般在初始化中将它置空就可以了,str=[])，ts是一个1×2的向量，ts(1)是采样周期，ts(2)是偏移量。 下面结合sf untmpl.m中的代码来讲具体的结构： switch f lag, %判断f lag，看当前处于哪个状态 case 0, [sys,x0,str,ts]=mdlInitializeSizes; f l ag=0表示处于初始化状态，此时用函数mdlInitializeSizes进行初始化，此函数在sf untmpl.m的149行 我们找到他，在初始化状态下，sys是一个结构体，用它来设置模块的一些参数，各个参数详细说明如下 size = simsizes;%用于设置模块参数的结构体用simsizes来生成 sizes.NumContStates = 0;%模块连续状态变量的个数 sizes.NumDiscStates = 0;%模块离散状态变量的个数 sizes.NumOutputs = 0;%模块输出变量的个数 sizes.NumInputs = 0;%模块输入变量的个数 sizes.DirFeedthrough = 1;%模块是否存在直接贯通（直接贯通我的理解是输入能%直接控制输出） sizes.NumSample Times = 1;%模块的采样时间个数，至少是一个 sys = simsizes(sizes); %设置完后赋给sys输出 举个例子，考虑如下模型： dx/dt=f c(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表示是连续系统，如果是离散系统在下面的mdlGet %TimeOfNextV arHit函数中具体介绍 嘻嘻，总算讲完了初始化，后面的应该快了 在sf untmpl的106行继续往下看： case 1, sys=mdlDerivatives(t,x,u); f l ag=1表示此时要计算连续状态的微分，即上面提到的dx/dt=f c(t,x,u)中的dx/dt,找到 mdlDerivatives函数（在193行）如果设置连续状态变量个数为0，此处只需sys=[]; 就可以了（如sf untmpl中一样），按我们上述讨论的那个模型，此处改成 sys=f c(t,x(1),u)或sys=A*x(1)+B*u %我们这儿x(1)是连续状态变量，而x(2)是离散的，这儿只用到连续的，此时的输出sys就是微分 继续，在sf untmpl的112行： case 2, sys=mdlU pdate(t,x,u);
函数是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 是仿真过程中的状态标志（以它来判断当前是初始化还是运行等）；s ys输出根据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;%模块连续状态变量的个数
期末作业 1.求解微分方程???? ?????=++++=++=8102432 1.31. 232221.1x x x x x x x x x x x ，初始条件1230x x x ===。要求交mdl 文件及仿真结果的截图。 解：用Simulink 连接的仿真图如下所示： 求解出的x1，x2，x3的解如下图（a ）所示的仿真结果： 2.当电源电压)30100 sin(10)(+=t t x πV 时，求解表达式y(t)的波形：
y x dt t dy 314)(-+=，其中y(0)=1。要求交mdl 文件及仿真结果的截图。 解：用Simulink 连接的仿真图如下所示： 仿真结果如下图（b ）所示： 3.使用S 函数实现y=5*x+3,建立仿真模型并得出仿真结果。要求交mdl 文件，S 函数及仿真结果的截图。 解：用Simulink 连接的仿真图如下所示： S 函数的代码如下： function [sys,x0,str,ts] = mysfun3 (t,x,u,flag) switch flag, case 0, [sys,x0,str,ts]=mdlInitializeSizes;
case 1, sys=mdlDerivatives(t,x,u); case 2, sys=mdlUpdate(t,x,u); case 3, sys=mdlOutputs(t,x,u); case 4, sys=mdlGetTimeOfNextVarHit(t,x,u); case 9, sys=mdlTerminate(t,x,u); otherwise error(['Unhandled flag = ',num2str(flag)]); end function [sys,x0,str,ts]=mdlInitializeSizes sizes = simsizes; sizes.NumContStates = 0; sizes.NumDiscStates = 0; sizes.NumOutputs = 1; sizes.NumInputs = 1; sizes.DirFeedthrough = 1; sizes.NumSampleTimes = 1; sys = simsizes(sizes); x0 = []; str = []; ts = [0 0]; function sys=mdlDerivatives(t,x,u) sys = []; function sys=mdlUpdate(t,x,u) sys = []; function sys=mdlOutputs(t,x,u) sys = 5*u+3; function sys=mdlGetTimeOfNextVarHit(t,x,u) sampleTime = 1; sys = t + sampleTime; function sys=mdlTerminate(t,x,u) sys = []; 仿真结果如下图（c）所示：
s函数是system Function的简称，用它来写自己的simulink模块。（够简单吧，^_^，详细的概念介绍大伙看帮助吧）可以用matlab、C、C++、Fortran、Ada等语言来写，这儿我只介绍怎样用matlab语言来写吧（主要是它比较简单）先讲讲为什么要用s函数，我觉得用s函数可以利用matlab的丰富资源，而不仅仅局限于simulink提供的模块，而用c 或c++等语言写的s函数还可以实现对硬件端口的操作，还可以操作windows API等的先介绍一下simulink的仿真过程（以便理解s函数），simulink的仿真有两个阶段： 1.初始化：这个阶段主要是设置一些参数，像系统的输入输出个数、状态初值、采样时间等； 2.运行阶段：这个阶段里要进行计算输出、更新离散状态、计算连续状态等等。这个阶段需要反复运行，直至结束。 在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.NumInp
程序代码：(代码标记[code]...[/code] ) 1-32是：图形应用篇 33-66是：界面设计篇 67-84是：图形处理篇 85-100是：数值分析篇 实例1：三角函数曲线（1） function shili01 h0=figure('toolbar','none',... 'position',[198 56 350 300],... 'name','实例01'); h1=axes('parent',h0,... 'visible','off'); x=-pi:0.05:pi; y=sin(x); plot(x,y); xlabel('自变量X'); ylabel('函数值Y'); title('SIN( )函数曲线'); grid on 实例2：三角函数曲线（2） function shili02 h0=figure('toolbar','none',... 'position',[200 150 450 350],... 'name','实例02'); x=-pi:0.05:pi; y=sin(x)+cos(x); plot(x,y,'-*r','linewidth',1); grid on xlabel('自变量X'); ylabel('函数值Y'); title('三角函数');
实例3：图形的叠加 function shili03 h0=figure('toolbar','none',... 'position',[200 150 450 350],... 'name','实例03'); x=-pi:0.05:pi; y1=sin(x); y2=cos(x); plot(x,y1,... '-*r',... x,y2,... '--og'); grid on xlabel('自变量X'); ylabel('函数值Y'); title('三角函数'); 实例4：双y轴图形的绘制 function shili04 h0=figure('toolbar','none',... 'position',[200 150 450 250],... 'name','实例04'); x=0:900;a=1000;b=0.005; y1=2*x; y2=cos(b*x); [haxes,hline1,hline2]=plotyy(x,y1,x,y2,'semilogy','plot'); axes(haxes(1)) ylabel('semilog plot'); axes(haxes(2)) ylabel('linear plot'); 实例5：单个轴窗口显示多个图形 function shili05 h0=figure('toolbar','none',... 'position',[200 150 450 250],... 'name','实例05'); t=0:pi/10:2*pi;
MATLAB S函数编写示范（MATLAB2010a） 2010-05-05 21:22:47| 分类：MATLAB之控制系统| 标签：matlab s-function 编写教程 s builder |字号大中小订阅 1.M语言S函数 在命令行输入命令“edit sfuntmpl”即可打开M语言s函数模板。M语言S函数是以flag为标志进行运算的。flag==0时，初始化参数；flag==1时，计算微分；flag==2时，状态更新；flag==3时，计算输出；flag==4时适用于变步长结算；其他参数不常用。 （1）flag==0时，调用mdlInitializeSizes函数，sizes结构及解释如下： sizes.NumContStates = 0; //连续变量个数 sizes.NumDiscStates = 0; //离散变量个数 sizes.NumOutputs = 0; //输出个数 sizes.NumInputs = 0; //输入个数 sizes.DirFeedthrough = 1; //是否直通，当在输出里包含了输入时，为直通 sizes.NumSampleTimes = 1; //采样时间 对状态变量进行初始化,如果初始条件是0初始条件，则 x0=[0 0] （2）flag==1时，调用mdlDerivatives(t,x,u)函数计算微分，如有一系统为状态方程为x'=2x+u，y=5x+u时，编写如下：sys=2*x+u; （3）flag==3时，计算输出，如上例系统，则编写如下：sys=5*x+u; ok,下面只需要将编写好的s函数名称写入s-function模块。 例子：对如下系统（零初始条件） x1'=-3x1-x2+u; x2'=2x1; y=x2;编写s函数如下： function [sys,x0,str,ts,simStateCompliance] = sfun_m_hyj(t,x,u,flag) switch flag, case 0, [sys,x0,str,ts,simStateCompliance]=mdlInitializeSizes; case 1, sys=mdlDerivatives(t,x,u); case 2, sys=mdlUpdate(t,x,u); case 3, sys=mdlOutputs(t,x,u); case 4, sys=mdlGetTimeOfNextVarHit(t,x,u); case 9, sys=mdlTerminate(t,x,u);