A bilinear pseudo-spectral method (BPSM) is proposed for solving two-dimensional parabolic optimal control problems (OCPs). Firstly, the OCP is converted to a partial differential equation system including the state equation of the main problem, the adjoint equation, and the gradient equation which should be solved. Secondly, the coupled system is discretized in the space domain by a BPSM based on Chebyshev polynomials and a coupled matrix differential equation is obtained. Lastly, by using the eigenvalue decomposition method, the equation is transformed to a set of independent coupled ordinary differential equation system in the new coordinate with opposite orientation and after another eigenvalue decomposition technique, these transformed equations are discretized in the time space by some higher-order formulas. The error estimate is investigated theoretically. Numerical results confirm the desirable accuracy and efficiency of the method.

Keywords