数学模型的相互转换

数学模型的相互转换第1页，课件共30页，创作于2023年2月2.1连续系统数学模型

2.1.1表达形式描述控制系统的主要模型有微分方程、状态空间表达式等形式的时域描述法和用传递函数描述的频域描述法。即对于一个连续的控制系统，数字仿真常用的数学模型一般有3种表示方式：◆直接用微分方程描述；◆用传递函数描述；

多项式形式零极点形式◆状态方程描述；这三种描述方式是可以相互转换的。

第2页，课件共30页，创作于2023年2月(1).微分方程设连续系统的输出量为y(t)，输入量为u(t)，采用微分方程的形式来表示的系统数学模型一般式可描述如下：（2.1）上式中，为常数。(2).传递函数第3页，课件共30页，创作于2023年2月对式（2.1）等号两边逐项进行拉氏变换，并考虑系统输出、输入及其各阶导数的初值均为零，可得到

（2.2）式中，-系统输出的拉氏变换；

-系统输入的拉氏变换；可得系统的传递函数为：(2.3)第4页，课件共30页，创作于2023年2月微分方程或传递函数是用系统的输入、输出之间的关系来描述系统的，表示了系统的外部特征，所以称其为外部模型。用微分方程表示的系统可以是非线性或线性系统，而对于传递函数表示的系统，只适用于单输入-单输出的线性定常系统，所以传递函数的模型表示有一定的局限性。

(3).状态空间表达式状态空间表达式可以由两个途径获得，由微分方程或系统结构方框图导出，这里对微分方程推导作简单说明。设系统由不含输入量导数项的n阶微分方程表示：(2.4)第5页，课件共30页，创作于2023年2月定义n个状态变量为，且令写出各个状态变量的一阶微分方程形式第6页，课件共30页，创作于2023年2月将上述n个一阶微分方程写成矩阵向量形式为(2.5)上式称为状态空间表达式，其中A、B、C为系数矩阵，x为状态变量。第7页，课件共30页，创作于2023年2月2.1.2.数学模型的相互转换由于要解决的控制问题所需的数学模型与所给定的已知数学模型往往是不一致的，不同的应用场合需要对控制系统的数学模型进行转换

（1）微分方程转化为传递函数和状态空间表达式例2.1已知某控制系统的微分方程为将其分别表示为传递函数、一阶微分方程组和状态空间描述。第8页，课件共30页，创作于2023年2月解：①将给定系统微分方程的两端取拉氏变换，并令初始值为零，则可用以下传递函数表示根据传递函数定义有第9页，课件共30页，创作于2023年2月③按照状态空间描述，将各变量和系数矩阵表达为②由于是二阶导数，可以引入两个状态变量，将给定的二阶微分方程写成一阶微分方程组形式第10页，课件共30页，创作于2023年2月（2）传递函数转换成状态空间表达式转换采用的方法是状态变量图法，用基本模拟单元替代系统的传递函数得到的图形式系统结构图，在系统结构图上标上状态变量的图形是状态变量图，然后再求出状态空间表达式。

图2-1积分器的系统结构图和状态变量图由状态变量图根据积分器的输入、输出关系写出：输出方程状态方程对于初始条件为零的积分器

第11页，课件共30页，创作于2023年2月对于带反馈的积分器，其传递函数为图2-2带反馈积分器的状态变量图由积分器输入、输出关系得到第12页，课件共30页，创作于2023年2月从上面得到由系统结构图到状态变量图并到处状态空间表达式的步骤如下：根据系统的传递函数，画出系统结构图，n阶系统有n个积分器；把积分器输出处定为状态变量x，积分器输入处为状态变量微分，并把状态变量x，和状态变量微分分别标在积分器输入和输出处，得到状态变量图；根据积分器输入、输出的方程写出系统的状态方程和输出方程。对于高阶、复杂系统采用级联法、并联法和串联法得到代表实际系统传递函数的系统结构图及相应的状态变量图，依据同样方法求得状态空间表达式。第13页，课件共30页，创作于2023年2月对于一个三阶系统的传递函数，

(2.6)（1）级联法

采用如下步骤进行：◆用传递函数的最高阶次除以传递函数分子分母多项式得到◆利用信号流图法画出该系统的信号流图及系统状态变量图

第14页，课件共30页，创作于2023年2月图2-3系统信号流图图2-4系统状态变量图第15页，课件共30页，创作于2023年2月根据积分器输入、输出关系得到如下方程◆写成矩阵表达式

A、B、C为系数矩阵第16页，课件共30页，创作于2023年2月其中，分析系数矩阵A、B、C可见：系数矩阵A是一个方阵，以I表示行号，J表示列号，最末一行元素和传递函数分母多项式系数按s0升幂排列的负值一一对应，其余各行的元素在J=I+1时为1，其他全部为0；系数矩阵B是一个单列矩阵，最后一行元素为1，其余为零；系数矩阵C是一个单行矩阵，各列元素与传递函数分子多项式系数按s0升幂排列值相同。第17页，课件共30页，创作于2023年2月推广到n阶方程，系数矩阵

A、B、C分别为

这种形式的矩阵称为可控标准型状态表达式。(2.7)第18页，课件共30页，创作于2023年2月例2.2已知某控制系统的传递函数为利用式（2.7）将系统模型转换为状态空间描述和一阶微分方程组描述。

解：①这是一个三阶系统，将给定的系统传递函数按照状态空间描述的系数矩阵A、B、C的对应关系，依据式（2.7）可得各系数矩阵如下第19页，课件共30页，创作于2023年2月组合为状态空间描述有②将上述矩阵展开即可得到系统模型的一阶微分方程组表示形式第20页，课件共30页，创作于2023年2月（Ⅱ）并联法并联法的思路是把高阶系统的传递函数转变成若干个一阶环节传递函数之和，如将式（2-6）表示成下式形式然后，对各个一阶环节的传递函数画出系统结构图，并标上状态变量后，得到如图2-5所示的状态变量图。图2.5并联法系统状态变量图第21页，课件共30页，创作于2023年2月写成矩阵表达式由图2-5可到状态空间表达式为第22页，课件共30页，创作于2023年2月其中利用级联法和并联法得到的状态空间表达式的系数矩阵A、B、C是不同的。可以用特征方程判断特征值是否相同，确定几个状态空间表达式是否属于同一个外部模型，因为同一系统有相同的特征值。

第23页，课件共30页，创作于2023年2月（3）状态空间表达式转换成传递函数已知系统状态空间表达式①对上式两边取拉氏变换②在零初始条件下整理上式并在等式两边乘以单位矩阵（2.9）

（2.10）第24页，课件共30页，创作于2023年2月③将X(s)代入式（2.10）得到

可得（2.11）

在系数矩阵已知的情况下，根据式（2.11）就可以求出状态空间表达式所对应的外部模型的传递函数。

在MATLAB控制系统工具箱中提供了大量的控制系统模型相互转换的函数，如表2-1所示。第25页，课件共30页，创作于2023年2月表2-1数学模型转换函数及其功能函数名函数功能ss2tf将系统状态空间模型转换为传递函数模型ss2zp将系统状态空间模型转换为零极点增益模型tf2ss将系统传递函数模型转换为状态空间模型tf2zp将系统传递函数模型转换为零极点增益模型zp2ss将系统零极点增益模型换为状态空间模型zp2tf零极点增益模型换为传递函数模型第26页，课件共30页，创作于2023年2月[z,p,k]=tf2zp(num,den)[A,B,C,D]=tf2ss(num,den)[num,den]=zp2tf(z,p,k)[A,B,C,D]=zp2ss(z,p,k)[num,den]=ss2tf(A,B,C,D)[z,p,k]=ss2zp(A,B,C,D)第27页，课件共30页，创作于2023年2月对线性定常系统，式中s的系数均为常数，且a1不等于零，这时系统在MATLAB中可以方便地由分子和分母系数构成的两个向量唯一地确定出来，这两个向量分别用num和den表示。

num=[b1,b2,…,bm,bm+1] den=[a1,a2,…,an,an+1]

注意：它们都是按s的降幂进行排列的。MATLAB表示数学模型举例连续系统的传递函数模型第28页，课件共30页，创作于2023年2月在MATLAB中零极点增益模型用[z,p,K]矢量组表示。即：z=[z1,z2,…,zm]p=[p1,p2,...,pn]K=[k]函数tf2zp()可以用来求传递函数的零极点和增益。零极点增益模型K为系统增益，zi为零点，pj为极点零极点模型实际上是传递函数模型的另一种表现形式，其原理是分别对原系统传递函数的分子、分母进行分解因式处理，以获得系统的零点和极点的表示形式。第29页，课件共30页，创作于2023年2月用法举例：1）已知系统状态空间模型为：》A=[01;-1-2];B=[0;1];》C=[