matlabppt课件-第3章_控制系统的数学模型及其转换.ppt

《matlabppt课件-第3章_控制系统的数学模型及其转换.ppt》由会员分享，可在线阅读，更多相关《matlabppt课件-第3章_控制系统的数学模型及其转换.ppt（97页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第三章第三章 控制系统的数学模型及其转换控制系统的数学模型及其转换6/8/20231本章内容本章内容u 利用利用MATLAB描述在控制系统中常见的几种数学模型；描述在控制系统中常见的几种数学模型；u 利用利用MATLAB实现任意数学模型之间的相互转换；实现任意数学模型之间的相互转换；u 利用利用MATLAB求解系统经过串联、并联和反馈连接后的系求解系统经过串联、并联和反馈连接后的系统模型；统模型；u 利用利用MATLAB获取一些典型系统的模型；获取一些典型系统的模型；u 利用利用MATLAB实现连续系统的离散化和离散系统的连续化；实现连续系统的离散化和离散系统的连续化；u 利用利用MATLAB

2、求取系统的特性函数。求取系统的特性函数。6/8/20232l控制系统计算机仿真是建立在控制系统数学模型基础之上的一门技术。对系统进行仿真，首先应该知道系统的数学模型，然后才可以在此基础上设计一个合适的控制器，使得原系统的响应达到预期的效果。在线性系统理论中，一般常用的数学模型形式有：传递函数模型传递函数模型（系统的外部模型）状态方程模型状态方程模型（系统的内部模型）零极点增益模型零极点增益模型（传递函数模型的一种）这些模型之间都有着内在的联系，可以相互进行转换。6/8/20233系统类型系统类型一一.连续和离散系统连续和离散系统根据系统自变量(时间)是连续变化还是离散变化，系统分为连续系统和离

3、散系统。(1)连续系统系统输入、输出信号都是连续时间信号。(2)离散系统系统输入、输出信号都是离散时间信号。(3)混合系统系统输入、输出信号包含连续信号和离散信号。连续时间的数学模型用微分方程描述。离散时间系统的数学模型用差分方程描述。例如：一般L、R、C电路都是连续时间系统。数字计算机是典型的离散时间系统。实际上离散时间系统经常与连续时间系统组合运用，此时称为混合系统(或采样系统)。如自动控制系统和数字通信系统。6/8/20234二二二二.线性和非线性系统线性和非线性系统线性和非线性系统线性和非线性系统 根据输入输出关系是否同时满足齐次性和叠加性，系统分为线性和非线性。假设系统在没有外界信号

4、作用之前处于静止状态，在输入信号为任意实数，和或和作用下，有式中，为输入输出之间函数关系。那么，该系统称为线性系统，否则是非线性系统。根据模型参数是否随时间变化，线性系统又可细分为线性定常系统和线性时变系统。参数不随时间变化的系统，称为时不变系统或定常系统，否则称为时变系统。例如，线性定常系统：6/8/20235线性时变系统：非线性定常系统：式中，分别为系统输入、输出。6/8/20236三三.确定和随机系统确定和随机系统 根据系统输入、输出和内部状态呈现的规律，系统分为确定性系统与随机性系统。输入输出之间函数关系能够用确定性模型描述的系统，称为确定性系统，否则称为随机系统(或不确定性系统)。在

5、控制系统中，随机因素可能作用在系统的入口，也可能作用在系统的出口，还可能影响系统模型本身。例如噪声输入的随机控制系统，其微分方程具有如下形式。式中，分别为状态变量和输出变量，为输入噪声，通常是正态分布的白噪声。6/8/20237控制系统常用数学模型控制系统常用数学模型 根据系统输入、输出与内部状态变量之间关系，控制系统模型可分为外部模型和内部模型。一般地，把建立系统输入输出关系的数学模型称为外部模型，包括时域模型和频域模型。把建立系统输入、输出与内部状态变量之间关系的数学模型称为内部模型，相应的数学模型称为系统的状态空间方程(连续状态空间方程或离散状态空间方程)。控制系统模型分类控制系统模型分

6、类6/8/20238一一一一.连续系统连续系统连续系统连续系统 1微分方程微分方程一个连续系统可以表示成高阶微分方程，即式中，分别为系统输入量、输出量，n为系统的阶次，为系统的结构参数，结构参数，为输入函数的及各阶导数的初始值为它们均为实常数。已知输出变量控制系统常用数学模型控制系统常用数学模型 6/8/20239稍加整理，并记 2.传递函数传递函数 若系统的初始条件为零，即系统在t=0 时已处于一个稳定状态，那么对式(3.1)两边取拉普拉斯变换后可得 称为系统的传递函数。6/8/2023103.状态空间描述状态空间描述 微分方程和传递函数都只描述了系统输入与输出之间的关系，而没有考虑系统内部

7、状态的动态运动，仅仅实现系统输入与输出之间的关系是不够的，还必须复现模型的内部变量即状态变量的动态变化规律。状态空间描述考虑了“输入-状态-输出”这一过程，引进经典控制理论所忽略的中间内部状态，因此状态空间表达式能够完全反映系统的全部独立变量的变化，而且能够方便处理初始条件。6/8/202311 在用状态空间法分析系统时，系统的动态特性是用由状态变量构成的一阶微分方程组来描述。状态空间表达式包括状态方程和输出方程。线性定常系统的状态空间描述为 =AX+BU 状态方程(3.3)Y=CX+DU 输出方程(3.4)维系统矩阵式中 维输入矩阵 维输出矩阵 维直接传递矩阵 6/8/202312 r维输入

8、向量，n维状态向量，m维输出向量 6/8/202313 对于线性时变系统，系数矩阵对于线性时变系统，系数矩阵A，B，C，D，均与时间，均与时间 t 有有关，状态空间描述为关，状态空间描述为 系统的状态空间分析法是时域内的一种矩阵运算方法，不再只局限于输入量、输出量、误差量，为提高系统性能提供了有力工具。状态空间分析法特别适合于用计算机来计算，有利于把工程技术人员从烦琐的计算中解脱出来。6/8/202314 二二二二.离散系统离散系统离散系统离散系统1差分方程差分方程设系统差分方程为(3.6)引进后移算子为 (3.7)控制系统常用数学模型控制系统常用数学模型 假设系统输入、输出 及其内部状态变量

9、 均是时间序列，其中T 为离散时间间隔，为书写简便，用 数学模型有3种形式。表示。与连续时间系统类似，离散时间系统6/8/202315式(3.6)可写为(3.8)令 则(3.9)6/8/2023162离散传递函数离散传递函数(Z传函传函)假设系统的初始条件为零，即 则得(3.10)系统传递函数为(3.11)在初始条件为零时，与等价。6/8/202317 3.离散状态空间模型离散状态空间模型类似在连续系统中，从微分方程或传递函数建立状态空间表达式，叫做系统的实现。在离散系统中，从差分方程或脉冲传递函数求取离散状态空间表达式，也是一种实现。多变量离散状态空间表达式(3.12)6/8/2023183

10、.1 3.1 控制系统的数学模型控制系统的数学模型 MATLAB的控制系统工具箱是提供自动控制系统建模、分析和设计方面函数的集合，提供传递函数模型、零极点增益模型、状态空间模型三种形式线性时不变(LTI)模型。有关模型表示的函数如表3-1所示。函数功能sys=tf(num，den)生成传递函数模型sys=zpk(z，p，k)生成零极点增益模型sys=ss(a，b，c，d)生成状态空间模型6/8/2023193.1.1 传递函数nnum=b1,b2,bm,bm+1nden=a1,a2,an,an+1n注意：它们都是按s的降幂进行排列的。6/8/202320MATLABMATLAB输入语句输入语句

11、6/8/202321 MATLAB Control 工具箱中，用命令tf()可以建立一个传递函数模型，或将零极点增益模型和状态空间模型变化为传递函数模型。tf()函数调用格式如下：n sys=tf(num，den)；%用于生成连续传递函数(S传递函数)；n sys=tf(num，den，Ts)；%用于生成离散传递函数(Z传递函数)；n sys=tf(num，den，Property1，Value1，.，PropertyN，ValueN)；%用于生成具有LTI模型属性的传递函数；n sys=tf(num，den，Ts，Property1，Value1，.，PropertyN，ValueN)；%用

12、于生成具有LTI模型属性的传递函数；n sys=tf(s)；%用于生成拉普拉斯变量s有理传递函数；n sys=tf(z，Ts)；%用于生成采样周期为Ts的z有理传递函数；传递函数模型命令tf()调用格式6/8/202322传递函数输入举例传递函数输入举例【例【例【例【例3.13.13.13.1】输入传递函数模型输入传递函数模型输入传递函数模型输入传递函数模型MATLABMATLAB输入语句输入语句 num=6 12 6 10;den=1 2 3 1 1;sys1=tf(num,den)Transfer function:6 s3+12 s2+6 s+10-s4+2 s3+3 s2+s+16/8

13、/202323 当传递函数的分子或分母由若干个多项式乘积表示时，它可由MATLAB 提供的多项式乘法运算函数conv()来处理，以便获得分子和分母多项式向量，此函数的调用格式为 c c=conv(=conv(a,ba,b)其中a和b分别为由两个多项式系数构成的向量，而c为a和b多项式的乘积多项式系数向量。conv()函数的调用是允许多级嵌套的。6/8/202324【例【例3-23-2】若给定系统的传递函数为【解】【解】则可以将其用下列MATLAB语句表示 num=4*conv(1 2,1 6 6)；den=conv(1 0,conv(1 1,conv(1 1,conv(1 1,1 3 2 5)

14、；G=tf(num,den)6/8/202325【例【例3-23-2】更简便的输入法更简便的输入法s=tf(s);%用于生成拉普拉斯变量用于生成拉普拉斯变量s有理传递函数有理传递函数G=4*(s+2)*(s2+6*s+6)/(s*(s+1)3*(s3+3*s2+2*s+5)6/8/202326【例【例3-33-3】对于单输入多输出系统 【解】【解】则可将其用下列MATLAB语句表示num=3 2;1 0 2 5;den=3 5 2 1;tf(num,den)Transfer function from input to output.3 s+2#1:-3 s3+5 s2+2 s+1 s3+2

15、s+5#2:-3 s3+5 s2+2 s+16/8/202327 【例3.4】给定SISO系统输入为“flow”，输出为“Temp”，传递函数为使用MATLAB表示该传递函数num=1.3 2 2.5;den=1 0.5 1.2 1;sys=tf(num,den,inputdelay,2,inputName,flow,OutputName,Temp)运行结果：Transfer function from input flow to output Temp:1.3 s2+2 s+2.5exp(-2*s)*-s3+0.5 s2+1.2 s+16/8/202328 【例3.5】给定一个多输入-多输出

16、MIMO系统结果为：Transfer function from input 1 to output.1#1:-s+1#2:2 Transfer function from input 2 to output.1#1:-s+2 s+1#2:-s+2MATLAB命令：num=1 1;2 1 1;den=1 1,1 2;1,1 2;sys3=tf(num,den)6/8/202329 【例3.6】若一采样周期为0.2s的离散MIMO传递函数为结果如下：Transfer function from input 1 to output .z+1#1:-z2+2 z+1 1#2:-2 z+1 Trans

17、fer function from input 2 to output.Transfer z#1:-z2+2 2#2:-z+1Sampling time:0.2MATLAB命令如下：num=1 1,1 0;1,2;den=1 2 1,1 0 2;2 1,1 1;sys4=tf(num,den,0.2)6/8/2023303.1.2 零极点增益形式n 零极点模型实际上是传递函数模型的另一种表现形式，其原理是分别对原系统传递函数的分子、分母进行分解因式处理，以获得系统的零点和极点的表示形式。MATLABMATLAB输入语句输入语句K为系统增益，zi为零点，pj为极点6/8/202331 在MATL

18、AB Control 工具箱中，用命令zpk()可以建立零极点增益模型，或将传递函数模型和状态空间模型变化为零极点增益模型。n sys=zpk(z，p，k)n sys=zpk(z，p，k，Ts)n sys=zpk(z，p，k，Property1，Value1，.，PropertyN，ValueN)n sys=zpk(z，p，k，Ts，Property1，Value1，.，PropertyN，ValueN)n sys=zpk(s)n sys=zpk(z)zpk()函数调用方法与tf()一致。零极点增益模型zpk函数调用格式6/8/202332 【例3.7】给定一零极点增益模型，使用MATLAB表

19、示该传递函数。MATLAB命令如下：z=;-0.5;p=0.3;0.1-j,0.1+j;k=1;2;sys=zpk(z,p,k,)结果如下：Zero/pole/gain from input to output.1#1:-(z-0.3)2(z+0.5)#2:-(z2 -0.2z+1.01)Sampling time:unspecified6/8/202333l MATLABMATLAB工具箱中的函数工具箱中的函数 poly()poly()和和roots()roots()可用来实现多项式和零极点间的转换，例如在命可用来实现多项式和零极点间的转换，例如在命令窗口中进行如下操作可实现互相转换。令窗口

20、中进行如下操作可实现互相转换。P=1 3 5 2;R=roots(P)R=-1.2267+1.4677i -1.2267-1.4677i -0.5466P1=poly(R)P1=1.0000 3.0000 5.0000 2.00006/8/202334例：建立下述传递函数模型的matlab表示 num=12,24,0,20;den=2 4 6 2 2;sys=tf(num,den)借助多项式乘法函数借助多项式乘法函数conv来处理：来处理：num=4*conv(1,2,conv(1,6,6,1,6,6);den=conv(1,0,conv(1,1,conv(1,1,conv(1,1,1,3,2

21、,5);sys=tf(num,den)2）1)6/8/2023353）零极点增益模型：4）零极点增益模型：z=;p=-1,-2,-3-4j,-3+4j;k=5;sys=zpk(z,p,k)z=-3,0;p=-1,50,-10;k=1;sys=zpk(z,p,k)6/8/2023363.1.3 部分分式形式传递函数表示成部分分式或留数形式 为极点，为各极点对应的留数，为分子除以分母的余子式。6/8/2023372.1.4 状态空间表达式 状态空间模型(state-space model:SS)线性定常状态空间模型描述为式中，为状态向量，为输入向量，是输出向量。6/8/202338 在MATLAB

22、 Control 工具箱中，用命令ss()可以建立状态空间模型，或将传递函数模型和零极点增益模型转化为状态空间模型。语法调用格式:n sys=ss(a，b，c，d)n sys=ss(a，b，c，d，Ts)n sys=ss(a，b，c，d，Property1，Value1，.，PropertyN，ValueN)n sys=ss(a，b，c，d，Ts，Property1，Value1，.，PropertyN，ValueN)a，b，c，d是状态空间模型系数矩阵。sys为状态空间模型对象。ss函数的调用方法与tf()、zpk()一致。6/8/202339【例3.8】设系统的状态空间表达式为A=0 0

23、1;-3/2 -2 -1/2;-3 0 -4;B=1 1;-1 -1;-1 -3;C=1 0 0;0 1 0;D=zeros(2,2);sys=ss(A,B,C,D)6/8/202340 【例3.9】状态空间模型结果为：a=x1 x2 x3 x1 0 1 0 x2 0 0 1 x3 -5 -20 -1 b=u1 x1 0 x2 0 x3 1 c=x1 x2 x3 y1 1 0 0 d=u1 y1 0 Continuous-time model。用MATLAB表示为A=0 1 0;0 0 1;-5-20-1;B=0;0;1;C=1 0 0;D=0;sys=ss(A,B,C,D)6/8/20234

24、1 传递函数模型、零极点增益模型和状态空间模型传递函数模型、零极点增益模型和状态空间模型之间模型转换函数如图所示。之间模型转换函数如图所示。3.2 3.2 系统数学模型间的相互转换系统数学模型间的相互转换6/8/202342 在不同的场合需要不同的模型，因此系统的数学模型有多种表示方式，这些数学模型需要相互转换。MATLAB control工具箱提供了一系列用于模型转换的函数如表3-2所示。函数功能sys=tf(ss_sys);sys=tf(zpk_sys);ss或zpk模型转换传递函数模型sys=zpk(ss_sys);sys=zpk(tf_sys)ss或tf模型生成零极点增益模型sys=s

25、s(tf_sys);sys=ss(zpk_sys)tf或zpk模型生成状态空间模型A，B，C，D=tf2ss(num，den)传递函数转换成状态空间模型num，den=ss2tf(A，B，C，D，iu)状态空间模型转换成传递函数z，p，k=tf2zp(num，den)传递函数转换成零极点增益模型num，den=zp2tf(z，p，k)零极点增益模型转换成传递函数z，p，k=ss2zp(A，B，C，D，iu)状态空间模型转换成零极点增益模型A，B，C，D=zp2ss(z，p，k)零极点增益模型转换成状态空间模型6/8/202343 MATLAB实现模型转换有两种不同的方式。方式方式1：简单的模型

26、转换：简单的模型转换 首先生成任一指定的模型对象(tf，ss，zpk)，然后将该模型对象类作为输入，调用欲转换的模型函数即可。例如：欲将传递函数转换为状态空间模型sys=tf(num，den)；%生成传递函数模型。a，b，c，d=ss(sys)；%将传递函数模型作为状态空间模型的输入。6/8/202344方式方式2：直接调用模型转换函数：直接调用模型转换函数 连续模型之间三种形式的数学模型相互转换函数类型包括tf2ss、ss2tf、zp2tf、tf2zp、ss2zp、zp2ss，共六个函数，如图3.2所示。函数的调用格式和功能见表3-2。传递函数模型、零极点增益模型和状态空间模型之间模型转换函

27、数如图3.2所示。6/8/202345 3.2.1 3.2.1 状态空间表达式到传递函数的转换状态空间表达式到传递函数的转换 在在MATLABMATLAB控制系统工具箱中，给出一个根据状态控制系统工具箱中，给出一个根据状态空间表达式求取系统传递函数的函数空间表达式求取系统传递函数的函数ss2tf()ss2tf()，其，其调用格式为调用格式为 num,den=ss2tf(A,B,C,D,iu)n 对多输入的系统，必须具体化iu。例如，如果系统有三个输入(u1，u2，u3)，则iu必须为1、2或3中的一个，其中1表示u1，2表示u2，3表示u3。n 如果系统只有一个输入，则可采用 num，den=

28、ss2tf(A，B，C，D)或 num，den=ss2tf(A，B，C，D，1)6/8/202346*【例3-10】对于例3-8中给出的多变量系统，可以由下面的命令分别对各个输入信号求取传递函数向量，然后求出这个传递函数阵。【解】利用下列MATLAB语句num1,den1=ss2tf(A,B,C,D,1)num1=0 1.0000 5.0000 6.0000 0-1.0000-5.0000-6.0000 den1=1 6 11 66/8/202347num2,den2=ss2tf(A,B,C,D,2)num2=0 1.0000 3.0000 2.0000 0 -1.0000 -4.0000 -

29、3.0000 den2=1 6 11 6则可得系统的传递函数阵6/8/2023483.2.2 3.2.2 状态空间形式到零极点形式的转换状态空间形式到零极点形式的转换MATLAB函数ss2zp()的调用格式为 Z,P,K=ss2zp(A,B,C,D,iu)iu值的设置与ss2tf()用法相同。6/8/202349 【例3-11】将状态空间模型转换为传递函数和零极点增益模型。MATLAB命令：A=0 1 0;0 0 1;-5-20-1;B=0;0;1;C=1 0 0;D=0;num,den=ss2tf(A,B,C,D);%状态空间模型转换成传递函数模型printsys(num,den)%输出传递

30、函数模型z,p,k=ss2zp(A,B,C,D);%状态空间模型转换成零极点增益模型zpk(z,p,k)%生成零极点增益模型对象类6/8/202350 【例3-13】多入单出系统状态空间表达式包括两个传递函数：Y1(s)/U1(s)、Y1(s)/U2(s)(当考虑输入u1时，可设u2为零。反之亦然)，试将其化成传递函数形式。MATLAB命令：A=0 1 0;0 1 1;-5.008 25.1026 -5.032471;B=0 1;25.04 2;121.005 3;C=1 0 0;D=0 0;num,den=ss2tf(A,B,C,D,1);sys1=tf(num,den)num,den=ss

31、2tf(A,B,C,D,2);sys2=tf(num,den)6/8/202351 传递函数Y1(s)/U1(s)为Transfer function:-3.553e-015 s2+25.04 s+247 -s3+4.032 s2 30.14 s+5.008传递函数Y1(s)/U2(s)为Transfer function:s2+6.032 s 17.07 -s3+4.032 s2 30.14 s+5.008以上就是下列两个传递函数的MATLAB表达式：6/8/202352 2.2.3 2.2.3 传递函数到状态空间表达式的转换传递函数到状态空间表达式的转换 如果已知系统的传递函数模型，求取系

32、统状态空间表达式的过程又称为系统的实现。由于状态变量可以任意地选取，所以实现的方法并不是唯一的，这里只介绍一种比较常用的实现方法。l在MATLAB控制系统工具箱中称这种方法为能控标准型实现方法，并给出了直接实现函数，该函数的调用格式为 A,B,C,D=tf2ss(num,den)6/8/202353 【例3-14】将传递函数化成状态空间表达式num=0 2 3;1 2 1;den=1 0.4 1;A,B,C,D=tf2ss(num,den)结果为：A=-0.4000 -1.0000 1.0000 0B=1 0C=2.0000 3.0000 1.6000 0D=0 16/8/2023543.2.

33、4 3.2.4 传递函数形式到零极点形式的转换传递函数形式到零极点形式的转换 MATLAB函数tf2zp()的调用格式为 Z,P,K=tf2zp(num,den)3.2.5 3.2.5 零极点形式到状态空间表达式的转换零极点形式到状态空间表达式的转换 MATLAB函数zp2ss()的调用格式为 A,B,C,D=zp2ss(Z,P,K)3.2.6 3.2.6 零极点形式到传递函数形式的转换零极点形式到传递函数形式的转换 MATLAB函数zp2tf()的调用格式为 num,den=zp2tf(Z,P,K)6/8/202355 【例3-15】将传递函数化为零极点增益模型MATLAB代码命令如下。nu

34、m=2 1;den=3 4 1;z,p,k=tf2zp(num,den)结果：z=-0.5000 p=-1.0000 -0.3333 k=0.6667因此，零点为-1/2，极点为-1/3和-1，增益为2/3，因此零极点增益模型为6/8/202356 【例3-16】零极点增益模型化成状态空间表达式MATLAB代码命令如下。z=-2;p=-1-3-3;k=2;A,B,C,D=zp2ss(z,p,k)6/8/2023573.2.7 3.2.7 部分分式与传递函数形式之间的转换部分分式与传递函数形式之间的转换 num,denresidue(P,R,H):部分分式模型转换为 传递函数模型 P,R,H=r

35、esidue(num,den):传递函数模型转换为 部分分式模型6/8/202358 sys=append(sys1，sys2，.，sysN)子系统合成对角形式子系统合成对角形式sys=parallel(sys1，sys2，inp1，inp2，out1，out2)并联连接并联连接sys=series(sys1，sys2，outputs1，inputs2)串联连接串联连接sys=feedback(sys1，sys2)反馈连接反馈连接sys=lft(sys1，sys2，nu，ny)模型连接模型连接sysc=connect(sys，Q，inputs，outputs)框图建模框图建模3.3 3.3 系

36、统模型的连接系统模型的连接 在一般情况下，控制系统常常由若干个环节通过串联、并联和反馈连接的方式而组成，对在各种连接模式下的系统能够进行分析就需要对系统的模型进行适当的处理,在MATLAB的控制系统工具箱中提供了大量的对控制系统的简单模型进行连接的函数。6/8/202359 ，3.3.1 3.3.1 串联连接串联连接 在在MATLABMATLAB的控制系统工具箱中提供了系统的串联连的控制系统工具箱中提供了系统的串联连接处理函数接处理函数series()series()，它既可处理由状态方程表示的，它既可处理由状态方程表示的系统，也可处理由传递函数阵表示的单输入多输出系统。系统，也可处理由传递函

37、数阵表示的单输入多输出系统。SISO系统模型的串连MIMO系统模型的串连6/8/202360n函数调用格式：sys=series(sys1,sys2)：求SISO系统sys1和sys2的串连模型sys=series(sys1,sys2,outputs1,inputs2):求MIMO系统的串连，其中outputs1为sys1的输出端口数，inputs2为sys2的输入端口数。若MIMO系统有多个连接，在outputs1和inputs2为同维数向量。a,b,c,d=series(a1,b1,c1,d1,a2,b2,c2,d2)a,b,c,d=series(a1,b1,c1,d1,a2,b2,c2,

38、d2,output1,input2)num,den=series(num1,den1,num2,den2)num,den=series(num1,den1,num2,den2,output1,input2)6/8/202361 【例【例3-173-17】求下列两系统串联后的系统模型 A1=2 3;-1 4;B1=1;0;C1=2 4;D1=1;A2=0 3;-3-1;B2=0;1;C2=1 3;D2=2;A,B,C,D=series(A1,B1,C1,D1,A2,B2,C2,D2)6/8/202362 对于单入单出对于单入单出SISOSISO系统，并联连接是指各子系统输入相系统，并联连接是指各

39、子系统输入相同，输出是各子系统输出之代数和。如图所示。同，输出是各子系统输出之代数和。如图所示。3.3.2 3.3.2 并联连接并联连接 在MATLAB的控制系统工具箱中提供了系统的并联连接处理函数parallel()。SISO系统模型的并连MIMO系统模型的并连6/8/202363na,b,c,d=parallel(a1,b1,c1,d1,a2,b2,c2,d2)na,b,c,d=parallel(a1,b1,c1,d1,a2,b2,c2,d2,in1,in2,out1,out2)nnum,den=parallel(num1,den1,num2,den2)其调用格式如下：nsys=paral

40、lel(sys1,sys2)：两个SISO系统模型的并联nsys=parallel(sys1,sys2,in1,in2,out1,out2)：生成两个MIMO系统的并联模型，并联输入端口由向量in1和in2定义，并联的输出端口由向量out1和out2定义。in1和out1对应sys1用于并联的输入、输出端口向量，in2和out2对应sys2用于并联的输入、输出端口向量。6/8/202364【例3-18】在上图所示的SISO系统并联模型结构中，已知sys1、sys2的传递函数分别为：s=tf(s);%定义拉普拉斯变量ssys1=(s+1)/(s*(s2+s+2);%定义SISO系统sys1的传递

41、函数模型sys2=5/(5*s+1);%定义SISO系统的sys2的传递函数模型disp(并联系统模型为)sys=parallel(sys1,sys2)%生成sys1与sys2的并联模型执行上述命令后在Command Windows得如下结果：并联系统模型为Transfer function:5 s3+10 s2+16 s+1-5 s4+6 s3+11 s2+2 s6/8/202365 【例3-19】求下列两系统并联后的系统模型 num1=3;den1=1 4;num2=2 4;den2=1 2 3;num,den=parallel(num1,den1,num2,den2);G=tf(num,

42、den)Transfer function:5 s2+18 s+25-s3+6 s2+11 s+126/8/202366 系统模型的反馈连接分SISO系统模型反馈连接和MIMO系统模型反馈连接。3.3.3 3.3.3 反馈连接反馈连接SISO系统模型的反馈MIMO系统模型的反馈6/8/202367nsys=feedback(sys1,sys2,sign)：生SISO系统sys1与sys2构成的反馈系统传递函数模型，sys2为反馈通道传递函数模型，反馈类型由sign指定，当sign=1为正反馈，当sign=-1为负反馈，此时sign可省略。nsys=feedback(sys1,sys2,feed

43、in,feedout,sign)：生成MIMO系统的部分反馈模型，feedin指定sys1中连入反馈回路的输入端口向量。feedout指定sys1中连入反馈回路中的输出端口向量。feedbackfeedback（）调用格式（）调用格式 na,b,c,d=feedback(a1,b1,c1,d1,a2,b2,c2,d2,sign)na,b,c,d=feedback(a1,b1,c1,d1,a2,b2,c2,d2,inp1,out1)nnum,den=feedback(num1,den1,num2,den2,sign)6/8/202368【例】给定一个5输入4输出的状态空间子系统sys1和一个3输

44、入2输出的状态空间反馈控制器sys2。sys1的输出1、3和4连接到控制器sys2的输入端，控制器sys2的输出端连接到子系统sys1的输入端4和2。则采用命令如下：feedin=4 2;feedout=1 3 4;cloop=feedback(sys1，sys2，feedin，feedout);6/8/202369 【例【例3-203-20】求下列两系统反馈连接后的系统模型numg=2 5 1;deng=1 2 3;numh=5 10;denh=1 10;num,den=feedback(numg,deng,numh,denh);printsys(num,den)6/8/2023703.3.

45、4 3.3.4 系统的组合系统的组合 MATLAB MATLAB的组合函数的组合函数append()append()的调用格式为的调用格式为 A,B,C,D=append(A1,B1,C1,D1,A2,B2,C2,D2)6/8/202371 3.3.5 3.3.5 根据框图建模根据框图建模 利用connect()函数，可以根据系统的方框图按指定方式求取系统模型，其函数调用格式为：A1,B1,C1,D1=connect(A,B,C,D,Q,inputs,outputs)其中（A,B,C,D）为由函数append()生成的无连接对角方块系统的状态空间模型系数矩阵；Q矩阵用于指定系统（A,B,C,D

46、）的内部连接关系，矩阵的每一行对应于一个有连接关系的输入，其第一个元素为输入编号，其后为连接该输入的输出编号，如采用负连接，则以负值表示；inputs和outputs用于指定系统（A1,B1,C1,D1）的输入和输出的编号；（A1,B1,C1,D1）为在指定输入和输出并按要求的内部连接关系下所生成的系统。6/8/202372 【例【例3-213-21】用框图建模法重做【例用框图建模法重做【例3-173-17】6/8/202373 【例【例3-213-21】用框图建模法求系统的传递函数。用框图建模法求系统的传递函数。6/8/202374 【例【例3-223-22】用框图建模法求系统的传递函数。用

47、框图建模法求系统的传递函数。6/8/202375 【例【例3-233-23】如图如图3-93-9所示的系统连接关系，建立该框图所示的系统连接关系，建立该框图的状态空间表达式和传递函数模型。的状态空间表达式和传递函数模型。图图3-9 3-9 用方框图表示系统的连接关系用方框图表示系统的连接关系6/8/202376Step 1:Step 1:定义子系统数学模型定义子系统数学模型(传递函数传递函数tftf、状态空间模型、状态空间模型ssss或零极点增益模型或零极点增益模型zpk)zpk)。6/8/202377Step 2Step 2：建立无连接的状态空间模型，形成一个无连接关系：建立无连接的状态空间

48、模型，形成一个无连接关系的对角块的对角块(a(a，b b，c c，d)d)。6/8/202378Step 3Step 3：指定框图间的连接关系。：指定框图间的连接关系。6/8/202379Step 4Step 4：选择输入选择输入/输出。输出。6/8/202380Step 5Step 5：调用：调用connectconnect函数，完成框图模型的状态空间模型。函数，完成框图模型的状态空间模型。sysc=connect(sys,Q,inputs,outputs);sysc=connect(sys,Q,inputs,outputs);6/8/202381Step 6Step 6：将状态空间模型转换

49、成传递函数模型。将状态空间模型转换成传递函数模型。6/8/202382 6/8/2023833.4 3.4 典型系统的生成典型系统的生成1.建立二阶系统模型 可利用 MATLAB所提供的函数 ord2()来建立其调用格式为 num,den=ord2(n,)或 A,B,C,D=ord2(n,)6/8/202384 【例【例3-243-24】已知=0.4，n=2.4rad/s,求二阶系统的传递函数。num,den=ord2(2.4,0.4);printsys(num,den)num/den=1 -s2+1.92 s+5.766/8/2023852.纯时延系统的Pad近似可利用 MATLAB所提供的

50、函数 pade()来建立其调用格式为 num,den=pade(tau，n)6/8/2023863.5 3.5 系统的特性值系统的特性值 在分析控制系统的时候，经常用到系统的一些在分析控制系统的时候，经常用到系统的一些特性函数，如系统的增益、阻尼系数和自然频率等特性函数，如系统的增益、阻尼系数和自然频率等等，等，MATLABMATLAB的控制系统工具箱中提供了相应的函数的控制系统工具箱中提供了相应的函数用来计算系统的特性函数，如下表所示。用来计算系统的特性函数，如下表所示。函 数 名功 能damp()求系统的阻尼系数和自然频率ddamp()求离散系统的阻尼系数和自然频率dcgain()求连续控