Nonlinear Identification for Control by Using NARMAX Models
Abstract
:1. Introduction
2. Theoretical Background
2.1. A NARMAX Model with 2 Inputs and 2 Outputs
2.2. Identification of NARMAX Parameters through Multi-Step Orthogonal Least Squares
2.3. Defining an Optimization Criterion to Find the Best Structural Indices
2.4. On the Stability of Identification Models
2.5. Optimization of Structural Indices through Cuckoo Search Algorithm
- (a)
- Each cuckoo can only lay one egg at a time and put it in a randomly chosen host nest. This means the number of cuckoos, host nests, and eggs is the same. The host nest represents a candidate solution to the optimization problem. Several host nests can simultaneously be searched for by a flock of cuckoos. This means a population of virtual cuckoos can be employed in the algorithm.
- (b)
- The best host nests, with the highest egg quality (i.e., accepted by the host bird), will be carried on to the next stage of the development (the next generation or iteration) based on an elitist survival mechanism. The cost function (such as ) measures the host nest quality, depending on its position in the search space. In this phase, the host nest is abandoned only if a host nest with a better quality was found. This strengthens the exploitation of the search space and contributes to the algorithm convergence towards an optimal point (host nest).
- (c)
- Cuckoo is a suspicious mind bird. Therefore, strangely enough, a second phase occurs. Normally, the worst host nest (in terms of cost function) should be replaced by other newly generated host nest. But this does not always happen. Sometimes, good host nests are abandoned. A probability of a nest to be agreed by the cuckoo and to become the host of its egg can be set. Thus, there is a probability of that the cuckoo chooses not to leave its egg in that nest, regardless the nest quality, in the hope of finding a better nest. The opposite number is referred to as renewal rate of the nests population. Usually, , which means the cuckoo is very selective with the future host nests of its eggs. Statistically, the majority of the nests playing the role of possible hosts are under question whether they will be preserved in the population or not. This means the exploitation of the search space is allowed, which increases the chances to find a global, or nearly global, solution.
3. Design and Implementation Details of Numerical Procedures
- identification through MS-EOLS;
- simulation of identification model by multi-step prediction (msp);
- CSA.
3.1. MS-EOLS Numerical Procedure
Algorithm 1 MS-EOLS Procedure to Identify a 2-MISO-NARMAX Model |
end for
|
3.2. Numerical Procedure to Simulate MISO-NARMAX Models
Algorithm 2 Simulation Procedure of a MISO-NARMAX Model, Based on msp Technique |
|
3.3. Numerical Procedure to Simulate MISO-NARXA Models
Algorithm 3 Simulation Procedure of a MISO-NARXA Model, Based on msp Technique |
|
3.4. CSA Numerical Procedure
Algorithm 4 Basic Numerical Procedure of CSA Metaheuristic |
|
4. Simulation Results and Discussion
4.1. Simulation Settings
4.2. Performance Results
5. Concluding Remarks
Author Contributions
Funding
Data Availability Statement
Conflicts of Interest
Acronyms
ARMAX | Auto-Regressive Moving Average with eXogenous input model |
ARMA | Auto-Regressive Moving Average model |
ARX | Auto-Regressive with eXogenous input model |
CSA | Cuckoo Search Algorithm |
ELS | Extended Least Squares |
ERR | Error Reduction Ratio |
MA | Moving Average model |
MISO | Multiple-Inputs-Single-Output (system) |
MS-EOLS | Multi-Step Extended OLS |
NARMAX | Nonlinear ARMAX model |
NARMA | Nonlinear ARMA model |
NARX | Nonlinear ARX model |
NARXA | NARX model combined with ARMA and AR models |
NSI | Nonlinear SI |
OLS | Orthogonal Least Squares |
PSO | Particle Swarm Optimization |
SI | System Identification |
SISO | Single-Input-Single-Output (system) |
SNR | Signal-to-Noise Ratio(s) |
dB | decibels (logarithmic scale) |
i/o | input/output (data, behavior, etc.) |
id | identification |
msa | multi-step ahead (prediction) |
osa | one-step ahead (prediction) |
ppd-prg | prescribed probability distribution prg |
prg | pseudo-random generator |
upd-prg | uniform probability distribution prg |
va | validation |
References
- Söderström, T.; Stoica, P. System Identification; Prentice Hall International Ltd.: London, UK, 1989; ISBN 0-I 3-881236-5. [Google Scholar]
- Ljung, L. System Identification—Theory for the User, 2nd ed.; Prentice Hall, Upper Saddle River: Bergen, NJ, USA, 1999; ISBN 978-0-13-656695-3. [Google Scholar]
- Alessandrini, M.; Falaschetti, L.; Turchetti, C. Nonlinear Dynamic System Identification in the Spectral Domain Using Particle-Bernstein Polynomials. Electronics 2022, 11, 3100. [Google Scholar] [CrossRef]
- Pal, P.S.; Kar, R.; Ghoshal, S.P. Identification of NARMAX Hammerstein Models with Performance Assessment Using Brain Storm Optimization Algorithm. Int. J. Adapt. Control Signal Process. 2016, 30, 1043–1070. [Google Scholar] [CrossRef]
- Rahrooh, A.; Shepard, S. Identification of Nonlinear Systems Using NARMAX Model. Nonlinear Anal. Theory Methods Appl. 2009, 71, E1198–E1202. [Google Scholar] [CrossRef]
- Billings, S.A.; Leontaritis, I.J. Identification of Nonlinear Systems Using Parameter Estimation Techniques. In Proceedings of the IEEE Conference Proceedings on Control and Its Applications, University of Warwick, Coventry, UK, 23–25 March 1981; pp. 183–187. [Google Scholar]
- Haber, R.; Keviczky, L. Nonlinear System Identification—Input-Output Modeling Approach; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1999; ISBN 978-0-7923-5858-9. [Google Scholar]
- Billings, S.A. Nonlinear System Identification: NARMAX Methods in the Time, Frequency, and Spatio-Temporal Domains; John Wiley & Sons: Hoboken, NJ, USA, 2013; ISBN 978-1-119-94359-4. [Google Scholar]
- Cao, Y.; Wang, Z.J.; Wang, W.D. Modeling of Weld Penetration Control System in GMAW-P Using NARMAX Methods. J. Manuf. Process. 2021, 65, 512–524. [Google Scholar] [CrossRef]
- Aggoune, L.; Chetouani, Y. Modeling of a distillation column based on NARMAX and Hammerstein models. Int. J. Model. Simul. Sci. Comput. 2017, 8, 1750034. [Google Scholar] [CrossRef]
- Guan, X.C.; Zhao, D.Y.; Zhu, Q.M. NARMAX Modelling and U-Model Control Design for Continuous Stirred Tank Reactor (CSTR). In Proceedings of the 35th Chinese Control Conference (CCC), Chengdu, China, 27–29 July 2016; pp. 1964–1969. [Google Scholar]
- Garces, H.O.; Rojas, A.J.; Arias, L.E. Selection of nonlinear structures for total radiation modeling. In Proceedings of the IEEE International Conference on Automatica (ICA-ACCA), Curico, Chile, 19–21 October 2016. [Google Scholar] [CrossRef]
- Fagundes, L.P.; Avelar, H.J.; Vincenzi, F. Improvements in Identification of Fuel Cell Temperature Model. In Proceedings of the 13th IEEE Brazilian Power Electronics Conference/1st Southern Power Electronics Conference (COBEP/SPEC), Fortaleza, Brazil, 29 November–2 December 2015. [Google Scholar] [CrossRef]
- Deng, Z.H.; Chen, Q.H.; Fu, Z.C. Data Driven NARMAX Modeling for PEMFC Air Compressor. Int. J. Hydrogen Energy 2020, 45, 20321–20328. [Google Scholar] [CrossRef]
- Jazayeri, P.; Rosehart, W.; Westwick, D.T. A Multistage Algorithm for Identification of Nonlinear Aggregate Power System Loads. IEEE Trans. Power Syst. 2007, 22, 1072–1079. [Google Scholar] [CrossRef]
- Fernandes, D.L.; Lopes, F.R.; Ayala, H.V.H. System Identification of an Elastomeric Series Elastic Actuator Using Black-Box Models. In Proceedings of the 31st Mediterranean Conference on Control and Automation (MED), Limassol, Cyprus, 26–29 June 2023; pp. 569–574. [Google Scholar]
- Pradhan, S.K.; Subudhi, B. NARMAX Modeling of a Two-Link Flexible Robot. In Proceedings of the Annual IEEE India Conference (INDICON): Engineering Sustainable Solutions, Hyderabad, India, 16–18 December 2011. [Google Scholar] [CrossRef]
- Iglesias, R.; Kyriacou, T.; Billings, S. Route Training in Mobile Robotics through System Identification. In Proceedings of the Conference of the World-Academy-of-Science-Engineering-and-Technology, Barcelona, Spain, 22–24 October 2006; Volume 15, pp. 181–186. [Google Scholar]
- Kelley, J.; Hagan, M.T. Comparison of Neural Network NARX and NARMAX Models for Multi-Step Prediction Using Simulated and Experimental Data. Expert Syst. Appl. 2024, 237, 121437. [Google Scholar] [CrossRef]
- Oruc, O.; Cook, A.; Mu, B.X. Nonlinear System Identification for Heading and Pitch Control of a Tethered Uncrewed Underwater Vehicle in Changing and Uncertain Environments. In Proceedings of the OCEANS Hampton Roads Conference, Hampton Roads, VA, USA, 17–20 October 2022. [Google Scholar] [CrossRef]
- Chiu, H.L. Identification Approach for Nonlinear MIMO Dynamics of Closed-Loop Active Magnetic Bearing System. Appl. Sci. 2022, 12, 8556. [Google Scholar] [CrossRef]
- Barbosa, M.P.S.; da Costa, D.P.; Ayala, H.V.H. Evaluation of Nonlinear System Identification to Model Piezoacoustic Transmission. In Proceedings of the 21st IFAC World Congress on Automatic Control—Meeting Societal Challenges, Berlin, Germany, 12–17 July 2020; Volume 53, pp. 8802–8807. [Google Scholar]
- Mohamad, M.S.A.; Yassin, I.M.; Adnan, R. Comparison Between PSO and OLS for NARX Parameter Estimation of a DC Motor. In Proceedings of the IEEE Symposium on Industrial Electronics and Applications (ISIEA), Kuching, China, 22–25 September 2013; pp. 27–32. [Google Scholar]
- Faieghi, M.R.; Azimi, S.M. Design an Optimized PID Controller for Brushless DC Motor by Using PSO and Based on NARMAX Identified Model with ANFIS. In Proceedings of the 12th International Conference on Computer Modelling and Simulation (UKSim), Cambridge, UK, 24–26 March 2010; pp. 16–21. [Google Scholar]
- Obeid, S.; Ahmadi, G.; Jha, R. NARMAX Identification Based Closed-Loop Control of Flow Separation over NACA 0015. Airfoil Fluids 2020, 5, 100. [Google Scholar] [CrossRef]
- Hall, R.J.; Wei, H.L.; Hanna, E. Complex Systems Modelling of UK Winter Wheat Yield. Comput. Electron. Agric. 2023, 209, 107855. [Google Scholar] [CrossRef]
- Krishnanathan, K.; Anderson, S.R.; Kadirkamanathan, V. A Data-Driven Framework for Identifying Nonlinear Dynamic Models of Genetic Parts. ACS Synth. Biol. 2012, 1, 375–384. [Google Scholar] [CrossRef] [PubMed]
- Gu, Y.L.; Yang, Y.; Wei, H.L. Nonlinear Modeling of Cortical Responses to Mechanical Wrist Perturbations Using the NARMAX Method. IEEE Trans. Biomed. Eng. 2021, 68, 948–958. [Google Scholar] [CrossRef] [PubMed]
- Piskulak, P.; Lewenstein, K. Modeling of Sleep Disordered Breathing Using NARMAX Methodology. In Proceedings of the International Mechatronics Conference—Recent Advances towards Industry 4.0, Warsaw, Poland, 16–19 September 2019; pp. 438–444. [Google Scholar]
- Sun, Y.M.; Simpson, I.; Wei, H.L.; Hanna, E. Probabilistic Seasonal Forecasts of North Atlantic Atmospheric Circulation Using Complex Systems Modelling and Comparison with Dynamical Models. Meteorol. Appl. 2024, 31, e2178. [Google Scholar] [CrossRef]
- Udaichi, K.; Nagappan, R.C.; Bhukya, S.N. Large-Scale System Identification Using Self-Adaptive Penguin Search Algorithm. IET Control Theory Appl. 2023, 17, 2292–2303. [Google Scholar] [CrossRef]
- Cheng, C.M.; Peng, Z.K.; Meng, G. Volterra-Series-Based Nonlinear System Modeling and Its Engineering Applications: A State-of-the-Art Review. Mech. Syst. Signal Process. 2017, 87, 340–364. [Google Scholar] [CrossRef]
- Culita, J.; Stefanoiu, D.; Nica, A.M. Nonlinear Identification for Control by Using HARMAX Models. In Proceedings of the 9th International Conference on Control, Decision and Information Technologies (CODIT 2023), Rome, Italy, 3–6 July 2023; pp. 1868–1875. [Google Scholar]
- Worden, K.; Becker, W.E.; Cross, E.J. On the Confidence Bounds of Gaussian Process NARX Models and Their Higher-Order Frequency Response Functions. Mech. Syst. Signal Process. 2018, 104, 188–223. [Google Scholar] [CrossRef]
- Li, P.; Wei, H.L.; Boynton, R. Nonlinear Model Identification from Multiple Data Sets Using an Orthogonal Forward Search Algorithm. J. Comput. Nonlinear Dyn. 2013, 8, 041001. [Google Scholar] [CrossRef]
- Zakaria, M.Z.; Mansor, Z.; Ahmad, R. NARMAX Model Identification Using Multi-Objective Optimization Differential Evolution. Int. J. Integr. Eng. 2018, 10, 188–203. [Google Scholar] [CrossRef]
- Yan, J.Y.; Deller, J.R. NARMAX Model Identification Using a Set-Theoretic Evolutionary Approach. Signal Process. 2016, 123, 30–41. [Google Scholar] [CrossRef]
- Zhu, Q.M.; Wang, Y.J.; Billings, S.A. Review of Rational (Total) Nonlinear Dynamic System Modelling, Identification and Control. Int. J. Syst. Sci. 2015, 46, 2122–2133. [Google Scholar] [CrossRef]
- Palanthandalam-Madapusi, H.J.; Edamana, B.; Ridley, A.J. NARMAX Identification for Space Weather Prediction Using Polynomial Radial Basis Functions. In Proceedings of the 46th IEEE Conference on Decision and Control, New Orleans, LA, USA, 12–14 December 2007; pp. 4774–4779. [Google Scholar]
- Cordova, J.; Yu, W. Two Types of Haar Wavelet Neural Networks for Nonlinear System Identification. Neural Process. Lett. 2012, 35, 283–300. [Google Scholar] [CrossRef]
- Xia, X.Z.; Han, L.J.; Cheng, L. A Compliant Elbow Exoskeleton with an SEA at Interaction Port. In Proceedings of the 30th International Conference on Neural Information Processing (ICONIP) of the Asia-Pacific-Neural-Network-Society (APNNS), Changsha, China, 20–23 November 2023; pp. 146–157. [Google Scholar]
- Bernat, J.; Kolota, J.; Superczynska, P. NARMAX Approach for the Identification of a Dielectric Electroactive Polymer Actuator. Int. J. Control Autom. Syst. 2023, 21, 3080–3090. [Google Scholar] [CrossRef]
- Watanabe, R.N.; Kohn, A.F. System Identification of a Motor Unit Pool Using a Realistic Neuromusculoskeletal Model. In Proceedings of the 5th IEEE RAS & EMBS International Conference on Biomedical Robotics and Biomechatronics (BioRob), Sao Paulo, Brazil, 12–15 August 2014; pp. 610–615. [Google Scholar]
- Falsone, A.; Piroddi, L.; Prandini, M. A Novel Randomized Approach to Nonlinear System Identification. In Proceedings of the 53rd IEEE Annual Conference on Decision and Control (CDC), Los Angeles, CA, USA, 15–17 December 2014; pp. 6516–6521. [Google Scholar]
- Yuan, X.L.; Bai, Y. Stochastic Nonlinear System Identification Using Multi-objective Multi-population Parallel Genetic Programming. In Proceedings of the 21st Chinese Control and Decision Conference, Guilin, China, 17–19 June 2009; pp. 1148–1153. [Google Scholar]
- Culita, J.; Stefanoiu, D.; Dumitrascu, A. ASTANK2: Analytical Modeling and Simulation. In Proceedings of the 20th International Conference on Control Systems and Computer Science (CSCS), Bucharest, Romania, 27–29 May 2015; pp. 141–148. [Google Scholar]
- Yang, X.S.; Deb, S. Cuckoo Search via Lévy Flights. In Proceedings of the World Congress on Nature & Biologically Inspired Computing, Coimbatore, India, 9–11 December 2009; pp. 210–214. [Google Scholar]
- Stefanoiu, D.; Borne, P.; Popescu, D.; Filip, F.G.; El Kamel, A. Optimization in Engineering Sciences—Metaheuristics, Stochastic Methods and Decision Support; John Wiley & Sons & ISTE Press: London, UK, 2014; ISBN 978-1-84821-498-9. [Google Scholar]
- Mantegna, R.N. Fast, Accurate Algorithm for Numerical Simulation of Lévy Stable Stochastic Processes. Phys. Rev. E 1994, 49, 4677–4683. [Google Scholar] [CrossRef]
- Stefanoiu, D.; Culita, J. Optimal Identification and Metaheuristic PID Control of a Two-Tank System. Electronics 2021, 10, 1101. [Google Scholar] [CrossRef]
- Oppenheim, A.V.; Schafer, R. Digital Signal Processing; Prentice Hall International Ltd.: London, UK, 1985; ISBN 0-13-214107-8. [Google Scholar]
- Proakis, J.G.; Manolakis, D.G. Digital Signal Processing. Principles, Algorithms and Applications; Prentice Hall Inc.: Upper Saddle River, NJ, USA, 1996; ISBN 0-13-394338-9. [Google Scholar]
Parameter | ||||||||
Value | 15 | 15 | 0.25 | 0.01 | 20 | 20 | 0.25 | 0.01 |
Model | 2-MISO-NARMAX | 2-MISO-NARXA | 2-MISO-ARMAX |
Runtime h | 355 | 208 | 97 |
Major Main | ||||||||||
2-MISO-NARMAX | 1 | 2 | 1 | 1 | 3 | – | 2 | 3 | 3 | 3 |
2-MISO-NARXA | 1 | 2 | – | – | – | 1 | 2 | 1 | – | – |
Major auxiliary | ||||||||||
2-MISO-NARX | 1 | 3 | 2 | – | – | – | 1 | – | – | – |
Minor Main | |||||||||||||||||
MISO-NARMAX1 | 5 | 26 | 30 | 30 | 30 | 26 | – | 26 | – | – | 10 | – | – | 4 | 30 | 8 | 8 |
MISO-NARXA1 | 10 | 15 | 9 | 24 | 10 | 8 | – | – | – | – | – | – | – | – | – | – | – |
MISO-ARMAX1 | 10 | – | 4 | 28 | – | – | – | – | – | – | – | – | – | – | – | – | 29 |
Minor auxiliary | |||||||||||||||||
MISO-NARX1 | 2 | 23 | 8 | 20 | 13 | 5 | 4 | 12 | 10 | – | – | – | – | – | – | – | – |
SISO-AR-ARMA1 | 100 | 41 | – | – | – | – | – | – | – | – | – | – | – | – | – | – | 5 |
MISO-ARX1 | 10 | – | 20 | 2 | – | – | – | – | – | – | – | – | – | – | – | – | – |
Minor main | |||||||||||||||||
MISO-NARMAX2 | 4 | – | 27 | 9 | 18 | 1 | – | 30 | 1 | 14 | 1 | 30 | 30 | 1 | 1 | 24 | 30 |
MISO-NARXA2 | 1 | 1 | 11 | 19 | 18 | 28 | – | 4 | – | – | – | – | – | – | – | – | – |
MISO-ARMAX2 | 1 | – | 17 | 30 | – | – | – | – | – | – | – | – | – | – | – | – | 30 |
Minor auxiliary | |||||||||||||||||
MISO-NARX2 | 7 | – | 27 | 15 | 22 | – | – | – | – | – | – | – | – | – | – | – | – |
SISO-AR-ARMA2 | 5 | 32 | – | – | – | – | – | – | – | – | – | – | – | – | – | – | 64 |
MISO-ARX2 | 1 | – | 1 | 29 | – | – | – | – | – | – | – | – | – | – | – | – | – |
Model | 2-MISO-NARMAX | 2-MISO-NARXA | 2-MISO-ARMAX |
Average fitness [%] | 84.07 | 85.71 | 81.38 |
Fitness [%] | 91.21 | 84.68 | 75.37 |
Model | 2-MISO-NARMAX | 2-MISO-NARXA | 2-MISO-ARMAX | |||
[%] | 75.78 | 74.74 | 76.37 | 74.57 | 78.66 | 75.07 |
[%] | 1.04 | 1.80 | 3.59 |
Model | ||||||||
2-MISO-NARMAX | 0.3908 | 0.5786 | 0.1827 | 0.2716 | 0.5286 | 0.9731 | 0.5372 | 0.7487 |
2-MISO-NARXA | 0.4181 | 1.4526 | 0.0984 | 0.1495 | 0.4050 | 0.6320 | 0.0848 | 5.5355 |
2-MISO-ARMAX | 1.0524 | 5.4345 | 0.2276 | 5.4742 | 0.4378 | 0.8409 | 0.3151 | 4.9893 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Stefanoiu, D.; Culita, J.; Voinea, A.-C.; Voinea, V. Nonlinear Identification for Control by Using NARMAX Models. Mathematics 2024, 12, 2252. https://doi.org/10.3390/math12142252
Stefanoiu D, Culita J, Voinea A-C, Voinea V. Nonlinear Identification for Control by Using NARMAX Models. Mathematics. 2024; 12(14):2252. https://doi.org/10.3390/math12142252
Chicago/Turabian StyleStefanoiu, Dan, Janetta Culita, Andreea-Cristina Voinea, and Vasilica Voinea. 2024. "Nonlinear Identification for Control by Using NARMAX Models" Mathematics 12, no. 14: 2252. https://doi.org/10.3390/math12142252
APA StyleStefanoiu, D., Culita, J., Voinea, A. -C., & Voinea, V. (2024). Nonlinear Identification for Control by Using NARMAX Models. Mathematics, 12(14), 2252. https://doi.org/10.3390/math12142252