简化了一种求取非线性常微分方程高阶谐波解的近似解析计算方法。对平方和立方非线性项的傅里叶展开过程进行改进和简化,使计算过程变为两次矩阵运算即可完成展开过程,且两次矩阵运算过程一致,易于编程。以Duffing方程为算例,计算结果与数值方法一致,运算效率有所提高。