2025-11-17T05:46:12.770787

Data-Driven Structured Controller Design Using the Matrix S-Procedure

Yang, Zhong, Yang et al.
This paper focuses on the data-driven optimal structured controller design for discrete-time linear time-invariant (LTI) systems, considering both the $H_2$ performance and the $H_\infty$ performance. Specifically, we consider three scenarios: (i) the model-based structured control, (ii) the data-driven unstructured control, and (iii) the data-driven structured control. For the $H_2$ performance, we primarily investigate cases (ii) and (iii), since case (i) has been extensively studied in the literature. For the $H_\infty$ performance, all three scenarios are considered. For the structured control, we introduce a linearization technique that transforms the original nonconvex problem into a semidefinite programming (SDP) problem. Based on this transformation, we develop an iterative linear matrix inequality (ILMI) algorithm. For the data-driven control, we describe the set of all possible system matrices that can generate the sequence of collected data. Additionally, we propose a sufficient condition to handle all possible system matrices using the matrix S-procedure. The data-driven structured control is followed by combining the previous two cases. We compare our methods with those in the existing literature and demonstrate our superiority via several numerical simulations.
academic

Data-Driven Structured Controller Design Using the Matrix S-Procedure

基本信息

  • 论文ID: 2503.14949
  • 标题: Data-Driven Structured Controller Design Using the Matrix S-Procedure
  • 作者: Zhaohua Yang, Yuxing Zhong, Nachuan Yang, Xiaoxu Lyu, Ling Shi
  • 分类: math.OC (数学优化与控制)
  • 发表时间: 2025年10月14日 (Preprint submitted to Automatica)
  • 论文链接: https://arxiv.org/abs/2503.14949

摘要

本文专注于离散时间线性时不变(LTI)系统的数据驱动最优结构化控制器设计,同时考虑H2H_2性能和HH_\infty性能。具体研究三种场景:(i)基于模型的结构化控制,(ii)数据驱动的非结构化控制,(iii)数据驱动的结构化控制。对于H2H_2性能,主要研究场景(ii)和(iii),因为场景(i)已在文献中得到广泛研究。对于HH_\infty性能,考虑所有三种场景。针对结构化控制,引入了一种线性化技术,将原始非凸问题转换为半定规划(SDP)问题,并基于此开发了迭代线性矩阵不等式(ILMI)算法。

研究背景与动机

问题背景

  1. 网络化控制系统的挑战:在大规模网络化控制系统中,集中式控制策略会带来巨大的通信负担,因此需要设计具有稀疏结构的控制器。
  2. 结构化控制器设计的困难:带有结构约束的控制器设计具有NP-hard的复杂性,现有方法主要分为基于LMI和基于梯度的方法。
  3. 数据驱动控制的需求:在实际应用中,许多系统的动力学模型往往无法获得,这促使了数据驱动控制的研究。

现有方法的局限性

  1. 传统方法依赖已知模型:现有的结构化控制器设计方法假设系统模型已知。
  2. 数据驱动方法的保守性:现有数据驱动方法如Miller等15的工作存在以下问题:
    • 将瞬时有界噪声过度近似为能量有界噪声序列,引入额外保守性
    • 仅在非凸可行域的凸子集内设计结构化控制器
    • 缺乏性能单调性保证

核心贡献

  1. 改进的数据驱动框架
    • 将所有可能系统矩阵的集合表示为每个数据点集合的交集,比15中的集合更小
    • 利用矩阵S-过程处理所有可能的系统矩阵,证明了可行域包含15中的可行域
    • 证明了性能相对于数据长度的单调非递减性
  2. 新颖的线性化技术
    • 开发了将原始非凸问题转换为SDP问题的线性化技术
    • 提出了ILMI算法来计算结构化控制器
    • 能够探索比15更大的可行域
  3. 全面的HH_\infty控制研究
    • 首次对基于模型的结构化控制、数据驱动非结构化控制和数据驱动结构化控制进行了全面研究

方法详解

任务定义

考虑离散时间LTI系统:

x_{k+1} = Ax_k + Bu_k + Gd_k
y_k = Cx_k + Du_k + Hd_k

其中xkRnxx_k \in \mathbb{R}^{n_x}是系统状态,ukRnuu_k \in \mathbb{R}^{n_u}是控制输入,ykRnyy_k \in \mathbb{R}^{n_y}是控制输出,dkRndd_k \in \mathbb{R}^{n_d}是外部扰动。

使用状态反馈控制器uk=Kxku_k = Kx_k,其中KRnu×nxK \in \mathbb{R}^{n_u \times n_x}需要满足结构约束KSK \in S

数据收集机制

假设矩阵A,BA, B未知但确定,矩阵C,D,G,HC, D, G, H已知。收集的数据包括:

X = [x_0, x_1, ..., x_T] ∈ ℝ^{n×(T+1)}
U = [u_0, u_1, ..., u_{T-1}] ∈ ℝ^{m×T}

关键假设:过程噪声有界,即wk2ϵ\|w_k\|_2 \leq \epsilon

核心技术创新

1. 改进的系统矩阵集合表示

基于第ii个时间步的数据(xi,ui,xi+1)(x_i, u_i, x_{i+1}),定义:

Σ_i = {(A,B) | x_{i+1} = Ax_i + Bu_i + w_i, w_i w_i^T ⪯ ε²I}

所有与数据一致的系统矩阵集合为:

Σ = ∩_{i=0}^{T-1} Σ_i

这比15中的集合更精确,减少了保守性。

2. 矩阵S-过程的应用

对于H2H_2控制,利用矩阵S-过程,条件(5a)对所有(A,B)Σ(A,B) \in Σ成立的充分条件是存在非负的α0,α1,...,αT1α_0, α_1, ..., α_{T-1}β>0β > 0使得:

P - GG^T - βI & 0 \\ 0 & -\begin{bmatrix} I \\ K \end{bmatrix} P \begin{bmatrix} I \\ K \end{bmatrix}^T \end{bmatrix} - \sum_{i=0}^{T-1} α_i Ψ_i ⪰ 0$$ #### 3. 线性化技术 对于结构化控制,关键挑战是约束$Y ⪯ P^{-1}$的非凸性。通过在给定点$\tilde{P}$附近线性化$P^{-1}$: ``` Y ⪯ P̃^{-1} - P̃^{-1}(P - P̃)P̃^{-1} ``` ### 算法框架 #### 算法1:非结构化最优控制器设计 ``` 输入:数据X, U,噪声界ε 输出:控制器K*, 性能界γ* 1. 根据问题类型选择相应的SDP问题 2. 求解SDP获得L*, P*, γ* 3. 返回K* = L*P*^{-1}, γ* ``` #### 算法2:结构化最优控制器设计 ``` 输入:数据X, U,结构约束I_S 输出:控制器K*, 性能界γ* 1. 初始化λ, μ, δ, ε, P_0 = I, k = 0 2. 重复: a. 用P̃ = P_k和λ求解SDP b. 更新K_{k+1}, P_{k+1}, γ_{k+1} c. 如果λ < δ则λ = μλ d. k = k + 1 3. 直到收敛 4. 返回K_k, γ_k ``` ## 实验设置 ### $H_2$控制案例 **系统参数**: ``` A = [-0.4095 0.4036 -0.0874 0.5154 -0.0815 0.1069 1.6715 0.7718 -0.3376] B = [ 0 0 -0.6359 -0.1098 -0.0325 2.2795] 结构约束:I_S = [1 1 0 0 1 1] ``` ### $H_\infty$控制案例 **系统参数**: ``` A = [0.8 0.2 0.1 0.1 0.7 -0.3 -0.3 0.5 0.9] B = [1 0 0 1 1 1] 结构约束:I_S = [1 1 0 1 1 0] ``` ### 评价指标 - $H_2$/$H_\infty$性能界 - 计算时间 - 与基准方法的性能对比 ## 实验结果 ### $H_2$控制结果 #### 不同噪声水平下的性能对比(T=20) | 设计方法 | (A,B)已知 | ε=0.05 | ε=0.1 | ε=0.2 | |----------|-----------|--------|-------|-------| | 非结构化[15] | 2.1537 | 2.3448 | 3.0939 | 5.2814 | | 非结构化(本文) | 2.1537 | 2.2663 | 2.5660 | 2.7889 | | 结构化[15] | 2.9794 | 3.5494 | 4.6806 | 11.3186 | | 结构化(本文) | 2.7165 | 2.9154 | 3.2249 | 4.0422 | #### 不同数据长度下的性能对比(ε=0.1) | 设计方法 | (A,B)已知 | T=6 | T=10 | T=15 | |----------|-----------|-----|------|------| | 非结构化[15] | 2.1537 | 2.9911 | 2.8156 | 2.6836 | | 非结构化(本文) | 2.1537 | 2.7494 | 2.5645 | 2.4374 | | 结构化[15] | 2.9794 | 4.4036 | 4.4323 | 4.4456 | | 结构化(本文) | 2.7165 | 3.5249 | 3.2429 | 3.1485 | ### $H_\infty$控制结果 #### 不同噪声水平下的性能对比(T=50) | 设计方法 | (A,B)已知 | ε=0.01 | ε=0.05 | ε=0.15 | |----------|-----------|--------|--------|--------| | 非结构化[15] | 0.7815 | 0.8035 | 0.9063 | 1.6483 | | 非结构化(本文) | 0.7815 | 0.7921 | 0.8207 | 1.0303 | | 结构化对角P | 3.7464 | 4.0950 | 6.6234 | 不可行 | | 结构化(本文) | 1.0580 | 1.0890 | 1.1826 | 1.5969 | ### 主要发现 1. **性能优势**:本文方法在所有测试场景下都显著优于基准方法[15] 2. **单调性**:本文方法保证性能界随数据长度T单调非递增,而[15]不具备此性质 3. **计算成本**:由于迭代过程和更多优化变量,计算时间有所增加 ## 相关工作 ### 结构化控制器设计 - **LMI方法**:Polyak等[2]、Fardad等[5]的工作 - **梯度方法**:Lin等[3]的增广拉格朗日方法 - **内点法**:Yang等[6]的新算法 ### 数据驱动控制 - **基础理论**:Willems等[7]的基本引理 - **鲁棒稳定化**:De Persis和Tesi[8]、Berberich等[9]的工作 - **矩阵S-引理应用**:van Waarde等[11]、Bisoffi等[12,13]的研究 ## 结论与讨论 ### 主要结论 1. 提出了更精确的系统矩阵集合表示,减少了保守性 2. 开发了有效的线性化技术处理结构约束 3. 证明了方法的理论优势并通过数值实验验证 ### 局限性 1. **计算复杂性**:迭代算法和额外变量增加了计算负担 2. **收敛性**:未提供算法收敛性的理论保证 3. **结构约束限制**:方法的有效性依赖于合理的结构约束 ### 未来方向 1. 算法收敛性的理论分析 2. 计算效率的进一步优化 3. 更复杂结构约束的处理 ## 深度评价 ### 优点 1. **理论贡献扎实**:改进的系统集合表示和矩阵S-过程应用具有理论价值 2. **方法创新性强**:线性化技术巧妙地处理了非凸约束 3. **实验验证充分**:在多种场景下验证了方法的有效性 4. **性能提升显著**:相比现有方法有明显的性能改进 ### 不足 1. **计算开销较大**:迭代算法的计算成本是实际应用的障碍 2. **算法收敛性缺乏保证**:理论分析不够完整 3. **适用范围限制**:方法对结构约束的合理性有要求 ### 影响力 1. **学术价值**:为数据驱动结构化控制提供了新的理论框架 2. **实用价值**:在网络化控制系统中有重要应用前景 3. **可复现性**:算法描述清晰,便于复现 ### 适用场景 - 大规模网络化控制系统 - 通信受限的分布式控制 - 系统模型未知但有历史数据的场景 - 对控制性能有较高要求的应用 ## 参考文献 论文引用了42篇相关文献,涵盖了结构化控制、数据驱动控制、鲁棒控制等多个领域的重要工作,为研究提供了坚实的理论基础。