Half-order optimally scaled Fourier expansion method for solving nonlinear dynamical system

In this study, an optimally scaled Fourier expansion method with half-order technique, referred to as OSFEM(H), is developed based on the combination of the optimal scaling and the half-order techniques. The optimal scaling technique is employed to reduce the ill-posedness which may arise from using high order Fourier expansion series to approximate periodic solutions. The half-order technique is a powerful tool with which one can use an m-order Fourier expansion to interpolate as many as up to 4m+1 points. In this paper, the best scales in the multi-scale Fourier expansion interpolation are derived based on the idea of equating the norm of each column of the interpolated matrix, such that the condition number of the preconditioned matrix is minimized. Then, the present half-order OSFEM, i.e. OSFEM(H), is used to solve the Duffing equation. It is shown that excellent results are achieved by comparing with the differential transformation method and the harmonic balance method. Numerical simulations verify the accuracy and effectiveness of the presently proposed method.