Comparative Simulation of Optimal Control and Modal Control for Three-Degree-of-Freedom Systems

In control engineering, designing efficient and stable control strategies for multi-degree-of-freedom systems (such as three-degree-of-freedom systems) is crucial. Optimal control and modal control are two classical methods. This article provides a comparative analysis from both theoretical principles and simulation performance perspectives.

System Modeling Fundamentals Three-degree-of-freedom systems typically consist of mass blocks, springs, and dampers. State-space models can be established using Lagrange's equations or Newtonian mechanics. State variables include displacements and velocities for each degree of freedom, with system matrices reflecting mass, stiffness, and damping coupling characteristics. Code Implementation Tip: In MATLAB, system matrices can be constructed using mass (M), damping (C), and stiffness (K) matrices, then converted to state-space form with `ss(A,B,C,D)` where A = [zeros(n) eye(n); -M\K -M\C] for n=3 degrees of freedom.

Optimal Control Implementation Optimal control based on Linear Quadratic Regulator (LQR) solves feedback gain matrices by minimizing a cost function that balances state errors and control energy consumption. The key lies in selecting appropriate weighting matrices—excessive control weighting may lead to sluggish response, while excessive state weighting can cause actuator saturation. Algorithm Insight: The LQR solution involves solving the algebraic Riccati equation using MATLAB's `lqr(A,B,Q,R)` function, where Q and R represent state and control weighting matrices respectively. Proper tuning requires iterative simulation with different Q/R ratios.

Modal Control Principle Modal control leverages the orthogonality of system modes to decouple differential equations through coordinate transformation. Independent PID or pole placement controllers are designed for dominant modes, making it suitable for suppressing vibrations in specific frequency bands. However, high-frequency mode truncation may introduce observation errors. Implementation Approach: Modal control implementation typically involves: 1) Eigenvalue decomposition using `[V,D]=eig(K,M)` to obtain mode shapes, 2) Designing SISO controllers for selected modes, 3) Transforming back to physical coordinates. Decoupling accuracy depends on proper mode selection.

Simulation Comparison Dimensions Dynamic Response: LQR provides globally optimized coordinated convergence across all degrees of freedom, while modal control offers more precise regulation for target modes. Robustness: LQR demonstrates stronger adaptability to parameter variations; modal control may fail due to modal spillover in strongly coupled systems. Computational Cost: LQR requires online solution of Riccati equations, whereas modal control necessitates pre-implementation modal analysis. Simulation Tip: Use MATLAB's `step()` and `bode()` functions to compare time-domain and frequency-domain responses, while robustness can be tested through parameter variation studies using Monte Carlo simulations.

Application Recommendations Optimal control suits scenarios with strict global performance requirements (e.g., spacecraft attitude control), while modal control is more appropriate for mechanical systems with known primary vibration sources (e.g., building seismic reduction). Practical designs can combine both advantages, such as hybrid strategies using LQR as the primary controller with modal control assisting low-frequency vibration suppression. Code Integration: Implement hybrid control by designing LQR for full-state regulation and adding modal control outputs through weighted summation in the control law calculation.

Simulation experiments should focus on the effects of non-ideal factors (such as sensor noise and time delays), which will be key directions for subsequent in-depth research. Consider adding noise injection and delay blocks in Simulink models to evaluate practical performance.