第六章(上) 系统模型.ppt

1、第六章 系统模型,本章主要内容如下： 6.1 系统仿真概述 6.2 控制系统建模 6.3 系统模型的MATLAB表示 6.4 系统模型的连接 6.5 系统建模举例,6.1 系统仿真概述,6.1.1 系统仿真及其分类 系统仿真的定义 “仿真”译自英文Simulation，意指在实际系统尚不存在的情况下，系统或活动本质的复现。在工程技术中是指通过对系统模型的实验，研究一个存在的或设计中的系统。 系统仿真的分类,6.1 系统仿真概述,6.1.2 仿真模型与仿真研究 仿真模型,6.1 系统仿真概述,6.1.2 仿真模型与仿真研究 计算机仿真过程 建模 所建立的计算机模型（仿真数学模型）应与对象的功能和

2、参数之间具有相似性和对应性 模型实现 利用优秀的算法将计算机模型编制成可运行的计算机程序（MATLAB软件）。 仿真分析 通过运行仿真程序，对仿真结果进行分析。,6.2 控制系统建模,6.2.1 控制系统的模型表示 微分方程形式,设线性定常系统输入、输出量是单变量，分别为u(t),y(t),模型参数形式为：,输出系统向量 ， n+1维,输入系统向量 ， m+1维,(1),6.2 控制系统建模,6.2.1 控制系统的模型表示 传递函数形式,在零初始条件下，将(1) 方程两边进行拉氏变换，则有,(2),模型参数可表示为,传递函数分母系数向量:,传递函数分子系数向量: B=b0 , b1, .bm,

3、用num=B,den=A分别表示分子，分母参数向量，则可简练的表示为(num,den)。,6.2 控制系统建模,6.2.1 控制系统的模型表示 零、极点增益形式,将(2)中的分子，分母分解为因式连乘形式，则有,(3),模型参数可表示为： 系统零点向量： 系统极点向量： 简记为(Z,P,K)形式，称为零极点增益三对组模型参数。,6.2 控制系统建模,6.2.1 控制系统的模型表示 状态空间形式,当控制系统输入、输出为多变量时，可用向量分别表示为U(t), Y(t),系统的内部状态变量为X(t).,(4),模型参数形式为： 系统系数矩阵 A，系统输入矩阵 B 系统输出矩阵 C，直接传输矩阵 D 简

4、记为(A,B,C,D)形式。,6.2 控制系统建模,6.2.2 控制系统建模的基本方法 机理模型法 采用由一般到特殊的推理演绎方法，对已知结构、参数的物理系统运用相应的物理定律或定理，经过合理分析简化而建立起来的描述系统各物理量动、静态变化性能的数学模型。 实验模型法 采用由特殊到一般的逻辑、归纳方法,根据一定数量的在系统运行过程中实测、观察的数据,运用统计规律、系统辨识等理论合理估计出反应实际系统各物理量相互制约关系的数学模型. 混合模型法 当对控制的内部结构和特性有部分了解，但又难以完全用机理模型的方法表述出来，这是需要结合一定的实验方法确定另外一部分不甚了解的结构特性，或是通过实际测定来

5、求取模型参数。这种方法是机理模型法和实验模型法的结合，故称为混合模型法。,6.3 系统模型的MATLAB表示,6.3.1 MATLAB建模 传递函数模型 MATLAB建模 其中，num为传递函数分子系数向量，den为传递函数分母系数向量。,sys = tf(num,den),6.3 系统模型的MATLAB表示,6.3.1 MATLAB建模 零极点增益模型 MATLAB建模 其中，Z、P、K分别为系统的零点向量、极点向量和增益。,sys = zpk(Z,P,K),6.3 系统模型的MATLAB表示,6.3.1 MATLAB建模 时间延迟系统模型 G(s)=G1(s)e-s 其中 G1(s)为系统

6、无时延时的模型传递函数，为延迟时间。 MATLAB建模 sys=tf(num,den, InputDelay,tao) sys=zpk(z,p,k, InputDelay,tao),【说明】 InputDelay为关键词，也可写成OuputDelay，对于线性SISO系统，二者是等价的。 tao为系统延迟时间的数值。,6.3 系统模型的MATLAB表示,6.3.2 系统模型的转换 newsys=tf(sys) 将非传递函数形式的系统模型sys转化成传递函数模型newsys. newsys=zpk(sys) 将非零极点增益形式的系统模型sys转化成零极点增益模型 newsys.,6.4 系统模型

7、的连接,6.4.1 模型串联 两个线性模型串联及其等效模型如图所示，且 sys = sys1sys2 MATLAB对串联模型的运算如下 sys=series( sys1,sys2) 上式可等价写成： sys=sys1*sys2,6.4 系统模型的连接,6.4.2 模型并联 两个线性模型并联及其等效模型如图所示，且 sys = sys1+sys2 MATLAB对并联模型的运算如下 sys=parallel( sys1,sys2) 上式可等价写成： sys=sys1+sys2,6.4 系统模型的连接,6.4.3 反馈连接 两个线性模型反馈连接及其等效模型如图所示. MATLAB对模型反馈连接的运算

8、如下 sys=feedback( sys1,sys2,sign) 说明：sign表示反馈连接符号：负反馈连接sign= -1，正反馈连接sign=1。 上式可等价写成： sys=minreal(sys1/(1+sys2*sys1) (负反馈) sys=minreal(sys1/(1-sys2*sys1) (正反馈) 说明：minreal()用来对消传递函数中相同的零极点。,【例5-12】反馈连接运算演示：G前向环节、H反馈环节,负反馈。,6.5 系统建模举例,6.5.1 机械加速度计建模(教材101页),加速度计的物理模型如图所示，其质量m的位移y近似与被测运动物体ms的加速度d2x/dt2成

9、正比，现求加速度计输出y与运动物体的作用力f之间的动力学关系。,加速度计的力平衡方程,2. 求传递函数,用step(Y_F)命令(见教材第6章),可得该加速度计的单位阶跃响应.由图可知,加速度计的稳态输出与输入力成比例,也即与ms的加速度成比例.,6.5.2 磁悬浮系统建模(教材103页),磁悬浮是通过调节磁场的磁性,使被控轴始终悬浮于轴承的中心,而不与轴承接触,可使这类轴承没有摩擦损失。,图示的基本磁悬浮系统中，电磁力大小可由电流i控制。浮球的位置由光探测器检测，e=kex为探测器的输出；V0为电磁力的预设值以平衡浮球重力mg；u为反馈控制信号；为作用在浮球上的外部扰动力。,1. 建立系统动力学方程,电磁力,线圈电流,控制电压,浮球的力平衡方程,系统动力学方程,(比例+微分控制),2. 求浮球位移对扰动的传递函数,设m =20 g，ki =0.5 N/A，kx =20 N/m，ke =100V/m，Kd=8，Kp=100，可得该系统的MATLAB模型:,m=20; ki=0.5; kx=20; ke=100; Kd=8; Kp=100; de