导入matlab后添加另一个gearsolve函数,使用ode45进行t为0.001:0.001:1000和刚度ky数组,即可获得齿轮组的全自由度运动图像