The reluctivity matrix is obtained by using the trigonometric function transformation and the form of coefficient matrix is simplified.

该方法适用于任意阶多述化曲线非线性谐波问 题求解，并提高了计算精度。