《状态变量分析导论》学习总结

导言

状态变量分析导论是于2021年10月购买的一本二手书籍，作者是P.F-布莱克曼。虽是一本很老的教材。但内容生动，形象，是一本学习状态空间法的不错教材。因此借隔离的机会把这部分的内容加深下理解。2021-12-25（圣诞节。）

第一章-状态变量分析原理

1. 状态变量分析法是通过向量–向量的变化关系从而描述完整的系统运行形式。该变化关系采用矩阵的形式描述。
2. 矩阵描述系统运行方程的合理性在于“向量”的广义性。向量可以是一组数，可以包含导数，可以是面积，可以是二阶导数，具体的含义依据实际的应用场景确定。
3. 对于一个三阶常系数微分方程而言：
   $$\begin{gathered} \frac{d^{3}y}{d{t}^3}+a_2\frac{d^{2}y}{d{t}^2}+a_1\frac{dy}{dt}+a_{0}y=u \\ \text{令：}{x}_1=y;x_2=\dot{y};x_3=\ddot{y} \\ \text{则：} \\ \left[\begin{array}{c}{\dot{x}}_1\\ {\dot{x}}_2\\ {\dot{x}}_3\end{array}\right]=\left[\begin{array}{ccc}0& 1& 0\\ 0& 0& 1\\ -a_0& -a_1& -a_2\end{array}\right]\left[\begin{array}{c}{x}_1\\ x_2\\ x_3\end{array}\right]+\left[\begin{array}{c}0\\ 0\\ 1\end{array}\right]u \end{gathered}$$
   对于主对角线以上为1的系统矩阵，表示原系统为串联形式。该形式称为同伴型矩阵；
4. 状态变量可以是系统中不存在的量，以保证系统内部运行关系的完善，但实际应用中通常选择为可测的，有明确物理意义的变量。
5. 任何系统都可以用一组一阶微分方程描述，或者叫向量微分方程。描述向量之间的变化关系。$$\{\begin{array}{l}\dot{x}=Ax+Bu\\ y=Cx+Du\end{array}$$
   矩阵A的主对角元素表述为自反馈项。主对角之外的元素表述该变量微分项与其他变量的关系。表述为耦合项。
6. 采用状态变量分析法的一个优点：状态与输入的关系不变，但系统描述可以多变，从而可以选择简单的系统描述形式对原始系统描述简化。
7. 系统矩阵A描述了固有运动形式，在输入u的作用下，系统运动轨迹发生变化。受到控制。
   8.状态变化的相轨迹可以这样理解：以不同状态作为坐标轴，如x1,x2—xn，当dx1/dt为正时，状态向量的分量x1朝正方向变化，当导数越大时，变化的速度越快。
8. **状态空间形式的稳态计算。**对于双状态变量，单输入，当系统稳定后，状态向量导数置0，状态将由输入决定，由于输入只有一个。因此相轨迹终点在一条直线上。表现为状态的控制程度。
9. 状态向量微分方程的变换：对于一阶微分方程为例，零输入响应相当于给状态导数施加一个冲击输入，该状态导数的积分即状态初值。因此两者是相似的。**对于一般的输入响应，可以等效为无限多的冲击信号的叠加。**因此可采用卷积积分的形式获得状态的解。$$x\left(t\right)={\int }_0^t{e}^{a\left(t-\tau \right)}bu\left(\tau \right)d\tau$$
10. 向量微分方程拉普拉斯变化法的求解：
    $$\begin{gathered} \mathbf{x}\left(s\right)={\left[\mathbf{s}\mathbf{I}-\mathbf{A}\right]}^{-1}\left[\mathbf{x}\left(0\right)+\mathbf{B}\mathbf{u}\left(s\right)\right] \\ \mathbf{x}\left(s\right)=\frac{adj{\left[\mathbf{s}\mathbf{I}-\mathbf{A}\right]}^{-1}}{\det \left[\mathbf{S}\mathbf{I}-\mathbf{A}\right]}\left[\mathbf{x}\left(0\right)+\mathbf{B}\mathbf{u}\left(s\right)\right] \end{gathered}$$
    det[SI-A]为系统矩阵的特征方程。系统矩阵的特征值表现为某些s导致系统矩阵特征方程为0 的值，由于特征方程出现在每一个状态变量的初值响应和输入响应，因此每一个输入和初始条件都将激发所有的状态变量出现该模态。
11. 状态变量传递函数的零点，也就是不同模态的大小由系统矩阵的伴随矩阵和B阵决定。
12. 转移矩阵的定义：系统从初始条件向量到状态向量变化的轨迹。或者叫状态的解。
13. 一阶微分方程的解：
    $$x\left(t\right)=e^{at}x\left(0\right)+{\int }_0^t{e}^{a\left(t-\tau \right)}bu\left(\tau \right)d\tau$$
14. 矩阵指数就是状态转移矩阵。描述了状态运动的轨迹。状态的解矩阵。
15. 对于输入的状态响应，需要计算状态转移矩阵的积分。
16. 由于状态转移矩阵如此重要，因此由大量的学者研究状态转移矩阵的计算和性质。
    状态转移矩阵（矩阵指数）的定义：
    $${\mathcal{L}}^{-1}{\left[\mathbf{s}\mathbf{I}-\mathbf{A}\right]}^{-1}=e^{\mathbf{A}t}$$
17. 关于矩阵指数计算的一个重要定理：凯莱哈密顿定理。
    每个矩阵满足自己的特征方程。即若矩阵A由特征方程：
    $$\det \left[\mathbf{s}\mathbf{I}-\mathbf{A}\right]=s^n+a_{n-1}{s}^{n-1}+...+a_{1}s+a_0=0$$
    则：$${\mathbf{A}}^n+a_{n-1}{\mathbf{A}}^{n-1}+...+a_1\mathbf{A}+a_0\mathbf{I}=0$$
    理解: 由于一般的特征方程都是有限阶数的，因此凯莱哈密顿定理提供了矩阵A的有限维关系式。

第二章-特征值、特征向量和轨迹

18. 系统的极点：使任一闭合回路的传递函数值为1的那些s数值。

19. 对于固有振型回路传输增益为1。即T(S)=1。对于单位负反馈，则开环环路增益=-1时，闭合回路符合增益为1特性。p27页。

20. 对于二阶系统，采用类似的方法处理得到闭环系统的极点条件。

21. 闭环极点的理解：相当于满足回路公式的方程的解的模态。即表述方程组的解的特征。

22. 典型系统：对角线系统，当选择主对角线元素接近原始系统矩阵特征值时，系统表现为彼此解耦的独立积分环节。典型系统具有重要意义。
    $$\begin{gathered} \left[\begin{array}{c}{\dot{d}}_1\\ {\dot{d}}_2\\ {\dot{d}}_3\end{array}\right]=\left[\begin{array}{ccc}{\lambda }_1& 0& 0\\ 0& {\lambda }_2& 0\\ 0& 0& {\lambda }_3\end{array}\right]\left[\begin{array}{c}{d}_1\\ d_2\\ d_3\end{array}\right] \\ \dot{\mathbf{d}}=\Lambda \mathbf{d} \\ \mathbf{s}\mathbf{d}\left(\mathbf{s}\right)-\mathbf{d}\left(0\right)=\Lambda \mathbf{d}\left(\mathbf{s}\right) \\ \mathbf{d}\left(\mathbf{s}\right)={\left[\mathbf{s}\mathbf{I}-\Lambda \right]}^{-1}\mathbf{d}\left(0\right) \end{gathered}$$
    根据拉普拉斯反变换可以解得：
    $$\left[\begin{array}{c}{d}_1\left(t\right)\\ d_2\left(t\right)\\ d_3\left(t\right)\end{array}\right]=\left[\begin{array}{ccc}{e}^{{\lambda }_{1}t}& 0& 0\\ 0& e^{{\lambda }_{2}t}& 0\\ 0& 0& e^{{\lambda }_{3}t}\end{array}\right]\left[\begin{array}{c}{d}_1\left(0\right)\\ d_2\left(0\right)\\ d_3\left(0\right)\end{array}\right]$$

23. 对于上三角矩阵或者下三角矩阵，由于未构成附加闭环，因此系统矩阵的特征根仍然是主对角线项。

24. 特征向量在状态变量的解的应用：
    [ x 1 x 2 x 3 ] = e λ 1 t [ w 11 w 21 w 31 ] { f 1 } + e λ 2 t [ w 12 w 22 w 32 ] { f 2 } + e λ 2 t [ w 13 w 23 w 33 ] { f 3 } \left[ \begin{array}{c} x_1\\ x_2\\ x_3\\ \end{array} \right] =e^{\lambda _1t}\left[ \begin{array}{c} w_{11}\\ w_{21}\\ w_{31}\\ \end{array} \right] \left\{ f_1 \right\} +e^{\lambda _2t}\left[ \begin{array}{c} w_{12}\\ w_{22}\\ w_{32}\\ \end{array} \right] \left\{ f_2 \right\} +e^{\lambda _2t}\left[ \begin{array}{c} w_{13}\\ w_{23}\\ w_{33}\\ \end{array} \right] \left\{ f_3 \right\} ⎣⎡​x1​x2​x3​​⎦⎤​=eλ1​t⎣⎡​w11​w21​w31​​⎦⎤​{f1​}+eλ2​t⎣⎡​w12​w22​w32​​⎦⎤​{f2​}+eλ2​t⎣⎡​w13​w23​w33​​⎦⎤​{f3​}
    特征向量是特征根对应的模态映射至状态向量的解的特定组合。特征向量的分量给出了固有振型在各个状态上的分布。

25. 对角线系统可以得到各个振型的信号，特征向量又描述了状态向量的解的不同模态的组合。因此可以通过特征向量来实现系统状态向量和对角线系统的输出向量的关系。
    $$\left[\begin{array}{c}{x}_1\\ x_2\\ x_3\end{array}\right]=\left[\begin{array}{ccc}{w}_{11}& w_{12}& w_{13}\\ w_{21}& w_{22}& w_{23}\\ w_{31}& w_{32}& w_{33}\end{array}\right]\left[\begin{array}{c}{d}_1\\ d_2\\ d_3\end{array}\right]$$
    其中：w阵是以特征向量为列向量的矩阵。
    w阵实现了对角线系统的向量d—>状态向量x的转化。
    
    由于w阵可以实现x状态向量和d的转化，因此可以通过w阵实现d对线线系统的初值确定，在对角线系统计算完模态后，通过w阵再反算x的状态解。如下图所示：
    

26. 转化为对角线系统的过程：
    $$\begin{gathered} \dot{\mathbf{x}}=\mathbf{A}\mathbf{x},\text{令}\mathbf{x}=\mathbf{W}\mathbf{d} \\ \text{则：}\mathbf{d}={\mathbf{W}}^{-1}\mathbf{x} \\ \text{则：}\dot{\mathbf{d}}={\mathbf{W}}^{-1}\mathbf{A}\mathbf{W}\mathbf{d} \end{gathered}$$
    上式表明;特征向量矩阵是这种转换的主要特征。

27. 相似矩阵：存在相同的特征根，但系统矩阵不同。

28. 逆向思路：状态向量x的分量每一个都包含不同的振型分量，而且由分量由特征向量决定。那么：$$\mathbf{d}={\mathbf{W}}^{-1}\mathbf{x}$$
    可以实现单一振型的提取，可以实现振型滤波器的效果。

29. 特征向量与状态轨迹
    
    首先根据：
    $$\mathbf{d}={\mathbf{W}}^{-1}\mathbf{x}$$
    1） 根据初始状态，得到d(0)的初始状态。
    2） 在以d为状态变量的对角系统，容易获得d向量各分量的时域解；
    3） 通过x = W*d获得原始状态向量x的时域解；
    由:$$\mathbf{x}=\mathbf{W}\mathbf{d}$$
    可以看到：x的最终结果是以W（特征向量）的列向量为坐标系，以d(t)分量的轨迹。也就是说特征向量是状态变量X变化的参考方向。d是该方向的分量，共同构成最终的状态变量x的解。

30. 存在共轭复特征根的系统矩阵对应的特征值不能组成对角阵，而是在转置的位置存在非0值；

第三章：传递函数的状态变量描述

**学习本章的目的：**建立经典控制理论的传递函数分析方法和状态变量分析方法之间的联系。
31. 以一个例子复习单输入单输出传递函数的同伴型矩阵表示法。
能控标准型的转化方式链接

$$\begin{gathered} \frac{y\left(s\right)}{u\left(s\right)}=\frac{s^2+4s+5}{s^3+6s^2+11s+6} \\ \text{取中间变量}x: \\ x=\frac{u\left(s\right)}{s^3+6s^2+11s+6} \\ \text{可得：}\frac{d^{3}x}{d{t}^3}=-6\frac{d^{2}x}{d{t}^2}-11\frac{dx}{dt}-6x+u\left(t\right) \\ y\left(s\right)=s^{2}x+4sx+5x \\ \left[\begin{array}{c}{\dot{x}}_1\\ {\dot{x}}_2\\ {\dot{x}}_3\end{array}\right]=\left[\begin{array}{ccc}0& 1& 0\\ 0& 0& 1\\ -6& -11& -6\end{array}\right]\left[\begin{array}{c}{x}_1\\ x_2\\ x_3\end{array}\right]+\left[\begin{array}{c}0\\ 0\\ 1\end{array}\right]u\left(t\right) \\ y\left(t\right)=\left[\begin{array}{ccc}5& 4& 1\end{array}\right]\left[\begin{array}{c}{x}_1\\ x_2\\ x_3\end{array}\right] \end{gathered}$$

32. 状态空间描述表达式：
    $$\frac{y\left(s\right)}{u\left(s\right)}=\mathbf{C}{\left[s\mathbf{I}-\mathbf{A}\right]}^{-1}\mathbf{B}$$
    转化为对角阵形式，若W为系统矩阵A对应的一组特征向量。
    $$\begin{gathered} \frac{y\left(s\right)}{u\left(s\right)}=\mathbf{C}{\left[s\mathbf{I}-{\mathbf{W}\Lambda \mathbf{W}}^{-1}\right]}^{-1}\mathbf{B} \\ \frac{y\left(s\right)}{u\left(s\right)}=\mathbf{C}{\left[\mathbf{W}s{\mathbf{I}\mathbf{W}}^{-1}-{\mathbf{W}\Lambda \mathbf{W}}^{-1}\right]}^{-1}\mathbf{B} \\ \frac{y\left(s\right)}{u\left(s\right)}=\mathbf{C}{\left[\mathbf{W}\left(s\mathbf{I}-\Lambda \right){\mathbf{W}}^{-1}\right]}^{-1}\mathbf{B} \end{gathered}$$
    根据逆矩阵的结合律可得：
    $$\frac{y\left(s\right)}{u\left(s\right)}=\mathbf{C}\mathbf{W}{\left[s\mathbf{I}-\Lambda \right]}^{-1}{\mathbf{W}}^{-1}\mathbf{B}$$
    
33. 闭环系统传递函数–状态空间模型的对应关系
    
    采用闭环传递函数描述的系统，等价于修改了状态空间描述的中的系统矩阵。也就是说在状态空间描述中采用状态反馈或输出反馈等效于传递函数方法中加入的反馈矫正。原理是统一的。从矩阵运算的基本原理可以很明显的得到。**闭环传递函数和状态反馈（输出反馈）之间的联系。

第四章 状态变量反馈

难点

1. 凯莱哈密顿定理的原理与应用
   状态转移矩阵的计算
2. 特征向量的计算
   
   特征值与特征向量
3. 范德蒙德矩阵的应用
   
   矩阵的转置复习：
   矩阵的转置和加减乘除一样，也是一种运算，且满足下列运算规律（假设运算都是可行的）：
   
   矩阵的逆复习：
   

后续计划：

1. 存在复特征根的电机原理学习，holtz教授论文；
2. 经典控制理论多重极点的问题分析；

阅读时长：

2021-12-26：1+0.5(13:31)+晚上2小时+