The Lyapunov inequality is an indispensable tool for stability analysis in linear control theory. It provides a necessary and sufficient condition for the stability of an autonomous linear-time invariant system in terms of the existence of a symmetric positive-definite Lyapunov matrix. This work proposes a new variant of this inequality in which the constituent Lyapunov matrix is allowed to be asymmetric. After analysing the properties of the proposed inequality for a class of matrices, we derive new results for the stabilisation of linear systems. Subsequently, we utilize the developed results to obtain sufficient conditions for the suboptimal linear quadratic control design problem, where addition to having an asymmetric Lyapunov matrix, which serves as a design matrix for this problem, we provide a characterization of the cost associated with the computed stabilizing suboptimal control laws by deriving an expression for the upper bound on cost in terms of the initial conditions of the system. We demonstrate the applicability of the proposed results using two numerical examples -- one for suboptimal control design for a linear time-invariant system and another for the consensus (state-agreement) protocol design for a multi-agent system where-in we see how the asymmetry of the design matrix emerges as an inherent requirement for the problem.
Lyapunov-like Stability Inequality with an Asymmetric Matrix and Application to Suboptimal LQ Control Design 论文ID : 2502.11556标题 : Lyapunov-like Stability Inequality with an Asymmetric Matrix and Application to Suboptimal LQ Control Design作者 : Avinash Kumar分类 : math.OC (Optimization and Control)发表时间 : 2025年2月 (arXiv预印本)论文链接 : https://arxiv.org/abs/2502.11556 Lyapunov不等式是线性控制理论中稳定性分析的重要工具,它通过对称正定Lyapunov矩阵的存在性为自治线性时不变系统的稳定性提供了充分必要条件。本文提出了该不等式的新变体,允许Lyapunov矩阵为非对称形式。在分析了该不等式对一类矩阵的性质后,作者推导出了线性系统镇定的新结果。随后,利用这些结果获得了次优线性二次控制设计问题的充分条件,其中非对称Lyapunov矩阵作为设计矩阵,并通过推导基于系统初始条件的成本上界表达式来刻画计算得到的镇定次优控制律的相关成本。
传统Lyapunov理论的局限性 : 经典Lyapunov不等式 A T P + P A ≺ 0 A^T P + PA \prec 0 A T P + P A ≺ 0 要求Lyapunov矩阵P必须是对称正定的,这一假设在所有基于Lyapunov不等式的结果中都普遍存在。对称假设的原因 :对称矩阵具有良好性质(实特征值、可对角化) 二次型 z T M z z^T Mz z T M z 满足 z T M z = z T M s y m z z^T Mz = z^T M_{sym}z z T M z = z T M sy m z ,其中 M s y m = 1 2 ( M + M T ) M_{sym} = \frac{1}{2}(M + M^T) M sy m = 2 1 ( M + M T ) 正定性概念容易从二次型转换到对称矩阵 次优LQ控制的需求 : 最优LQ控制解可能由于输入饱和约束、反馈增益矩阵的范数有界要求或结构化条件而无法实现,因此需要研究次优解。扩大设计空间 : 通过放松对称性假设,将设计矩阵P的搜索空间从对称正定矩阵扩展到更大的矩阵集合避免预设成本上界 : 现有方法需要预先指定成本上界γ,可能导致不可行性,特别是在网络化控制设计中满足结构化控制需求 : 某些问题(如多智能体一致性协议)本质上需要非对称设计矩阵提出新的Lyapunov-like稳定性不等式(LSI) : A P + P T A ≺ 0 AP + P^T A \prec 0 A P + P T A ≺ 0 ,其中P不要求对称建立非对称矩阵的稳定性理论 : 证明了当A为对称负定时,满足LSI的矩阵P具有正实部特征值开发次优LQ控制新方法 : 基于初始条件计算成本上界,而非预先指定提供成本量化结果 : 推导出显式的成本上界表达式 γ ( x 0 ) = trace ( P − P ˉ + P ^ ) x 0 T x 0 \gamma(x_0) = \text{trace}(P - \bar{P} + \hat{P})x_0^T x_0 γ ( x 0 ) = trace ( P − P ˉ + P ^ ) x 0 T x 0 考虑线性时不变系统:
x ˙ ( t ) = A x ( t ) + B u ( t ) \dot{x}(t) = Ax(t) + Bu(t) x ˙ ( t ) = A x ( t ) + B u ( t )
目标是设计控制律 u ( t ) = − R − 1 B T P x ( t ) u(t) = -R^{-1}B^T Px(t) u ( t ) = − R − 1 B T P x ( t ) 使系统稳定,并最小化二次成本:
J = ∫ 0 ∞ ( x T ( t ) Q x ( t ) + u T ( t ) R u ( t ) ) d t J = \int_0^{\infty} (x^T(t)Qx(t) + u^T(t)Ru(t))dt J = ∫ 0 ∞ ( x T ( t ) Q x ( t ) + u T ( t ) R u ( t )) d t
命题1 : 假设 A ∈ R n × n A \in \mathbb{R}^{n \times n} A ∈ R n × n 是对称负定矩阵。若存在矩阵 P ∈ R n × n P \in \mathbb{R}^{n \times n} P ∈ R n × n 满足LSI:
A P + P T A ≺ 0 AP + P^T A \prec 0 A P + P T A ≺ 0
则 Re ( λ i ( P ) ) > 0 , ∀ i ∈ { 1 , 2 , … , n } \text{Re}(\lambda_i(P)) > 0, \forall i \in \{1,2,\ldots,n\} Re ( λ i ( P )) > 0 , ∀ i ∈ { 1 , 2 , … , n } ,因此 trace ( P ) > 0 \text{trace}(P) > 0 trace ( P ) > 0 。
证明思路 : 利用任意特征值λ和对应特征向量v,通过 v H ( A P + P T A ) v < 0 v^H(AP + P^T A)v < 0 v H ( A P + P T A ) v < 0 推导出 Re ( λ ) > 0 \text{Re}(\lambda) > 0 Re ( λ ) > 0 。
命题2 : 若存在矩阵 X ≻ 0 X \succ 0 X ≻ 0 和 Y Y Y 满足:
X A T + A X − Y T B R − 1 B T − B R − 1 B T Y ≺ 0 XA^T + AX - Y^T BR^{-1}B^T - BR^{-1}B^T Y \prec 0 X A T + A X − Y T B R − 1 B T − B R − 1 B T Y ≺ 0
则矩阵 A − B R − 1 B T P A - BR^{-1}B^T P A − B R − 1 B T P 是Hurwitz的,其中 P = Y X − 1 P = YX^{-1} P = Y X − 1 。
定理3 : 假设存在矩阵 X ≻ 0 X \succ 0 X ≻ 0 , W ≻ 0 W \succ 0 W ≻ 0 , Y Y Y , P P P 满足:
X A T + A X − Y T B R − 1 B T − B R − 1 B T Y ≺ 0 XA^T + AX - Y^T BR^{-1}B^T - BR^{-1}B^T Y \prec 0 X A T + A X − Y T B R − 1 B T − B R − 1 B T Y ≺ 0 [ A T P + P T A − Q + W ( I + Y T A ) T ( P T B ) ( P T B ) T R ] ≻ 0 \begin{bmatrix} A^T P + P^T A - Q + W & (I + Y^T A)^T(P^T B) \\ (P^T B)^T & R \end{bmatrix} \succ 0 [ A T P + P T A − Q + W ( P T B ) T ( I + Y T A ) T ( P T B ) R ] ≻ 0
若闭环矩阵 A − B R − 1 B T P A - BR^{-1}B^T P A − B R − 1 B T P 对称且存在 P ^ \hat{P} P ^ 使得:
( A − B R − 1 B T P ) P ^ + P ^ T ( A − B R − 1 B T P ) + W = 0 (A - BR^{-1}B^T P)\hat{P} + \hat{P}^T(A - BR^{-1}B^T P) + W = 0 ( A − B R − 1 B T P ) P ^ + P ^ T ( A − B R − 1 B T P ) + W = 0
则控制律解决Problem 2,成本上界为:
γ ( x 0 ) = trace ( P − P ˉ + P ^ ) x 0 T x 0 \gamma(x_0) = \text{trace}(P - \bar{P} + \hat{P})x_0^T x_0 γ ( x 0 ) = trace ( P − P ˉ + P ^ ) x 0 T x 0
非对称设计矩阵 : 首次系统性地研究了非对称Lyapunov矩阵在控制设计中的应用基于初始条件的成本界 : 直接根据系统初始状态计算成本上界,避免预设参数结构化控制适用性 : 特别适合需要非对称设计矩阵的结构化控制问题系统矩阵 : A = [ 1 2 0 2 ] A = \begin{bmatrix} 1 & 2 \\ 0 & 2 \end{bmatrix} A = [ 1 0 2 2 ] , B = [ 4 2 0 2 ] B = \begin{bmatrix} 4 & 2 \\ 0 & 2 \end{bmatrix} B = [ 4 0 2 2 ] 初始条件 : x 0 = [ 0.1 , − 0.2 ] T x_0 = [0.1, -0.2]^T x 0 = [ 0.1 , − 0.2 ] T 权重矩阵 : Q = 10 I 2 × 2 Q = 10I_{2 \times 2} Q = 10 I 2 × 2 , R = 0.05 I 2 × 2 R = 0.05I_{2 \times 2} R = 0.05 I 2 × 2 工具 : CVX工具箱求解LMI智能体数量 : 4个标量不稳定智能体动态方程 : x ˙ i ( t ) = x i ( t ) + u i ( t ) \dot{x}_i(t) = x_i(t) + u_i(t) x ˙ i ( t ) = x i ( t ) + u i ( t ) 通信拓扑 : 有向图,信息交换受限目标 : 实现状态一致性 ∣ x i ( t ) − x j ( t ) ∣ → 0 |x_i(t) - x_j(t)| \to 0 ∣ x i ( t ) − x j ( t ) ∣ → 0 稳定性 : 闭环系统的Hurwitz性质成本性能 : 实际成本J与理论上界 γ ( x 0 ) \gamma(x_0) γ ( x 0 ) 的比较可行性 : LMI约束的可解性最优成本 : J ∗ = 0.0207 J^* = 0.0207 J ∗ = 0.0207 次优控制律 : u ( t ) = − [ 1.3604 − 1.0302 0.0034 3.0638 ] x ( t ) u(t) = -\begin{bmatrix} 1.3604 & -1.0302 \\ 0.0034 & 3.0638 \end{bmatrix}x(t) u ( t ) = − [ 1.3604 0.0034 − 1.0302 3.0638 ] x ( t ) 设计矩阵 : P = [ 0.0170 − 0.0129 − 0.0169 0.0895 ] P = \begin{bmatrix} 0.0170 & -0.0129 \\ -0.0169 & 0.0895 \end{bmatrix} P = [ 0.0170 − 0.0169 − 0.0129 0.0895 ] (非对称)成本上界 : γ ( x 0 ) = 1.1374 \gamma(x_0) = 1.1374 γ ( x 0 ) = 1.1374 实际成本 : J = 0.0627 J = 0.0627 J = 0.0627 结构化增益矩阵 : 由通信拓扑决定的稀疏结构非对称设计矩阵 : P = [ 0.0350 0.0350 0 0 0.1049 0 0 0 0.1049 ] P = \begin{bmatrix} 0.0350 & 0.0350 & 0 \\ 0 & 0.1049 & 0 \\ 0 & 0 & 0.1049 \end{bmatrix} P = 0.0350 0 0 0.0350 0.1049 0 0 0 0.1049 成本上界 : γ ( x 0 ) = 1.8426 \gamma(x_0) = 1.8426 γ ( x 0 ) = 1.8426 实际成本 : J = 0.1407 J = 0.1407 J = 0.1407 非对称性的必要性 : 例子2展示了某些问题本质上要求非对称设计矩阵成本界的有效性 : 理论上界能够有效约束实际成本方法的实用性 : 能够处理结构化控制约束问题经典结果 : Lyapunov不等式 A T P + P A ≺ 0 A^T P + PA \prec 0 A T P + P A ≺ 0 提供稳定性的充分必要条件LMI发展 : 线性矩阵不等式在控制理论中的广泛应用现有方法 : 预先指定成本上界γ,推导充分条件局限性 : 可能导致不可行性,特别是在不知道最优成本的情况下多智能体系统 : 一致性协议设计中的结构化约束网络化控制 : 通信拓扑对控制器结构的限制理论贡献 : 提出了允许非对称Lyapunov矩阵的新稳定性不等式实用价值 : 为次优LQ控制提供了新的设计方法应用前景 : 特别适用于需要结构化控制的网络化系统适用条件 : 要求闭环矩阵为对称负定,限制了方法的普适性计算复杂性 : 需要求解多个耦合的LMI约束保守性 : 提供的是充分条件而非必要条件扩展到更一般的系统 : 研究非对称系统矩阵的情况网络化控制应用 : 探索在更复杂网络拓扑中的应用鲁棒性分析 : 考虑参数不确定性和外部干扰创新性强 : 首次系统性地放松了Lyapunov矩阵的对称性假设理论严谨 : 提供了完整的数学证明和理论分析实用性好 : 通过数值例子验证了方法的有效性应用价值 : 特别适合结构化控制问题适用范围受限 : 对系统矩阵的对称性要求限制了方法的普适性实验不够充分 : 仅提供了两个相对简单的数值例子与现有方法比较不足 : 缺乏与其他次优控制方法的定量比较计算效率分析缺失 : 未讨论LMI求解的计算复杂性理论贡献 : 为Lyapunov理论提供了新的视角和工具实用价值 : 为结构化控制设计提供了新思路研究启发 : 可能开启非对称设计矩阵研究的新方向多智能体系统 : 一致性协议设计网络化控制 : 具有通信约束的分布式控制结构化控制 : 需要特定增益矩阵结构的控制问题次优控制 : 当最优解不可实现时的替代方案论文引用了控制理论的经典文献,包括:
Lyapunov稳定性理论的基础文献 线性矩阵不等式的相关研究 次优LQ控制的现有方法 多智能体系统一致性控制的相关工作 总体评价 : 这是一篇具有创新性的控制理论论文,通过放松传统Lyapunov矩阵的对称性假设,为次优控制设计提供了新的理论工具。虽然在适用范围和实验验证方面还有改进空间,但其理论贡献和实用价值值得认可,特别是在结构化控制问题中展现出的优势。