第2章控制系统的状态空间描述.ppt2

第二章 控制系统的状态空间描述2.1 2.1 引引 言言2.2 2.2 状态空间模型状态空间模型2.3 2.3 状态空间表达式的建立状态空间表达式的建立2.4 2.4 系统状态方程的线性变换系统状态方程的线性变换2.5 2.5 由状态空间表达式求传递函数阵由状态空间表达式求传递函数阵2.6 2.6 离散时间系统的状态空间表达式离散时间系统的状态空间表达式2.7 2.7 利用利用MATLABMATLAB进行系统数学模型的转换进行系统数学模型的转换小小 结结 2.1 2.1 引引 言言 20 20世纪世纪6060年代年代, ,人们将状态空间的概念引入控人们将状态空间的概念引入控制理论，产生了以状态空间描述为基础制理论，产生了以状态空间描述为基础, ,最优控制最优控制为核心的现代控制理论为核心的现代控制理论系统动态特性的状态空间系统动态特性的状态空间描述由两个数学方程组成，描述由两个数学方程组成，一个是反映系统内部状一个是反映系统内部状态变量和输入变量间因果关系的态变量和输入变量间因果关系的状态方程状态方程; ;另一个另一个是表征系统内部状态变量及输入变量与输出变量转是表征系统内部状态变量及输入变量与输出变量转换关系的换关系的输出方程输出方程。

状态空间法具备如下优点： （1）在数字计算机上求解一阶微分方程组或者差分方程组，比求解与它相当的高阶微分方程或差分方程要容易 （2）状态空间法引入了向量矩阵，大大简化了一阶微分方程组的数学表示法 （3）在控制系统的分析中，系统的初始条件对经典法感到困难的问题，采用状态空间法就迎刃而解了 （4）状态空间法能同时给出系统的全部独立变量的响应，不但反映了系统的输入输出外部特性，而且揭示了系统内部的结构特性，既适用单输入单输出系统又适用多输入多输出系统5）状态空间法可利用计算机进行分析设计以及实时控制，所以可应用求解大量的非线性系统、时变系统、随机过程和采样系统6）利用现代空间法进行系统综合时，是非常有利的 建立动态系统的状态空间模型是状态空间分建立动态系统的状态空间模型是状态空间分析和综合的基本问题和前提析和综合的基本问题和前提, ,本章本章2.32.3节节在介绍状在介绍状态空间分析法基本概念的基础上，讨论动态系统态空间分析法基本概念的基础上，讨论动态系统状态空间表达式建立问题状态空间表达式建立问题; ; 2.42.4节节介绍动态系统介绍动态系统数学模型的等效变换数学模型的等效变换, ,包括状态向量的线性变换包括状态向量的线性变换与状态空间表达式标准型、系统的高阶微分方程与状态空间表达式标准型、系统的高阶微分方程描述化为状态空间描述、系统的传递函数描述化描述化为状态空间描述、系统的传递函数描述化为状态空间描述、由系统状态空间表达式求传递为状态空间描述、由系统状态空间表达式求传递函数阵函数阵; 2.2~2.5; 2.2~2.5节以连续系统为研究对象节以连续系统为研究对象, , 2.62.6节节讨论离散系统的状态空间模型讨论离散系统的状态空间模型; ; 2.72.7节节介绍应介绍应用用MATLABMATLAB进行系统模型变换进行系统模型变换。

1.系统的基本概念 2. 动态系统的两类数学描述 3. 状态的基本概念 2.2 2.2 状态空间模型状态空间模型2.2.12.2.1状态空间的基本概念状态空间的基本概念 1.1.系统的基本概念系统的基本概念 ■■系统系统: :是由相互制约的各个部分有机结合，且具有一定是由相互制约的各个部分有机结合，且具有一定功能的整体功能的整体■■静态系统静态系统: :对于任意时刻对于任意时刻t t，系统的输出惟一地取决于，系统的输出惟一地取决于同一时刻的输入，这类系统称为静态系统静态系统同一时刻的输入，这类系统称为静态系统静态系统亦称为无记忆系统静态系统的输入、输出关系为代亦称为无记忆系统静态系统的输入、输出关系为代数方程 ■■动态系统动态系统: :对任意时刻，系统的输出不仅与对任意时刻，系统的输出不仅与t t时刻的输时刻的输入有关，而且与入有关，而且与t t时刻以前的累积有关时刻以前的累积有关( (这种累积在这种累积在t t0 0(t(t0 0

动态系统的输入、输出关系为微亦称为有记忆系统动态系统的输入、输出关系为微分方程 2. 2. 动态系统的两类数学描述动态系统的两类数学描述(1)(1)外部描述外部描述 外部描述通常称为输入、输出描述，这种外部描述通常称为输入、输出描述，这种描述把系统的输出取为系统外部输入的直接响描述把系统的输出取为系统外部输入的直接响应，显然这种描述回避了表征系统内部的动态应，显然这种描述回避了表征系统内部的动态过程即把系统当成一个过程即把系统当成一个““黑匣”，认为系统，认为系统的内部结构和内部信息全然不知，系统描述直的内部结构和内部信息全然不知，系统描述直接反映了输出变量与输入变量间的动态因果关接反映了输出变量与输入变量间的动态因果关系系  考察图考察图2 2－－1 1所示的所示的n n级级RCRC网络图中虚线框内为网络图中虚线框内为具有放大器隔离的具有放大器隔离的n n级级RCRC电路电路, ,设放大器的输入阻抗设放大器的输入阻抗为无穷大，输出阻抗为零，放大倍数为为无穷大，输出阻抗为零，放大倍数为1 1 图图图图2 2 2 2－－－－1 1 n n级级级级RCRC网络网络网络网络 （2-1） 系统以输入系统以输入u u、输出、输出y y作为变量的外部描述为作为变量的外部描述为式式(2-1)(2-1)所示的高阶线性常系数微分方程所示的高阶线性常系数微分方程, ,即即(2)(2)内部描述内部描述 状态空间描述是内部描述的基本形式，这种描述是基于状态空间描述是内部描述的基本形式，这种描述是基于系统内部结构分析的一类数学模型。

其由两个数学方程组成系统内部结构分析的一类数学模型其由两个数学方程组成: : 一个是反映系统内部状态变量一个是反映系统内部状态变量x1,x2,x1,x2,……,xn ,xn 和输入变量和输入变量u1,u2,u1,u2,……,ur,ur间因果关系的数学表达式，称为状态方程，其数间因果关系的数学表达式，称为状态方程，其数学表达式的形式对于连续时间系统为一阶微分方程组，对于学表达式的形式对于连续时间系统为一阶微分方程组，对于离散时间系统为一阶差分方程组离散时间系统为一阶差分方程组; ; 另一个是表征系统内部状态变量另一个是表征系统内部状态变量x1,x2,x1,x2,……,xn,xn及输入变量及输入变量u1,u2,u1,u2,……,ur,ur与输出变量与输出变量y1,y2,y1,y2,……,ym ,ym 转换关系的数学表达式，转换关系的数学表达式，称为输出方程称为输出方程, ,其数学表达式的形式为代数方程其数学表达式的形式为代数方程 重新考察图重新考察图2-12-1的电网络，利用电路知识容易得到如下一的电网络，利用电路知识容易得到如下一阶微分方程组阶微分方程组 （2-2） 及及 (2-3) 在已知输入在已知输入u u的情况下，解方程式（的情况下，解方程式（2-22-2）、式）、式（（2-32-3））, ,不仅可求出输出响应不仅可求出输出响应y y，而且能得知系统内，而且能得知系统内部电容上电压随时间变化的动态过程信息。

因此部电容上电压随时间变化的动态过程信息因此, ,式式（（2-22-2）、式）、式(2-3)(2-3)是图是图2-12-1所示电网络系统的一种完所示电网络系统的一种完全描述 3.状态的基本概念 (1) (1) 状态状态 状态是完全地描述动态系统运动状况的信息，状态是完全地描述动态系统运动状况的信息，系统在某一时刻的运动状况可以用该时刻系统系统在某一时刻的运动状况可以用该时刻系统运动的一组信息表征，定义系统运动信息的集运动的一组信息表征，定义系统运动信息的集合为状态合为状态2)(2)状态变量状态变量 定义完全表征动态系统时间域运动行为的定义完全表征动态系统时间域运动行为的信息组中的元素为状态变量状态变量组常用信息组中的元素为状态变量状态变量组常用符号符号x x x x1 1 1 1( (t t t t),),x x x x2 2 2 2( (t t t t),),……, ,x x x xn n n n( (t t t t) )表示，且它们相互独表示，且它们相互独立（即变量的数目最小立（即变量的数目最小）） 【例2】确定图确定图2 2－－2 2所示电路的状态变量。

所示电路的状态变量 图图图图2 2 2 2－－－－2 2 2 2 RLCRLCRLCRLC电路电路电路电路 要惟一地确定要惟一地确定要惟一地确定要惟一地确定t t t t时刻电路的运动行为，除了要知道输时刻电路的运动行为，除了要知道输时刻电路的运动行为，除了要知道输时刻电路的运动行为，除了要知道输入电压入电压入电压入电压u u u u( ( ( (t t t t) ) ) )外，还必须给出流过电感上的初始电流外，还必须给出流过电感上的初始电流外，还必须给出流过电感上的初始电流外，还必须给出流过电感上的初始电流i i i i( ( ( (t t t t0 0 0 0) ) ) )和电容上的初始电压和电容上的初始电压和电容上的初始电压和电容上的初始电压u u u uC C C C ( ( ( (t t t t0 0 0 0) ) ) ) ，或者说，或者说，或者说，或者说u u u uC C C C ( ( ( (t t t t) ) ) )和和和和i i i i( ( ( (t t t t) ) ) )这两个这两个这两个这两个变量可用来完全地描述该电路的运动行为，且它们之间变量可用来完全地描述该电路的运动行为，且它们之间变量可用来完全地描述该电路的运动行为，且它们之间变量可用来完全地描述该电路的运动行为，且它们之间是独立的，故是独立的，故是独立的，故是独立的，故u u u uC C C C ( ( ( (t t t t) ) ) )和和和和i i i i( ( ( (t t t t) ) ) )是该电路的状态变量。

是该电路的状态变量是该电路的状态变量是该电路的状态变量 (3)状态向量 设x1(t),x2(t),…,xn(t)是系统的一组状态变量，把这些状态变量看做向量x x(t)的分量，则x x(t)就称为状态向量，记为 (4)状态空间 以x1(t),x2(t),…,xn(t)为坐标轴构成的一个n维欧氏空间，称为状态空间 (5)(5)状态轨迹状态轨迹 状态向量的端点在状态空间中的位置代表了状态向量的端点在状态空间中的位置代表了某一特定时刻系统的状态系统的状态是时间某一特定时刻系统的状态系统的状态是时间t t t t的的函数在不同时刻，系统状态不同，则随着函数在不同时刻，系统状态不同，则随着t t t t的变的变化化, ,状态向量的端点不断移动，其移动的路径就称状态向量的端点不断移动，其移动的路径就称为系统的状态轨迹。

为系统的状态轨迹 ( (6)6)状态方程状态方程 描述系统状态变量间或状态变量与系统输入描述系统状态变量间或状态变量与系统输入变量间关系的一个一阶微分方程组变量间关系的一个一阶微分方程组( (连续系统连续系统) )或或一阶差分方程组（离散系统）一阶差分方程组（离散系统）, ,称为状态方程称为状态方程 【例3】建立图2－2所示RLC电路的状态方程 取电容上的电压uC (t)和电感中的电流i(t)作为状态变量，根据电路原理有 将上式中状态变量的一阶导数放在方程左边，其余项移至方程右边，整理得一阶微分方程组为 上式即为图1所示电路的状态方程,并将其写成向量-矩阵形式，即 (2－4) 式（2-4）可简写为 令，记，， (2－5) 式中, ,(7)状态空间表达式 状态方程和输出方程合起来构成对一个动态系统完整的描述，称为动态系统的状态空间表达式 图2－1所示电路, 若uC (t)为输出,取x1=uC (t),x2=i(t)作为状态变量,则其状态空间表达式为 (2-6) 2.2.2 2.2.2 2.2.2 2.2.2 系统的状态空间表达式的一般形式系统的状态空间表达式的一般形式系统的状态空间表达式的一般形式系统的状态空间表达式的一般形式 一、状态方程：描述系统状态与输入之间关系一、状态方程：描述系统状态与输入之间关系的、一阶微分方程（组）：的、一阶微分方程（组）：、输出方程：描述系统输出与状态、输入之间关、输出方程：描述系统输出与状态、输入之间关系的数学表达式：系的数学表达式：二二 状态空间表达式——描述系统u(t)、X(t)、Y(t)之间关系的状态方程和输出方程总合。

构成了对系统动态行为的完整描述1.1.非线性系统非线性系统 用状态空间表达式描述非线性系统的动态特性，其状态方程是一组一阶非线性微分方程，输出方程是一组非线性代数方程，即 (2-7) 2. 2. 线性系统的状态空间描述线性系统的状态空间描述 若向量方程中 和 的所有组成元都是变量 和 的线性函数，则称相应的系统为线性系统而线性系统的状态空间描述可表示为如下形式： （（2 2－－8 8）） 式中，各个系数矩阵分别为 （2－9） 3 3．．线性时变系统的状态空间描述线性时变系统的状态空间描述 （2－10）一个动态系统的状态向量、输入向量和输出向量自然是时间的函数，而矩阵 、 、 和 的各个元素如果与时间有关，则称这种系统是线性时变系统。

4 4．．线性定常系统的状态空间描述线性定常系统的状态空间描述 矩阵 ， ， 和 的各个元素如果与时间无关，则称这种系统是线性定常系统 式中的各个系数矩阵为常数矩阵  当系统的输出与输入无直接关系（即 ）时，称为惯性系统；相反，系统的输出与输入有直接关系（即 ）时，称为非惯性系统非惯性系统大多数控制系统为惯性系统，所以，它们的动态方程为 （2－11）5 5．．离散系统的状态空间描述离散系统的状态空间描述 当系统的各个变量只在离散的时刻取值时，这种系统称为离散时间系统简称离散系统其状态空间描述只反映离散时刻的变量组之间的因果关系和转换关系是用 来表示离散的时刻，那么离散系统状态空间描述的最一般形式为：（2－12） 对于线性离散时间系统，则上述状态空间描述还可进一步化为如下形式 ：（2－13） 6.6.单输入单输出线性定常连续系统单输入单输出线性定常连续系统 设单输入单输出线性定常n阶连续系统， n个状态变量为x1(t),x2(t),…,xn(t),其状态方程的一般形式为 输出方程的一般形式为 则其向量- 矩阵方程形式的状态空间表达式为 上式简记为 式中 ,为n维状态向量; 称为系统矩阵或状态矩阵; 称为输入矩阵或控制矩阵； 称为输出矩阵或观测矩阵； D是标量，反映输出与输入的直接关联。

7.7.多输入多输出系统多输入多输出系统系统状态空间表达式：系数矩阵中当r>l，m>1时，系统为多输人多输出系统(multi—input and multi output，MIMO)这种系统也称为多变量系统它有r个输入变量和m个输出变量，输入变量u和输出变量y都是向量，为n维状态向量，所以各个矩阵相应的维数为 是 n×n方阵， 是n×r矩阵， 是m×n矩阵，而 是一个m×r矩阵【例】考察图2-10电路，取电压源e为输入变量，R1上的电压为输出变量,建立该电网络的状态空间表达式, 电压和电流为关联参考方向 图2-10四、系统的状态空间描述列写举例网络中只含有电容C、电感L两个独立储能元件，选电容端电压uC、流经电感的电流iL作为状态变量 解 （1）选取状态变量（2）利用电路基本定理列原始方程 回路Ⅰ: （2-14） 回路Ⅱ: (2-15) 代入式(2-14),得 将(2-16) （3）导出状态变量的一阶微分方程组 (2-17) （4）导出状态方程和输出方程 将状态变量的一阶导数看成待定量，用解代数方程方法求解式(2-17)即可求出状态方程。

将式(2-17)写成向量-矩阵形式的方程,即 (2-18) 解之,得向量-矩阵形式的状态方程 (2-19)输出方程为 (2-20) (5) 列写状态空间表达式 将式(2-19)和式(2-20)合起来即为状态空间表达式,若令 则可得状态空间表达式的一般式,即 (2-21) 例例2.2 2.2 系统如图系统如图取状态变量取状态变量: :得：得：系统输出方程为：系统输出方程为：写成矩阵形式的状态空间表达式为：写成矩阵形式的状态空间表达式为：五、状态变量的选取 1．状态变量的选取具有非惟一性2．动态方程或状态空间描述具有非惟一性 3．完全描述一个动态系统所需状态变量的个数由 系统的阶次决定，状态变量必须是相互独立的 4．一般来说，状态变量不一定是有实际物理意义或可以测量的量，但是从工程实际的角度出发，总是选择物理上有意义或可测量的量作为状态变量，如电感中的电流、电容上的电压、电机的转速等列写状态空间表达式的一般步骤： (1)确定系统的状态变量、输入变量、输出变量; (2)根据变量应遵循的物理、化学定理,列出描述系统动态特性或运动规律的微分方程; (3)消去中间变量,得出状态变量的导数与各状态变量、输入变量的关系及输出变量与各状态变量、输入变量的关系; (4)将方程整理成状态方程、输出方程的表准形式。

二 说明： ②系统的状态空间分析法是时域内的一种矩阵运算方法，特别适合于用计算机来计算 ④确定最小的状态变量组以及与之对应的状态空间描述的形式、特点、它们之间的联系与转换等问题，需要进一步分析解决一 步骤： ③状态变量的选择不是唯一的2.2.3 2.2.3 状态空间模型的图示状态空间模型的图示一、结构图一、结构图 线性系统状态空间表达式可用结构图来表示不仅适用于多线性系统状态空间表达式可用结构图来表示不仅适用于多输入多输出系统，当然也适用于单输入单输出系统这种表输入多输出系统，当然也适用于单输入单输出系统这种表示法的实质是把系统分成两部分，如图示法的实质是把系统分成两部分，如图2－－2所示与古典控所示与古典控制理论类似，状态空间表达式也可用图制理论类似，状态空间表达式也可用图2－－4所示的方框结构所示的方框结构图来表示值得注意的是：图中的信号传输线一般是表示列图来表示值得注意的是：图中的信号传输线一般是表示列向量，方框中的字母代表矩阵，每一方框的输入输出关系规向量，方框中的字母代表矩阵，每一方框的输入输出关系规定为：定为： 输出向量输出向量=(方块所示矩阵方块所示矩阵)×(输入向量输入向量) D ∫ A B Cu(t)Y(t)X(t)二、状态变量图二、状态变量图 在状态空间分析中，常以状态变量图来表示系统各变量之间的关系，其来源出自模拟计算机的模拟结构图，这种图为系统提供了一种物理图像，有助于加深对状态空间概念的理解。

所谓状态变量图是由积分器、加法器和放大器构成的图形 绘制步骤绘制步骤：（：（1 1）） 绘制积分器绘制积分器 （（2 2）） 画出加法器和放大器画出加法器和放大器 （（3 3）） 用线连接各元件，并用箭头用线连接各元件，并用箭头 示出信号传递的方向示出信号传递的方向例例 设一阶系统状态方程为设一阶系统状态方程为则其状态图为则其状态图为例例 设三阶系统状态空间表达式为设三阶系统状态空间表达式为则其状态图为则其状态图为2.3 状态空间表达式的建立2.3.1.2.3.1.由物理机理直接建立状态空间表达式由物理机理直接建立状态空间表达式：：例例 系统如图所示系统如图所示选择状态变量：选择状态变量：整理得：整理得：状态方程为：状态方程为：输出方程为：输出方程为： 写成矩阵形式写成矩阵形式例例 系统如图系统如图取状态变量取状态变量: :得：得：系统输出方程为：系统输出方程为：写成矩阵形式的状态空间表达式为：写成矩阵形式的状态空间表达式为：2.3.2 2.3.2 由系统微分方程建立状态空间表达式由系统微分方程建立状态空间表达式 的情形：的情形： 化为能控标准型化为能控标准型取状态变量：取状态变量：则有：则有：即即写成矩阵形式：写成矩阵形式：其中：其中：称为友矩阵称为友矩阵。

能控标准型能控标准型取状态变量：取状态变量：化为能观测标准型化为能观测标准型整理得：整理得：则得能观标准型状态空间表达式：则得能观标准型状态空间表达式：的情形：的情形：计算：计算：定义状态变量定义状态变量: :写成矩阵形式的状态空间表达式写成矩阵形式的状态空间表达式2.3.3. 2.3.3. 由传递函数建立状态空间表达式：由传递函数建立状态空间表达式：(1) (1) 直接分解法直接分解法单输入单输出线性定常系统传递函数单输入单输出线性定常系统传递函数：输出为：输出为：令：令：则有：则有：的的L L氏反变换，氏反变换，则系统的状态空间则系统的状态空间表达式为表达式为：：：：令：令：分别表示分别表示例例 考虑系统考虑系统试写出其能控标准型状态空间表达式试写出其能控标准型状态空间表达式则状态空间表达式为：则状态空间表达式为：选择状态变量：选择状态变量：(2) (2) 并联分解法并联分解法极点两两相异时极点两两相异时其中：其中：令：令：则有：则有：则有：则有：系统的矩阵式表达：系统的矩阵式表达：二、传递函数含重实极点时二、传递函数含重实极点时 设n阶严格有理真分式传递函数为 当传递函数含重实极点时，不失一般性，假设 其中，q重极点 所对应的部分分式系数 按式 计算对于单极点 ，对应的部分分式的系数则按下式计算 选择系统状态变量的拉氏变换为 （2－22）整理 得整理 得（2－23）取拉氏反变换，得输出方程为 式（2－23）取拉氏反变换，得状态方程为 系统的向量-矩阵形式的状态空间表达式为 系统状态变量图 2.4 系统状态方程的线性变换2.4.1 2.4.1 系统状态的线性变换系统状态的线性变换考虑系统：考虑系统：取线性非奇异变换取线性非奇异变换：, ,矩阵矩阵P P 非奇异非奇异整理得：整理得：其中：其中：例例 考虑系统考虑系统取变换：取变换：状态空间表达式变为：状态空间表达式变为：2.4.2 2.4.2 对角标准型对角标准型， 充要条件充要条件：：n n 阶系统矩阵阶系统矩阵 A A 有有n n 个线性无关的个线性无关的特征向量。

特征向量化对角标准型的步骤：化对角标准型的步骤：求取系统矩阵求取系统矩阵 的的 个特征根个特征根和对应的特征向量和对应的特征向量令令 解：解：1)1) 求系统特征根求系统特征根例例 将下系统化为对角标准型将下系统化为对角标准型2)2)求特征向量求特征向量对对由由得得对对由由得得对对由由得得3) 3) 新的状态方程为：新的状态方程为：构成状态转移矩阵构成状态转移矩阵 2.4.3 2.4.3 约当标准型约当标准型重特征根重特征根设矩阵设矩阵具有具有满足满足是是所对应的特征向量若所对应的特征向量若变换化为约当标准型变换化为约当标准型可通过可通过则则称为广义特征向量矩阵称为广义特征向量矩阵线性线性A求约当标准型的步骤：求约当标准型的步骤：求解求解 令令 2.4.42.4.4线性变换的基本性质线性变换的基本性质一、系统特征方程和特征值的不变性一、系统特征方程和特征值的不变性 系统的状态空问表达式为系统的状态空问表达式为 系统的特征方程为系统的特征方程为 上式的根就是系统的特征值而同一系统经线性非奇异变上式的根就是系统的特征值而同一系统经线性非奇异变换之后为换之后为 它的特征方程为 二、传递函数矩阵的不变性二、传递函数矩阵的不变性 传递函数矩阵是系统的输入输出描述，系统状态空间的坐标变换，即内部描述的改变显然不会影响传递函数矩阵。

2.5 2.5 由状态空间表达式求传递函数阵由状态空间表达式求传递函数阵2.5.1 SISO2.5.1 SISO系统系统取取L L氏变换得：氏变换得：A A的特征值即为系统的极点的特征值即为系统的极点结论：结论： （1）系统矩阵A A的特征多项式等同于传递函数的分母多项式 （2）传递函数的极点就是系统矩阵A A的特征值 （3）多项式与之和即为传递函数的分子多项式 （4）由于状态变量选择的不同，同一系统的状态空间描述不是惟一的，但从表征系统状态空间描述的不同的A A，B B，C C和D变换到表征系统输入输出描述的传均函数G(s)是惟一的选择不同的状态变量可得至两个不同的动态方程式，所求得的传递函数却是相同的，这称为传递函数的不变性2.5.2 MIMO2.5.2 MIMO系统系统其中：其中：2.5.2 2.5.2 组合系统的传递函数阵组合系统的传递函数阵一、一、 并联并联: :系统如图，二子系统并联连接系统如图，二子系统并联连接特点：特点：二、二、 反馈反馈系统如图，二子系统并联连接系统如图，二子系统并联连接特点特点：：(1) (1) 动态反馈动态反馈传递矩阵：传递矩阵：(2) (2) 静态反馈静态反馈闭环系统状态空间描述为：闭环系统状态空间描述为：闭环系统传递矩阵为：闭环系统传递矩阵为：2.6 离散时间系统状态空间表达式离散时间系统差分方程表示：离散时间系统差分方程表示：其对应脉冲传函为：其对应脉冲传函为：定义定义：：取：取：对其进行对其进行Z 反变换得：反变换得： 写成矩阵形式：写成矩阵形式：例例2.6.1 2.6.1 考虑离散系统考虑离散系统试写出其状态空间表达式。

试写出其状态空间表达式得状态空间表达式为：得状态空间表达式为：解：取解：取Z变换得变换得2. 7 MATLAB性系统动态分析中的应用2. 7. 1 2. 7. 1 应用应用MATLABMATLAB计算线性定常系统的矩阵计算线性定常系统的矩阵 指数指数( (状态转移矩阵状态转移矩阵) )2. 7. 2 2. 7. 2 应用应用MATLAB MATLAB 求定常系统时间响应求定常系统时间响应 2. 7.3 2. 7.3 应用应用MATLAB MATLAB 变连续状态空间模型为离变连续状态空间模型为离 散状态空间模型散状态空间模型 2. 7. 1 应用MATLAB计算矩阵指数 1. 应用MATLAB 符号数学工具箱求矩阵指数闭合解析式 基于2.3节矩阵指数的拉普拉斯变换求解法，可调用MATLAB 符号数学工具箱（Symbolic Math Toolbox）中的符号运算函数先算出“预解矩阵” ，再对“预解矩阵”进行拉普拉斯反变换即求得 。

另外，MATLAB 符号数学工具箱中有专用于计算矩阵指数的指令expm（）可调用 【例】 已知 ,应用MATLAB求 %MATLAB Program 2_1a syms s t %定义基本符号变量s 和tA=[4,0,0;0,3,1;0,1,3];FS=inv(s*eye(3)-A); %求预解矩阵 eAt=ilaplace(FS,s,t); %求 eAt=simplify(eAt) %化简 的表达式 解 MATLAB Program 2_1a给出了基于拉普拉斯变换求 的MATLAB 程序 2.应用数值矩阵的指数运算函数expm（）求 对应于 （ 为某一常数）的值 MATLAB Program 2_2给出了调用expm（）求例 中矩阵A的矩阵指数 对应于 的值 的MATLAB 程序 %MATLAB Program 2_2 A=[4,0,0;0,3,1;0,1,3]; T=0.1;eAT=expm(A*T) 3.应用MATLAB 符号数学工具箱求离散系统状态转移矩阵解析式 【例2-14】 已知离散系统状态方程为 应用MATLAB求其状态转移矩阵 的解析式 解 例2-10中已采用四种方法求出了系统的 ， MATLAB Program 2_3给出了基于Z变换求 的MATLAB 程序。

%MATLAB Program 2_3 syms z k %定义基本符号变量z和k G=[0,1;-0.2,-0.9]; Fz=(inv(z*eye(2)-G))*z; %求 Fk=iztrans(Fz,z,k) %调用Z反变换指令求 Fk=simple(Fk) %将符号运算结果表达式转换为最简形式 与例2-10求解结果一致，MATLAB Program 2_3 程序运行结果如下： Fk = [ 5*(-2/5)^k-4*(-1/2)^k, 10*(-2/5)^k-10*(-1/2)^k] [ -2*(-2/5)^k+2*(-1/2)^k, -4*(-2/5)^k+5*(-1/2)^k] 2. 7. 2 应用MATLAB 求定常系统时间响应 1.1.状态方程的数值解状态方程的数值解 常微分方程数值解一般使用逐步积分的方法实现，Runge–Kutta法是应用最多的一种微分方程数值解法MATLAB提供的ode23（）、ode45（）是分别采用2/3阶、4/5阶Runge–Kutta法的常微分方程数值求解的函数，一般ode45（）较ode23（）运算速度快，两者调用格式相同，即 其中,xfun为由m函数定义的一阶微分方程组的m函数名，该m函数必须以状态向量x的一阶导数为输出。

若原方程为高阶微分方程，，应通过第1章的“实现”方法将其转换为一阶微分方程组，即状态空间表达式； t0和tf分别为积分的起始和终止时间，单位为秒；x0 0为状态向量的初始值；t和x均为返回值，其中t为离散时间列向量；x为解向量构成的矩阵，其第j列为第j个状态变量与t相对应的解向量，，j=1,2,…n 例 已知系统状态空间表达式为 设x(0)=0，，试求u(t)=1(t)为单位阶跃函数时系统时间响应的数值解 解 MATLAB Program 2_4a建立了描述系统状态方程的m函数ode_example.m %MATLAB Program 2_4a %ode_example.m function sx=ode_example(t,x) % sx为状态列向量x的导数的导数 sx(1,1)=-10*x(1)-35*x(2)-50*x(3)-24*x(4)+1; %sx应按状态方程编写 sx(2,1)=x(1); % sx是与x同维的列向量 sx(3,1)=x(2); sx(4,1)=x(3); 将MATLAB Program 2_4a保存为名为ode_example.m的m文件，且将保存ode_example.m的路径设置成当前路径。

MATLAB Program 2_4b为调用求解函数求状态方程和输出响应数值解的程序，图2-4所示为状态方程和输出响应数值解曲线 %MATLAB Program 2_4b x0=[0;0;0;0]; %设置初值条件t0=0;tf=6;tspan=[t0,tf] %设置积分起始和终止时间[t,x]=ode45('ode_example',tspan,x0);%调用求解函数求状态方程数值解y=24*x(:,4); %据输出方程求输出响应的数值解subplot(1,2,1) plot(t,x(:,1),' k',t,x(:,2),'-.r',t,x(:,3),':b',t,x(:,4),'-k') %绘状态方程数值解曲线gtext('x1')gtext('x2') gtext('x3')gtext(‘x4’)subplot(1,2,2)plot(t,y,‘k’) %绘输出响应数值解曲线gtext('y') 图2-4状态方程和输出响应数值解曲线 2 2．状态方程的解析解．状态方程的解析解 MATLAB Symbolic Math Toolbox提供的dsolve（）为求常微分方程解析解的指令，其调用格式为 S=dsolve（‘eqn1’，‘eqn2’，…） 其中，‘eqn1’,‘eqn2’,…为输入参数，其为描述常微分方程、初始条件及独立变量的字符表达式。

微分方程是必不可少的输入参数，多个方程或初始条件可在一个输入参数内联立输入，且以逗号分隔；若独立变量默认,则小写字母t为独立变量；若要定义其它独立变量，则由全部输入参数‘eqn1’,‘eqn2’,…中的最后一个参数定义  在输入参数中，描述常微分方程规定用字符D代表对独立变量的导数（因此，用户所定义的字符变量不应含有字符D），例如若t为独立变量，y为t的函数，则Dy代表dy/dt，D2y代表 ，D3y代表 ,…；初始条件可采用形如‘y(a)=b’或’Dy(a)=b’的字符（串）表达式给出S为返回的存放符号微分方程解的构架数组 3.3.基于状态空间模型的时域响应分析基于状态空间模型的时域响应分析 MATLAB Control System Toolbox 提供了连续系统单位阶跃响应计算函数step( )、单位脉冲响应计算函数impulse()、零输入响应计算函数initial()、任意输入（包括系统初始状态）响应计算函数lsim()，与此对应，dstep( ) 、 dimpulse()、dinitial()、dlsim()分别为计算离散系统单位阶跃响应、单位脉冲响应、零输入响应、任意输入（包括系统初始状态）响应的函数。

 例如，若给定线性定常连续系统、离散系统分别如式（2-83）、式（2-84）所示，则 执行step(A，，B，，C，，D)指令，可得一组单位阶跃响应曲线，每条曲线对应于式（2-83）所示连续系统的输入/输出组合即在某一输入端单独施加单位阶跃信号作用下的某一输出响应，时间向量t的范围自动设定；； 执行step(A，，B，，C，，D，，t)指令与执行step(A，，B，，C，，D)指令一样，可得一组单位阶跃响应曲线，但时间向量t是由用户设定的；； 执行step(A，，B，，C，，D，iu)指令，可得式（2-83）所示连续系统从第iu个输入到所有输出的单位阶跃响应曲线；  执行[y,x,t]=step(A，，B，，C，，D，iu)指令，可得式（2-83）所示连续系统从第iu个输入到所有输出y及状态x的单位阶跃响应数据，且返回函数自动设定的时间向量t，但不绘制响应曲线； 执行dinitial (G，，H，，C，，D，，x0)指令可得式（2-84）所示离散系统每一个输出的零输入响应曲线，取样点数由函数自动设定；  执行lsim（A，，B，，C，，D，，u，，t，，x0）指令可针对系统初始状态x0和输入u绘制系统所有输出（全）响应曲线，其中t为用户设定的线性等间距的时间向量；对多输入系统,u为数值矩阵，其列数等于输入信号数，第j个输入信号对应于t的离散序列构成u的第j列，行数等于时间向量t的维数。

…… 【例】设双输入双输出系统状态空间表达式为 且设 、 ，系统初始状态为零 1)分别求 、 单独作用下系统的输出响应；2)求 和 共同作用下系统的输出响应 解 1）MATLAB Program 2_6a为调用step（）函数求 、 单独作用下系统输出响应曲线的程序，图2-5为程序运行结果 %MATLAB Program 2_6a A=[-1,-1;25,-2];B=[1,1;0,2];C=[1,0;0,1];D=[0,0;0,0];step(A,B,C,D)grid 图2-5 、 单独作用下系统输出响应 2）MATLAB Program 2_6b为求 和 共同作用下系统输出响应的MATLAB程序，图2-6为程序运行结果 %MATLAB Program 2_6b A=[-1,-1;25,-2];B=[1,1;0,2];C=[1,0;0,1];D=[0,0;0,0];t=0:0.01:4; %生成时间向量t LT=length(t); %求时间向量t的维数（长度） u1=ones(1,LT); u2=ones(1,LT); %生成单位阶跃信号对应于向量t的离散 序列，u1和u2均为与向量t同维的向量 u=[u1;u2]'; % u1和u2的转置分别构成u的第1和第2列 lsim(A,B,C,D,u,t) grid 图2-6 和 共同作用下系统输出响应 2. 7. 3 应用MATLAB 变连续状态空间模型为离散状态空间模型 MATLAB Control System Toolbox提供的c2d（）函数可简化线性定常连续状态方程离散化系数矩阵的求解，若要将式（2-83）所示的线性定常连续系统变换为（2-84）所示的离散系统，且设控制输入端采用零阶保持器，T为采样周期，其调用格式为 [G，H]=c2d(A,B,T) 。