第2章数学模型的相互转换

《第2章数学模型的相互转换》由会员分享，可在线阅读，更多相关《第2章数学模型的相互转换（33页珍藏版）》请在文档大全上搜索。

1、第第2章章 连续系统数字仿真的连续系统数字仿真的 基本算法基本算法2.1 连续系统数学模型连续系统数学模型2.2 数值积分算法数值积分算法 2.3 数值积分算法的基本分析数值积分算法的基本分析 2.4 连续系统仿真的离散相似算法连续系统仿真的离散相似算法2.5 常用快速数字仿真算法常用快速数字仿真算法2.6 实时数字仿真算法实时数字仿真算法 小结小结2.1连续系统数学模型连续系统数学模型2.1 .1表达形式表达形式描述控制系统的主要模型有微分方程、状态空间表达式等形式的时域描述法和用传递函数描述的频域描述法。即对于一个连续的控制系统，数字仿真常用的数学模型一般有3种表示方式：直接用微分方程描述

2、；用传递函数描述； 多项式形式 零极点形式状态方程描述；这三种描述方式是可以相互转换的。(1). 微分方程微分方程设连续系统的输出量为y(t)，输入量为u(t) ，采用微分方程的形式来表示的系统数学模型一般式可描述如下：121211201.nnnnnnnnmmmmmd ydydydyaaaa ydtdtdtdtd ud uccc udtdt（2.1）上式中，12101,nnma aaa c ccLL为常数。(2). 传递函数传递函数对式（2.1）等号两边逐项进行拉氏变换，并考虑系统输出、输入及其各阶导数的初值均为零，可得到121212012( )( )( )( )( )( )( )( )nnn

3、nmmmns Y sa sY sa sY sa Y sc s U sc sU sc sU sc U sLL（2.2）( )Y s( )U s式中，-系统输出的拉氏变换； -系统输入的拉氏变换；可得系统的传递函数为：120121212( )( )( )mmmnnnnnc sc sc scY sG sU ssa sa saLL(2.3)微分方程或传递函数是用系统的输入、输出之间的关系来描述系统的，表示了系统的外部特征，所以称其为外部模型。用微分方程表示的系统可以是非线性或线性系统，而对于传递函数表示的系统，只适用于单输入-单输出的线性定常系统，所以传递函数的模型表示有一定的局限性。(3). 状态空

4、间表达式状态空间表达式状态空间表达式可以由两个途径获得，由微分方程或系统结构方框图导出，这里对微分方程推导作简单说明。设系统由不含输入量导数项的n阶微分方程表示：1212112nnnnnnnnd ydydydyaaaa yudtdtdtdtL(2.4)定义n个状态变量为12,nx xxL，且令1212322111nnnnxydyxxdtd yxxdtdyxxdt&M&写出各个状态变量的一阶微分方程形式1223341121121nnnnnnnxxxxxxxxxa xa xaxa xu &M&L将上述n个一阶微分方程写成矩阵向量形式为& xAxBuyCx(2.5)上式称为状态空间表达式，其中122

5、10100000010000001001nnnAaaaaa121nnBBBBB121nnxxxxxA、B、C为系数矩阵，x为状态变量。2.1.2. 数学模型的相互转换数学模型的相互转换由于要解决的控制问题所需的数学模型与所给定的已知数学模型往往是不一致的，不同的应用场合需要对控制系统的数学模型进行转换（1）微分方程转化为传递函数和状态空间表达式）微分方程转化为传递函数和状态空间表达式例2.1 已知某控制系统的微分方程为222.56210d ydyduyudtdtdt将其分别表示为传递函数、一阶微分方程组和状态空间描述。解：解：将给定系统微分方程的两端取拉氏变换，并令初始值为零，则可用以下传递函

6、数表示2( )2.5( )6 ( )2( ) 10 ( )s Y ssY sY ssU sU s2(2.56) ( )(210) ( )ssY ssU s根据传递函数定义有2( )210( )( )2.56Y ssG sU sss122121262.5102xxxxxuyxx &按照状态空间描述，将各变量和系数矩阵表达为112201062.51xxuxx &12102xyx由于是二阶导数，可以引入两个状态变量，将给定的二阶微分方程写成一阶微分方程组形式（2）传递函数转换成状态空间表达式）传递函数转换成状态空间表达式转换采用的方法是状态变量图法 ，用基本模拟单元替代系统的传递函数得到的图形式系统

7、结构图，在系统结构图上标上状态变量的图形是状态变量图，然后再求出状态空间表达式 图2-1 积分器的系统结构图和状态变量图xu&yx由状态变量图根据积分器的输入、输出关系写出：输出方程状态方程对于初始条件为零的积分器对于带反馈的积分器，其传递函数为1( )G ssa图2-2 带反馈积分器的状态变量图由积分器输入、输出关系得到xaxuyx &从上面得到由系统结构图到状态变量图并到处状态空间表达式的步骤如下：根据系统的传递函数，画出系统结构图，n阶系统有n个积分器；把积分器输出处定为状态变量x，积分器输入处为状态变量微分 ，并把状态变量x，和状态变量 微分分别标在积分器输入和输出处，得到状态变量图；

8、根据积分器输入、输出的方程写出系统的状态方程和输出方程。x &x &对于高阶、复杂系统采用级联法、并联法和串联法得到代表实际系统传递函数的系统结构图及相应的状态变量图，依据同样方法求得状态空间表达式。对于一个三阶系统的传递函数，23232( )712ssG ssss(2.6)（1）级联法）级联法采用如下步骤进行： 用传递函数的最高阶次除以传递函数分子分母多项式得到232132( )7121sssG sss利用信号流图法画出该系统的信号流图及系统状态变量图图2-3 系统信号流图图2-4 系统状态变量图232132( )7121sssG sss根据积分器输入、输出关系得到如下方程122332312

9、7xxxxxxxu &12323yxxx写成矩阵表达式yxAxBu= Cx&A、B、C为系数矩阵0100010127A1001B23 1C分析系数矩阵A、B、C可见：可见：系数矩阵A是一个方阵，以是一个方阵，以I表示行号，J表示列号，最末一行元素和传递函数分母多项式系数按s0升幂排列的负值一一对应，其余各行的元素在J=I+1时为1，其他全部为0；系数矩阵B是一个单列矩阵，最后一行元素为是一个单列矩阵，最后一行元素为1，其余为零；系数矩阵C是一个单行矩阵，各列元素与传递函数分子多项是一个单行矩阵，各列元素与传递函数分子多项式系数按式系数按s0升幂排列值相同。推广到n阶方程，系数矩阵 A、B、C分