本文专注于离散时间线性时不变(LTI)系统的数据驱动最优结构化控制器设计,同时考虑性能和性能。具体研究三种场景:(i)基于模型的结构化控制,(ii)数据驱动的非结构化控制,(iii)数据驱动的结构化控制。对于性能,主要研究场景(ii)和(iii),因为场景(i)已在文献中得到广泛研究。对于性能,考虑所有三种场景。针对结构化控制,引入了一种线性化技术,将原始非凸问题转换为半定规划(SDP)问题,并基于此开发了迭代线性矩阵不等式(ILMI)算法。
考虑离散时间LTI系统:
x_{k+1} = Ax_k + Bu_k + Gd_k
y_k = Cx_k + Du_k + Hd_k
其中是系统状态,是控制输入,是控制输出,是外部扰动。
使用状态反馈控制器,其中需要满足结构约束。
假设矩阵未知但确定,矩阵已知。收集的数据包括:
X = [x_0, x_1, ..., x_T] ∈ ℝ^{n×(T+1)}
U = [u_0, u_1, ..., u_{T-1}] ∈ ℝ^{m×T}
关键假设:过程噪声有界,即。
基于第个时间步的数据,定义:
Σ_i = {(A,B) | x_{i+1} = Ax_i + Bu_i + w_i, w_i w_i^T ⪯ ε²I}
所有与数据一致的系统矩阵集合为:
Σ = ∩_{i=0}^{T-1} Σ_i
这比15中的集合更精确,减少了保守性。
对于控制,利用矩阵S-过程,条件(5a)对所有成立的充分条件是存在非负的和使得:
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篇相关文献,涵盖了结构化控制、数据驱动控制、鲁棒控制等多个领域的重要工作,为研究提供了坚实的理论基础。