On the Solutions of Optimal Control Problems Constrained by Ordinary Differential Equations with Vector-Matrix Coefficients Using FICO Xpress Mosel

A. S.  Afolabi; A. A.  Oyewale

A. S. Afolabi Department of Mathematical Sciences, Federal University of Technology, Akure, Ondo State, Nigeria
A. A. Oyewale Department of Mathematical Sciences, Federal University of Technology, Akure, Ondo State, Nigeria

Keywords: Hamiltonian Function, Vector-Matrix Coefficient, Discretization, Augmented Lagrangian Method, Conjugate Gradient Method, Fico Xpress Mosel

Abstract

This study addresses a general class of quadratic optimal control problems (OCPs) constrained by ordinary differential equations (ODEs) with vector-matrix coefficients. Due to the intractability of analytical solutions for complex dynamic systems, the focus is on developing and comparing efficient numerical methods. An analytical framework is first established by applying first-order optimality conditions to the Hamiltonian, yielding a system of first-order ODEs. The associated Riccati differential equation is then solved using a state transformation approach. For numerical solutions, the objective functional is discretized using Simpson's 1331 rule, and the system dynamics are approximated using a fifth-order implicit integration scheme. The discretized problem is reformulated as an unconstrained optimization problem via the Augmented Lagrangian Method and solved using both the CGM and FICO Xpress Mosel. Comparative results reveal that FICO Xpress Mosel provides faster convergence and greater numerical stability, especially for high-dimensional problems. These findings underscore the effectiveness of commercial solvers like FICO Xpress Mosel in solving large-scale quadratic OCPs with enhanced accuracy and efficiency.

References

[1] Afolabi, A. S. (2022). On the constrained optimal control problems with mixed constraints. International Journal of Mathematical Analysis and Modelling. Vol. 5: 25-35.
[2] Afolabi, A.S. and Olotu, O. (2020), Analytical Solutions of Optimal Control Problems with Mixed Constraints. Iconic Research and Engineering Journals. Vol. 3, Issue 12: 1–9.
[3] Akpanibah, C. E., and Ini, U. O. (2021). Refund Clause of Contributions with Predetermined Interest under CEV Model. International Journal of Mathematical Sciences and Optimization: Theory and Applications, 7(1): 1-16.
[4] Belloti, P., Berthold,T. and Neves, K. (2020). Algorithms for discrete nonlinear optimization in FICO Xpress.Technische Universität Berlin. Berlin, Germany: 10-20.
[5] Brandenberg, R. and Silbernagl, M. (2014). Implementing a unit commitment power market model in FICO Xpress Mosel. FICO Community: 1-37.
[6] Betts, J. T. (2001). Practical Methods for Optimal Control Using Nonlinear Programming. Philadelphia. PA SIAM: 23-30.
[7] Belloti, P. (2010). Dash Optimization. Xpress-MP Essentials, 2nd ed. Birmingham. UK: Dash Associates: 1-10.
[8] Forsgren, A., Gill, P. E. and Wright, M. H. (2002). Interior-Point Methods for Nonlinear Optimization. SIAM Review. Vol. 44: 525-597.
[9] Hassan, M.K., Ahamed, M.S., Alam, M.S. and Hossain, M.B. (2013). An implicit method for numerical solution of singular and stiff initial value problems. Applied Mathematics and Nonlinear Sciences. Vol. 18: 23-45.
[10] FICO® Analytics. How Mosel Fits into FICO® Xpress Optimization. Analytic Cloud Youtube channel,https://www.youtube.com/watch?v=iGMhjLi09uw. Accessed: Nov. 26, 2019
[11] Kogut, P. I. and Leugering, G.(2022). Matrix-valued
L
−
1
L
−1
-optimal controls in the coefficients of linear elliptic problems. DFG-Cluster CE315: Engineering of Advanced Materials: 1-20.
[12] Kogut,P. I., Kupenko, O. P. and Leugering, G. (2020). Optimal Control in Matrix-Valued Coefficients for Nonlinear Monotone Problems: Optimality Conditions. Part I: 1-20.
[13] Lenhart, S. and Workman, J.T. (2007). Optimal Control Applied to Biological Models. Boca Raton. FL: Chapman & Hall/CRC: 25-33.
[14] Miquel, T. (2020). Introduction to Optimal Control. Master's course material. Journal of Optimization Theory and Applications. Vol. 187: 99-124.
[15] MIT OpenCourseWare. (2024).MIT Lecture Notes on Lebesgue Integration, 2024. [Online]. Available: https://math.mit.edu/ rbm/18-102-S14/Chapter2. pdf. Accessed: Nov. 26, 2024.
[16] Mattheij, R. and Molenaar, J. (2002). Ordinary Differential Equations in Theory and Practice. Philadelphia, SIAM: 231-237.
[17] Rao, S.S. (2009). Engineering Optimization: Theory and Practice 4th ed. John Wiley & Sons, Hoboken, NJ: 2-9.
[18] Naidu, D. S. (2003). Optimal Control Systems. Boca Raton FL: CRC Press: 25-65.
[19] Olaiya, O. O. (2024). One-step fourth derivative block integrator for the numerical solution of singularly perturbed third-order boundary value problems. International Journal of Mathematical Sciences and Optimization: Theory and Applications, 10(1): 60-71.
[20] Olotu, O and Adekunle, A. I. (2010). Analytic and numeric solutions of discretized constrained optimal control problem with vector-matrix coefficients. Advanced Modeling and Optimization. Vol. 12, No. 1: 119-131.
[21] Olotu, O and Olorunsola, S. A. (2008). Convergence Ratio Profile for Optimal Control Problems Governed by Ordinary Differential Equations with Matrix Coefficients. American Journal of Applied Sciences. No. 2: 89-92.
[22] Orji, J.C., Miswato, M. and Afolabi, A.S. (2024). Optimal Control Problems Constrained by Non-Linear Ordinary Differential Equations. International Research and Innovation Journal. 1(1):134-149
[23] Oyewale, A.A. and Afolabi, A.S. (2024). On the Optimal Control Problems constrained by Ordinary Differential Equations Using Conjugate Gradient Method and FICO Xpress Mosel. Journal of Mathematical and Computational Science. ISSN:1927-5307, 14:18.
[24] Rahaman, M. H., Hasan, M. K., Li, M. A. and Alam, A. S. (2020). Implicit methods for numerical solution of singular initial value problems. Applied Mathematics and Nonlinear Sciences. Vol. 5: 18-30.
[25] Ross, I. M. and Fahroo, F. (2004). Pseudospectral Methods for Optimal Control. Journal of Guidance, Control, and Dynamics. Vol. 27, No. 3: 397-405.
[26] Walker, P. Worcester Polytechnic Institute, Handout on Lipschitz Continuity, [Online]. Available:https://users.wpi.edu/ walker/MA500/HANDOUTS/LipschitzContinuity.pdf. Accessed: Nov. 26, 2024.