This paper focuses on data-driven optimal structured controller design for discrete-time linear time-invariant (LTI) systems considering both and performance. Three scenarios are investigated: (i) model-based structured control, (ii) data-driven unstructured control, and (iii) data-driven structured control. For performance, scenarios (ii) and (iii) are primarily studied, as scenario (i) has been extensively researched in the literature. For performance, all three scenarios are considered. For structured control, a linearization technique is introduced that transforms the original non-convex problem into a semidefinite programming (SDP) problem, based on which an iterative linear matrix inequality (ILMI) algorithm is developed.
Consider a discrete-time LTI system:
x_{k+1} = Ax_k + Bu_k + Gd_k
y_k = Cx_k + Du_k + Hd_k
where is the system state, is the control input, is the controlled output, and is the external disturbance.
A state feedback controller is employed, where must satisfy structural constraints .
Assuming matrices are unknown but deterministic, while matrices are known. The collected data includes:
X = [x_0, x_1, ..., x_T] ∈ ℝ^{n×(T+1)}
U = [u_0, u_1, ..., u_{T-1}] ∈ ℝ^{m×T}
Key Assumption: Process noise is bounded, i.e., .
Based on data at the -th time step, define:
Σ_i = {(A,B) | x_{i+1} = Ax_i + Bu_i + w_i, w_i w_i^T ⪯ ε²I}
The set of all system matrices consistent with the data is:
Σ = ∩_{i=0}^{T-1} Σ_i
This is more precise than the set in 15, reducing conservatism.
For control, utilizing the matrix S-procedure, a sufficient condition for condition (5a) to hold for all is the existence of non-negative and such that:
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. Linearization Technique For structured control, the key challenge is the non-convexity of the constraint $Y ⪯ P^{-1}$. By linearizing $P^{-1}$ in the neighborhood of a given point $\tilde{P}$: ``` Y ⪯ P̃^{-1} - P̃^{-1}(P - P̃)P̃^{-1} ``` ### Algorithm Framework #### Algorithm 1: Unstructured Optimal Controller Design ``` Input: Data X, U, noise bound ε Output: Controller K*, performance bound γ* 1. Select the appropriate SDP problem based on problem type 2. Solve SDP to obtain L*, P*, γ* 3. Return K* = L*P*^{-1}, γ* ``` #### Algorithm 2: Structured Optimal Controller Design ``` Input: Data X, U, structural constraints I_S Output: Controller K*, performance bound γ* 1. Initialize λ, μ, δ, ε, P_0 = I, k = 0 2. Repeat: a. Solve SDP with P̃ = P_k and λ b. Update K_{k+1}, P_{k+1}, γ_{k+1} c. If λ < δ then λ = μλ d. k = k + 1 3. Until convergence 4. Return K_k, γ_k ``` ## Experimental Setup ### $H_2$ Control Case **System Parameters**: ``` 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] Structural Constraint: I_S = [1 1 0 0 1 1] ``` ### $H_\infty$ Control Case **System Parameters**: ``` 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] Structural Constraint: I_S = [1 1 0 1 1 0] ``` ### Evaluation Metrics - $H_2$/$H_\infty$ performance bounds - Computational time - Performance comparison with benchmark methods ## Experimental Results ### $H_2$ Control Results #### Performance Comparison Under Different Noise Levels (T=20) | Design Method | (A,B) Known | ε=0.05 | ε=0.1 | ε=0.2 | |---------------|-------------|--------|-------|-------| | Unstructured [15] | 2.1537 | 2.3448 | 3.0939 | 5.2814 | | Unstructured (This Work) | 2.1537 | 2.2663 | 2.5660 | 2.7889 | | Structured [15] | 2.9794 | 3.5494 | 4.6806 | 11.3186 | | Structured (This Work) | 2.7165 | 2.9154 | 3.2249 | 4.0422 | #### Performance Comparison Under Different Data Lengths (ε=0.1) | Design Method | (A,B) Known | T=6 | T=10 | T=15 | |---------------|-------------|-----|------|------| | Unstructured [15] | 2.1537 | 2.9911 | 2.8156 | 2.6836 | | Unstructured (This Work) | 2.1537 | 2.7494 | 2.5645 | 2.4374 | | Structured [15] | 2.9794 | 4.4036 | 4.4323 | 4.4456 | | Structured (This Work) | 2.7165 | 3.5249 | 3.2429 | 3.1485 | ### $H_\infty$ Control Results #### Performance Comparison Under Different Noise Levels (T=50) | Design Method | (A,B) Known | ε=0.01 | ε=0.05 | ε=0.15 | |---------------|-------------|--------|--------|--------| | Unstructured [15] | 0.7815 | 0.8035 | 0.9063 | 1.6483 | | Unstructured (This Work) | 0.7815 | 0.7921 | 0.8207 | 1.0303 | | Structured Diagonal P | 3.7464 | 4.0950 | 6.6234 | Infeasible | | Structured (This Work) | 1.0580 | 1.0890 | 1.1826 | 1.5969 | ### Key Findings 1. **Performance Advantage**: The proposed method significantly outperforms the benchmark method [15] across all test scenarios 2. **Monotonicity**: The proposed method guarantees monotonic non-increasing performance bounds with respect to data length T, which [15] does not possess 3. **Computational Cost**: Computational time increases due to the iterative process and additional optimization variables ## Related Work ### Structured Controller Design - **LMI Methods**: Works by Polyak et al. [2], Fardad et al. [5] - **Gradient Methods**: Augmented Lagrangian method by Lin et al. [3] - **Interior Point Methods**: Novel algorithm by Yang et al. [6] ### Data-Driven Control - **Foundational Theory**: Fundamental lemma by Willems et al. [7] - **Robust Stabilization**: Works by De Persis and Tesi [8], Berberich et al. [9] - **Matrix S-Lemma Applications**: Research by van Waarde et al. [11], Bisoffi et al. [12,13] ## Conclusions and Discussion ### Main Conclusions 1. A more precise system matrix set representation is proposed, reducing conservatism 2. An effective linearization technique is developed to handle structural constraints 3. Theoretical advantages of the method are established and validated through numerical experiments ### Limitations 1. **Computational Complexity**: The iterative algorithm and additional variables increase computational burden 2. **Convergence Guarantees**: Theoretical convergence guarantees for the algorithm are not provided 3. **Structural Constraint Limitations**: The method's effectiveness depends on reasonable structural constraints ### Future Directions 1. Theoretical analysis of algorithm convergence 2. Further optimization of computational efficiency 3. Handling of more complex structural constraints ## In-Depth Evaluation ### Strengths 1. **Solid Theoretical Contributions**: The improved system set representation and matrix S-procedure application have theoretical value 2. **Strong Method Innovation**: The linearization technique cleverly addresses non-convex constraints 3. **Sufficient Experimental Validation**: The method's effectiveness is verified across multiple scenarios 4. **Significant Performance Improvement**: Notable performance gains compared to existing methods ### Weaknesses 1. **Large Computational Overhead**: The computational cost of the iterative algorithm is an obstacle to practical applications 2. **Lack of Convergence Guarantees**: Theoretical analysis is incomplete 3. **Limited Applicability**: The method requires reasonable structural constraints ### Impact 1. **Academic Value**: Provides a new theoretical framework for data-driven structured control 2. **Practical Value**: Significant application prospects in networked control systems 3. **Reproducibility**: Clear algorithm descriptions facilitate reproduction ### Applicable Scenarios - Large-scale networked control systems - Distributed control with communication constraints - Scenarios with unknown system models but available historical data - Applications with high control performance requirements ## References The paper cites 42 relevant references covering important works in structured control, data-driven control, robust control, and other related fields, providing a solid theoretical foundation for the research.