816 控制理论攻坚

8.5.5 能控性与能观性的综合应用

在 816 考研专业课现代控制理论大题中,能控性与能观性常与参数辨识、坐标变换、极点配置、零极点对消等知识点综合考察。核心考题模型通常分为三类: 1. 参数待定与条件约束型:已知系统含未定参数,求解使系统完全能控、完全能观(或不可控/不可观)的参数充要条件; 2. 结构配置与多输入多输出设计型:

8.5.5 能控性与能观性的综合应用

一、核心考点与理论精要

在 816 考研专业课现代控制理论大题中,能控性与能观性常与参数辨识、坐标变换、极点配置、零极点对消等知识点综合考察。核心考题模型通常分为三类:

  1. 参数待定与条件约束型:已知系统含未定参数,求解使系统完全能控、完全能观(或不可控/不可观)的参数充要条件;
  2. 结构配置与多输入多输出设计型:根据约旦块特征值重数,论证单输入/单输出配置的局限性;
  3. 模态变换与对角化求解型:通过 Vandermonde 矩阵等变换化为对角型,求解待定参数。

二、标准解题模型与计算步骤

💡 【METHOD】 含参系统综合分析解题模板

  1. 代数判据法(低阶系统首选)
    • 构造能控性矩阵 M(p)=[BAB]M(p) = \begin{bmatrix} B & AB & \cdots \end{bmatrix} 与能观性矩阵 N(q)=[CCA]N(q) = \begin{bmatrix} C \\ CA \\ \cdots \end{bmatrix}
    • 计算行列式 det(M)\det(M)det(N)\det(N)
    • 建立方程/不等式:完全能控/能观     det0\iff \det \neq 0;不完全能控/能观     det=0\iff \det = 0
  2. 模态变换法(标准型或高阶系统首选)
    • 求出系统特征值 λi\lambda_i 及变换矩阵 PP
    • 计算变换后输入矩阵 Bˉ=P1B\bar{B} = P^{-1}B 与输出矩阵 Cˉ=CP\bar{C} = CP
    • 根据模态判据直接对 Bˉ\bar{B} 各行与 Cˉ\bar{C} 各列列出代数方程。

三、经典例题与真题精析

例题 1:二阶四参数系统可控可观条件确定

💡 【PROBLEM】

设系统的状态空间表达式为:

x˙(t)=[acdb]x(t)+[11]u(t),y(t)=[10]x(t)\dot{x}(t) = \begin{bmatrix} -a & c \\ -d & -b \end{bmatrix} x(t) + \begin{bmatrix} 1 \\ 1 \end{bmatrix} u(t), \quad y(t) = \begin{bmatrix} 1 & 0 \end{bmatrix} x(t)

试确定系统既状态完全能控、又状态完全能观时,参数 a,b,c,da, b, c, d 应满足的充要条件。

💡 【SOLUTION】 规范解答

解:

第一步:分析能控性条件 系统阶数 n=2n = 2。计算 AbAb

Ab=[acdb][11]=[cabd]Ab = \begin{bmatrix} -a & c \\ -d & -b \end{bmatrix} \begin{bmatrix} 1 \\ 1 \end{bmatrix} = \begin{bmatrix} c - a \\ -b - d \end{bmatrix}

构造能控性矩阵 QcQ_c

Qc=[bAb]=[1ca1bd]Q_c = \begin{bmatrix} b & Ab \end{bmatrix} = \begin{bmatrix} 1 & c - a \\ 1 & -b - d \end{bmatrix}

计算其行列式:

det(Qc)=1×(bd)1×(ca)=abcd\det(Q_c) = 1 \times (-b - d) - 1 \times (c - a) = a - b - c - d

系统状态完全能控的充要条件为:

det(Qc)0    abcd0\det(Q_c) \neq 0 \iff a - b - c - d \neq 0

第二步:分析能观性条件 计算 cAcA

cA=[10][acdb]=[ac]cA = \begin{bmatrix} 1 & 0 \end{bmatrix} \begin{bmatrix} -a & c \\ -d & -b \end{bmatrix} = \begin{bmatrix} -a & c \end{bmatrix}

构造能观性矩阵 QoQ_o

Qo=[ccA]=[10ac]Q_o = \begin{bmatrix} c \\ cA \end{bmatrix} = \begin{bmatrix} 1 & 0 \\ -a & c \end{bmatrix}

计算其行列式:

det(Qo)=1×c0×(a)=c\det(Q_o) = 1 \times c - 0 \times (-a) = c

系统状态完全能观的充要条件为:

det(Qo)0    c0\det(Q_o) \neq 0 \iff c \neq 0

结论: 当且仅当参数满足 abcd0a - b - c - d \neq 0c0c \neq 0 时,系统既完全能控又完全能观。


例题 2:重特征值约旦系统可控可观性论证

💡 【PROBLEM】

设系统状态空间表达式为:

x˙=[λ1000λ1000λ3]x+[abc]u,y=[def]x(λ1λ3)\dot{x} = \begin{bmatrix} \lambda_1 & 0 & 0 \\ 0 & \lambda_1 & 0 \\ 0 & 0 & \lambda_3 \end{bmatrix} x + \begin{bmatrix} a \\ b \\ c \end{bmatrix} u, \quad y = \begin{bmatrix} d & e & f \end{bmatrix} x \quad (\lambda_1 \neq \lambda_3)

试分析: (1) 能否通过选择参数 a,b,ca, b, c 使系统状态完全能控? (2) 能否通过选择参数 d,e,fd, e, f 使系统状态完全能观?

💡 【SOLUTION】 规范解答

解:

  1. 能控性分析
    • 矩阵 AA 中特征值 λ1\lambda_1 为二重特征根,且对应于两个独立的 1 阶约旦块(几何重数为 2)。
    • 根据 Gilbert 模态判据,对于相同特征值的多个独立约旦块,输入矩阵中对应各块末行的行向量组必须线性无关
    • 在本题中,输入为单输入(p=1p=1),对应 λ1\lambda_1 两个约旦块末行的行标量为 aabb。标量 aabb 作为一维向量必线性相关(即总存在常数使 k1a+k2b=0k_1 a + k_2 b = 0)。
    • 结论: 无论如何选择参数 a,b,ca, b, c,单输入系统绝不可能状态完全能控(至少需要 2 个独立控制输入)。
  2. 能观性分析
    • 同理,二重特征根 λ1\lambda_1 对应两个独立约旦块,要求 CC 中对应的列向量组线性无关。
    • 本题为单输出系统(q=1q=1),测量矩阵对应列为标量 ddee,一维标量必线性相关。
    • 结论: 无论如何选择参数 d,e,fd, e, f,单输出系统绝不可能状态完全能观

例题 3:二阶含参系统不可控与不可观参数求解

💡 【PROBLEM】

已知系统状态空间表达式为:

x˙=[11210]x+[p1]u,y=[q1]x\dot{x} = \begin{bmatrix} 1 & 12 \\ 1 & 0 \end{bmatrix} x + \begin{bmatrix} p \\ -1 \end{bmatrix} u, \quad y = \begin{bmatrix} q & 1 \end{bmatrix} x

试确定参数 ppqq 取何值时,系统分别处于不能控与不能观状态。

💡 【SOLUTION】 规范解答

解:

1. 求解不可控参数 pp: 计算 AbAb

Ab=[11210][p1]=[p12p]Ab = \begin{bmatrix} 1 & 12 \\ 1 & 0 \end{bmatrix} \begin{bmatrix} p \\ -1 \end{bmatrix} = \begin{bmatrix} p - 12 \\ p \end{bmatrix}

构造能控性矩阵:

M=[bAb]=[pp121p]M = \begin{bmatrix} b & Ab \end{bmatrix} = \begin{bmatrix} p & p - 12 \\ -1 & p \end{bmatrix}

计算行列式:

det(M)=p2(1)(p12)=p2+p12=(p+4)(p3)\det(M) = p^2 - (-1)(p - 12) = p^2 + p - 12 = (p + 4)(p - 3)

系统不完全能控     det(M)=0    p=4\iff \det(M) = 0 \iff p = -4p=3p = 3

2. 求解不可观参数 qq: 计算 cAcA

cA=[q1][11210]=[q+112q]cA = \begin{bmatrix} q & 1 \end{bmatrix} \begin{bmatrix} 1 & 12 \\ 1 & 0 \end{bmatrix} = \begin{bmatrix} q + 1 & 12q \end{bmatrix}

构造能观性矩阵:

N=[ccA]=[q1q+112q]N = \begin{bmatrix} c \\ cA \end{bmatrix} = \begin{bmatrix} q & 1 \\ q + 1 & 12q \end{bmatrix}

计算行列式:

det(N)=q(12q)1(q+1)=12q2q1=(3q1)(4q+1)\det(N) = q(12q) - 1(q + 1) = 12q^2 - q - 1 = (3q - 1)(4q + 1)

系统不完全能观     det(N)=0    q=13\iff \det(N) = 0 \iff q = \frac{1}{3}q=14q = -\frac{1}{4}


例题 4:能控标准型通过对角化设计不可观测量参数

💡 【PROBLEM】

已知线性系统状态方程为能控标准型:

x˙=[0100016116]x+[001]u,y=[c1c2c3]x\dot{x} = \begin{bmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -6 & -11 & -6 \end{bmatrix} x + \begin{bmatrix} 0 \\ 0 \\ 1 \end{bmatrix} u, \quad y = \begin{bmatrix} c_1 & c_2 & c_3 \end{bmatrix} x

试通过对角化方法,求出使系统不可观测的一组不全为 0 的参数 c1,c2,c3c_1, c_2, c_3

💡 【SOLUTION】 规范解答

解:

第一步:求特征值 矩阵 AA 为能控标准型,直接由最后一行读出特征多项式:

λIA=λ3+6λ2+11λ+6=(λ+1)(λ+2)(λ+3)=0|\lambda I - A| = \lambda^3 + 6\lambda^2 + 11\lambda + 6 = (\lambda + 1)(\lambda + 2)(\lambda + 3) = 0

解得互异特征值:λ1=1,λ2=2,λ3=3\lambda_1 = -1, \lambda_2 = -2, \lambda_3 = -3

第二步:构造对角化变换矩阵 PP(Vandermonde 阵) 对应特征向量构成的变换矩阵为:

P=[111λ1λ2λ3λ12λ22λ32]=[111123149]P = \begin{bmatrix} 1 & 1 & 1 \\ \lambda_1 & \lambda_2 & \lambda_3 \\ \lambda_1^2 & \lambda_2^2 & \lambda_3^2 \end{bmatrix} = \begin{bmatrix} 1 & 1 & 1 \\ -1 & -2 & -3 \\ 1 & 4 & 9 \end{bmatrix}

x=Pzx = Pz,则变换后对角阵为 Λ=diag(1,2,3)\Lambda = \text{diag}(-1, -2, -3)

第三步:计算模态测量向量 cˉ=cP\bar{c} = cP

cˉ=[c1c2c3][111123149]\bar{c} = \begin{bmatrix} c_1 & c_2 & c_3 \end{bmatrix} \begin{bmatrix} 1 & 1 & 1 \\ -1 & -2 & -3 \\ 1 & 4 & 9 \end{bmatrix} =[c1c2+c3c12c2+4c3c13c2+9c3]= \begin{bmatrix} c_1 - c_2 + c_3 & c_1 - 2c_2 + 4c_3 & c_1 - 3c_2 + 9c_3 \end{bmatrix}

第四步:确定不可观参数条件 根据对角型模态判据,系统不可观     cˉ\iff \bar{c} 中至少有一个元素为 00

c1c2+c3=0c12c2+4c3=0c13c2+9c3=0c_1 - c_2 + c_3 = 0 \quad \text{或} \quad c_1 - 2c_2 + 4c_3 = 0 \quad \text{或} \quad c_1 - 3c_2 + 9c_3 = 0

取特解:令 c1c2+c3=0c_1 - c_2 + c_3 = 0,可取 c1=1,c2=2,c3=1c_1 = 1, c_2 = 2, c_3 = 1(此时第 1 个模态不可观)。 验算:取 c=[121]c = \begin{bmatrix} 1 & 2 & 1 \end{bmatrix} 时,系统状态不完全能观。


例题 5:按指定模态配置参数使状态分别满足可控/不可观

💡 【PROBLEM】

设系统状态空间表达式为:

x˙(t)=[0123]x(t)+[b1b2]u(t),y(t)=[c1c2]x(t)\dot{x}(t) = \begin{bmatrix} 0 & 1 \\ -2 & -3 \end{bmatrix} x(t) + \begin{bmatrix} b_1 \\ b_2 \end{bmatrix} u(t), \quad y(t) = \begin{bmatrix} c_1 & c_2 \end{bmatrix} x(t)

试确定参数 b1,b2,c1,c2b_1, b_2, c_1, c_2 满足的条件,使得系统经变换后的一个模态既能控又能观,而另一个模态既不能控又不能观。

💡 【SOLUTION】 规范解答

解:

第一步:求系统特征值与对角化矩阵

sIA=det[s12s+3]=s(s+3)+2=s2+3s+2=(s+1)(s+2)=0|sI - A| = \det \begin{bmatrix} s & -1 \\ 2 & s + 3 \end{bmatrix} = s(s + 3) + 2 = s^2 + 3s + 2 = (s + 1)(s + 2) = 0

得特征值 s1=1,s2=2s_1 = -1, s_2 = -2。 对应的特征向量为 v1=[11],v2=[12]v_1 = \begin{bmatrix} 1 \\ -1 \end{bmatrix}, v_2 = \begin{bmatrix} 1 \\ -2 \end{bmatrix}。 变换矩阵及其逆矩阵为:

P=[1112],P1=[2111]P = \begin{bmatrix} 1 & 1 \\ -1 & -2 \end{bmatrix}, \quad P^{-1} = \begin{bmatrix} 2 & 1 \\ -1 & -1 \end{bmatrix}

第二步:计算模态矩阵 Bˉ\bar{B}Cˉ\bar{C}x=Pzx = Pz,变换后的对角系统为:

Aˉ=[1002]\bar{A} = \begin{bmatrix} -1 & 0 \\ 0 & -2 \end{bmatrix} Bˉ=P1B=[2111][b1b2]=[2b1+b2b1b2]\bar{B} = P^{-1}B = \begin{bmatrix} 2 & 1 \\ -1 & -1 \end{bmatrix} \begin{bmatrix} b_1 \\ b_2 \end{bmatrix} = \begin{bmatrix} 2b_1 + b_2 \\ -b_1 - b_2 \end{bmatrix} Cˉ=CP=[c1c2][1112]=[c1c2c12c2]\bar{C} = CP = \begin{bmatrix} c_1 & c_2 \end{bmatrix} \begin{bmatrix} 1 & 1 \\ -1 & -2 \end{bmatrix} = \begin{bmatrix} c_1 - c_2 & c_1 - 2c_2 \end{bmatrix}

第三步:按要求列出参数约束方程

  • 情形 1:设模态 z1z_1 能控且能观,模态 z2z_2 不能控且不能观

    {2b1+b20c1c20b1b2=0    b1=b2c12c2=0    c1=2c2\begin{cases} 2b_1 + b_2 \neq 0 \\ c_1 - c_2 \neq 0 \\ -b_1 - b_2 = 0 \implies b_1 = -b_2 \\ c_1 - 2c_2 = 0 \implies c_1 = 2c_2 \end{cases}

    代入验证:当 b1=b20b_1 = -b_2 \neq 0c1=2c20c_1 = 2c_2 \neq 0 时,2b1+b2=b102b_1+b_2 = b_1 \neq 0c1c2=c20c_1-c_2 = c_2 \neq 0 自然满足。

  • 情形 2:设模态 z2z_2 能控且能观,模态 z1z_1 不能控且不能观

    {2b1+b2=0    b2=2b1c1c2=0    c1=c2b1b20c12c20\begin{cases} 2b_1 + b_2 = 0 \implies b_2 = -2b_1 \\ c_1 - c_2 = 0 \implies c_1 = c_2 \\ -b_1 - b_2 \neq 0 \\ c_1 - 2c_2 \neq 0 \end{cases}

结论: 满足要求的参数条件为:

{b1=b20c1=2c20{b2=2b10c1=c20\begin{cases} b_1 = -b_2 \neq 0 \\ c_1 = 2c_2 \neq 0 \end{cases} \quad \text{或} \quad \begin{cases} b_2 = -2b_1 \neq 0 \\ c_1 = c_2 \neq 0 \end{cases}

四、易错陷阱与避坑指南

💡 【PITFALL】 考场易错点

  1. 变换矩阵 PPP1P^{-1} 的乘法次序混淆
    • 模态输入矩阵是 Bˉ=P1B\bar{B} = P^{-1}B(左乘逆矩阵);
    • 模态输出矩阵是 Cˉ=CP\bar{C} = CP(右乘原矩阵)。千万不可颠倒!
  2. 非零条件漏判
    • 题目要求“使一个状态可控可观”,解出等式关系后(如 b1=b2b_1 = -b_2),必须补充非零条件(b10b_1 \neq 0),否则若取 b1=b2=0b_1 = b_2 = 0 会导致整个系统全不可控。

Discussion

Comments

Post

Thoughts, corrections, and follow-up notes are welcome here.