M. Aliane, Mohand Bentobache, N. Moussouni, P. Marthon
In this work, we have proposed a new approach for solving the linear-quadratic optimal control problem, where the quality criterion is a quadratic function, which can be convex or non-convex. In this approach, we transform the continuous optimal control problem into a quadratic optimization problem using the Cauchy discretization technique, then we solve it with the active-set method. In order to study the efficiency and the accuracy of the proposed approach, we developed an implementation with MATLAB, and we performed numerical experiments on several convex and non-convex linear-quadratic optimal control problems. The obtained simulation results show that our method is more accurate and more efficient than the method using the classical Euler discretization technique. Furthermore, it was shown that our method fastly converges to the optimal control of the continuous problem found analytically using the Pontryagin's maximum principle.
{"title":"Direct method to solve linear-quadratic optimal control problems","authors":"M. Aliane, Mohand Bentobache, N. Moussouni, P. Marthon","doi":"10.3934/NACO.2021002","DOIUrl":"https://doi.org/10.3934/NACO.2021002","url":null,"abstract":"In this work, we have proposed a new approach for solving the linear-quadratic optimal control problem, where the quality criterion is a quadratic function, which can be convex or non-convex. In this approach, we transform the continuous optimal control problem into a quadratic optimization problem using the Cauchy discretization technique, then we solve it with the active-set method. In order to study the efficiency and the accuracy of the proposed approach, we developed an implementation with MATLAB, and we performed numerical experiments on several convex and non-convex linear-quadratic optimal control problems. The obtained simulation results show that our method is more accurate and more efficient than the method using the classical Euler discretization technique. Furthermore, it was shown that our method fastly converges to the optimal control of the continuous problem found analytically using the Pontryagin's maximum principle.","PeriodicalId":44957,"journal":{"name":"Numerical Algebra Control and Optimization","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"80217227","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This paper discusses exact (approximate) controllability and exact (approximate) observability of stochastic implicit systems in Banach spaces. Firstly, we introduce the stochastic GE-evolution operator in Banach space and discuss existence and uniqueness of the mild solution to stochastic implicit systems by stochastic GE-evolution operator in Banach space. Secondly, we discuss conditions for exact (approximate) controllability and exact (approximate) observability of the systems considered in terms of stochastic GE-evolution operator and the dual principle. Finally, an illustrative example is given.
{"title":"Controllability and observability of stochastic implicit systems and stochastic GE-evolution operator","authors":"Zhao-qiang Ge","doi":"10.3934/NACO.2021009","DOIUrl":"https://doi.org/10.3934/NACO.2021009","url":null,"abstract":"This paper discusses exact (approximate) controllability and exact (approximate) observability of stochastic implicit systems in Banach spaces. Firstly, we introduce the stochastic GE-evolution operator in Banach space and discuss existence and uniqueness of the mild solution to stochastic implicit systems by stochastic GE-evolution operator in Banach space. Secondly, we discuss conditions for exact (approximate) controllability and exact (approximate) observability of the systems considered in terms of stochastic GE-evolution operator and the dual principle. Finally, an illustrative example is given.","PeriodicalId":44957,"journal":{"name":"Numerical Algebra Control and Optimization","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"79128356","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this article, we define a new class of functions on Riemannian manifolds, called geodesic begin{document}$ mathcal{E} $end{document}-prequasi-invex functions. By a suitable example it has been shown that it is more generalized class of convex functions. Some of its characteristics are studied on a nonlinear programming problem. We also define a new class of sets, named geodesic slack invex set. Furthermore, a sufficient optimality condition is obtained for a nonlinear programming problem defined on a geodesic local begin{document}$ mathcal{E} $end{document}-invex set.
In this article, we define a new class of functions on Riemannian manifolds, called geodesic begin{document}$ mathcal{E} $end{document}-prequasi-invex functions. By a suitable example it has been shown that it is more generalized class of convex functions. Some of its characteristics are studied on a nonlinear programming problem. We also define a new class of sets, named geodesic slack invex set. Furthermore, a sufficient optimality condition is obtained for a nonlinear programming problem defined on a geodesic local begin{document}$ mathcal{E} $end{document}-invex set.
{"title":"Geodesic $ mathcal{E} $-prequasi-invex function and its applications to non-linear programming problems","authors":"Akhlad Iqbal, P. Kumar","doi":"10.3934/naco.2021040","DOIUrl":"https://doi.org/10.3934/naco.2021040","url":null,"abstract":"<p style='text-indent:20px;'>In this article, we define a new class of functions on Riemannian manifolds, called geodesic <inline-formula><tex-math id=\"M2\">begin{document}$ mathcal{E} $end{document}</tex-math></inline-formula>-prequasi-invex functions. By a suitable example it has been shown that it is more generalized class of convex functions. Some of its characteristics are studied on a nonlinear programming problem. We also define a new class of sets, named geodesic slack invex set. Furthermore, a sufficient optimality condition is obtained for a nonlinear programming problem defined on a geodesic local <inline-formula><tex-math id=\"M3\">begin{document}$ mathcal{E} $end{document}</tex-math></inline-formula>-invex set.</p>","PeriodicalId":44957,"journal":{"name":"Numerical Algebra Control and Optimization","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"79765985","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The current work aims at finding the approximate solution to solve the nonlinear fractional type Volterra integro-differential equation begin{document}$ begin{equation*} sumlimits_{k = 1}^{m}F_{k}(x)D^{(kalpha )}y(x)+lambda int_{0}^{x}K(x, t)D^{(alpha )}y(t)dt = g(x)y^{2}(x)+h(x)y(x)+P(x). end{equation*} $end{document} In order to solve the aforementioned equation, the researchers relied on the Bernstein polynomials besides the fractional Caputo derivatives through applying the collocation method. So, the equation becomes nonlinear system of equations. By solving the former nonlinear system equation, we get the approximate solution in form of Bernstein's fractional series. Besides, we will present some examples with the estimate of the error.
The current work aims at finding the approximate solution to solve the nonlinear fractional type Volterra integro-differential equation begin{document}$ begin{equation*} sumlimits_{k = 1}^{m}F_{k}(x)D^{(kalpha )}y(x)+lambda int_{0}^{x}K(x, t)D^{(alpha )}y(t)dt = g(x)y^{2}(x)+h(x)y(x)+P(x). end{equation*} $end{document} In order to solve the aforementioned equation, the researchers relied on the Bernstein polynomials besides the fractional Caputo derivatives through applying the collocation method. So, the equation becomes nonlinear system of equations. By solving the former nonlinear system equation, we get the approximate solution in form of Bernstein's fractional series. Besides, we will present some examples with the estimate of the error.
{"title":"Application of the bernstein polynomials for solving the nonlinear fractional type Volterra integro-differential equation with caputo fractional derivatives","authors":"Miloud Moussai","doi":"10.3934/NACO.2021021","DOIUrl":"https://doi.org/10.3934/NACO.2021021","url":null,"abstract":"The current work aims at finding the approximate solution to solve the nonlinear fractional type Volterra integro-differential equation begin{document}$ begin{equation*} sumlimits_{k = 1}^{m}F_{k}(x)D^{(kalpha )}y(x)+lambda int_{0}^{x}K(x, t)D^{(alpha )}y(t)dt = g(x)y^{2}(x)+h(x)y(x)+P(x). end{equation*} $end{document} In order to solve the aforementioned equation, the researchers relied on the Bernstein polynomials besides the fractional Caputo derivatives through applying the collocation method. So, the equation becomes nonlinear system of equations. By solving the former nonlinear system equation, we get the approximate solution in form of Bernstein's fractional series. Besides, we will present some examples with the estimate of the error.","PeriodicalId":44957,"journal":{"name":"Numerical Algebra Control and Optimization","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78215987","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this paper, we introduce and study a modified extragradient algorithm for approximating solutions of a certain class of split pseudo-monotone variational inequality problem in real Hilbert spaces. Using our proposed algorithm, we established a strong convergent result for approximating solutions of the aforementioned problem. Our strong convergent result is obtained without prior knowledge of the Lipschitz constant of the pseudo-monotone operator used in this paper, and with minimized number of projections per iteration compared to other results on split variational inequality problem in the literature. Furthermore, numerical examples are given to show the performance and advantage of our method as well as comparing it with related methods in the literature.
{"title":"A modified extragradient algorithm for a certain class of split pseudo-monotone variational inequality problem","authors":"G. Ogwo, C. Izuchukwu, O. Mewomo","doi":"10.3934/NACO.2021011","DOIUrl":"https://doi.org/10.3934/NACO.2021011","url":null,"abstract":"In this paper, we introduce and study a modified extragradient algorithm for approximating solutions of a certain class of split pseudo-monotone variational inequality problem in real Hilbert spaces. Using our proposed algorithm, we established a strong convergent result for approximating solutions of the aforementioned problem. Our strong convergent result is obtained without prior knowledge of the Lipschitz constant of the pseudo-monotone operator used in this paper, and with minimized number of projections per iteration compared to other results on split variational inequality problem in the literature. Furthermore, numerical examples are given to show the performance and advantage of our method as well as comparing it with related methods in the literature.","PeriodicalId":44957,"journal":{"name":"Numerical Algebra Control and Optimization","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81110145","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this paper, we study excess demand equilibrium problems in Euclidean spaces. Applying the Glicksberg's fixed point theorem, sufficient conditions for the existence of solutions for the reference problems are established. We introduce a concept of well-posedness, say Levitin–Polyak well-posedness in the sense of Painlevé–Kuratowski, and investigate sufficient conditions for such kind of well-posedness.
{"title":"Existence and well-posedness for excess demand equilibrium problems","authors":"L. Q. Anh, Pham Thanh Duoc, T. Q. Duy","doi":"10.3934/naco.2021043","DOIUrl":"https://doi.org/10.3934/naco.2021043","url":null,"abstract":"In this paper, we study excess demand equilibrium problems in Euclidean spaces. Applying the Glicksberg's fixed point theorem, sufficient conditions for the existence of solutions for the reference problems are established. We introduce a concept of well-posedness, say Levitin–Polyak well-posedness in the sense of Painlevé–Kuratowski, and investigate sufficient conditions for such kind of well-posedness.","PeriodicalId":44957,"journal":{"name":"Numerical Algebra Control and Optimization","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"89002903","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The present manuscript deals with a 3-D food chain ecological model incorporating three species phytoplankton, zooplankton, and fish. To make the model more realistic, we include predation delay in the fish population due to the vertical migration of zooplankton species. We have assumed that additional food is available for both the predator population, viz., zooplankton, and fish. The main motive of the present study is to analyze the impact of available additional food and predation delay on the plankton-fish dynamics. The positivity and boundedness (with and without delay) are proved to make the system biologically valid. The steady states are determined to discuss the stability behavior of non-delayed dynamics under certain conditions. Considering available additional food as a control parameter, we have estimated ranges of alternative food for maintaining the sustainability and stability of the plankton-fish ecosystem. The Hopf-bifurcation analysis is carried out by considering time delay as a bifurcation parameter. The predation delay includes complexity in the system dynamics as it passes through its critical value. The direction of Hopf-bifurcation and stability of bifurcating periodic orbits are also determined using the centre manifold theorem. Numerical simulation is executed to validate theoretical results.
{"title":"Dynamical complexity in a delayed Plankton-Fish model with alternative food for predators","authors":"R. Kaur, Amit Sharma, A. Sharma","doi":"10.3934/naco.2021036","DOIUrl":"https://doi.org/10.3934/naco.2021036","url":null,"abstract":"The present manuscript deals with a 3-D food chain ecological model incorporating three species phytoplankton, zooplankton, and fish. To make the model more realistic, we include predation delay in the fish population due to the vertical migration of zooplankton species. We have assumed that additional food is available for both the predator population, viz., zooplankton, and fish. The main motive of the present study is to analyze the impact of available additional food and predation delay on the plankton-fish dynamics. The positivity and boundedness (with and without delay) are proved to make the system biologically valid. The steady states are determined to discuss the stability behavior of non-delayed dynamics under certain conditions. Considering available additional food as a control parameter, we have estimated ranges of alternative food for maintaining the sustainability and stability of the plankton-fish ecosystem. The Hopf-bifurcation analysis is carried out by considering time delay as a bifurcation parameter. The predation delay includes complexity in the system dynamics as it passes through its critical value. The direction of Hopf-bifurcation and stability of bifurcating periodic orbits are also determined using the centre manifold theorem. Numerical simulation is executed to validate theoretical results.","PeriodicalId":44957,"journal":{"name":"Numerical Algebra Control and Optimization","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87020706","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The aim of this work is to apply the results of R. Gabasov et al. [ 4 , 14 ] to an extended class of optimal control problems in the Bolza form, with intermediate phase constraints and multivariate control. In this paper, the developed iterative numerical method avoids the discretization of the dynamical system. Indeed, by using a piecewise constant control, the problem is reduced for each iteration to a linear programming problem, this auxiliary task allows to improve the value of the quality criterion. The process is repeated until the optimal or the suboptimal control is obtained. As an application, we use this method to solve an extension of the deterministic optimal cash management model of S.P. Sethi [ 31 , 32 ]. In this extension, we assume that the bank overdrafts and short selling of stock are allowed, but within the authorized time limit. The results of the numerical example show that the optimal decision for the firm depends closely on the intermediate moment, the optimal decision for the firm is to purchase until a certain date the stocks at their authorized maximum value in order to take advantage of the returns derived from stock. After that, it sales the stocks at their authorized maximum value in order to satisfy the constraint at the intermediate moment.
{"title":"Optimal control of a dynamical system with intermediate phase constraints and applications in cash management","authors":"Mourad Azi, Mohand Ouamer Bibi","doi":"10.3934/NACO.2021005","DOIUrl":"https://doi.org/10.3934/NACO.2021005","url":null,"abstract":"The aim of this work is to apply the results of R. Gabasov et al. [ 4 , 14 ] to an extended class of optimal control problems in the Bolza form, with intermediate phase constraints and multivariate control. In this paper, the developed iterative numerical method avoids the discretization of the dynamical system. Indeed, by using a piecewise constant control, the problem is reduced for each iteration to a linear programming problem, this auxiliary task allows to improve the value of the quality criterion. The process is repeated until the optimal or the suboptimal control is obtained. As an application, we use this method to solve an extension of the deterministic optimal cash management model of S.P. Sethi [ 31 , 32 ]. In this extension, we assume that the bank overdrafts and short selling of stock are allowed, but within the authorized time limit. The results of the numerical example show that the optimal decision for the firm depends closely on the intermediate moment, the optimal decision for the firm is to purchase until a certain date the stocks at their authorized maximum value in order to take advantage of the returns derived from stock. After that, it sales the stocks at their authorized maximum value in order to satisfy the constraint at the intermediate moment.","PeriodicalId":44957,"journal":{"name":"Numerical Algebra Control and Optimization","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"82189027","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In power systems, Economic Power dispatch Problem (EPP) is an influential optimization problem which is a highly non-convex and non-linear optimization problem. In the current study, a novel version of Differential Evolution (NDE) is used to solve this particular problem. NDE algorithm enhances local and global search capability along with efficient utilization of time and space by making use of two elite features: selfadaptive control parameter and single population structure. The combined effect of these concepts improves the performance of Differential Evolution (DE) without compromising on quality of the solution and balances the exploitation and exploration capabilities of DE. The efficiency of NDE is validated by evaluating on three benchmark cases of the power system problem having constraints such as power balance and power generation along with nonsmooth cost function and is compared with other optimization algorithms. The Numerical outcomes uncovered that NDE performed well for all the benchmark cases and maintained a trade-off between convergence rate and efficiency.
{"title":"A novel differential evolution algorithm for economic power dispatch problem","authors":"Pooja","doi":"10.3934/naco.2021042","DOIUrl":"https://doi.org/10.3934/naco.2021042","url":null,"abstract":"In power systems, Economic Power dispatch Problem (EPP) is an influential optimization problem which is a highly non-convex and non-linear optimization problem. In the current study, a novel version of Differential Evolution (NDE) is used to solve this particular problem. NDE algorithm enhances local and global search capability along with efficient utilization of time and space by making use of two elite features: selfadaptive control parameter and single population structure. The combined effect of these concepts improves the performance of Differential Evolution (DE) without compromising on quality of the solution and balances the exploitation and exploration capabilities of DE. The efficiency of NDE is validated by evaluating on three benchmark cases of the power system problem having constraints such as power balance and power generation along with nonsmooth cost function and is compared with other optimization algorithms. The Numerical outcomes uncovered that NDE performed well for all the benchmark cases and maintained a trade-off between convergence rate and efficiency.","PeriodicalId":44957,"journal":{"name":"Numerical Algebra Control and Optimization","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"75737841","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this paper, we propose a new alternate gradient (AG) method to solve a class of optimization problems with orthogonal constraints. In particular, our AG method alternately takes several gradient reflection steps followed by one gradient projection step. It is proved that any accumulation point of the iterations generated by the AG method satisfies the first-order optimal condition. Numerical experiments show that our method is efficient.
{"title":"An alternate gradient method for optimization problems with orthogonality constraints","authors":"Yanmei Sun, Yakui Huang","doi":"10.3934/NACO.2021003","DOIUrl":"https://doi.org/10.3934/NACO.2021003","url":null,"abstract":"In this paper, we propose a new alternate gradient (AG) method to solve a class of optimization problems with orthogonal constraints. In particular, our AG method alternately takes several gradient reflection steps followed by one gradient projection step. It is proved that any accumulation point of the iterations generated by the AG method satisfies the first-order optimal condition. Numerical experiments show that our method is efficient.","PeriodicalId":44957,"journal":{"name":"Numerical Algebra Control and Optimization","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"77617137","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}