“…Regarding that, a more concise and elegant high dimensional harmonic balance (HDHB) method was proposed [22,33] using collocation, where the collocation number equals to the number of unknown Fourier coefficients, which is much less than the sampling rate of the AFT method. The HDHB method is computationally fast [34,35], and a series of modified versions, e.g., the Chebyshev-based Time-spectral Method (C-TSM) [36], the Supplemental-frequency Harmonic Balance (SF-HB) method [37], etc., have been developed for specific problems. However, both the AFT method and the family of HDHB methods are impaired by aliasing when dealing with non-polynomial nonlinear problems [21,32].…”