关于matlab中s-function参数设置的两个例子讲解详尽受益匪浅

《关于matlab中s-function参数设置的两个例子讲解详尽受益匪浅》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、例10.1设系统的微分方程为：'x1x2t'0.5txxe22试建立系统仿真模型。2例10.2利用Simulink构建函数曲线y=5t+16。1例10.3利用Simulink仿真求I=xln(1x)dx。0例10.4有初始状态为0的二阶微分方程x"+0.2x'+0.4x=0.2u(t),其中u(t)是单位阶跃函数，试建立系统模型并仿真。例10.5PID控制器是在自动控制中经常使用的模块，在工程应用中其标准的数学模型为1TdsU(s)K(1)E(s)pTsTs/Nid其中采用

2、了一阶环节来近似纯微分动作，为保证有良好的微分近似的效果，一般选N≥10。试建立PID控制器的模型并建立子系统。例10.6利用使能子系统构成一个正弦半波整流器。例10.7利用触发子系统将一锯齿波转换成方波。例10.8采用S函数实现模块y=nx，即模块的功能是把一个输入信号n倍以后再输出。(1)利用MATLAB语言编写S函数。程序如下：%*************************************************%S函数timesn.m，其输出是输入的n倍%***********

3、**************************************function[sys,x0,str,ts]=timesn(t,x,u,flag,n)switchflag,case0%初始化[sys,x0,str,ts]=mdlInitializeSizes;case3%计算输出量sys=mdlOutputs(t,x,u,n);case{1,2,4,9}%未使用的flag值sys=[];otherwise%出错处理error(['Unhandleflag=',num2str(flag)

4、]);end%*************************************************%mdlInitializeSizes：当flag为0时进行整个系统的初始化%*************************************************function[sys,x0,str,ts]=mdlInitializeSizes(T)%调用函数simsizes以创建结构体sizessizes=simsizes;%用初始化信息填充结构体sizessizes.N

5、umContStates=0;%无连续状态sizes.NumDiscStates=0;%无离散状态sizes.NumOutputs=1;%有一个输出量sizes.NumInputs=1;%有一个输入信号sizes.DirFeedthrough=1;%输出量中含有输入量sizes.NumSampleTimes=1;%单个采样周期%根据上面的设置设定系统初始化参数sys=simsizes(sizes);%给其他返回参数赋值。x0=[];%设置初始状态为零状态str=[];%将str变量设置为空字符串ts

6、=[-1,0];%假定继承输入信号的采样周期%初始化子程序结束%*************************************************%mdlOutputs：当flag值为3时，计算输出量%*************************************************functionsys=mdlOutputs(t,x,u,n)sys=n*u;%输出量计算子程序结束。(2)模块的封装与测试。例10.9采用S函数来构造非线性分段函数。31xx31

7、x33(xx3)234y24x52(xx5)25616x(1)利用MATLAB语言编写S函数。程序如下：function[sys,x0,str,ts]=sfunction(t,x,u,flag)switchflag,case0,[sys,x0,str,ts]=mdlInitializeSizes;case3,sys=mdlOutputs(t,x,u);case{1,2,4,9}2sys=[];otherwiseerror(['Unhandledflag='

8、,num2str(flag)]);endfunction[sys,x0,str,ts]=mdlInitializeSizessizes=simsizes;sizes.NumContStates=0;sizes.NumDiscStates=0;sizes.NumOutputs=1;sizes.NumInputs=1;sizes.DirFeedthrough=1;sizes.NumSampleTimes=1;sys=simsizes(sizes);x0=[];str=[]