حل مسائل کنترل بهینه با معادلات انتگرالی یا معادلات انتگرال- دیفرانسیل به‌کمک توابع مقیاس و موجک‌های بی‌اسپلاین مکعبی

نویسندگان
دانشگاه تربیت دبیر شهید رجایی، گروه ریاضی
چکیده
در این مقاله روشی عددی مبتنی بر توابع مقیاس و موجک‌های بی‌اسپلاین مکعبی برای حل مسائل کنترل بهینه با سیستم دینامیکی معادله انتگرالی یا معادله انتگرال-دیفرانسیل بحث می‌شود. ماتریس‌های عملیاتی مشتق و انتگرال حاصل‌ضرب دو بردار موجک‌های بی‌اسپلاین مکعبی، روش هم‌محلی و قاعده انتگرال‌گیری گاوس-لژاندر برای گسسته‌سازی مسئلۀ کنترل بهینه پیوسته و تبدیل آن به یک مسئلۀ برنامه‌ریزی غیر‌خطی به‌کار گرفته می‌شود. همگرایی توابع کنترل و حالت و تابعک معیار بهینه تقریبی حاصل از روش پیشنهادی و هم‌چنین کران بالای خطای آنها به‌دست آورده می‌شوند. مثال‌های عددی کارایی، دقت و مفید بودن ایدۀ پیشنهادی را نشان می‌دهند.
کلیدواژه‌ها

عنوان مقاله English

Solving optimal control problems with integral equations or integral equations - differential with the help of cubic B-spline scaling functions and wavelets

نویسندگان English

Hamid Mesgarani
Hamid Safdari
Abolfazl Ghasemian
Shahid Rajaee Teacher Training University,
چکیده English

Introduction

Optimal control problems (OCPs) appear in a wide class of applications. In the classical control problems, the state-space equations are expressed as differential equations. Many physical systems, technology, biology, viscoelastic, electrochemical, economic, and generally the systems that have a memory effect cannot properly be described as ordinary differential equations. Hence, the equation of these systems expresses as integral equations, integro-differential equations, fractional differential equations and fractional integro-differential equations. Almost every system of controlled ordinary differential equations or controlled integro-differential equations can be modeled by a class of systems of controlled Volterra integral equations. There are many methods for solving optimal control problems with the state space of the system in the form of ordinary, fractional, and integral equations; can be mentioned the Euler-Lagrange method, the method of using Pontryagin’s maximum principle, the numerical methods based on finite difference, finite element methods, conjugate gradient method, spectral methods, the methods of continuous orthogonal functions, the operational matrices of integrals and embedding method. The method which we used in this paper is based on using the operational matrix of cubic B-spline scaling functions and wavelets with collocation method to reduce the optimal control problem governed by the nonlinear integral equation and integro-differential equation system with quadratic performance index to a nonlinear programming. The semi-orthogonal B-spline scaling functions and wavelets and their dual functions used in this paper have compact support, vanishing moments. These properties make many of the operational matrix elements be very small compared with the largest ones. These scaling functions and wavelets can be represented in a closed form so working with them is easy. The convergence of control and state functions and the performance index of the optimal approximation of the proposed method and also the upper bound of the error are given.

Material and methods

In this paper, a numerical method based on cubic B-spline scaling functions and wavelets for solving optimal control problems with the dynamical system of the integral equation or the differential-integral equation is discussed. The Operational matrices of derivative and integration of the product of two cubic B-spline wavelet vectors, collocation method and Gauss-Legendre integration rule for the discretization of the continuous optimal control problem and its transformation into a problem of non-linear programming is used.

Results and discussion

We solve two examples of optimal control with the dynamical system of integral equation and two examples with the dynamical system of the integro-differential equation by using present method to demonstrate validity, applicability and the simplicity of the new technique, then compare the present method with hybrid pseudo-spectral and Legendre wavelets method. These results illustrate that the accuracy of our numerical solutions are a few better than the numerical solutions obtained in the other method and there is a good agreement between the approximate solution and exact solution. Also, the numerical results reported in the tables and convergence analysis demonstrate that the accuracy improve by increasing the. Therefore, to get more accurate results, using the larger is recommended.

Conclusion

The following conclusions were drawn from this research.

The operation matrix can be simply obtained for any basis of the approximation space and it is always available, therefore it can be applied to obtain the numerical solution of various kind of optimal control problems.
Numerical results and convergence analysis indicates that the approximation solution fairly matches with the exact solution and the upper bound of error exponentially decreases by growing of approximation space.
Due to the characteristics of the B-spline wavelet and dual of them, a nonlinear objective function can be obtained without calculating the integral.
The semi-orthogonal B-spline scaling functions and wavelets used in the present paper have the properties of compact support, vanishing moments, smoothness function and the representation by a closed-form expression. With these assumptions, time is reduced, computer memory is less occupied and the operation matrix is always available../files/site1/files/61/12.pdf

کلیدواژه‌ها English

 Operational matrix
Optimal control problems
Integral equations
Cubic scale functions and B-spline wavelets
Collocation methods
Numerical integration
1. Kuo B. C., "Automatic control systems", Prentice Hall PTR (1987). ## 2. Peyghami M. R., Hadizadeh M., Ebrahimzadeh A., "Some explicit class of hybrid methods for optimal control of Volterra integral equations", J. Inform. Comput. Sci, 7(4) (2012) 253-266. ## 3. Kochetkov Y. A., Tomshin V. K., "Optimal control of deterministic systems described by integro-differential equations", Avtomatika i Telemekhanika (1) (1978) 5-11. ## 4. Bagley R. L., Calico R. A., "Fractional order state equations for the control of viscoelastically damped structures", Journal of Guidance, Control, and Dynamics, 14(2) (1991) 304-311. ## 5. Hartley T. T., Lorenzo C. F., "Dynamics and control of initialized fractional-order systems", Nonlinear Dynamics, 29 (1) (2002), 201-233. ## 6. Neustadt L. W., Warga J., "Comments on the Paper” Optimal Control of Processes Described by Integral Equations", I” by VR Vinokurov. SIAM Journal on Control and Optimization, 8 (4) (1970) 572. ## 7. Bakke V. L., "A maximum principle for an optimal control problem with integral constraints", Journal of Optimization Theory and Applications, 13 (1) (1974) 32-55. ## 8. Belbas S. A., "A reduction method for optimal control of Volterra integral equations", Applied Mathematics and Computation, 197(2)(2008) 880-890. ## 9. Belbas S. A., "A new method for optimal control of Volterra integral equations. Applied Mathematics and Computation, 189 (2) (2007) 1902-1915. ## 10. Belbas S. A., "Iterative schemes for optimal control of Volterra integral equations", Nonlinear Analysis: Theory, Methods Applications, 37 (1) (1999) 57-79. ## 11. Tohidi E., Samadi O. R. N., "Optimal control of nonlinear Volterra integral equations via Legendre polynomials", IMA Journal of Mathematical Control and Information, 30 (1) (2012) 67-83. ## 12. Maleknejad K., Ebrahimzadeh A., "Optimal control of Volterra integro-differential systems based on Legendre wavelets and collocation method. World Academy of Science", Engineering and Technology, International Journal of Mathematical, Computational, Physical, Electrical and Computer Engineering, 8 (7) (2014) 1040-1044. ## 13. Maleknejad K., Ebrahimzadeh A., "An efficient hybrid pseudo-spectral method for solving optimal control of Volterra integral systems", Mathematical Communications, 19 (2) (2014), 417-435. ## 14. Maleknejad K., Ebrahimzadeh A., "The use of rationalized Haar wavelet collocation method for solving optimal control of Volterra integral equations", Journal of Vibration and Control, 21(10) (2015) 1958-1967. ## 15. Peygham M. R., Hadizadeh M., Ebrahimzadeh A., "Some explicit class of hybrid methods for optimal control of Volterra integral equations", J. Inform. Comput. Sci, 7 (4) (2012) 253-266. ## 16. Shamsi M., "A modified pseudospectral scheme for accurate solution of Bang‐Bang optimal control problems", Optimal Control Applications and Methods, 32 (6) (2011) 668-680. ## 17. Maleknejad K., Almasieh H., "Optimal control of Volterra integral equations via triangular functions", J Mathematical and Computer Modelling, 53 (9) (2011) 1902-1909. ## 18. Xinxiu L., "Operational method for solving fractional differential equations using cubic B-spline approximation", International Journal of Computer Mathematics, 91 (12) (2014) 2584-2602. ## 19. Maleknejad K., Nosrati Sahlan M., "The method of moments for solution of second kind Fredholm integral equations based on B-spline wavelets. International Journal of Computer Mathematics, 87 (7) (2010) 1602-1616. ## 20. Çağlar H., Çağlar N., Özer M., "B-spline solution of non-linear singular boundary value problems arising in physiology", Chaos, Solitons Fractals, 39 (3) (2009) 1232-1237. ## 21. Xinxiu L., "Numerical solution of fractional differential equations using cubic B-spline wavelet collocation method. Communications in Nonlinear Science and Numerical Simulation, 17 (10) (2012) 3934-3946. ## 22. Xinxiu L., "Numerical solution of fractional partial differential equations using cubic Bspline wavelet collocation method", Advances in Computational Mathematics and its Applications, 1 (3) (2012) 159-164. ## 23. Maleknejad K., Mollapourasl R., Shahabi M., "On the solution of a nonlinear integral equation on the basis of a fixed point technique and cubic B-spline scaling functions", Journal of Computational and Applied Mathematics, 239 (2013) 346-358. ## 24. Edrisi-Tabri Y., Lakestani M., Heydari A., "Two numerical methods for nonlinear constrained quadratic optimal control problems using linear B-spline functions", Iranian Journal of Numerical Analysis and Optimization 6.2 (2016) 17-38. ## 25. Loghmani G. B., "B-Spline Solution of Boundary Value Problems of Fractional Order Based on Optimal Control Strategy", Journal of Sciences, Islamic Republic of Iran, 23 (1) (2012) 59-65. ## 26. Chui CK., "Wavelets: a mathematical tool for signal analysis. Society for Industrial and Applied Mathematics; (1997 Jan 1). ## 27. Rahimkhani P., Ordokhani Y., Babolian E., "An efficient approximate method for solving delay fractional optimal control problems", Nonlinear Dynamics, 86(3) (2016) 1649-1661. ## 28. Rahimkhani P., Ordokhani Y., Babolian E., "Müntz-Legendre wavelet operational matrix of fractional-order integration and its applications for solving the fractional pantograph differential equations", Numerical Algorithms, 77 (4) (2018) 1283-1305. ## 29. Rahimkhani P., Ordokhani Y., "Generalized fractional-order Bernoulli–Legendre functions: an effective tool for solving two-dimensional fractional optimal control problems", IMA Journal of Mathematical Control and Information (2017). ## 30. Ashpazzadeh E., Lakestani M., "Biorthogonal cubic Hermite spline multiwavelets on the interval for solving the fractional optimal control problems", Computational Methods for Differential Equations, 4 (2) (2016) 99-115. ## 31. Edrisi Tabriz Y., Heydari A., "Generalized B-spline functions‎ method‎‎ for solving optimal control problems", Computational Methods for Differential Equations, 2 (4) (2014) 243-255. ## 32. Davie A. M., Stothers A. J., "Improved bound for the complexity of matrix multiplication", Proceedings of the Royal Society of Edinburgh Section A: Mathematics, 143 (2) (2013) 351-369. ##