《建模与仿真及其医学应用》

《《建模与仿真及其医学应用》》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、《建模与仿真及其医学应用》摘要：sys为原系统，sysb(sysh)为平衡实现系统，g为平衡对角线矩阵，T为状态变换矩阵，Ti是前者...ode23Runge-Kutta法三阶积分算法、二阶误差估计、变积分步长的低阶算法...关键词：矩阵,算法类别：专题技术来源：牛档搜索（Niudown.COM）　　本文系牛档搜索（Niudown.COM）根据用户的指令自动搜索的结果，文中内涉及到的资料均来自互联网，用于学习交流经验，作品其著作权归原作者所有。不代表牛档搜索（Niudown.COM）赞成本文的内容或立场，牛档搜索（Ni

2、udown.COM）不对其付相应的法律责任！26《建模与仿真及其医学应用》实验讲义天津医科大学生物医学工程系2004年26实验一系统建模的MATLAB实现一、实验目的：1．学习MATLAB基本知识。2．掌握数学模型的MATLAB实现：时域模型、状态空间模型和零极点模型。3．学习用MATLAB实现系统外部模型到内部模型的转换。4．学习用MATLAB实现系统模型的连接：串联、并联、反馈连接。5．了解模型降阶的MATLAB实现。二、实验内容1．系统的实现、外部模型到内部模型的转换(1)给定连续系统的传递函数，利用MATLAB

3、建立传递函数模型，微分方程，并转换为状态空间模型。(2)已知某系统的状态方程的系数矩阵为：利用MATLAB建立状态空间模型，并将其转换为传递函数模型和零极点模型。(3)已知系统的零极点传递函数为，利用MATLAB转换为传递函数模型和状态空间模型。2．系统的离散、连接、降阶（1）给定连续系统的传递函数26，将该连续系统的传递函数用零阶重构器和一阶重构器转换为离散型传递函数，抽样时间T=1秒。（2）该系统与系统分别①串联②并联③负反馈连接，求出组成的新系统的传递函数模型。（3）将串联组成的新系统进行降阶处理，求出降阶后系统

4、的模型，并用plot图形比较降阶前后系统的阶跃响应。要求：将以上过程用MATLAB编程（M文件）实现，运行输出结果。三、实验说明—关于系统建模的主要MATLAB函数1．建立传递函数模型：tf函数：格式：sys=tf(num,den)num=[bm,bm-1,……,b0]分子多项式系数den=[an,an-1,……,a0]分母多项式系数2．建立状态空间模型：ss函数：格式：sys=ss(a,b,c,d)%a,b,c,d为状态方程系数矩阵sys=ss(a,b,c,d,T)%产生离散时间状态空间模型3．建立零极点模型的函数：

5、zpk格式：sys=zpk(z,p,k)4．模型转换函数：tf2sstf2zpss2tfss2zpzp2tfzp2ss%2为to的意思格式：[a,b,c,d]=tf2ss(num,den)[z,p,k]=tf2zp(num,den)[num,den]=ss2tf(a,b,c,d,iu)%iu指定是哪个输入[z,p,k]=ss2zp(a,b,c,d,iu)]26[num,den]=zp2tf(z,p,k)[a,b,c,d]=zp2ss(z,p,k)5．模型的连接串联：sys=series(sys1,sys2)并联：sys

6、=parallel(sys1,sys2)反馈连接：sys=feedback(sys1,sys2,sign)%负反馈时sign可忽略；正反馈时为1。6．系统扩展：把若干个子系统组成系统组。格式：sys=append(sys1,sys2,…)7．模型降阶（1）基于平衡的状态空间实现--balreal格式：sysb=balreal(sys)[sysh,g,T,Ti]=balreal(sys)sys为原系统，sysb(sysh)为平衡实现系统，g为平衡对角线矩阵，T为状态变换矩阵，Ti是前者的逆矩阵。两种格式的区别：前者只给出

7、原系统的一个平衡的状态空间实现，而后者还给出平衡实现的对角线矩阵g，从中可以看出哪个状态变量该保留，哪个状态变量该删去，从而实现降阶。（2）降阶的实现—modred格式：rsys=modred(sys,elim)rsys=modred(sys,elim,’mde’)rsys=modred(sys,elim,’del’)强调：这里的sys应是函数balreal()变换的模型，elim为待消去的状态，’mde’指降阶中保持增益匹配，’del’指降阶中不保持增益匹配。268．连续系统模型离散化函数：C2DMConversio

8、nofcontinuousLTIsystemstoiscrete-time.格式：①[Ad,Bd,Cd,Dd]=C2DM(A,B,C,D,Ts,'method')将连续系统状态空间—离散系统状态空间'method':'zoh'零阶重构器zeroorderhold'foh'一阶重构器firstorderhold②[NUMd,DENd