《建模与仿真及其医学应用》
实 验 讲 义
天津医科大学生物医学工程系
2004年
实验一 系统建模的MATLAB实现
一、实验目的:
1.学习MATLAB基本知识。
2.掌握数学模型的MATLAB实现:时域模型、状态空间模型和零极点模型。
3.学习用MATLAB实现系统外部模型到内部模型的转换。 4.学习用MATLAB实现系统模型的连接:串联、并联、反馈连接。 5.了解模型降阶的MATLAB实现。 二、实验内容
1.系统的实现、外部模型到内部模型的转换
(s?8)(s2?2s?5)(1)给定连续系统的传递函数G(s)?,利用2(2s?3)(3s?4s?13)MATLAB建立传递函数模型,微分方程,并转换为状态空间模型。 (2)已知某系统的状态方程的系数矩阵为:
?10??00?1??0?10??11? d??10? c? a??b???11???????2?3??????02???01??利用MATLAB建立状态空间模型,并将其转换为传递函数模型和零极点模型。
(3)已知系统的零极点传递函数为G(s)?2(s?1),利用
(s?2)(s?3)(s?4)MATLAB转换为传递函数模型和状态空间模型。 2.系统的离散、连接、降阶
(s?8)(s2?2s?5)(1)给定连续系统的传递函数G(s)?,将该连2(2s?3)(3s?4s?13)续系统的传递函数用零阶重构器和一阶重构器转换为离散型传递函数,抽样时间T=1秒。 (2)该系统与系统H(s)?1分别①串联②并联③负反馈连2s?6s?5接,求出组成的新系统的传递函数模型。
(3)将串联组成的新系统进行降阶处理,求出降阶后系统的模型,并用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.建立零极点模型的函数:zpk 格式:sys=zpk(z,p,k) 4.模型转换函数:
tf2ss tf2zp ss2tf ss2zp zp2tf zp2ss %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)]