行星轮系固有频率及动力学求解matlab程序