摘要

A new, reliable algorithm for nonnegative, stabilizing solutions for the periodic Riccati differential equation (PRDE) is proposed based on Fourier series expansion and the precise integration method (PIM). Taking full advantages of periodicity, we expand coefficient matrices of the underlying linear time-varying periodic Hamiltonian system associated with the PRDE in Fourier series, and a novel extended PIM for the transition matrix of linear time-varying periodic systems is developed by combining the doubling algorithm with the increment-storage technique. This method needs to compute the matrix exponential and its related integrals only once for all evenly divided subintervals, which greatly improves the computational efficiency. Further, by introducing the Riccati transformation, a fast recursive formula for the PRDE is derived based on the block form of the transition matrix computed by the extended PIM. Finally, two numerical examples are presented to verify the numerical accuracy and efficiency of the proposed algorithm with compared results.

  • 出版日期2017-12
  • 单位大连理工大学; 工业装备结构分析国家重点实验室