线性系统数学模型.pptx

《线性系统数学模型.pptx》由会员分享，可在线阅读，更多相关《线性系统数学模型.pptx（125页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、知 识 要 点 线性系统的数学模型，拉普拉斯变换，传递函数的定义，非线性特性的线性化处理，方框图的简化，梅逊公式的含义和应用。第1页/共125页描述控制系统输入、输出变量以及内部各变量之间关系的数学表达式，称为系统的数学模型。常用的数学模型有微分方程、差分方程、传递函数、脉冲传递函数和状态空间表达式等。建立合理的数学模型，对于系统的分析研究是至关重要的。系统数学模型的建立，一般采用解析法或实验法。第2页/共125页目 录v2.1 线性系统的微分方程v2.2 微分方程的线性化v2.3 传递函数v2.4 方框图v2.5 信号流图v2.6 在MATLAB中数学模型的表示 v小 结第3页/共125页

2、2.1 线性系统的微分方程（1）分析系统工作原理，将系统划分为若干环节，确定系统和环节的输入、输出变量，每个环节可考虑列写一个方程；（2）根据各变量所遵循的基本定律(物理定律、化学定律)或通过实验等方法得出的基本规律，列写各环节的原始方程式，并考虑适当简化和线性化；（3）将各环节方程式联立，消去中间变量，最后得出只含输入、输出变量及其导数的微分方程；（4）将输出变量及各阶导数放在等号左边，将输入变量及各阶导数放在等号右边，并按降幂排列，最后将系统归化为具有一定物理意义的形式，成为标准化微分方程。第4页/共125页例2-1 试列写图中所示RC无源网络的微分方程。输入为ui(t)，输出为u0(t)

3、。解 根据基尔霍夫定理，可列出以下式子：第5页/共125页整理得：令T1=R1C1，T2=R2C2，T3=R1C2则得该网络的数学模型是一个二阶线性常微分方程。第6页/共125页例2-2 图为一弹簧阻尼系统，当外力F(t)作用于系统时，系统将产生运动。试列写外力F(t)与位移y(t)之间的微分方程。第7页/共125页解 弹簧和阻尼器有相应的弹簧阻力F1(t)和粘性摩擦阻力F2(t)，根据牛顿第二定律有：其中F1(t)和F2(t)可由弹簧、阻尼器特性写出式中k 弹簧系数f 阻尼系数第8页/共125页整理且标准化令称为时间常数；称为阻尼比；称为放大系数。得第9页/共125页例2-3 电枢控制的它激

4、直流电动机如图所示，电枢输入电压u0(t)，电动机输出转角为。Ra、La、ia(t)分别为电枢电路的电阻、电感和电流，if为恒定激磁电流，eb为反电势，f为电动机轴上的粘性摩擦系数，G为电枢质量，D为电枢直径，ML为负载力矩。第10页/共125页解电枢回路电压平衡方程为ce为电动机的反电势系数力矩平衡方程为式中为电动机电枢的转动惯量为电动机的力矩系数第11页/共125页整理得无量纲放大系数电机转速电磁时间常数机电时间常数时间常数电机传递系数第12页/共125页无量纲放大系数。时间常数电机传递系数第13页/共125页例2-4热水电加热系统，如图所示，为减小周围空气的热损耗，槽壁是绝热的，控温元件

5、是电动控温开关。第14页/共125页根据能量守恒定律其中Qh加热器供给的热量；QC贮槽内水吸收的热量；Q0热水流出槽所带走的热量:Qi冷水进入槽带入的热量:Ql隔热壁逸散的热量:C贮槽水的热容量；V流出槽水的流量；H水的比热；R热阻；Ti进入槽水的温度；T槽内水的温度；Te槽周围空气温度。第15页/共125页整理得 一般情况下，描述线性定常系统输入与输出关系的微分方程为:或返回第16页/共125页2.2 微分方程的线性化 实际的物理系统往往有间隙、死区、饱和等各类非线性现象。严格地讲，几乎所有实际物理和化学系统都是非线性的。目前，线性系统的理论已经相当成熟，但非线性系统的理论还远不完善。因此，

6、在工程允许范围内，尽量对所研究的系统进行线性化处理，然后用线性理论进行分析不失为一种有效的方法。第17页/共125页 当非线性因素对系统影响较小时，一般可直接将系统当作线性系统处理。另外，如果系统的变量只发生微小的偏移，则可通过切线法进行线性化，以求得其增量方程式。第18页/共125页 非线性函数的线性化，是指将非线性函数在工作点附近展开成泰勒级数，忽略掉高阶无穷小量及余项，得到近似的线性化方程，来替代原来的非线性函数。第19页/共125页 假如元件的输出与输入之间关系x2=f(x1)的曲线如图，元件的工作点为(x10，x20)。将非线性函数x2=f(x1)在工作点(x10，x20)附近展开成

7、泰勒级数第20页/共125页当(x1x10)为微小增量时，可略去二阶以上各项，写成 其中为工作点(x10，x20)处的斜率，即此时以工作点处的切线代替曲线，得到变量在工作点的增量方程，经上述处理后，输出与输入之间就成为线性关系。第21页/共125页图2-8为一铁芯线圈，输入为ui(t)，输出为i(t)。线圈的微分方程为第22页/共125页当工作过程中线圈的电压和电流只在工作点（u0,i0）附近变化时，即有线圈中的磁通对也有增量变化，假如在i0附近连续可微，将在i0附近展开成泰勒级数，即因是微小增量，将高阶无穷小量略去，得近似式第23页/共125页这就是铁芯线圈的增量化方程，为简便起见，常略去增

8、量符号而写成返回第24页/共125页2.3 传递函数 2.2.1 传递函数 在零初始条件下，线性定常系统输出量的拉普拉斯变换与输入量的拉普拉斯变换之比，定义为线性定常系统的传递函数。即，第25页/共125页若已知线性定常系统的微分方程为式中c(t)为输出量，r(t)为输入量。设c(t)和r(t)及其各阶导数初始值均为零，对式(2-47)取拉氏变换，得第26页/共125页则系统的传递函数为或写为传递函数与输入、输出之间的关系，可用图表示。第27页/共125页2.2.2 传递函数的特点 1.作为一种数学模型，传递函数只适用于线性定常系统，这是由于传递函数是经拉普拉斯变换导出的，而拉氏变换是一种线性

9、积分运算。2.传递函数是以系统本身的参数描述的线性定常系统输入量与输出量的关系式，它表达了系统内在的固有特性，只与系统的结构、参数有关，而与输入量或输入函数的形式无关。第28页/共125页3.传递函数可以是无量纲的，也可以是有量纲的，视系统的输入、输出量而定，它包含着联系输入量与输出量所必须的单位，它不能表明系统的物理特性和物理结构。许多物理性质不同的系统，有着相同的传递函数，正如一些不同的物理现象可以用相同的微分方程描述一样。4.传递函数只表示单输入和单输出(SISO)之间的关系，对多输入多输出(MIMO)系统，可用传递函数阵表示。第29页/共125页5.传递函数式(2-49)可表示成式中p

10、1，p2pn为分母多项式的根，称为传递函数的极点；z1、z2、zn为分子多项式的根，称为传递函数的零点；第30页/共125页6.传递函数分母多项式称为特征多项式，记为而D(s)=0称为特征方程。传递函数分母多项式的阶次总是大于或等于分子多项式的阶次，即nm。这是由于实际系统的惯性所造成的。第31页/共125页2.2.3 典型环节的传递函数 控制系统由许多元件组合而成，这些元件的物理结构和作用原理是多种多样的，但抛开具体结构和物理特点，从传递函数的数学模型来看，可以划分成几种典型环节，常用的典型环节有比例环节、惯性环节、积分环节、微分环节、振荡环节、延迟环节等。第32页/共125页1.比例环节环

11、节输出量与输入量成正比，不失真也无时间滞后的环节称为比例环节，也称无惯性环节。输入量与输出量之间的表达式为c(t)=Kr(t)比例环节的传递函数为式中K为常数，称为比例环节的放大系数或增益。第33页/共125页2.惯性环节(非周期环节)惯性环节的动态方程是一个一阶微分方程其传递函数为式中T惯性环节的时间常数 K惯性环节的增益或放大系数第34页/共125页当输入为单位阶跃函数时，其单位阶跃响应为单位阶跃响应曲线第35页/共125页惯性环节实例很多，如图所示的R-L网络，输入为电压u，输出为电感电流i，其传递函数式中第36页/共125页2.积分环节输出量正比于输入量的积分的环节称为积分环节，其动态

12、特性方程其传递函数式中Ti为积分时间常数。第37页/共125页积分环节的单位阶跃响应为它随时间直线增长，当输入突然消失，积分停止，输出维持不变，故积分环节具有记忆功能，如图所示。第38页/共125页上图为运算放大器构成的积分环节，输入ui(t)，输出u0(t)，其传递函数为式中Ti=RC第39页/共125页4.微分环节理想微分环节的特征输出量正比于输入量的微分，其动态方程其传递函数式中Td称微分时间常数它的单位阶跃响应曲线第40页/共125页如图所示，理想微分环节实际上难以实现，因此我们常采用带有惯性的微分环节，其传递函数其单位阶跃响应为第41页/共125页曲线如下图所示，实际微分环节的阶跃响

13、应是按指数规律下降，若K值很大而Td值很小时，实际微分环节就愈接近于理想微分环节。第42页/共125页5.二阶振荡环节(二阶惯性环节)二阶振荡环节的动态方程为其传递函数式中为无阻尼自然振荡角频率，为阻尼比，在后面时域分析中将详细讨论。第43页/共125页图中所示为RLC网络，输入为ui(t)、输出u0(t)，其动态特性方程其传递函数式中第44页/共125页6.延迟环节(时滞环节)延迟环节是输入信号加入后，输出信号要延迟一段时间后才重现输入信号，其动态方程为其传递函数是一个超越函数式中称延迟时间第45页/共125页需要指出，在实际生产中，有很多场合是存在迟延的，比如皮带或管道输送过程、管道反应和

14、管道混合过程，多个设备串联以及测量装置系统等。迟延过大往往会使控制效果恶化，甚至使系统失去稳定。返回第46页/共125页 2.4 方框图 在控制工程中，为了便于对系统进行分析和设计，常将各元件在系统中的功能及各部分之间的联系用图形来表示，即方框图和信号流图。第47页/共125页2.4.1方框图 方框图也称方块图或结构图，具有形象和直观的特点。系统方框图是系统中各元件功能和信号流向的图解，它清楚地表明了系统中各个环节间的相互关系。构成方框图的基本符号有四种，即信号线、比较点、传递环节的方框和引出点。第48页/共125页第49页/共125页2.4.2系统方框图的构成对于一个系统在清楚系统工作原理及

15、信号传递情况下，可按方框图的基本连接形式，把各个环节的方框图，连接成系统方框图。例2-5图中为一无源RC网络。选取变量如图所示，根据电路定律，写出其微分方程组为第50页/共125页第51页/共125页零初始条件下，对等式两边取拉氏变换，得第52页/共125页RC网络方框图各环节方框图第53页/共125页例2-6图中为电枢电压控制的直流他励电动机,描述其运动方程为第54页/共125页零初始条件下，对式中两边取拉氏变换第55页/共125页 将同一变量的信号线连接起来，将输入Ua(s)放在左端，输出(s)放在图形右端，得系统方框图如图所示。第56页/共125页2.4.3环节间的连接 环节的连接有串联

16、、并联和反馈三种基本形式。1.串联:在单向的信号传递中，若前一个环节的输出就是后一个环节的输入，并依次串接如图2-32所示，这种联接方式称为串联。n个环节串联后总的传递函数:第57页/共125页即环节串联后总的传递函数等于串联的各个环节传递函数的乘积。环节的串联RC网络第58页/共125页2.并联:若各个环节接受同一输入信号而输出信号又汇合在一点时，称为并联。如图2-34所示。由图可知总的传递函数为环节的并联第59页/共125页3.反馈：若将系统或环节的输出信号反馈到输入端，与输入信号相比较，就构成了反馈连接，如图所示。如果反馈信号与给定信号极性相反，则称负反馈连接。反之，则为正反馈连接，若反

17、馈环节H(s)=1称为单位反馈。反馈连接第60页/共125页反馈连接后，信号的传递形成了闭合回路。通常把由信号输入点到信号输出点的通道称为前向通道；把输出信号反馈到输入点的通道称为反馈通道。对于负反馈连接，给定信号r(t)和反馈信号b(t)之差，称为偏差信号e(t)即通常将反馈信号B(s)与误差信号E(s)之比，定义为开环传递函数，即开环传递函数=第61页/共125页输出信号C(s)与偏差信号E(s)之比，称为前向通道传递函数，即前向通道传递函数=而系统输出信号C(s)与输入信号R(s)之比称为闭环传递函数，记为(s)或GB(s)。第62页/共125页得闭环传递函数为对于正反馈连接，则闭环传递

18、函数为第63页/共125页2.4.4方框图的变换和简化 有了系统的方框图以后，为了对系统进行进一步的分析研究，需要对方框图作一定的变换，以便求出系统的闭环传递函数。方框图的变换应按等效原则进行。所谓等效，即对方框图的任一部分进行变换时，变换前、后输入输出总的数学关系式应保持不变。除了前面介绍的串联、并联和反馈连接可以简化为一个等效环节外，还有信号引出点及比较点前后移动的规则。第64页/共125页例例 2-7化 简图(a)所 示系统方框图，并求系统传递函数第65页/共125页第66页/共125页 图2-37(a)是一个交错反馈多路系统，采用引出点后移或前移，比较点前移等，逐步变换简化，可求得系统

19、的闭环传递函数为例2-8 试化简如图2-37(a)所示系统的方框图，并求闭环传递函数。第67页/共125页图2-37方框图的变换与简化第68页/共125页返回第69页/共125页 2.5 信号流图 信号流图是表示线性方程组变量间关系的一种图示方法，将信号流图用于控制理论中，可不必求解方程就得到各变量之间的关系，既直观又形象。当系统方框图比较复杂时，可以将它转化为信号流图，并可据此采用梅逊(Mason)公式求出系统的传递函数。第70页/共125页考虑如下简单等式这里变量xi和xj可以是时间函数、复变函数，aij是变量xj变换(映射)到变量xi的数学运算，称作传输函数，如果xi和xj是复变量s的函

20、数，称aij为传递函数Aij(s)，即上式写为2.5.1信号流图的定义第71页/共125页变量xi和xj用节点“”来表示，传输函数用一有向有权的线段(称为支路)来表示，支路上箭头表示信号的流向，信号只能单方向流动。信号流图第72页/共125页2.5.2系统的信号流图 在线性系统信号流图的绘制中应包括以下步骤：（1）将描述系统的微分方程转换为以s为变量的代数方程。（2）按因果关系将代数方程写成如下形式：第73页/共125页（3）用节点“”表示n个变量或信号，用支路表示变量与变量之间的关系。通常把输入变量放在图形左端，输出变量放在图形右端。例2-9如上图所示的电阻网络，v1为输入、v3为输出。选5

21、个变量v1、i1、v2、i2、v3，由电压、电流定律可写出四个独立方程第74页/共125页 将变量V1(s)、I1(s)、V2(s)、I2(s)、V3(s)作节点表示，由因果关系用支路把节点与节点联接，得信号流图。第75页/共125页2.5.3信号流图的定义和术语 节点：表示变量或信号的点，用“”表示。支路：连接两个节点之间的有向有权线段，方向用箭头表示，权值用传输函数表示。输入支路：指向节点的支路。输出支路：离开节点的支路。源节点：只有输出支路的节点，也称输入节点，如图中节点X1。汇节点：只有输入支路的节点，如图节点X7。第76页/共125页信号流图定义与术语混合节点：既有输入支路、又有输出

22、支路的节点，如图中的X2、X3、X4、X5、X6。通道(路径)：沿着支路箭头方向通过各个相连支路的路径，并且每个节点仅通过一次。如X1到X2到X3到X4或X2到X3又反馈回X2。第77页/共125页前向通道：从输入节点(源节点)到汇节点的通道。如图X1到X2到X3到X4到X5到X6到X7为一条前向通道，又如X1到X2到X3到X5到X6到X7也为另一条前向通道。闭通道(反馈通道或回环)：通道的起点就是通道的终点，如图X2到X3又反馈到X2；X4到X5又反馈到X4。自回环：单一支路的闭通道，如图中的-H3构成自回环。第78页/共125页通道传输或通道增益：沿着通道的各支路传输的乘积。如从X1到X7

23、前向通道的增益G1G2G3G4G5G6。不接触回环：如果一些回环没有任何公共的节点，称它们为不接触回环。如G2H1与G4H2。第79页/共125页2.5.4信号流图的性质（1）信号流图只适用于线性系统；（2）信号流图所依据的方程式，一定为因果函数形式的代数方程；（3）信号只能按箭头表示的方向沿支路传递；（4）节点上可把所有输入支路的信号叠加，并把总和信号传送到所有输出支路；（5）具有输入和输出支路的混合节点，通过增加一个具有单位传输的支路，可把其变为输出节点，即汇节点；（6）对于给定的系统，其信号流图不是唯一的。第80页/共125页2.5.5信号流图的简化（1）加法规则：n个同方向并联支路的总

24、传输，等于各个支路传输之和，如图(a)所示：（2）乘法规则：n个同方向串联支路的总传输，等于各个支路传输之积，如图（b）。第81页/共125页（3）混合节点可以通过移动支路的方法消去，如图（c）。（4）回环可根据反馈连接的规则化为等效支路，如图（d）。第82页/共125页例2-10 2-10 将图2-43所示系统方框图化为信号流图并化简求出系统的闭环传递函数第83页/共125页解：信号流图如图(a)所示。化G1与G2串联等效为G1G2支路，G3与G4并联等效为G3+G4支路，第84页/共125页如图(b)，G1G2与-H1反馈简化为 支路，又与G3+G4串联，等效为 如图(c)第85页/共12

25、5页进而求得闭环传递函数为第86页/共125页2.5.6信号流图的增益公式 给定系统信号流图之后，常常希望确定信号流图中输入变量与输出变量之间的关系，即两个节点之间的总增益或总传输。上节采用信号流图简化规则，逐渐简化，最后得到总增益或总传输。但是，这样很费时又麻烦，而梅逊(Mason)公式可以对复杂的信号流图直接求出系统输出与输入之间的总增益，或传递函数，使用起来更为方便。第87页/共125页梅逊增益公式可表示为式中，T 输出和输入之间的增益或传递函数；Pk 第k条前向通道的增益或传输函数；信号流图的特征值，Lj1所有不同回环增益之和；Lj2所有两两互不接触回环增益乘积之和；Lj3所有三个互不

26、接触回环增益乘积之和k 与第k条前向通道不接触的那部分信号流图的，称为第k条前向通道特征式的余子式。第88页/共125页例2-11利用梅逊公式求图中所示系统的传递函数 C(s)/R(s)。第89页/共125页解：输入量R(s)与输出量C(s)之间有三条前向通道，对应Pk与k为P1=G1G2G3G4G51=1P2=G1G6G4G52=1P3=G1G2G7G53=1P4=-G1G6G2G7G54=1图中有五个单回环，其增益为：L1=-G3H2，L2=-G5H1，L3=-G2G3G4G5H3，L4=-G6G4G5H3，L5=-G2G7G5H3，其中L1与L2是互不接触的，其增益之积L1L2=G3G5

27、H1H2第90页/共125页系统的特征式为系统的传递函数为第91页/共125页例2-12求图示信号流图的闭环传递函数解：系统单回环有：L1=G1，L2=G2，L3=G1G2，L4=G1G2，L5=G1G2系统的特征式为：第92页/共125页前向通道有四条：P1=-G11=1 P2=G22=1 P3=G1G23=1 P4=G1G24=1系统的传递函数为返回第93页/共125页2.6 在MATLAB中数学模型的表示 控制系统的数学模型在系统分析和设计中是相当重要的,在线性系统理论中常用的数学模型有微分方程、传递函数、状态空间表达式等,而这些模型之间又有着某些内在的等效关系。MATLAB主要使用传递

28、函数和状态空间表达式来描述线性时不变系统(LinearTimeInvariant简记为LTI)。第94页/共125页2.6.1传递函数 单输入单输出线性连续系统的传递函数为其中mn。G(s)的分子多项式的根称为系统的零点,分母多项式的根称为系统的极点。令分母多项式等于零,得系统的特征方程:D(s)=a0sn+a1sn1+an1s+an=0第95页/共125页 因传递函数为多项式之比,所以我们先研究MATLAB是如何处理多项式的。MATLAB中多项式用行向量表示,行向量元素依次为降幂排列的多项式各项的系数,例如多项式P(s)=s3+2s+4,其输入为P=1024注意尽管s2项系数为0,但输入P(

29、s)时不可缺省0。MATLAB下多项式乘法处理函数调用格式为C=conv(A,B)第96页/共125页例 如 给 定 两 个 多 项 式 A(s)=s+3和B(s)=10s2+20s+3,求C(s)=A(s)B(s),则应先构造多项式A(s)和B(s),然后再调用conv()函数来求C(s)A=1,3;B=10,20,3；C=conv(A,B)C=1050639即得出的C(s)多项式为10s3+50s2+63s+9第97页/共125页MATLAB提供的conv()函数的调用允许多级嵌套,例如G(s)=4(s+2)(s+3)(s+4)可由下列的语句来输入G=4*conv(1,2,conv(1,3

30、,1,4)第98页/共125页 有 了 多 项 式 的 输 入,系 统 的 传 递 函 数 在MATLAB下可由其分子和分母多项式唯一地确定出来,其格式为sys=tf(num,den)其中num为分子多项式，den为分母多项式num=b0,b1,b2,bm；den=a0,a1,a2,an；第99页/共125页对于其它复杂的表达式,如可由下列语句来输入num=conv(1,1,conv(1,2,6,1,2,6);den=conv(1,0,0,conv(1,3,1,2,3,4);G=tf(num,den)Transferfunction:第100页/共125页2.6.2传递函数的特征根及零极点图

31、传递函数G(s)输入之后,分别对分子和分母多项式作因式分解,则可求出系统的零极点,MATLAB提供了多项式求根函数roots(),其调用格式为roots(p)其中p为多项式。第101页/共125页例如,多项式p(s)=s3+3s2+4p=1,3,0,4;%p(s)=s3+3s2+4r=roots(p)%p(s)=0的根r=-3.35330.1777+1.0773i0.1777-1.0773i 反过来,若已知特征多项式的特征根,可调用MATLAB中的poly()函数,来求得多项式降幂排列时各项的系数,如上例poly(r)p=1.00003.00000.00004.0000第102页/共125页

32、而polyval函数用来求取给定变量值时多项式的值,其调用格式为polyval(p,a)其中p为多项式;a为给定变量值例如,求n(s)=(3s2+2s+1)(s+4)在s=5时值：n=conv(3,2,1,1,4);value=polyval(n,-5)value=66第103页/共125页p,z=pzmap(num,den)其中,p传递函数G(s)=numden的极点z传递函数G(s)=numden的零点例如,传递函数传递函数在复平面上的零极点图,采用pzmap()函数来完成,零极点图上,零点用“。”表示,极点用“”表示。其调用格式为第104页/共125页 用MATLAB求出G(s)的零极点

33、,H(s)的多项式形式,及G(s)H(s)的零极点图numg=6,0,1;deng=1,3,3,1;z=roots(numg)z=0+0.4082i00.4082i%G(s)的零点p=roots(deng)p=1.0000+0.0000i1.0000+0.0000i%G(s)的极点1.0000+0.0000i第105页/共125页n1=1,1;n2=1,2;d1=1,2*i;d2=1,-2*i;d3=1,3;numh=conv(n1,n2);denh=conv(d1,conv(d2,d3);printsys(numh,denh)numh/denh=%H(s)表达式pzmap(num,den)%

34、零极点图title(pole-zeroMap)第106页/共125页零极点图如图所示：第107页/共125页2.6.3 控制系统的方框图模型 若已知控制系统的方框图,使用MATLAB函数可实现方框图转换。1.串联串联如图所示G1(s)和G2(s)相串联,在MATLAB中可用串联函数series()来求G1(s)G2(s),其调用格式为num,den=series(num1,den1,num2,den2)其中：第108页/共125页2.并联并联如图所示G1(s)和G2(s)相并联,可由MATLAB的并联函数parallel()来实现,其调用格式为 num,den=parallel(num1,de

35、n1,num2,den2)其中：第109页/共125页3.反馈反馈反 馈 连 接 如 图 所 示。使 用 MATLAB中 的feedback()函数来实现反馈连接,其调用格式为 num,den=feedback(numg,deng,numh,denh,sign)式中：sign为反馈极性,若为正反馈其为1,若为负反馈其为1或缺省。第110页/共125页例如G(s)=,H(s)=,负反馈连接。numg=1,1;deng=1,2;numh=1;denh=1,0;num,den=feedback(numg,deng,numh,denh,1);printsys(num,den)num/den=第111页

36、/共125页MATLAB中 的 函 数 series,parallel和feedback可用来简化多回路方框图。另外,对于单位反馈系统,MATLAB可调用cloop()函数求闭环传递函数,其调用格式为num,den=cloop(num1,den1,sign)第112页/共125页2.6.4 控制系统的零极点模型 传递函数可以是时间常数形式,也可以是零极点形式,零极点形式是分别对原系统传递函数的分子和分母进行因式分解得到的。MATLAB控制系统工具箱提供了零极点模型与时间常数模型之间的转换函数,其调用格式分别为z,p,k=tf2zp(num,den)num,den=zp2tf(z,p,k)其中第

37、一个函数可将传递函数模型转换成零极点表示形式,而第二个函数可将零极点表示方式转换成传递函数模型。第113页/共125页例 如 G(s)=用MATLAB语句表示：num=12241220;den=24622;z,p,k=tf2zp(num,den)z=1.92940.03530.9287i0.03530.9287i第114页/共125页p=0.95671.2272i0.95671.2272i0.04330.6412i0.04330.6412ik=6即变换后的零极点模型为G(s)=第115页/共125页 可以验证MATLAB的转换函数,调用zp2tf()函数将得到原传递函数模型。num,den=z

38、p2tf(z,p,k)num=06.000012.00006.000010.0000den=1.00002.00003.00001.00001.0000即第116页/共125页2.6.5状态空间表达式 状态空间表达式是描述系统特性的又一种数学模型,它由状态方程和输出方程构成,即x(t)=Ax(t)+Bu(t)y(t)=Cx(t)+Du(t)式中x(t)Rn称为状态向量,n为系统阶次;ARnn称为系统矩阵;BRnp称为控制矩阵,p为输入量个数;CRqn称为输出矩阵;DRqp称为连接矩阵,q为输出量个数。第117页/共125页在一般情况下,控制系统的状态空间表达式项简记为(A,B,C,D)。例如：

39、设一个双输入双输出系统的状态空间表达式为第118页/共125页系统模型可由MATLAB命令直观地表示：A=1,2,4;3,2,6;0,1,5B=4,6;2,2;0,2C=0,0,1;0,2,0D=zeros(2,2)MATLAB的控制系统工具箱提供了由状态空间表达式转换成传递函数或由传递函数转换成状态空间表达式的转换函数ss2tf()和tf2ss()。其调用格式为num,den=ss2tf(A,B,C,D,iu)第119页/共125页 反过来,若已知系统的传递函数,求取系统状态空间表达式的调用格式为 A,B,C,D=tf2ss(num,den)例如系统的传递函数为系统的状态空间表达式为num=

40、1,2,3;den=1,3,6,1;A,B,C,D=tf2ss(num,den)第120页/共125页A=-3-6-1100010B=100C=123D=0返回第121页/共125页小 结 本章要求学生熟练掌握系统数学模型的建立和拉氏变换方法。对于线性定常系统，能够列写其微分方程，会求传递函数，会画方框图和信号流图，并掌握方框图的变换及化简方法。1.数学模型是描述元件或系统动态特性的数学表达式，是对系统进行理论分析研究的主要依据。用解析法建立实际系统的数学模型时，分析系统的工作原理，忽略一些次要因素，运用基本物理、化学定律，获得一个既简单又能足够精确地反映系统动态特性的数学模型。第122页/共

41、125页2.实际系统均不同程度地存在非线性，但许多系统在一定条件下可近似为线性系统，故我们尽量对所研究的系统进行线性化处理(如增量化法)，然后用线性理论进行分析。但应注意，不是任何非线性特性均可进行线性化处理。3.传递函数是经典控制理论中的一种重要的数学模型。其定义为：在零初始条件下，系统输出的拉普拉斯与输入的拉普拉斯变换之比。4.根据运动规律和数学模型的共性，任何复杂系统都可划分为几种典型环节的组合，再利用传递函数和图解法能较方便地建立系统的数学模型。第123页/共125页5.方框图是研究控制系统的一种图解模型，它直观形象地表示出系统中信号的传递特性。应用梅逊公式不经任何结构变换，可求出源节点和汇节点之间的传递函数。信号流图的应用更为广泛。6利用MATLAB来进行多项式运算，传递函数零点和极点的计算，闭环传递函数的求取，方框图模型的化简等。返回第124页/共125页第2章 线性系统的数学模型 感谢您的观看！第125页/共125页