广东海洋大学：《自动控制原理》课程教学资源（教学内容）第八章 根轨迹法 §8.1 引言

第8章控制系统的状态空间分析与综合 第1~7章涉及的内容属于经典控制理论的范畴，系统的数学模型是线性定常微分方程和 传递函数，主要的分析与综合方法是时域法、根轨迹法和频域法。经典控制理论通常用于单 输入一单输出线性定常系统，其缺点是只能反映输入一输出间的外部特性，难以揭示系统内 部的结构和运行状态，不能有效处理多输入一多输出系统、非线性系统、时变系统等复杂系 统的控制问题。 随着科学技术的发展，对控制系统速度、精度、适应能力的要求越来越高，经典控制理 论已不能满足要求。1960年前后，在航天技术和计算机技术的推动下，现代控制理论开始 发展，一个重要的标志就是美国学者卡尔曼引入了状态空间的概念。它是以系统内部状态为 基础进行分析与综合的控制理论，两个重要的内容如下。 ()最优控制：在给定的限制条件和评价函数下，寻求使系统性能指标最优的控制规 律。 (2)最优估计与滤波：在有随机干扰的情况下，根据测量数据对系统的状态进行最优 估计。 本章讨论控制系统的状态空间分析与综合，它是现代控制理论的基础。 8.1控制系统的状态空间描述 81.1系统数学描述的两种基本方法 控制 被控过程 控制输入 观测y 控制器 执行器 被控对象 +传感器 反馈控制 图8-1典型控制系统方块图 典型控制系统如图&1所示，由被控对象、传感器、执行器和控制器组成。被控过程 327

327 第 8 章 控制系统的状态空间分析与综合 第 1~7 章涉及的内容属于经典控制理论的范畴，系统的数学模型是线性定常微分方程和 传递函数，主要的分析与综合方法是时域法、根轨迹法和频域法。经典控制理论通常用于单 输入－单输出线性定常系统，其缺点是只能反映输入－输出间的外部特性，难以揭示系统内 部的结构和运行状态，不能有效处理多输入－多输出系统、非线性系统、时变系统等复杂系 统的控制问题。 随着科学技术的发展，对控制系统速度、精度、适应能力的要求越来越高，经典控制理 论已不能满足要求。1960 年前后，在航天技术和计算机技术的推动下，现代控制理论开始 发展，一个重要的标志就是美国学者卡尔曼引入了状态空间的概念。它是以系统内部状态为 基础进行分析与综合的控制理论，两个重要的内容如下。 （1）最优控制：在给定的限制条件和评价函数下，寻求使系统性能指标最优的控制规 律。 （2）最优估计与滤波：在有随机干扰的情况下，根据测量数据对系统的状态进行最优 估计。 本章讨论控制系统的状态空间分析与综合，它是现代控制理论的基础。 8.1 控制系统的状态空间描述 8.1.1 系统数学描述的两种基本方法 图 8-1 典型控制系统方块图 典型控制系统如图 8-1 所示，由被控对象、传感器、执行器和控制器组成。被控过程 控制 u 控制器 执行器 被控对象 传感器 控制输入 观测 y 被控过程 x 反馈控制

(见图8-2)具有若干输入端和输出端。数学描述通常有两种基本方法：一种是输入、输出 描述（外部描述），它将系统看成为“黑箱”，只是反映输入与输出间的关系，而不去表征系 统的内部结构和内部变量，如传递函数：另一种是状态空间描述（内部描述），它是基于系 统内部结构的一种数学模型，由两个方程组成。一个反映系统内部变量x和输入变量间的 关系，具有一阶微分方程组或一阶差分方程组的形式：另一个是表征系统输出向量y与内部 变量及输入变量间的关系，具有代数方程的形 式。外部描述虽能反映系统的外部特性，却不 能反映系统内部的结构与运行过程，内部结构 不同的两个系统也可能具有相同的外部特性， X1,2,…x。 因此外部描述通常是不完整的：内部描述则能 y 全面完整地反映出系统的动力学特征。 图8-2被控过程 8.1.2状态空间描述常用的基本概念 1输入和输出 由外部施加到系统上的激励称为输入，若输入是按需要人为施加的，又称为控制：系统 的被控量或从外部测量到的系统信息称为输出，若输出是由传感器测量得到的，又称为观测。 2.状态、状态变量和状态向量 能完整描述和惟一确定系统时域行为或运行过程的一组独立（数目最小）的变量称为系 统的状态，其中的各个变量称为状态变量。当状态表示成以各状态变量为分量组成的向量时， 称为状态向量。系统的状态x0由1=1时的初始状态x(。)及1之的输入4)惟一确 定。 对n阶微分方程描述的系统，当n个初始条件x(。,(t。),,x-(t)及1≥。的输入 )给定时，可惟一确定方程的解，故x,元，，m-这n个独立变量可选作状态变量。状 态对于确定系统的行为既是必要的，也是充分的。阶系统状态变量所含独立变量的个数为 m,当变量个数小于时，便不能完全确定系统的状态，而当变量个数大于n时，则存在多 余的变量，这些多余的变量就不是独立变量。判断变量是否独立的基本方法是看它们之间是 否存在代数约束。 状态变量的选取并不惟一，一个系统通常有多种不同的选取方法。但应尽量选取能测 328

328 （见图 8-2）具有若干输入端和输出端。数学描述通常有两种基本方法：一种是输入、输出 描述（外部描述），它将系统看成为“黑箱”，只是反映输入与输出间的关系，而不去表征系 统的内部结构和内部变量，如传递函数；另一种是状态空间描述（内部描述），它是基于系 统内部结构的一种数学模型，由两个方程组成。一个反映系统内部变量 x 和输入变量 u 间的 关系，具有一阶微分方程组或一阶差分方程组的形式；另一个是表征系统输出向量 y 与内部 变量及输入变量间的关系，具有代数方程的形 式。外部描述虽能反映系统的外部特性，却不 能反映系统内部的结构与运行过程，内部结构 不同的两个系统也可能具有相同的外部特性， 因此外部描述通常是不完整的；内部描述则能 全面完整地反映出系统的动力学特征。 8.1.2 状态空间描述常用的基本概念 1.输入和输出 由外部施加到系统上的激励称为输入，若输入是按需要人为施加的，又称为控制；系统 的被控量或从外部测量到的系统信息称为输出，若输出是由传感器测量得到的，又称为观测。 2.状态、状态变量和状态向量 能完整描述和惟一确定系统时域行为或运行过程的一组独立（数目最小）的变量称为系 统的状态，其中的各个变量称为状态变量。当状态表示成以各状态变量为分量组成的向量时， 称为状态向量。系统的状态 x(t) 由 0 t = t 时的初始状态 x ( 0 t ) 及 0 t  t 的输入 u(t) 惟一确 定。 对 n 阶微分方程描述的系统，当 n 个初始条件 ( ), ( ), , ( ) 0 ( 1) 0 0 x t x t x t  n−  及 0 t  t 的输入 u(t) 给定时，可惟一确定方程的解，故 ( 1) , , , n− x x   x 这 n 个独立变量可选作状态变量。状 态对于确定系统的行为既是必要的，也是充分的。n 阶系统状态变量所含独立变量的个数为 n，当变量个数小于 n 时，便不能完全确定系统的状态，而当变量个数大于 n 时，则存在多 余的变量，这些多余的变量就不是独立变量。判断变量是否独立的基本方法是看它们之间是 否存在代数约束。 状态变量的选取并不惟一，一个系统通常有多种不同的选取方法。但应尽量选取能测 n x , x , x 1 2 p u u u  2 1 q y y y  2 1 图 8-2 被控过程

量的物理量或独立贮能元件的贮能变量作为状态变量，以便实现系统设计。在机械系统中， 常选取位移和速度作为变量：在RL-C网络中，常选电感电流和电容电压作为状态变量： 在由传递函数绘制的方块图中，常取积分器的输出作为状态变量。 3,状态空间 以状态向量的n个分量作为坐标轴所组成的n维空间称为状态空间 4状态轨迹 系统在某个时刻的状态，可以看作是状态空间的一个点。随着时间的推移，系统状态不 断变化，便在状态空间中描绘出一条轨迹，该轨迹称为状态轨迹。 5.状态方程 描述系统状态变量与输入变量之间关系的一阶向量微分方程或差分方程称为系统的状 态方程，它不含输入的微积分项。状态方程表征了系统由输入所引起的状态变化，一般情况 下，状态方程既是非线性的，又是时变的，它可以表示为 ()fx(t).u(t).t (8-1) 6输出方程 描述系统输出变量与系统状态变量和输入变量之间函数关系的代数方程称为输出方程， 当输出由传感器得到时，又称为观测方程。输出方程的一般形式为 y0)=gt),), (8-2) 输出方程表征了系统状态和输入的变化所引起的系统输出变化。 7动态方程 状态方程与输出方程的组合称为动态方程，又称为状态空间表达式，其一般形式为 )=fx0,0, (8.3a) )=gx(t),t), 或离散形式 x(t)=f().u). (8-3b) ）=gx)t), 8线性系统： 线性系统的状态方程是一阶向量线性微分方程或差分方程，输出方程是向量代数方程。 线性连续时间系统动态方程的一般形式为 39

329 量的物理量或独立贮能元件的贮能变量作为状态变量，以便实现系统设计。在机械系统中， 常选取位移和速度作为变量；在 R-L-C 网络中，常选电感电流和电容电压作为状态变量； 在由传递函数绘制的方块图中，常取积分器的输出作为状态变量。 3.状态空间 以状态向量的 n 个分量作为坐标轴所组成的 n 维空间称为状态空间。 4.状态轨迹 系统在某个时刻的状态，可以看作是状态空间的一个点。随着时间的推移，系统状态不 断变化，便在状态空间中描绘出一条轨迹，该轨迹称为状态轨迹。 5.状态方程 描述系统状态变量与输入变量之间关系的一阶向量微分方程或差分方程称为系统的状 态方程，它不含输入的微积分项。状态方程表征了系统由输入所引起的状态变化，一般情况 下，状态方程既是非线性的，又是时变的，它可以表示为 x (t) = fx(t),u(t),t （8-1） 6.输出方程 描述系统输出变量与系统状态变量和输入变量之间函数关系的代数方程称为输出方程， 当输出由传感器得到时，又称为观测方程。输出方程的一般形式为 y(t) = gx(t),u(t),t （8-2） 输出方程表征了系统状态和输入的变化所引起的系统输出变化。 7.动态方程 状态方程与输出方程的组合称为动态方程，又称为状态空间表达式，其一般形式为       = = y t g x t u t t x t f x t u t t ( ) ( ), ( ), ( ) ( ), ( ), （8-3a） 或离散形式        = + = k k k k k k k k y t g x t u t t x t f x t u t t ( ) ( ), ( ), ( ) ( ), ( ), 1 （8-3b） 8.线性系统： 线性系统的状态方程是一阶向量线性微分方程或差分方程，输出方程是向量代数方程。 线性连续时间系统动态方程的一般形式为

()A(t)x(t)+B(t)u(t) (8-4) y(t)C(t)x(t)+D(t)u(t) 设状态x、输入、输出y的维数分别为n,P,g,称n×n矩阵A()为系统矩阵或状态矩 阵，称n×p矩阵B()为控制矩阵或输入矩阵，称q×n矩阵C()为输出矩阵或观测矩阵， 称g×D矩阵D)为前馈矩阵或输入输出矩阵。 9.线性定常系统 线性系统的A,B,C,D中的各元素全部是常数。即 (t)=Ax(t)+Bu(t) (8-5a) y(t)=Cx(t)+Du(t) 对应的离散形式为 x(k+1)=Gx(k)+Hu(k y(k)=Cx(k)+Du(k) (8-5b) 「x 47 「y yJ aia2…a1n A= …a B= b2b22 「dd C= D= d …Cg d …dp」 为书写方便，常把系统(&-5a)和系统(8-5b)分别简记为S(4,B,C,D)和S(G,HCD) 10.线性系统的结构图 线性系统的动态方程常用结构图表示。图83为连续系统的结构图：图84为离散系统 的结构图。图中，I为n×)单位矩阵，s是拉普拉斯算子，：为单位延时算子

330    = + = + y(t) C(t)x(t) D(t)u(t) x(t) A(t)x(t) B(t)u(t) （8-4） 设状态 x、输入 u、输出 y 的维数分别为 n, p,q ，称 nn 矩阵 A(t)为系统矩阵或状态矩 阵，称 n p 矩阵 Bt() 为控制矩阵或输入矩阵，称 q  n 矩阵 C(t)为输出矩阵或观测矩阵， 称 q  p 矩阵 D(t)为前馈矩阵或输入输出矩阵。 9.线性定常系统 线性系统的 A，B，C，D 中的各元素全部是常数。即    = + = + y (t) Cx(t) Du(t) x(t) Ax(t) Bu(t) （8-5a） 对应的离散形式为 ( ) ( ) ( ) ( 1) ( ) ( ) y k Cx k Du k x k Gx k Hu k = + + = + （8-5b）             = n x x x x  2 1               = p u u u u  2 1               = q y y y y  2 1             = n n nn n n a a a a a a a a a A       1 2 21 22 2 11 12 1               = n n np p p b b b b b b b b b B       1 2 21 22 2 11 12 1               = q q qn n n c c c c c c c c c C       1 2 21 22 2 11 12 1 11 12 1 21 22 2 1 2 p p q q qp d d d d d d D d d d       =         为书写方便，常把系统(8-5a)和系统(8-5b)分别简记为 S(A,B,C,D)和 S(G,H,C,D)。 10.线性系统的结构图 线性系统的动态方程常用结构图表示。图 8-3 为连续系统的结构图；图 8-4 为离散系统 的结构图。 图中，I 为( nn )单位矩阵，s 是拉普拉斯算子，z 为单位延时算子

D 图8-3线性连续时间系统结构图 图84线性离散时间系统结构图 由于状态变量的选取不是惟一的，因此状态方程、输出方程、动态方程也都不是惟一的。 但是，用独立变量所描述的系统的维数应该是惟一的，与状态变量的选取方法无关。 动态方程对于系统的描述是充分的和完整的，即系统中的任何一个变量均可用状态方程 和输出方程来描述。状态方程着眼于系统动态演变过程的描述，反映状态变量间的微积分约 束：而输出方程则反映系统中变量之间的静态关系，者眼于建立系统中输出变量与状态变量 间的代数约束，这也是非独立变量不能作为状态变量的原因之一。动态方程描述的优点是便 于采用向量、矩阵记号简化数学描述，便于在计算机上求解，便于考虑初始条件，便于了解 系统内部状态的变化特征，便于应用现代设计方法实现最优控制和最优估计，适用于时变、 非线性、连续、离散、随机、多变量等各类控制系统。 Ka) (b) 图85电路的独立变量 例8-1试确定图85中(a、(b)所示电路的独立状态变量。图中4、1分别是输入电 压和输入电流，y为输出电压，x,-1,2,3,为电容器电压或电感器电流。 解并非所有电路中的电容器电压和电感器电流都是独立变量。对图&5(a)所示电路， 不失一般性，假定电容器初始电压值均为0，有 C2 6+ 331

331 图 8-3 线性连续时间系统结构图 图 8-4 线性离散时间系统结构图 由于状态变量的选取不是惟一的，因此状态方程、输出方程、动态方程也都不是惟一的。 但是，用独立变量所描述的系统的维数应该是惟一的，与状态变量的选取方法无关。 动态方程对于系统的描述是充分的和完整的，即系统中的任何一个变量均可用状态方程 和输出方程来描述。状态方程着眼于系统动态演变过程的描述，反映状态变量间的微积分约 束；而输出方程则反映系统中变量之间的静态关系，着眼于建立系统中输出变量与状态变量 间的代数约束，这也是非独立变量不能作为状态变量的原因之一。动态方程描述的优点是便 于采用向量、矩阵记号简化数学描述，便于在计算机上求解，便于考虑初始条件，便于了解 系统内部状态的变化特征，便于应用现代设计方法实现最优控制和最优估计，适用于时变、 非线性、连续、离散、随机、多变量等各类控制系统。 （a） （b） 图 8-5 电路的独立变量 例 8-1 试确定图 8-5 中（a）、（b）所示电路的独立状态变量。图中 u、i 分别是输入电 压和输入电流，y 为输出电压， i x ，i=1，2，3，为电容器电压或电感器电流。 解 并非所有电路中的电容器电压和电感器电流都是独立变量。对图 8-5（a）所示电路， 不失一般性，假定电容器初始电压值均为 0，有 1 2 3 3 2 x c c c x + = 1 2 3 2 3 x c c c x + =

因此，只有一个变量是独立的，状态变量只能选其中一个，即用其中的任意一个变量作 为状态变量便可以确定该电路的行为。实际上，三个串并联的电容可以等效为一个电容。 对图8-5(b)所示电路，1=，因此两者相关，电路只有两个变量是独立的，即(x 和)或和，可以任用其中一组变量如(，)作为状态变量。 8.13系统的传递函数矩阵 设初始条件为零，对线性定常系统的动态方程进行拉氏变换，可以得到 X(s)=(sI-A)BU(s) (8-6) Y(s)=[C(sI-A)B+DWU(s) 系统的传递函数矩阵（简称传递矩阵）定义为 G(s)=C(sI-A)B+D (8-7) 例8-2已知系统动态方程为 [-86] 试求系统的传递函数矩阵。 "w4-a-0c-0-0 Gs)=C(1-A0B=

332 因此，只有一个变量是独立的，状态变量只能选其中一个，即用其中的任意一个变量作 为状态变量便可以确定该电路的行为。实际上，三个串并联的电容可以等效为一个电容。 对图 8-5（b）所示电路，x1 = x2，因此两者相关，电路只有两个变量是独立的，即（x1 和 x3）或(x2 和 x3)，可以任用其中一组变量如（x2，x3）作为状态变量。 8.1.3 系统的传递函数矩阵 设初始条件为零，对线性定常系统的动态方程进行拉氏变换，可以得到 1 1 ( ) ( ) ( ) ( ) [ ( ) ] ( ) X s sI A BU s Y s C sI A B D U s − − = − = − + (8-6) 系统的传递函数矩阵（简称传递矩阵）定义为 G s = C sI − A B + D −1 ( ) ( ) （8-7） 例 8-2 已知系统动态方程为             =                  +            − =      2 1 2 1 2 1 2 1 2 1 0 1 1 0 0 1 1 0 0 2 0 1 x x y y u u x x x x   试求系统的传递函数矩阵。 解 已知 0 1 1 0 1 0 , , , 0 0 2 0 1 0 1 A B C D       = = = =             − 故             + + =       + − − = − − 2 1 0 ( 2) 1 1 0 2 1 ( ) 1 1 s s s s s s sI A G(s) = C             + + =                  + +       − = − 2 1 0 ( 2) 1 1 0 1 1 0 2 1 0 ( 2) 1 1 0 1 1 0 ( ) 1 s s s s s s s s sI A B

8.1.4线性定常系统动态方程的建立 1.根据系统物理模型建立动态方程 例8-3试列写如图8-6所示的R-LC电路方程，选择几组状态变量并建立相应的动 态方程，并就所选状态变量间的关系进行讨论。 e 0 0 图8-6R-LC电路 解有明确物理意义的常用变量主要有：电流、电阻器电压、电容器的电压与电荷、电 感器的电压与磁通。 根据回路电压定律 8+会+= 电路输出量为 》设状变量为电器电流和电容器电压。即=1名=己山，则状态方程为 名=--+ 输出方程为 y=X2 其向量矩阵形式为

333 8.1.4 线性定常系统动态方程的建立 1.根据系统物理模型建立动态方程 例 8-3 试列写如图 8-6 所示的 R-L-C 电路方程，选择几组状态变量并建立相应的动 态方程，并就所选状态变量间的关系进行讨论。 图 8-6 R-L-C 电路 解 有明确物理意义的常用变量主要有：电流、电阻器电压、电容器的电压与电荷、电 感器的电压与磁通。 根据回路电压定律 idt e dt C di Ri + L + =  1 电路输出量为 1 y u C = =  idt 1) 设状态变量为电感器电流和电容器电压，即 x = i 1 ，  = idt C x 1 2 ，则状态方程为 e L x L x L R x 1 1  1 = − 1 − 2 + 2 1 1 x C x = 输出方程为 2 y = x 其向量-矩阵形式为              =         +                − − =      2 1 2 1 2 1 0 1 0 1 0 1 1 x x y L e x x C L L R x x  

简记为 ∫=Ax+b 式中， 2)设状态变量为电容器电流和电荷，即，=，x=∫1d山，则有 卧中 》设状态变量名=己仙+名=己仙，其中，高无明确意义的物理，可 以推出 高-名+定低-+ y=X 其向量矩阵形式为 可见对同一系统，状态变量的选择不具有惟一性， 图87双输入一三输出机械位移系统 动态方程也不是惟一的。 例8-4由质量块、弹簧、阻尼器组成的双输入三输出机械位移系统如图87所示，具 有力F和阻尼器气缸速度V两种外作用，输出量为质量块的位移，速度和加速度。试列写

334 图 8-7 双输入－三输出机械位移系统 简记为    y cx x Ax be =  = + 式中， , 0 1 0 1 , 0 1 1 , , 2 1 2 1 =         =           − − =       =       = b L c C C L R A x x x x x x    2) 设状态变量为电容器电流和电荷，即 1 2 x i x idt = = ,  ，则有 1 1 1 2 2 2 1 1 1 , 0 1 0 0 R x x x L LC L e y x x x C           − −       = + =                         3) 设状态变量 1 2 1 1 x idt Ri x idt , C C = + =   ，其中， 1 x 无明确意义的物理量，可 以推出 ( ) ( ) 1 1 2 1 2 1 x e L R x x dt RC di x = x + R = − + − + 2 2 1 2 ( ) 1 1 y x x x RC i C x =  = = − 其向量-矩阵形式为         =         +                − − − =      2 1 2 1 2 1 0 1 0 1 1 1 1 x x y L e R x x RC RC L RC R RC x x   可见对同一系统，状态变量的选择不具有惟一性， 动态方程也不是惟一的。 例 8-4 由质量块、弹簧、阻尼器组成的双输入三输出机械位移系统如图 8-7 所示，具 有力 F 和阻尼器气缸速度 V 两种外作用，输出量为质量块的位移，速度和加速度。试列写

该系统的动态方程。m,k,∫分别为质量、弹簧刚度、阻尼系数：x为质量块位移。 解根据牛顿力学可知，系统所受外力F与惯性力m求、阻尼力元一门和弹簧恢复力kx 构成平衡关系，系统微分方程如下 m成+f(-+kx=F 这是一个二阶系统，若已知质量块的初始位移和初始速度，系统在输入作用下的解便可 惟一确定，故选择质量块的位移和速度作为状态变量。设x=x,X=产。由题意知系统有 三个输出量，设 y=X1=x,当=元=x2,y3=求。 于是由系统微分方程可以导出系统状态方程 元= 名==上c-n-+月 其向量矩阵形式为 「101「001 m mm 当5 例日-5对于图8-8所示的机械系统，若不考虑重 m. 力对系统的作用，试列写该系统以拉力F为输入，以质 量块m和m的位移片和之为输出的动态方程。 中 解根据牛顿定律，系统微分方程为 m=ky2-)+f2(2-y)-ky-f m2j3=F-k,(y2-)-f(2-j) 图8-8双质量块机械系统 式中，k,k为弹簧刚度，不，方为阻尼系数 该系统有4个独立的储能元件，即弹簧k,k和质量m,m,故应选择其中4个相互独 立的变量作为系统的状态变量，现选择x=片，为2=2，X=头，x=2,经过整 理，可得到系统的动态方程

335 图 8-8 双质量块机械系统 该系统的动态方程。 m,k, f 分别为质量、弹簧刚度、阻尼系数；x 为质量块位移。 解 根据牛顿力学可知，系统所受外力 F 与惯性力 m  x  、阻尼力 f( x  －V)和弹簧恢复力 kx 构成平衡关系，系统微分方程如下: m  x  + f (x  −V) + kx = F 这是一个二阶系统，若已知质量块的初始位移和初始速度，系统在输入作用下的解便可 惟一确定，故选择质量块的位移和速度作为状态变量。设 x = x x = x  1 2 , 。由题意知系统有 三个输出量，设 y = x = x y = x  = x y =  x  1 1 2 2 3 , , 。 于是由系统微分方程可以导出系统状态方程  f x V k x F m x x x x = = − − − + = 2 2 1 1 2 ( ) 1    其向量-矩阵形式为               +              − − =      V F m f m x x m f m k x x 1 0 1 0 0 2 1 2 1   1 1 2 2 3 1 0 0 0 0 1 0 0 1 y x F y k f f x V y m m m m                 = +                       − −         例 8-5 对于图 8-8 所示的机械系统，若不考虑重 力对系统的作用，试列写该系统以拉力Ｆ为输入，以质 量块 m1 和 m2 的位移 y1 和 y2 为输出的动态方程。 解 根据牛顿定律，系统微分方程为 ( ) ( ) ( ) ( ) 2 2 2 2 1 2 2 1 1 1 1 2 1 2 2 1 1 1 1 1 m y F k y y f y y m y k y y f y y k y f y        = − − − − = − + − − − 式中， 1 2 k ,k 为弹簧刚度， 1 2 f , f 为阻尼系数。 该系统有 4 个独立的储能元件，即弹簧 1 2 k ,k 和质量 1 2 m ,m ，故应选择其中 4 个相互独 立的变量作为系统的状态变量，现选择 1 1 2 2 3 1 4 2 x = y , x = y , x = y  , x = y  ，经过整 理，可得到系统的动态方程

0 0 1 0 7 「0 叫上 1五 0 + 01 m m、 「]00x ,0100x3 L. 2.由高阶微分方程建立动态方程 (1)微分方程不含输入量的导数项 ym+a+a++ay+doy=Bou (8-8) 选n个状态变量为x=，x2=立，…，x。=y-,有 元=为 元2=x (8-9) 元n1=x 元=-a-a2--an1+Bu y=x 得到动态方程 =Ax+bu (8-10) y=cx 式中 「0 1 0 … 01 0 0 1 0 A= b c=0.0] 0 0 0 0 L-ao -%-4… -4- B. 36

336                   =                     +                             − − + − + = −             4 3 2 1 2 1 2 4 3 2 1 2 2 2 2 2 2 2 2 1 2 1 1 2 1 2 1 1 2 4 3 2 1 0 1 0 0 1 0 0 0 1 0 0 0 0 0 0 1 0 0 1 0 x x x x y y F m x x x x m f m f m k m k m f m f f m k m k k x x x x     2.由高阶微分方程建立动态方程 （1) 微分方程不含输入量的导数项 y a y a y a y a y u n n n n n 1 0 0 ( 2) 2 ( 1) 1 ( ) + + + + + =  − − − −   (8-8) 选 n 个状态变量为 ( 1) 1 2 , , , − = = = n n x y x y   x y ，有 1 2 2 3 1 0 1 1 2 1 0 1 n n n n x x x x x x x a x a x a u y x − − =   =    =  = − − − − +   =   （8-9） 得到动态方程 x Ax bu y cx = +   =  （8-10） 式中 1 2 1 0 1 2 1 0 0 1 0 0 0 0 0 1 0 0 , , 0 0 0 1 0 n n n x x x A b x x a a a a  − −                   = = =                               − − − − , c = 1 0  0