Search

A new reliable structure-preserving algorithm for computingH∞norm of linear continuous-time periodic systems is proposed in this paper. In the computation of theH∞norm, no Riccati differential equations are needed to solve and only eigenvalues of a monodromy matrix of the associated periodic Hamiltonian system will be evaluated. First, the monodromy matrix is expressed as the product of state transition matrices of the Hamiltonian system. Second, these state transition matrices, which have been proved to be symplectic matrices, are evaluated by a structure-preserving Magnus series method. Then, in order to preserve the standard symplectic form of the monodromy matrix, the structure-preserving matrices obtained by state transition matrices are employed to compute the monodromy matrix. At last, the effectiveness and the high accuracy of the proposed structure-preserving algorithm are demonstrated by numerical examples. [ABSTRACT FROM PUBLISHER]

A semianalytical spectral element method (SEM) is proposed for electromagnetic simulations of 3-D layered structures. 2-D spectral elements are employed to discretize the cross section of a layered structure, and the Legendre transformation is then used to cast the semidiscretized problem from the Lagrangian system into the Hamiltonian system. A Riccati equation-based high precision integration method is utilized to perform integration along the longitudinal direction, which is the undiscretized direction, to generate the stiffness matrix of the whole layered structure. The final system of equations by the semianalytical SEM will take the form of a set of linear equations with a block tri-diagonal matrix, which can be solved efficiently by the block Thomas algorithm. Numerical examples demonstrate the high efficiency and accuracy of the proposed method. [ABSTRACT FROM AUTHOR]