行星轮系特征频率计算