This article has multiple issues. Please help improve it or discuss these issues on the talk page. ( Learn how and when to remove these template messages) This article may be confusing or unclear to readers. Please help clarify the article. There might be a discussion about this on the talk page. (July 2013) ( Learn how and when to remove this message) This article may contain an excessive amount of intricate detail that may interest only a particular audience. Please help by spinning off or relocating any relevant information, and removing excessive detail that may be against Wikipedia's inclusion policy. (July 2013) ( Learn how and when to remove this message) This article may be too technical for most readers to understand. Please help improve it to make it understandable to non-experts, without removing the technical details. (February 2011) ( Learn how and when to remove this message) ( Learn how and when to remove this message)

Pseudospectral optimal control is a joint theoretical-computational method for solving optimal control problems. ^{[1] [2] [3] [4]} It combines pseudospectral (PS) theory with optimal control theory to produce a PS optimal control theory. PS optimal control theory has been used in ground and flight systems ^[1] in military and industrial applications. ^[5] The techniques have been extensively used to solve a wide range of problems such as those arising in UAV trajectory generation, missile guidance, control of robotic arms, vibration damping, lunar guidance, magnetic control, swing-up and stabilization of an inverted pendulum, orbit transfers, tether libration control, ascent guidance and quantum control. ^{[5] [6]}

There are a very large number of ideas that fall under the general banner of pseudospectral optimal control. ^[7] Examples of these are the Legendre pseudospectral method, the Chebyshev pseudospectral method, the Gauss pseudospectral method, the Ross-Fahroo pseudospectral method, the Bellman pseudospectral method, the flat pseudospectral method and many others. ^{[1] [3]} Solving an optimal control problem requires the approximation of three types of mathematical objects: the integration in the cost function, the differential equation of the control system, and the state-control constraints. ^[3] An ideal approximation method should be efficient for all three approximation tasks. A method that is efficient for one of them, for instance an efficient ODE solver, may not be an efficient method for the other two objects. These requirements make PS methods ideal because they are efficient for the approximation of all three mathematical objects. ^{[8] [9] [10]} In a pseudospectral method, the continuous functions are approximated at a set of carefully selected quadrature nodes. The quadrature nodes are determined by the corresponding orthogonal polynomial basis used for the approximation. In PS optimal control, Legendre and Chebyshev polynomials are commonly used. Mathematically, quadrature nodes are able to achieve high accuracy with a small number of points. For instance, the interpolating polynomial of any smooth function (C^{${\displaystyle \infty }$}) at Legendre–Gauss–Lobatto nodes converges in L² sense at the so-called spectral rate, faster than any polynomial rate. ^[9]

A basic pseudospectral method for optimal control is based on the covector mapping principle. ^[2] Other pseudospectral optimal control techniques, such as the Bellman pseudospectral method, rely on node-clustering at the initial time to produce optimal controls. The node clusterings occur at all Gaussian points. ^{[8] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20]}

Moreover, their structure can be highly exploited to make them more computationally efficient, as ad-hoc scaling ^[21] and Jacobian computation methods, involving dual number theory ^[22] have been developed. ^[19]

In pseudospectral methods, integration is approximated by quadrature rules, which provide the best numerical integration result. For example, with just N nodes, a Legendre-Gauss quadrature integration achieves zero error for any polynomial integrand of degree less than or equal to ${\displaystyle 2N-1}$. In the PS discretization of the ODE involved in optimal control problems, a simple but highly accurate differentiation matrix is used for the derivatives. Because a PS method enforces the system at the selected nodes, the state-control constraints can be discretized straightforwardly. All these mathematical advantages make pseudospectral methods a straightforward discretization tool for continuous optimal control problems.^{[ citation needed]}

• Bellman pseudospectral method
• Chebyshev pseudospectral method
• Covector mapping principle
• Flat pseudospectral methods
• Gauss pseudospectral method
• Legendre pseudospectral method
• Pseudospectral knotting method
• Ross–Fahroo lemma
• Ross–Fahroo pseudospectral methods
• Ross' π lemma

• How Stuff Works
• Pseudospectral optimal control: Part 1
• Pseudospectral optimal control: Part 2

• DIDO – MATLAB tool for optimal control named after Dido, the first queen of Carthage.
• GPOPS-II: General Purpose Optimal Control Software
• GESOP – Graphical Environment for Simulation and OPtimization
• OpenOCL – Open Optimal Control Library Archived 20 April 2019 at the Wayback Machine
• PROPT – MATLAB Optimal Control Software
• PSOPT – Open Source Pseudospectral Optimal Control Solver in C++ Archived 12 April 2016 at the Wayback Machine
• SPARTAN: Simple Pseudospectral Algorithm for Rapid Trajectory ANalysis
• OpenGoddard – Python Open Source Pseudospectral Optimal Control Software