A Mesh Refinement Algorithm for Solving Optimal Control Problems Using Pseudospectral Methods

A state approximation-mesh refinement algorithm is presented for determining solutions to optimal control problems using pseudospectral methods. In the method presented in this paper the polynomial approximation of the state is used to assess the difference between consecutive selections of the number of segments and the number of collocation point within each segment to determine better combinations of segments and collocation points on the subsequent grid. This process of computing the difference between polynomial approximations is continued until the difference lies below a user-specified threshold. Because the method developed in this paper combines dividing the problem into segments and determining the best number of collocation points within each segment, the approach conceived here is termed a hybrid global/local collocation method. The user-specified parameters of the method are provided and the approach is demonstrated successfully on four optimal control problems of varying complexity. It is found that the hybrid approach developed in this paper leads to a greater accuracy with lower overall computational time as compared to using a purely global approach.

[1]  Anil V. Rao,et al.  Algorithm 902: GPOPS, A MATLAB software for solving multiple-phase optimal control problems using the gauss pseudospectral method , 2010, TOMS.

[2]  Victor M. Zavala,et al.  Large-scale nonlinear programming using IPOPT: An integrating framework for enterprise-wide dynamic optimization , 2009, Comput. Chem. Eng..

[3]  P. Tsiotras,et al.  Trajectory Optimization Using Multiresolution Techniques , 2008 .

[4]  Lorenz T. Biegler,et al.  Convergence rates for direct transcription of optimal control problems using collocation at Radau points , 2008, Comput. Optim. Appl..

[5]  I. Michael Ross,et al.  Spectral Algorithm for Pseudospectral Methods in Optimal Control , 2008, Journal of Guidance, Control, and Dynamics.

[6]  Anil V. Rao,et al.  Direct Trajectory Optimization and Costate Estimation via an Orthogonal Collocation Method , 2006 .

[7]  Jorge Nocedal,et al.  Knitro: An Integrated Package for Nonlinear Optimization , 2006 .

[8]  Michael A. Saunders,et al.  SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization , 2002, SIAM J. Optim..

[9]  Arnold Neumaier,et al.  Introduction to Numerical Analysis , 2001 .

[10]  Anil V. Rao,et al.  Eigenvector approximate dichotomic basis method for solving hyper‐sensitive optimal control problems , 1999 .

[11]  J. Betts,et al.  MESH REFINEMENT IN DIRECT TRANSCRIPTION METHODS FOR OPTIMAL CONTROL , 1998 .

[12]  Hans Seywald,et al.  Finite difference scheme for automatic costate calculation , 1996 .

[13]  Gamal N. Elnagar,et al.  The pseudospectral Legendre method for discretizing optimal control problems , 1995, IEEE Trans. Autom. Control..

[14]  Anil V. Rao,et al.  Practical Methods for Optimal Control Using Nonlinear Programming , 1987 .

[15]  P. Davis,et al.  Methods of Numerical Integration , 1985 .