现代控制原理精华总结

机电工程系1.状态空间表达式：状态方程和输出方程的总和即称为状态空间表达式。它构成对一个系统动态行为的完整描述。Y:输出向量u:输入向量A:系数矩阵B:控制矩阵（输入矩阵）C:输出矩阵D:直接矩阵状态方程：由系统的状态变量与输入变量之间的关系构成的一阶微分方程组。状态变量：能完全表示系统运动状态的最小个数的一组变量，系统状态变量的数目是惟一的如果矩阵A,B,C,D中的所有元素都是实常数时，则称这样的系统为线性定常（LTI）系统。如果这些元素中有些是时间t的函数，则称系统为线性时变系统。2.由系统微分方程求其状态空间表达式步骤1）假设初始条件为零，将系统微分方程的拉氏变换2）移项变为系统的传递函数3）拼凑分开输入输出变量4）反拉氏变换，写出输入输出微分方程5）用X替换Z，写出X^和Y的含X的微分方程6）将微分方程改写为矩阵形式7）画模拟结构图。传递函数——系统初始松弛（即：初始条件为零）时，输出量的拉氏变换式与输入量的拉氏变换式之比。1.3.3传递函数（矩阵）描述和状态空间描述的比较1）传递函数是系统在初始松弛的假定下输入-输出间的关系描述，非初始松弛系统，不能应用这种描述；状态空间表达式即可以描述初始松弛系统，也可以描述非初始松弛系统。2）传递函数仅适用于线性定常系统；而状态空间表达式可以在定常系统中应用，也可以在时变系统中应用。3）对于数学模型不明的线性定常系统，难以建立状态空间表达式；用实验法获得频率特性，进而可以获得传递函数。综上所示，传递函数（矩阵）和状态空间表达式这两种描述各有所长，在系统分析和设计中都得到广泛应用。机电工程系2.化矩阵A为约当形如果矩阵A有重特征值，并且独立特征向量的个数小于n,这时不能化为对角阵，只能化为约当形。机电工程系例1-7将矩阵化为对角阵3210A解0)2)(1(321det]det[λλλλλAI11λ22λ解出111q212q211121qqQ变换矩阵1112211111QP20012111321011121PAPΛ求线性变换的目的：将系统矩阵变成为标准形，便于求解状态方程。机电工程系系统输出向量对输入向量的脉冲传递函数矩阵DHGICGyu1][)(zz例1-6已知线性定常离散系统方程为)(10)(3.04.010)1(kukkxx)(1011)(kkxy求其脉冲传递函数矩阵解103.04.011011][)(11zzzzHGICGyu)5.0)(8.0()5.0)(8.0(1zzzzzz脉冲传递函数矩阵机电工程系状态空间表达式的模拟结构图buaxxya+u+bxx一阶标量微分方程的结构模拟图1从状态空间表达式画系统模拟框图机电工程系例1-4系统状态方程式为u105610xxx11y求系统传递函数。解：1056111)(11ssssgbAIC22151adj00656111111115656det65sssssssssss机电工程系例1-5线性定常系统状态空间表达式为uxx100100211340010xy100001求系统的传递函数矩阵。解10010021134001100001)(11sssssBAICGyu)4()1(323116123sssssss机电工程系机电工程系机电工程系机电工程系例2-1线性定常系统的齐次状态方程为21213210xxxx求其状态转移矩阵ttAe)(解2211221221112112213)2)(1(1321][11sssssssssssssssAI于是ttAe)(Ltttttttts222211e2ee2e2eeee2][AI机电工程系机电工程系机电工程系机电工程系机电工程系机电工程系机电工程系机电工程系机电工程系机电工程系机电工程系机电工程系例3-16系统方程如下，要求按能控性进行结构分解。u112112xxx10y解211111rankrankranknAbbQC系统不能控111pCQ由于的秩为1。说明中线性独立的列向量只有一列。选择，再补充一个列向量，且与其线性无关，CQ102p110111011121ppPC1CCAPPABPBC1CCPC经过线性变换后u013011CCCCxxxxCCxx11y机电工程系机电工程系机电工程系机电工程系例4-2系统的状态方程如下，判别系统稳定性。)(21221xxxxx解而221121212))(()(xxxxxxxxVx将状态方程代入上式，化简后得)()(2221xxVx222122121)(21)(xxxxVx选取Lyapunov函数，显然是正定的，即满足00)(00)(xxxxVV可见，是负定的，即满足)(xV00)(00)(xxxxVV因此，是一致渐进稳定的。0ex当，有，故系统是一致大范围渐进稳定的。0exx)(xV机电工程系例4-3系统的状态方程为1222221)1(xxxaxxx其中，a为大于零的实数。判别系统的稳定性。解系统的平衡状态为0ex选取Lyapunov函数：2221)(xxVx00)(00)(xxxxVV显然它是正定的，即满足而221122)(xxxxVx将状态方程代入上式，化简后得2222)1(2)(xxaVx21xx1x可见，当和任意的时，有，而和任意时，。又因为，只要变化就不为零，因此在整条状态轨线上不会有。02x1x0)(xV02x1x0)(xV21xx0)(xV因此，是一致渐进稳定的。0ex当，有，故系统是一致大范围渐进稳定的。0exx)(xV机电工程系例4-4系统的状态方程为1221xxkxx其中，k为大于零的实数。分析系统平衡状态的稳定性。解系统的平衡状态为0ex选取Lyapunov函数：2221)(kxxVx00)(00)(xxxxVV显然它是正定的，即满足而02222)(21212211xkxkxxxkxxxVx由定理4-3可知，为Lyapunov意义下一致稳定。0ex机电工程系定理4-4设系统状态方程为)(xfx0ex)(xV)(xV0ex在的某邻域内，标量函数具有连续一阶偏导数，并且满足：1）为正定；2）为正定或半正定；则为不稳定的。)(xV例4-5系统的状态方程为21221xxxxx分析系统平衡状态的稳定性。解系统的平衡状态为0ex选取Lyapunov函数：2221)(xxVx00)(00)(xxxxVV显然它是正定的，即满足而222221212211222222)(xxxxxxxxxxVx由定理4-4可知，是不稳定的。0ex机电工程系例4-6线性定常系统的状态方程为xx1110--判别系统的稳定性。解系统的平衡状态为0ex为简单起见，可以令Q阵为单位矩阵I。IPAPAT1001111011102221121122211211PPPPPPPP解得121212322211211PPPP022211211PPPP011P有可见，P为正定的矩阵，故为大范围一致渐近稳定的。0ex机电工程系例4-7线性定常离散系统的状态方程如下，试判别其稳定性。)(02110)1(kkxx解系统的平衡状态为0ex为简单起见，可以令Q阵为单位矩阵I。IPPGGT解得100102110012102221121122211211PPPPPPPP380035P035P的各阶主子式均大于零，即0380035可见，P为正定的矩阵，故为大范围一致渐近稳定的。0ex机电工程系机电工程系机电工程系机电工程系机电工程系例5-1某位置控制系统（伺服系统）简化线路如下DiiiKu为了实现全状态反馈，电动机轴上安装了测速发电机TG，通过霍尔电流传感器测得电枢电流，即。已知折算到电动机轴上的粘性摩擦系数、转动惯量；电动机电枢回路电阻；电枢回路电感；电动势系数为、电动机转矩系数为。选择、、作为状态变量。将系统极点配置到和，求K阵。TGTGKuDim/(rad/s)N1f2mkg1DJΩ1DRH1.0DLV/(rad/s)1.0eKm/AN1mKoDi31j10机电工程系机电工程系