1、2022-12-1第第八八章章 系统状态变量分析系统状态变量分析 前面的分析方法称为前面的分析方法称为外部法外部法，它强调用，它强调用系统的输系统的输入、输出之间的关系来描述系统的特性入、输出之间的关系来描述系统的特性。其特点：。其特点：（1）只适用于单输入单输出系统，对于多输入多输出）只适用于单输入单输出系统，对于多输入多输出系统，将增加复杂性；系统，将增加复杂性；（2）只研究系统输出与输入的外部特性，而对系统的）只研究系统输出与输入的外部特性，而对系统的内部情况一无所知，也无法控制内部情况一无所知，也无法控制。本章将介绍的本章将介绍的内部法内部法状态变量法状态变量法是用是用n个状态个状态变
2、量的变量的一阶微分或差分方程组（状态方程）一阶微分或差分方程组（状态方程）来描述系来描述系统。优点有：统。优点有：（1）提供系统的内部特性以便研究。）提供系统的内部特性以便研究。（2）便于分析多输入多输出系统；）便于分析多输入多输出系统；（3）一阶方程组便于计算机数值求解。并容易推广用）一阶方程组便于计算机数值求解。并容易推广用于时变系统和非线性系统。于时变系统和非线性系统。2022-12-1第第八八章章 系统状态变量分析系统状态变量分析 8.1 8.1 状态变量与状态方程状态变量与状态方程一、状态变量与状态方程一、状态变量与状态方程二二、状状态方态方程的一般形式程的一般形式8.2 8.2 状
3、态方程的建立状态方程的建立一、电路状态方程的列写一、电路状态方程的列写二、二、由输入由输入-输出方程建立状态方程输出方程建立状态方程 8.3 8.3 离散系统状态方程的建立离散系统状态方程的建立8.4 8.4 连续系统状态方程的解连续系统状态方程的解8.5 8.5 离散系统状态方程的解离散系统状态方程的解点击目录点击目录 ，进入相关章节，进入相关章节2022-12-18.1 8.1 状态变量与状态方程状态变量与状态方程8.1 8.1 状态变量与状态方程状态变量与状态方程一、状态与状态变量的概念一、状态与状态变量的概念从一个电路系统实例引入从一个电路系统实例引入R1R2L1L2iL1iL2iCu
4、Cus1us2au以以u(t)和和iC(t)为输出为输出 若还想了解内部三个若还想了解内部三个变量变量uC(t),iL1(t),iL2(t)的变化情况。的变化情况。这时可列出方程这时可列出方程0dd12LLCiituCa0dd11111SCLLuutiLiR0dd22222CSLLuuiRtiL222222211111112111dd11dd11ddSLCLSLCLLLCuLiLRuLtiuLiLRuLtiiCiCtu2022-12-18.1 8.1 状态变量与状态方程状态变量与状态方程222222211111112111dd11dd11ddSLCLSLCLLLCuLiLRuLtiuLiLRu
5、LtiiCiCtuR1R2L1L2iL1iL2iCuCus1us2au 这是由三个内部变量这是由三个内部变量uC(t)、iL1(t)和和iL2(t)构成的一构成的一阶微分方程组。阶微分方程组。若初始值若初始值uC(t0)、iL1(t0)和和iL2(t0)已知，则根据已知，则根据tt0时时的给定激励的给定激励uS1(t)和和uS2(t)就可惟一地确定在就可惟一地确定在tt0时的解时的解uC(t)、iL1(t)和和iL2(t)。)()()()()()(21222titititutiRtuLLCSL 系统的输出容易地由系统的输出容易地由三个内部变量和激励求三个内部变量和激励求出：出：一组代数方程一组
6、代数方程 2022-12-18.1 8.1 状态变量与状态方程状态变量与状态方程状态与状态变量的定义状态与状态变量的定义 动态系动态系统在某一时刻统在某一时刻t0的的状态状态是指表示该系统是指表示该系统所必所必需最少需最少的一组数值，已知这组数值和的一组数值，已知这组数值和tt0时系统的激时系统的激励，就能完全确定励，就能完全确定tt0时系统的全部工作情况。时系统的全部工作情况。状态变量状态变量是描述状态随时间是描述状态随时间t 变化的一组变量，变化的一组变量，它们在某时刻的值就组成了系统在该时刻的它们在某时刻的值就组成了系统在该时刻的状态状态。对对n阶动态系统需有阶动态系统需有n个独立的状态
7、变量，通常用个独立的状态变量，通常用x1(t)、x2(t)、xn(t)表示。表示。说明说明（1）系统中任何响应均可表示成状态变量及）系统中任何响应均可表示成状态变量及输入的线性组合；输入的线性组合；（2）状态变量应线性独立；）状态变量应线性独立；（3）状态变量的选择并不是唯一的）状态变量的选择并不是唯一的。在初始时刻的值称为在初始时刻的值称为初始状态初始状态。2022-12-18.1 8.1 状态变量与状态方程状态变量与状态方程二、状态方程和输出方程二、状态方程和输出方程在选定状态变量的情况下在选定状态变量的情况下，用状态变量分析系统时，用状态变量分析系统时，一般分一般分两步两步进行：进行：（
8、1）第一步是根据系统的初始状态求出状态变量；）第一步是根据系统的初始状态求出状态变量；（2）第二步是用这些状态变量来确定初始时刻以后的）第二步是用这些状态变量来确定初始时刻以后的系统输出。系统输出。状态变量是通过求解由状态变量构成的一阶微分方状态变量是通过求解由状态变量构成的一阶微分方程组来得到，该一阶微分方程组称为程组来得到，该一阶微分方程组称为状态方程状态方程。状态方程状态方程描述了描述了状态变量的一阶导数与状态变量和状态变量的一阶导数与状态变量和激励激励之间的关系之间的关系。而描述而描述输出输出与状态变量和激励之与状态变量和激励之间关系的一组间关系的一组代数方程代数方程称为称为输出方程输
9、出方程。通常将状态方程和输出方程总称为通常将状态方程和输出方程总称为动态方程动态方程或或系统方程系统方程。2022-12-18.1 8.1 状态变量与状态方程状态变量与状态方程对于一般的对于一般的n阶多输入阶多输入-多多输出输出LTI连续系统，如图连续系统，如图。xi(t0)f1(t)f2(t)fp(t)y1(t)y2(t)yq(t)其状态方程和输出方程为其状态方程和输出方程为 pnpnnnnnnnnppnnppnnfbfbfbxaxaxaxfbfbfbxaxaxaxfbfbfbxaxaxax22112211222212122221212121211112121111pqpqqnqnqqqpp
10、nnppnnfdfdfdxcxcxcyfdfdfdxcxcxcyfdfdfdxcxcxcy221122112222121222212121212111121211112022-12-18.1 8.1 状态变量与状态方程状态变量与状态方程写成矩阵形式：写成矩阵形式：状态方程状态方程)()()(tttBfAxx输出方程输出方程)()()(tttDfCxy其中其中A为为nn方阵，称为方阵，称为系统矩阵系统矩阵，B为为np矩阵，称为矩阵，称为控制矩阵控制矩阵，C为为qn矩阵，称为矩阵，称为输出矩阵输出矩阵，D为为qp矩阵矩阵 对对离散系统离散系统，类似，类似状态方程状态方程)()()1(kkkBfAx
11、x输出方程输出方程)()()(kkkDfCxy状态变量分析的状态变量分析的关键在于状态变量的选取以及状态方程的建立。关键在于状态变量的选取以及状态方程的建立。2022-12-18.2 8.2 连续系统状态方程的建立连续系统状态方程的建立8.2 8.2 连续系统状态方程的建立连续系统状态方程的建立 一、由电路图直接建立状态方程一、由电路图直接建立状态方程 首先选择状态变量首先选择状态变量。通常选通常选电容电压电容电压和和电电感电流感电流为状态变量。为状态变量。必须保证所选状态变必须保证所选状态变量为量为独立的电容电压独立的电容电压和独立的电感电流和独立的电感电流。(a)任选两个电容电压独立(b)
12、任选一个电容电压独立(c)任选两个电感电流独立(d)任选一个电感电流独立uC1uC2uC3uC1uC2usiL1iL2iL3iL2iL1is四种非独立的电路结构四种非独立的电路结构2022-12-18.2 8.2 连续系统状态方程的建立连续系统状态方程的建立状态方程状态方程的建立：的建立：根据电路列出各状态变量的一阶微分方程。根据电路列出各状态变量的一阶微分方程。由于由于tuCiCCddtiLuLLdd为使方程中含有状态变量为使方程中含有状态变量uC的一阶导数的一阶导数，可对接有可对接有该电容的独立结点该电容的独立结点列写列写KCL电流方程；电流方程；为使方程中含有状态变量为使方程中含有状态变
13、量iL的一阶导数的一阶导数，可对含有可对含有该电感的独立回路该电感的独立回路列写列写KVL电压方程。电压方程。对列出的方程，只对列出的方程，只保留状态变量和输入激励保留状态变量和输入激励，设法，设法消消去其它中间的变量去其它中间的变量，经整理即可给出，经整理即可给出标准的状态方程标准的状态方程。对于对于输出方程输出方程，通常可用，通常可用观察法观察法由电路直接列出。由电路直接列出。2022-12-18.2 8.2 连续系统状态方程的建立连续系统状态方程的建立由电路图直接列写状态方程和输出方程的步骤：由电路图直接列写状态方程和输出方程的步骤：（1）选电路中所有）选电路中所有独立的电容电压和电感电
14、流作为独立的电容电压和电感电流作为状态变量状态变量；（2）对）对接有所选电容的独立结点列出接有所选电容的独立结点列出KCL电流方程，电流方程，对对含有所选电感的独立回路列写含有所选电感的独立回路列写KVL电压方程；电压方程；（3）若上一步所列的方程中含有除激励以外的非状）若上一步所列的方程中含有除激励以外的非状态变量，则利用适当的态变量，则利用适当的KCL、KVL方程方程将它们消去将它们消去，然后整理给出然后整理给出标准的状态方程形式标准的状态方程形式；（4）用观察法由电路或前面已推导出的一些关系直）用观察法由电路或前面已推导出的一些关系直接列写输出方程，并整理成标准形式。接列写输出方程，并整
15、理成标准形式。2022-12-18.2 8.2 连续系统状态方程的建立连续系统状态方程的建立例例：电路如图，以电阻电路如图，以电阻R1上的电压上的电压uR1和电阻和电阻R2上的电上的电流流iR2为输出，列写电路的状态方程和输出方程。为输出，列写电路的状态方程和输出方程。uCiLuR1iR2uS1uS2LCR1R2a解解 选状态变量选状态变量x1(t)=iL(t),x2(t)=uC(t)L 1(t)+R1x1(t)+x2(t)=uS1(t)x aC 2(t)+iR2(t)=x1(t)x 消去消去 iR2(t),列右网孔列右网孔KVL方程：方程：R2iR2(t)+uS2(t)-x2(t)=0 代入
16、整理得代入整理得2222()()()sRx tutitR1112121222211()()()()111()()()+()ssRx tx tx tutLLLx tx tx tutCR CR C 111122222110()()()111()()()0ssRx tx tutLLLx tx tutR CCR C整理出矩阵形式整理出矩阵形式：2022-12-1)()(1000)()(100)()(212212121tutuRtxtxRRtitussRR输出方程：输出方程：uR1(t)=R1x1(t)2222()()()sRx tutitR2022-12-18.2 8.2 连续系统状态方程的建立连续系
17、统状态方程的建立二、由输入二、由输入-输出方程建立状态方程输出方程建立状态方程 这里需要解决的问题是这里需要解决的问题是：已知系统的外部描述（已知系统的外部描述（输入输入-输出方程、系统函数、输出方程、系统函数、模拟框图、信号流图模拟框图、信号流图等）；如何写出其状态方程及输等）；如何写出其状态方程及输出方程。出方程。具体方法：具体方法：（1）由系统的）由系统的输入输入-输出方程输出方程或或系统函数系统函数，首先首先画出画出其其信号流图信号流图或或框图框图；（2）选）选一阶子系统一阶子系统(积分器）的积分器）的输出输出作为作为状态变量状态变量；（3）根据每个）根据每个一阶子系统一阶子系统的的输
18、入输出关系输入输出关系列状态方列状态方程；程；（4）在）在系统的输出端系统的输出端列输出方程。列输出方程。2022-12-18.2 8.2 连续系统状态方程的建立连续系统状态方程的建立例例1 某系统的微分方程为某系统的微分方程为 y(t)+3 y (t)+2y(t)=2 f (t)+8 f(t)试求该系统的状态方程和输出方程。试求该系统的状态方程和输出方程。解解由微分方程不难写出其系统函数由微分方程不难写出其系统函数 23)4(2)(2ssssH方法一方法一：画出直接形式的信号流图：画出直接形式的信号流图1-3-228f(t)y(t)设状态变量设状态变量x1(t)、x2(t)x1x2由后一个积
19、分器，有由后一个积分器，有21xx fxxx21232由前一个积分器，有由前一个积分器，有系统输出端，有系统输出端，有 y(t)=8 x1+2 x22022-12-18.2 8.2 连续系统状态方程的建立连续系统状态方程的建立方法二方法二：221423)4(2)(2sssssssH画出串联形式的信号流图画出串联形式的信号流图-1f(t)114y(t)-221设状态变量设状态变量x1(t)、x2(t)x2x1fxx11设中间变量设中间变量 y1(t)y1fxxxy1111341x 2x fxxxyx21212232系统输出端，有系统输出端，有 y(t)=2 x21123012121fxxxx20
20、22-12-18.2 8.2 连续系统状态方程的建立连续系统状态方程的建立方法三方法三：241623)4(2)(2ssssssH画出并联形式的信号流图画出并联形式的信号流图-116-2-41f(t)y(t)设状态变量设状态变量x1(t)、x2(t)x1x21x fxx112x fxx2221120012121fxxxx系统输出端，有系统输出端，有 y(t)=6x1-4 x2可见可见H(s)相同的系统，相同的系统，状态变量的选择并不状态变量的选择并不唯一。唯一。2022-12-18.2 8.2 连续系统状态方程的建立连续系统状态方程的建立例例2 某系统框图如图，状态变量如图标示，试列出其某系统框
21、图如图，状态变量如图标示，试列出其状态方程和输出方程。状态方程和输出方程。f(t)y1(t)y2(t)x2(t)x1(t)x3(t)解解 对三个一阶系统对三个一阶系统211yxx其中，其中，y2=f-x3fxxx311112242xxxxfxx313fxxxx3212232333xxx3233xxx输出方程输出方程 y1(t)=x2y2(t)=-x3+f2022-12-18.3 8.3 离散离散系统状态方程的建立系统状态方程的建立8.3 8.3 离散离散系统状态方程的建立系统状态方程的建立 与连续系统类似，具体方法为：与连续系统类似，具体方法为：（1）由系统的）由系统的输入输入-输出方程输出方
22、程或或系统函数系统函数，首先首先画出其画出其信号流图信号流图或或框图框图；（2）选）选一阶子系统一阶子系统(迟延器）的迟延器）的输出输出作为作为状态变量状态变量；（3）根据每个）根据每个一阶子系统一阶子系统的的输入输出关系输入输出关系列状态方程；列状态方程；（4）在）在系统的输出端系统的输出端列输出方程。列输出方程。2022-12-18.3 8.3 离散离散系统状态方程的建立系统状态方程的建立例例1：某离散系统的差分方程为某离散系统的差分方程为 y(k)+2y(k 1)y(k 2)=f(k 1)f(k 2)列出其动态方程。列出其动态方程。解解：不难写出系统函数不难写出系统函数 212121)(
23、zzzzzH画信号流图：画信号流图：1-21-1y(k)1f(k)设状态变量设状态变量x1(k)，x2(k)：x1x2x1(k+1)=x2(k)：x2(k+1)x2(k+1)=x1(k)2x2(k)+f(k)：输出方程输出方程y(k)=x1(k)+x2(k)2022-12-18.3 8.3 离散离散系统状态方程的建立系统状态方程的建立例例2 某离散系统有两个输入某离散系统有两个输入f1(k)、f2(k)和两个输出和两个输出y1(k)、y2(k)，其信号流图如图示。列写该系统的状态方程和，其信号流图如图示。列写该系统的状态方程和输出方程。输出方程。z-1x1(k)x2(k)x3(k)p1(k)p
24、2(k)-1-2-1111122-13abcdy1(k)y2(k)f1(k)f2(k)解解p1(k)=2x1(k)+2x3(k)p2(k)=3p1(k)-x3(k)+f2(k)=6x1(k)+5x3(k)+f2(k)112:(1)3()()ax kx kx k332132:(1)2()()6()5()()bx kx kp kx kx kf k2212112312:(1)2()()()()7()2()5()()()cx kx kx kp kf kx kx kx kf kf k 2022-12-1112(1)3()()x kx kx k3132(1)6()7()()x kx kx kf k2123
25、12(1)7()2()5()()()x kx kx kx kf kf k)()(101100)()()(706527013)1()1()1(21321321kfkfkxkxkxkxkxkx11()()y kx k213()2()2()y kx kx k)()()(202001)()(32121kxkxkxkyky2022-12-18.4 8.4 连续状态方程的求解连续状态方程的求解8.4 8.4 连续系统状态方程的求解连续系统状态方程的求解 状态方程和输出方程的一般形式为状态方程和输出方程的一般形式为)()()(tttBfAxx)()()(tttDfCxy用拉普拉斯变换法求解状态方程用拉普拉斯
26、变换法求解状态方程 sX(s)-x(0-)=A X(s)+BF(s)(sI-A)X(s)=x(0-)+BF(s)X(s)=(sI-A)-1x(0-)+(sI-A)-1BF(s)=(s)x(0-)+(s)BF(s)式中式中(s)=(sI-A)-1常称为预解矩阵常称为预解矩阵。Y(s)=CX(s)+DF(s)Yx(s)=C(s)x(0-)Yf(s)=C(s)B+D F(s)H(s)=C(s)B+D H(s)称为系统的系统函数矩阵或转移函数矩阵称为系统的系统函数矩阵或转移函数矩阵=C(s)x(0-)+C(s)B+D F(s)2022-12-1所以，所以，(s)的极点就是的极点就是H(s)的极点的极点
27、.即即|sI-A|=0的的根（系根（系统特征根）。统特征根）。1()()()det()adj sIAssIAsIA()det()H()()det()Cadj sIA BDsIAsCs BDsIA1.判判断断|sI-A|=0的根的根是否处于是否处于S平面的左半平面可以平面的左半平面可以判断因果系统是否稳定。判断因果系统是否稳定。2.系系统是否稳定只与状态方程的系统矩阵统是否稳定只与状态方程的系统矩阵A有关。有关。2022-12-18.4 8.4 连续状态方程的求解连续状态方程的求解例例1 描述描述LTI因果系统的状态方程和输出方程为因果系统的状态方程和输出方程为)(10)()(4121)()(2
28、121tftxtxtxtx)(1)()(11)(21tftxtxty解解412141211001)(ssssAI)det()adj()()(1AIAIAIssss1124)3)(2(1ssssX(s)=(s)x(0-)+BF(s)1 10231124)3)(2(1ssss起始状态起始状态x1(0-)=3，x2(0-)=2，输入，输入f(t)=(t)。求状态变量。求状态变量和输出。并判断该系统是否稳定。和输出。并判断该系统是否稳定。263939212)3)(2(3)3)(2()6(3ssssssssss2022-12-18.4 8.4 连续状态方程的求解连续状态方程的求解y(t)=1 1x(t)
29、+f(t)=)(e6e9e9e12)(2332ttttttx)()(e6e9e9e12112332tttttt=(t)+6e-2t(t)由于由于H(s)的极点均在左半平面，故该因果系统稳定。的极点均在左半平面，故该因果系统稳定。H(s)的极点就是的极点就是|sI-A|=0的根。的根。|sI-A|=(s+2)(s+3)2022-12-12022-12-12022-12-12022-12-1()()zizsxtxt2022-12-1X(s)=(sI-A)-1x(0-)+(sI-A)-1BF(s)=(s)x(0-)+(s)BF(s)比较比较 可知可知 1()()AtessIA 记记 为()Atet称
30、称为状态转移矩阵为状态转移矩阵 ()()(0)()*()()y tCt xCt Bf tDf t()0()()*()tA teBfdt Bf t()()zizsytyt2022-12-1()()*()()zsytCt Bf tDf t()()*()Ct BDtf t()*()h tf t比比较：较：Yzs(s)=C(s)B+D F(s)知：知：()()()()()h tCt BDtH sCs BD2022-12-12022-12-12022-12-12022-12-12022-12-1已知（t），求系统矩阵A=？（t）=e，由矩阵指数函数性质1得At（t）0dAdtt2022-12-18.5
31、8.5 离散状态方程的求解离散状态方程的求解8.5 8.5 离散系统状态方程的求解离散系统状态方程的求解 )()()1(kkkBfAxx)()()(kkkDfCxy用用Z变换法求解状态方程变换法求解状态方程 取单边取单边z变换变换:zX(z)-zx(0)=AX(z)+BF(z)Y(z)=CX(z)+DF(z)X(z)=(zI-A)-1zx(0)+(zI-A)-1BF(z)定定义：义：(z)=(zI-A)-1 z，称为预解矩阵。，称为预解矩阵。X(z)=(z)x(0)+z-1(z)BF(z)Y(z)=C(z)x(0)+Cz-1(z)B+DF(z)类似于连续系统，状态转移矩阵(k)(z)2022-
32、12-18.5 8.5 离散状态方程的求解离散状态方程的求解yzi(k)=Z-1C(z)x(0)，yzs(k)=Z-1(Cz-1(z)B+D)F(z)H(z)=Cz-1(z)B+D(z)的极点就是的极点就是H(z)的极点的极点.即即|zI-A|=0的根的根。(1)判判断断特特征根是征根是否在否在z平面的单位圆内可以判断因果系平面的单位圆内可以判断因果系统是否稳定。统是否稳定。(2)系统是否稳定只与系统矩阵系统是否稳定只与系统矩阵A有关。有关。11()(),()()det()adj zIAzIAzzIAzzIA11H()()()()det()det()zCzz BDC zIABDCadj zIA
33、DzIAzIA2022-12-18.5 8.5 离散状态方程的求解离散状态方程的求解例例 已知某离散因果系统的状态方程和输出方程分别为已知某离散因果系统的状态方程和输出方程分别为)(10)()(5610)1()1(2121kfkxkxkxkx)()(1211)()(2121kxkxkyky初始状态为初始状态为 21)0()0(21xx，激励，激励f(k)=(k)。求状态方程的解。求状态方程的解和系统的输出。和系统的输出。解解 zI-A=)3)(2()3)(2(6)3)(2()3)(2(522zzzzzzzzzzzzz(z)=zI-A-1z=zI-A-1=165zz511(2)(3)6zzzz2
34、022-12-18.5 8.5 离散状态方程的求解离散状态方程的求解X(z)=(z)x(0)+z-1BF(z)=11(2)22(1)(3)13(23)1322(1)(3)13zzz zzzzzzzzzzzzz)()3(31 21)3(1 21)(kkkkx212510(2)(3)(2)(3)1216(2)(3)(2)(3)zzzzzzzzzzzzzzzz 2251(2)(3)(2)(3)2161(2)(3)(2)(3)zzzzzzzzzzzzzzz2022-12-18.5 8.5 离散状态方程的求解离散状态方程的求解)()3(31 21)3(1 211211)()(1211)()(2121kkxkxkykykk)()3(1 21)3(21kkk
