【现代控制理论笔记】——第五章：能控、能观和传递函数

现代控制理论

第五章

说明：

对于最小实现问题，先判断系统是否是严格真的，如果不是则用D换成严格真的；如果是，则判断是不是可简约的，如果是则化成不可简约的，如果不是，则写出能控标准型实现即为最小实现

注意，其中可能涉及到能控/能观分解，原理是针对能控性/能观性矩阵，找到能控列/能观行部分作为T矩阵的前半部分，再找到不能控/不能观部分对应的与其线性无关的部分作为后半部分，然后能观用$TA{T}^{-1},TB,C{T}^{-1}$，能控用$T^{-1}AT,T^{-1}B,CT$得到分解后的矩阵

1. 不完全能控不完全能观系统的结构分解

2. 不完全能控不完全能观系统的传递函数

这从另一个角度说明了为什么系统的传递函数并不一定能完全展示系统的结构：

只有对完全能观能控的系统,传递函数才足以表征系统的结构描述才是完全的。

3. 能控能观系统的传递函数

对于系统：

$$\begin{gathered} \dot{x}=Ax+Bu \\ y=Cx \end{gathered}$$

其系统传递函数为：

$G(s)=C(sI-A{)}^{-1}B=C\frac{adj(sI-A)}{det(sI-A)}B$

其也可以写作：

$G(s)=C\frac{P(s)}{\varphi (s)}B$

的形式，其中$P(s)$是矩阵多项式，$\varphi (s)$是最小多项式。

说明：

① A的最小多项式是其特征多项式的因式，必满足$\varphi (A)=0$

② $P(s)A=AP(s)$

4. 能控能观系统的传递函数不可简约性

在复平面无根即P(s)B=0无解，反证法。

→

可以推出：

说明：

上式说明：若$(A,B,C)$能控能观，则传递函数是不可简约的（既约的）。

也就是说，

但是若传递函数不可简约，并不一定系统就能控能观。

5. 最小实现

如果可以找到一个关系式：

$$\begin{gathered} \dot{x}=Ax+Bu \\ y=Cx \end{gathered}$$

使$C(sI-A{)}^{-1}B+D=G(s)$

则称其为传递函数矩阵$G(s)$的一个实现。

那么给定传递函数是否一定可以找到实现呢？

① 对于严格真的传递函数矩阵：p为输入维数

与能控标准型和能观标准型类似。

② 对于真的传递函数矩阵：

可以先令$D=li{m}_{s\to \infty }G(s)$，用$\bar{G}(s)=G(s)-D$求出$(A,B,C)$，则其实现为$(A,B,C,D)$

其中，A的维数为实现的维数。$G(s)$中维数最小的实现，称为最小实现。

这说明对于严格真的传递函数，能控/能观实现为系统的一个最小实现。

且有$G(s)$的任意两个最小实现代数等价。

5.1 如何求最小实现

注意：

一定要把矩阵化简成首一式：$s^n$的系数为1！！！！

1）对于单输入单输出的系统，传递函数

$G(s)=\frac{c_{n-1}{s}^{n-1}+\cdots +c_{1}s+c_0}{s^n+a_{n-1}{s}^{n-1}+\cdots +a_{1}s+a_0}$

不可简约，其可控性实现和可观性实现即为最小实现。

$A=\left[\begin{array}{cccc}0& 1& & \\ & \ddots & \ddots & \\ & & 0& 1\\ -a_0& -a_1& \cdots & -a_{n-1}\end{array}\right]$

$B=\left[\begin{array}{c}0\\ ⋮\\ 0\\ 0\\ 1\end{array}\right]$

$C=\left[\begin{array}{cccc}{c}_0& c_1& \cdots & c_{n-1}\end{array}\right]$

2）对单输入多输出的系统，传递函数：

$G(s)=\frac{C_{n-1}{s}^{n-1}+\cdots +C_{1}s+C_0}{s^n+a_{n-1}{s}^{n-1}+\cdots +a_{1}s+a_0}$

不可简约，因此，其可控性实现/可观性实现即为最小实现。

$A=\left[\begin{array}{cccc}0& 1& & \\ & \ddots & \ddots & \\ & & 0& 1\\ -a_0& -a_1& \cdots & -a_{n-1}\end{array}\right]$

$B=\left[\begin{array}{c}0\\ ⋮\\ 0\\ 1\end{array}\right]$

$C=\left[\begin{array}{cccc}{C}_0& C_1& \cdots & C_{n-1}\end{array}\right]$

3）给定$q\times p$严格真的传递函数矩阵$G(s)$有n个单的实数极点$s_1,...s_n$，则传递函数矩阵可写为：

$G(s)=\frac{P_1}{s-s_1}+\frac{P_2}{s-s_2}+\cdots +\frac{P_n}{s-s_n},$

其中，P可以写成：

$P_i={\lim }_{s\to s_i}(s-s_i)G(s)$

对P进行满秩分解，则有：

$P_i=C_i{B}_i$

满足：$rank{B}_i=rank{C}_i=rank{P}_i$

则系统的最小实现为：

$A=\left[\begin{array}{cccc}{s}_1{I}_{r_1}& & & \\ & s_2{I}_{r_2}& & \\ & & \ddots & \\ & & & s_n{I}_{r_n}\end{array}\right]$

$B=\left[\begin{array}{c}{B}_1\\ B_2\\ ⋮\\ B_n\end{array}\right]$

$C=\left[\begin{array}{cccc}{C}_1& C_2& \cdots & C_n\end{array}\right]$

4）对于更一般的情况，可以先写出能控性实现或能观性实现，再进行能控/能观分解，得到能控能观子系统，即为最小实现。

5.2 满秩分解

对于一个矩阵A：

$A=[\begin{array}{rrrr}1& 2& 3& 0\\ 0& 2& 1& -1\\ 1& 0& 2& 1\end{array}]$

其经过变换矩阵P进行初等行变换得到B：

$[A,I]=[\begin{array}{rrrrrrr}1& 2& 3& 0& 1& 0& 0\\ 0& 2& 1& -1& 0& 1& 0\\ 0& 0& 0& 0& -1& 1& 1\end{array}]=\left[\begin{array}{c}B,P\end{array}\right]$

阶梯矩阵B可以写成：

$B=[D;0]=[\begin{array}{cccc}1& 2& 3& 0\\ 0& 2& 1& -1\\ 0& 0& 0& 0\end{array}]$

的形式，且$P^{-1}$可以写成分块的形式：

$P^{-1}=[C,M]=[\begin{array}{rrr}1& 0& 0\\ 0& 1& 0\\ 1& -1& 1\end{array}]$

进而可以得到满秩分解：

$A=[C,M]\cdot [D;0]=C\cdot D$

6. Jordan标准型实现

6.1 由状态空间求Jordan标准型：

对于系统：

$$\begin{gathered} \dot{\mathbf{x}}=\left(\begin{array}{ccc}0& 1& 0\\ 0& 0& 1\\ 2& 3& 0\end{array}\right)\mathbf{x}+\left(\begin{array}{c}0\\ 0\\ 1\end{array}\right)\mathbf{u} \\ \mathbf{y}=(1,0,0)\mathbf{x} \end{gathered}$$

其特征值为：$\lambda =-1,-1,2$

求特征值对应的特征向量：

1）对于单根，由$A{P}_i={\lambda }_i{P}_i$：

$[\begin{array}{ccc}0& 1& 0\\ 0& 0& 1\\ 2& 3& 0\end{array}]\left[\begin{array}{c}{p}_{11}\\ p_{21}\\ p_{31}\end{array}\right]=-\left[\begin{array}{c}{p}_{11}\\ p_{21}\\ p_{31}\end{array}\right]$

得：

$P_1=\left[\begin{array}{c}1\\ -1\\ 1\end{array}\right]$

进一步可以得到变换矩阵：

$T=\left[\begin{array}{ccc}{p}_1& p_2& p_3\end{array}\right]$

则Jordan标准形实现为：

$A=\left[\begin{array}{ccc}-1& 1& 0\\ 0& -1& 0\\ 0& 0& 2\end{array}\right]$

$B=T^{-1}B$

$C=CT$

6.2 由传递函数求Jordan标准型

对于传递函数：

$W(s)=\frac{6(s+1)}{s(s+2)(s+3{)}^2}$

使用留数法进行因式分解：

$K_2=-\frac{10}{3}=li{m}_{s\to -3}\frac{1}{(2-1)!}\frac{d}{d{s}^{2-1}}[G(s)(s+3{)}^2]$

可得：

$W(S)=\frac{-4}{(S+3{)}^2}+\frac{-\frac{10}{3}}{S+3}+\frac{3}{S+2}+\frac{\frac{1}{3}}{S}$

则可以直接写出Jordan标准型：

$A=\left[\begin{array}{cccc}-3& 1& 0& 0\\ 0& -3& 0& 0\\ 0& 0& -2& 0\\ 0& 0& 0& 0\end{array}\right]$

$B=\left[\begin{array}{c}0\\ 1\\ 1\\ 1\end{array}\right]$

$C=\left[\begin{array}{cccc}-4& -\frac{10}{3}& 3& \frac{1}{3}\end{array}\right]$

6.3例题

对于系统：

$\begin{array}{rl}G(s)& =\frac{3s^2+17s+25}{{\left(s+2\right)}^3\left(s+3\right)}\end{array}$

其Jordan标准型实现为：

$A=\left[\begin{array}{cccc}-2& 1& 0& 0\\ 0& -2& 1& 0\\ 0& 0& -2& 0\\ 0& 0& 0& -3\end{array}\right]$

$B=\left[\begin{array}{c}0\\ 0\\ 1\\ 1\end{array}\right]$

$C=\left[\begin{array}{cccc}3& 2& 1& -1\end{array}\right]$